Laserfiche WebLink
<br />the two different approaches with respect to stability, <br />convergence, accuracy and practicality; 'and finally, (3) to <br />compare relative run times and simulation costs tetween the <br />two methods for similar problems. The follo~ling paragraphs <br />present the fundamental equatic.ns used by the two models. <br /> <br />Governing Equations <br />-li.olh models incorporate similar fOr1'lS of the so-called phenom- <br />enological equations for monl2ntum, along V1ith the continuity <br />equation and a form of the convective-diffusion equation for <br />thermal or material transport in the vertical. Note, however, <br />that the FEM model retains the vertical n,omentum equation, <br />which is replaced by the hydrostatic pressure distribution in <br />the FDM model. Both models utilize a Cartesian coordinate <br />system with the longitudinal x dimension positive downstream. <br />The vertical z dimension is referenced positive upward from <br />the x-axis in the FEM model, while it is positive downward <br />from the x-axis in LARM. Both models allow for a variable <br />width in the lateral y direction. <br /> <br />F_EM Hydrodynamic Model <br />Momentum Equation: <br /> <br />b (ll!+ ll!+ au) + <br />p at uax waz <br /> <br />d.(PEl _ Exx.L (b-"J:!.) _ EXZ .L (b~) <br />ax ax ax az az <br /> <br />2 2 <br />+ pgC- uJulAb - ,Va Cos ~ As = D <br /> <br />(4) <br /> <br />bp (~~ + ~ + w~~) + df!J- + bpg - EXZ ~x (b~~) <br /> <br />E. a (aw) <br />- zz - b- <br />az az <br /> <br />= 0 <br /> <br />(5) <br /> <br />Continuity Equation: <br />fx-{bU) + ~z(bW) = 0 (6) <br /> <br />Convective-Diffusion Equation for Density: <br />.a.~ + b (ul.e. + wap) - D _L (b"E.) - D ~- (b~-)= 0 (7) <br />at ax az x ax ax z. az az <br /> <br />whe re <br />u, w = <br /> b = <br /> P = <br />Dx' Dz = <br /> Ab = <br /> As = <br /> <br />fluid velocity in the x and z directions respectively <br />breadth <br />pressure <br />eddy diffusion coefficients in the x and z directions <br />respectively <br />area over which bottom stress is effective <br />the area over \',hich the wind stress is effective <br /> <br />Other variables have previously been defined. <br /> <br />1;1, <br />