CMS-Flow:Hydro Eqs

Governing Equations

The depth-averaged 2-D continuity equation may be written as

$\frac{\partial h}{\partial t}+\mathrm{\nabla }\cdot (h\mathbf{𝐔})=S$

(1)

where $h$ is the total water depth $h=\zeta +\eta$, $\eta$ is the water surface elevation, $\zeta$ is the still water depth, $\mathbf{𝐔}=(U_1,U_2)$ is the depth-averaged current velocity, $S$ is a source term due to precipitation and evaporation, and $\mathrm{\nabla }=({\mathrm{\nabla }}_1,{\mathrm{\nabla }}_2)$ is the divergence operator.

The momentum equation can be written as

$\frac{\partial (h{U}_i)}{\partial t}+\mathrm{\nabla }\cdot (h\mathbf{𝐔}{U}_i)-\mathbf{𝐁}\mathbf{𝐔}=-gh{\mathrm{\nabla }}_i\eta +\mathrm{\nabla }\cdot \left({\nu }_{t}h\mathrm{\nabla }{U}_i\right)+\frac{1}{\rho }({\tau }_{wi}+{\tau }_{Si}-{\tau }_{bi})$

(2)

where $g$ is the gravitational constant, $\mathbf{𝐁}=\left(\begin{array}{cc}0& f_c\\ -f_c& 0\\ \end{array}\right)$ where $f_c$ is the Coriolis parameter, ${\nu }_t$ is the eddy viscosity, ${\tau }_{wi}$ is the wind stress, ${\tau }_{Si}$ is the wave stresses, and ${\tau }_{bi}$ is the combined wave-current mean bed shear stress.

Numerical Methods

General Transport Equation: Discretization

All of the governing equations may be written in general form

$\underset{\text{Temporal Term}}{\underbrace{\frac{\partial (h\varphi )}{\partial t}}}+\underset{\text{Advection Term}}{\underbrace{\mathrm{\nabla }\cdot (h\mathbf{𝐔}\varphi )}}=\underset{\text{Diffusion Term}}{\underbrace{\mathrm{\nabla }\cdot \left(\mathrm{\Gamma }h\mathrm{\nabla }\varphi \right)}}+\underset{\text{Source Term}}{\underbrace{S}}$

(3)

where $\varphi$ is a general scalar, $t$ is time, $h$ is the total water depth, $\mathbf{𝐔}$ is the depth averaged current velocity, $\mathrm{\Gamma }$ is the diffusion coefficient for $\varphi$, $\mathrm{\nabla }=({\mathrm{\nabla }}_1,{\mathrm{\nabla }}_2)$ is the gradient operator, and $S$ includes all other terms. Note that in the case of the continuity and momentum equations $\varphi$ is equal to 1 and $U_i$ respectively.

Temporal Term

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

$\underset{A}{\int }\frac{\partial (h\varphi )}{\partial t}\text{d}A=\frac{\partial }{\partial t}\underset{A}{\int }h\varphi \text{d}A=\frac{h^{n+1}{\varphi }_{}^{n+1}-h^n{\varphi }_{}^n}{\mathrm{\Delta }t}\mathrm{\Delta }A$

(3)

where $\mathrm{\Delta }A$ is the cell area, and $\mathrm{\Delta }t$ is the hydrodynamic time step.

Advection Term

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

$\underset{A}{\int }\mathrm{\nabla }\cdot (h\mathbf{𝐔}\varphi )\text{d}A=\underset{L}{\oint }h\varphi (\mathbf{𝐔}\cdot \mathbf{𝐧})\text{d}L=\sum _f^{}{\overline{h}}_f\mathrm{\Delta }{l}_f{\left({\hat{n}}_i{U}_i\right)}_f{\tilde{\varphi }}_f$

(4)

where $\mathbf{𝐧}={\hat{n}}_i=({\hat{n}}_1,{\hat{n}}_2)$ is the outward unit normal on cell face f, $\mathrm{\Delta }{l}_f$ is the cell face length and ${\overline{h}}_f$ 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 $U_f=U_i\in f\mathrm{\perp }i$ the above equation simplifies to

$\sum _f^{}{\overline{h}}_f\mathrm{\Delta }{l}_f{\left({\hat{n}}_i{U}_i\right)}_f{\tilde{\varphi }}_f=\sum _f^{}{n}_f{F}_f{\tilde{\varphi }}_f$

(5)

