MEMO Model - Numerical Solution of The Equation System

Numerical Solution of The Equation System

The discretized equations are solved numerically on a staggered grid, i.e. the scalar quantities, and are defined at the cell centre while the velocity components, and are defined at the centre of the appropriate interfaces.

Temporal discretization of the prognostic equations is based on the explicit second order Adams-Bashforth scheme. There are two deviations from the Adams-Bashforth scheme: The first refers to the implicit treatment of the nonhydrostatic part of the mesoscale pressure perturbation . To ensure non-divergence of the flow field, an elliptic equation is solved. The elliptic equation is derived from the continuity equation wherein velocity components are expressed in terms of . Since the elliptic equation is derived from the discrete form of the continuity equation and the discrete form of the pressure gradient, conservativity is guaranteed (Flassak and Moussiopoulos, 1988). The discrete pressure equation is solved numerically with a fast elliptic solver in conjunction with a generalized conjugate gradient method. The fast elliptic solver is based on fast Fourier analysis in both horizontal directions and Gaussian elimination in the vertical direction (Moussiopoulos and Flassak, 1989).

The second deviation from the explicit treatment is related to the turbulent diffusion in vertical direction. In case of an explicit treatment of this term, the stability requirement may necessitate an unacceptable abridgement of the time increment. To avoid this, vertical turbulent diffusion is treated using the second order Crank–Nicolson method.

On principle, advective terms can be computed using any suitable advection scheme. In the present version of MEMO, a 3-D second-order total-variation-diminishing (TVD) scheme is implemented which is based on the 1-D scheme proposed by Harten (1986). It achieves a fair (but not any) reduction of numerical diffusion, the solution being independent of the magnitude of the scalar (i.e., preserving transportivity).

Read more about this topic:  MEMO Model

Famous quotes containing the words numerical, solution, equation and/or system:

    The terrible tabulation of the French statists brings every piece of whim and humor to be reducible also to exact numerical ratios. If one man in twenty thousand, or in thirty thousand, eats shoes, or marries his grandmother, then, in every twenty thousand, or thirty thousand, is found one man who eats shoes, or marries his grandmother.
    Ralph Waldo Emerson (1803–1882)

    Coming out, all the way out, is offered more and more as the political solution to our oppression. The argument goes that, if people could see just how many of us there are, some in very important places, the negative stereotype would vanish overnight. ...It is far more realistic to suppose that, if the tenth of the population that is gay became visible tomorrow, the panic of the majority of people would inspire repressive legislation of a sort that would shock even the pessimists among us.
    Jane Rule (b. 1931)

    A nation fights well in proportion to the amount of men and materials it has. And the other equation is that the individual soldier in that army is a more effective soldier the poorer his standard of living has been in the past.
    Norman Mailer (b. 1923)

    A religion so cheerless, a philosophy so sorrowful, could never have succeeded with the masses of mankind if presented only as a system of metaphysics. Buddhism owed its success to its catholic spirit and its beautiful morality.
    W. Winwood Reade (1838–1875)