CMS-Flow Numerical Methods:General Transport Equation: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
No edit summary
No edit summary
Line 9: Line 9:
where <math>\phi</math> is a general scalar, <math>t</math> is time, <math>h</math> is the total water depth, <math>bf{U}</math> is the depth averaged current velocity, <math> \Gamma </math> is the diffusion coefficient for <math> \phi </math>, <math> \nabla =({{\nabla }_{1}},{{\nabla }_{2}}) </math> is the gradient operator, and <math> S </math> includes all other terms. Note that in the case of the continuity and momentum equations <math> \phi </math> is equal to 1 and <math> U_i </math> respectively.  
where <math>\phi</math> is a general scalar, <math>t</math> is time, <math>h</math> is the total water depth, <math>bf{U}</math> is the depth averaged current velocity, <math> \Gamma </math> is the diffusion coefficient for <math> \phi </math>, <math> \nabla =({{\nabla }_{1}},{{\nabla }_{2}}) </math> is the gradient operator, and <math> S </math> includes all other terms. Note that in the case of the continuity and momentum equations <math> \phi </math> is equal to 1 and <math> U_i </math> respectively.  


= General Transport Equation =
== Spatial Discretization ==
== Spatial Discretization ==



Revision as of 19:40, 31 July 2014

General Transport Equation

All of the governing equations may be written in general form

 

(8)

where is a general scalar, is time, is the total water depth, is the depth averaged current velocity, is the diffusion coefficient for , is the gradient operator, and includes all other terms. Note that in the case of the continuity and momentum equations is equal to 1 and respectively.

Spatial Discretization

Temporal Discretization

The temporal term of the momentum equations is discretized using a first order implicit Euler scheme

 

(9)

where is the cell area, and is the hydrodynamic time step.

Cell-face interpolation operator

The general formula for estimating the cell-face value of is given by

 

(12)

where is a linear interpolation factor given by and is the gradient operator in the direction parallel to face f. By definition . Note that for neighboring cells without any refinement and are zero and thus the above equation is consistent with non-refined cell faces.

Advection Schemes

The advection scheme obtained using the divergence theorem as where is the outward unit normal on cell face f, is the cell face length and is the total water depth linearly interpolated to the cell face. Here the overbar indicates a cell face interpolation operator described in the following section. For Cartesian grids the cell face unit vector is always aligned with one of the Cartesian coordinates which simplifies the calculation. Defining the cell face normal velocity as the above equation simplifies to

 

(10)

where is the outward unit normal on cell face f, is the cell face length and is the total water depth linearly interpolated to the cell face. Here the overbar indicates a cell face interpolation operator described in the following section. For Cartesian grids the cell face unit vector is always aligned with one of the Cartesian coordinates which simplifies the calculation. Defining the cell face normal velocity as the above equation simplifies to

 

(11)

where , , with being the basis vector. is equal to 1 for West and South faces and equal to -1 for North and East cell faces. Lastly, is the advective value of on cell face f, and is calculated using either the Hybrid, Exponential, HLPA (Zhu 1991) schemes. The cell face velocities are calculated using the momentum interpolation method of Rhie and Chow (1983) described in the subsequent section. The advection value is calculated as , where the superscripts and indicate low and high order approximations and the superscripts and indicate either explicit and implicit treatment. The explicit term is solved directly while the implicit term is implemented through a deferred correction in which the terms are approximated using the values from the previous iteration step.


Diffusion term

The diffusion term is discretized in general form using the divergence theorem

 

(13)

The discritization of the cell-face gradient is described in the next section. On a Cartesian grid the above expression may be further simplified as

 

(14)

where is gradient in the direction perpendicular to the cell face and .