where $F_f={\overline{h}}_f\mathrm{\Delta }{l}_f{U}_f$, $n_f=n_{\mathrm{\perp }}={\left({\hat{e}}_i{\hat{n}}_i\right)}_f$, with ${\hat{e}}^i=({\hat{e}}_1,{\hat{e}}_2)$ being the basis vector. $n_f$ is equal to 1 for West and South faces and equal to -1 for North and East cell faces. Lastly, ${\tilde{\varphi }}_f^{}$ is the advective value of $\varphi$ on cell face f, and is calculated using either the Hybrid, Exponential, HLPA (Zhu 1991) schemes. The cell face velocities $U_f$ are calculated using the momentum interpolation method of Rhie and Chow (1983) described in the subsequent section.

Cell-face interpolation operator

The general formula for estimating the cell-face value of ${\tilde{\varphi }}_f^{}$ is given by

${\overline{\varphi }}_f=L_{\mathrm{\perp }}{\varphi }_N+(1-L_{\mathrm{\perp }}){\varphi }_P+{\left({\mathrm{\nabla }}_{\parallel }\varphi \right)}_N{L}_{\mathrm{\perp }}(x_{\parallel ,O}-x_{\parallel ,N})+{\left({\mathrm{\nabla }}_{\parallel }\varphi \right)}_P(1-L_{\mathrm{\perp }})(x_{\parallel ,O}-x_{\parallel ,P})$

(6)

where $L_{\mathrm{\perp }}$ is a linear interpolation factor given by $L_{\mathrm{\perp }}=\mathrm{\Delta }{x}_{\mathrm{\perp },P}/(\mathrm{\Delta }{x}_{\mathrm{\perp },P}+\mathrm{\Delta }{x}_{\mathrm{\perp },N})$ and ${\mathrm{\nabla }}_{\parallel }$ is the gradient operator in the direction parallel to face f. By definition $\parallel =2\left|{\hat{n}}_1\right|+1\left|{\hat{n}}_2\right|$. Note that for neighboring cells without any refinement <math< {{x}_{\parallel ,O}}-{{x}_{\parallel ,P}} </math> and $x_{\parallel ,O}-x_{\parallel ,N}$ are zero and thus the above equation is consistent with non-refined cell faces.

Diffusion term

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

$\underset{A}{\int }\mathrm{\nabla }\cdot \left(\mathrm{\Gamma }h\mathrm{\nabla }\varphi \right)\text{d}A=\underset{S}{\oint }\mathrm{\Gamma }h(\mathrm{\nabla }\varphi \cdot \mathbf{𝐧})\text{d}S=\sum _f^{}{\overline{\mathrm{\Gamma }}}_f^{}{\overline{h}}_f\mathrm{\Delta }{l}_f{\left({\hat{n}}_i{\mathrm{\nabla }}_i\varphi \right)}_f$

(7)

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

$\sum _f^{}{n}_f{\overline{\mathrm{\Gamma }}}_f^{}{\overline{h}}_f\mathrm{\Delta }{l}_f{\left({\mathrm{\nabla }}_{\mathrm{\perp }}\varphi \right)}_f=\sum _f^{}{D}_f[{\varphi }_N-{\varphi }_P+{\left({\mathrm{\nabla }}_{\parallel }\varphi \right)}_N(x_{\parallel ,O}-x_{\parallel ,N})-{\left({\mathrm{\nabla }}_{\parallel }\varphi \right)}_P(x_{\parallel ,O}-x_{\parallel ,P})]$

(8)

where ${\mathrm{\nabla }}_{\mathrm{\perp }}\varphi$ is gradient in the direction perpendicular to the cell face and $D_f=\frac{{\overline{\mathrm{\Gamma }}}_f^{}{\overline{h}}_f\mathrm{\Delta }{l}_f}{\left|\delta x_{\mathrm{\perp }}\right|}$.

Cell-centered node-based gradient operator

Cell-centered face-based gradient operator

Cell-face gradient operator

Source terms

Hydrodynamic Solver

---

Symbol	Description	Units
$t$	Time	sec
$h$	Total water depth $h=\zeta +\eta$	m
$\zeta$	Still water depth	m
$\eta$	Water surface elevation with respect to the still water elevation	m
$U_j$	Current velocity in the jth direction	m/sec
$S$	Sum of Precipitation and evaporation per unit area	m/sec
$g$	Gravitational constant	m/sec2
$\rho$	Water density	kg/m3
$p_a$	Atmospheric pressure	Pa
${\nu }_t$	Turbulent eddy viscosity	m2/sec

---

Documentation Portal