TITLE

Stability of Leapfrog Constant-Coefficients Semi-Implicit Schemes for the Fully Elastic System of Euler Equations: Flat-Terrain Case

AUTHOR(S)
Bénard, P.; Laprise, R.; Vivoda, J.; Smoíková, P.
PUB. DATE
May 2004
SOURCE
Monthly Weather Review;May2004, Vol. 132 Issue 5, p1306
SOURCE TYPE
Academic Journal
DOC. TYPE
Article
ABSTRACT
The stability of semi-implicit schemes for the hydrostatic primitive equations system has been studied extensively over the past 20 yr, since this temporal scheme and this system represented a standard for NWP. However, with the increase of computational power, the relaxation of the hydrostatic approximation through the use of nonhydrostatic fully elastic systems is now emerging for future NWP as an attractive solution valid at any scale. In this context, several models employing the so-called Euler equations together with a constant- coefficients semi-implicit time discretization have already been developed, but no solid justification for the suitability of this algorithmic combination has been presented so far, especially from the point of view of robustness. The aim of this paper is to investigate the response of this system/scheme in terms of stability in presence of explicitly treated residual terms, as it inevitably occurs in the reality of NWP. This study is restricted to the impact of thermal and baric residual terms (metric residual terms linked to the orography are not considered here). It is shown that, conversely to what occurs with hydrostatic primitive equations, the choice of the prognostic variables used to solve the system in time is of primary importance for the robustness with Euler equations. For an optimal choice of prognostic variables, unconditionally stable schemes can be obtained (with respect to the length of the time step), but only for a smaller range of reference states than in the case of hydrostatic primitive equations. This study also indicates that (i) vertical coordinates based on geometrical height and on mass behave similarly in terms of stability for the problems examined here, and (ii) hybrid coordinates induce an intrinsic instability, the practical importance of which is, however, not completely elucidated in the theoretical context of this paper.
ACCESSION #
13019644

 

Related Articles

  • Exact Pressure Evolution Equation for Incompressible Fluids. Tessarotto, M.; Ellero, M.; Aslan, N.; Mond, M.; Nicolini, P. // AIP Conference Proceedings;12/31/2008, Vol. 1084 Issue 1, p224 

    An important aspect of computational fluid dynamics is related to the determination of the fluid pressure in isothermal incompressible fluids. In particular this concerns the construction of an exact evolution equation for the fluid pressure which replaces the Poisson equation and yields an...

  • Advanced multiple time scale molecular dynamics. Omelyan, Igor P. // Journal of Chemical Physics;9/14/2009, Vol. 131 Issue 10, p104101 

    A novel method for integrating the equations of motion in molecular dynamics simulations of many-particle systems is proposed. It is based on canonical transformations of the phase space and high-accuracy reversible decompositions of the time evolution propagator into terms with different time...

  • Gravitational-thermocapillary convection of fluid in the horizontal layer in co-current gas flow. Goncharova, O. N.; Kabov, O. A. // Doklady Physics;May2009, Vol. 54 Issue 5, p242 

    The article presents a study that investigates on the solution for gravitational-thermocapillary convection of fluids in co-current gas flows. The study was conducted using Oberbeck-Boussinesq equations wherein interaction mechanisms, fluid conviction and motions were assessed obtain fluid-gas...

  • Reply. Kasahara, Akira // Journal of the Atmospheric Sciences;8/1/2004, Vol. 61 Issue 15, p1987 

    Recently, two scientists analyzed the normal modes of a compressible, stratified, nonhydrostatic tangent-plane atmospheric model, bounded by top and bottom rigid boundaries. As seen from the case of spherical geometry, the presence of cos Coriolis terms prevents the separation of basic dynamical...

  • UNIQUE MOMENT SET FROM THE ORDER OF MAGNITUDE METHOD. Struchtrup, Henning // Kinetic & Related Models;Jun2012, Vol. 5 Issue 2, p417 

    The order of magnitude method [Struchtrup, Phys. Fluids 16, 3921-3934 (2004)] is used to construct a unique moment set for 1-D transport with scattering. Simply speaking, the method uses a series of leading order Chapman-Enskog expansions in the Knudsen number to construct the moments such that...

  • ANALYTIC DESCRIPTION OF THE VISCOUS FINGERING INTERFACE IN A ROTATING HELE-SHAW CELL. DJONDJOROV, PETER A.; VASSILEV, VASSIL M.; HADZHELAZOVA, MARIANA TS.; IVAÏLO M. MLADENOV, IVAÏLO M. // Proceedings of the Annual International Conferences on Geometry,;2012, Vol. 13, p107 

    The determination of the interface shape between two fluids dur-ing the process of viscous fingering (Saffman-Taylor instability) in a Hele-Shaw cell is addressed here. The parametric equations describing this inter-face are obtained in an explicit analytic form. With these results, a number of...

  • Conjugate flows and flat solitary waves for a continuously stratified fluid. Lamb, Kevin G.; Bangiun Wan // Physics of Fluids;Aug98, Vol. 10 Issue 8, p2061 

    Examines the use of conjugate flows in determining the vertical structure of solitary waves. Comparison between solutions gathered with and without Boussinesq approximation; Guideline for stratifications with single or double pycnocline; Measurement of non-Boussinesq effect.

  • Quasi-hydrostatic primitive equations for ocean global circulation models. Lucas, Carine; Petcu, Madalina; Rousseau, Antoine // Chinese Annals of Mathematics;Nov2010, Vol. 31 Issue 6, p939 

    Global existence of weak and strong solutions to the quasi-hydrostatic primitive equations is studied in this paper. This model, that derives from the full non-hydrostatic model for geophysical fluid dynamics in the zero-limit of the aspect ratio, is more realistic than the classical hydrostatic...

  • On the Relation between the AINV and the FAPINV Algorithms. Salkuyeh, Davod Khojasteh; Roohani, Hadi // International Journal of Mathematics & Mathematical Sciences;2009, p1 

    The approximate inverse (AINV) and the factored approximate inverse (FAPINV) are two known algorithms in the field of preconditioning of linear systems of equations. Both of these algorithms compute a sparse approximate inverse of matrix A in the factored form and are based on computing two sets...

Share

Read the Article

Courtesy of

Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics