A Modified Characteristic Finite Element Method for a Fully Nonlinear Formulation of the Semigeostrophic Flow Equations

This paper develops a fully discrete modified characteristic finite element method for a coupled system consisting of the fully nonlinear Monge-Ampere equation and a transport equation. The system is the Eulerian formulation in the dual space for B. J. Hoskins' semigeostrophic flow equations, which are widely used in meteorology to model frontogenesis. To overcome the difficulty caused by the strong nonlinearity, we first formulate (at the differential level) a vanishing moment approximation of the semigeostrophic flow equations, a methodology recently proposed by the authors [X. Feng and M. Neilan, J. Sci. Comput., 38 (2009), pp. 74-98], which involves approximating the fully nonlinear Monge-Ampere equation by a family of fourth-order quasilinear equations. We then construct a fully discrete modified characteristic finite element method for the regularized problem. It is shown that under certain mesh constraint, the proposed numerical method converges with an optimal order rate of convergence. In particular, the obtained error bounds show explicit dependence on the regularization parameter $\varepsilon$. Numerical tests are also presented to validate the theoretical results and to gauge the efficiency of the proposed fully discrete modified characteristic finite element method.

[1]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[2]  Adam M. Oberman,et al.  Exact semi-geostrophic flows in an elliptical ocean basin , 2004 .

[3]  Michael Cullen,et al.  Lagrangian Solutions of Semigeostrophic Equations in Physical Space , 2006, SIAM J. Math. Anal..

[4]  Xiaobing Feng,et al.  Vanishing Moment Method and Moment Solutions for Fully Nonlinear Second Order Partial Differential Equations , 2009, J. Sci. Comput..

[5]  Xiaobing Feng,et al.  Vanishing moment method and moment solutions for second order fully nonlinear partial differential equations , 2007, 0708.1758.

[6]  Brian J. Hoskins,et al.  The Geostrophic Momentum Approximation and the Semi-Geostrophic Equations. , 1975 .

[7]  L. Caffarelli,et al.  Fully Nonlinear Elliptic Equations , 1995 .

[8]  Xiaobing Feng,et al.  Galerkin Methods for the Fully Nonlinear Monge-Ampère Equation , 2007 .

[9]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[10]  Mario Milman,et al.  Monge Ampère equation : applications to geometry and optimization : NSF-CBMS Conference on the Monge Ampère Equation : Applications to Geometry and Optimization, July 9-13, 1997, Florida Atlantic University , 1999 .

[11]  Roland Glowinski,et al.  Numerical methods for fully nonlinear elliptic equations of the Monge-Ampère type , 2006 .

[12]  Adam M. Oberman.  Wide stencil finite difference schemes for the elliptic Monge-Ampère equation and functions of the eigenvalues of the Hessian , 2008 .

[13]  G. Barles,et al.  Convergence of approximation schemes for fully nonlinear second order equations , 1991 .

[14]  Michael Neilan,et al.  A nonconforming Morley finite element method for the fully nonlinear Monge-Ampère equation , 2010, Numerische Mathematik.

[15]  Y. Brenier.  Polar Factorization and Monotone Rearrangement of Vector-Valued Functions , 1991 .

[16]  L. D. Prussner,et al.  On the numerical solution of the equation ∂2z/∂x2 ∂2z/∂y2−(∂2z/∂x∂y)2=f and its discretizations. I , 1988 .

[17]  W. D. Evans,et al.  PARTIAL DIFFERENTIAL EQUATIONS , 1941 .

[18]  Shirley Dex,et al.  JR 旅客販売総合システム（マルス）における運用及び管理について , 1991 .

[19]  G. Loeper,et al.  A Fully Nonlinear Version of the Incompressible Euler Equations: The Semigeostrophic System , 2006, SIAM J. Math. Anal..

[20]  G. Barles,et al.  Convergence of approximation schemes for fully nonlinear second order equations , 1990, 29th IEEE Conference on Decision and Control.

[21]  Yann Brenier,et al.  Weak Existence for the Semigeostrophic Equations Formulated as a Coupled Monge-Ampère/Transport Problem , 1998, SIAM J. Appl. Math..

[22]  R. Salmon,et al.  Geophysical Fluid Dynamics , 2019, Classical Mechanics in Geophysical Fluid Dynamics.

[23]  John Norbury,et al.  Generalised Lagrangian solutions for atmospheric and oceanic flows , 1991 .

[24]  Roland Glowinski,et al.  Numerical methods for fully nonlinear elliptic equations , 2009 .

[25]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[26]  G. Loeper.  A pr 2 00 5 A fully non-linear version of the incompressible Euler equations : the semi-geostrophic system , 2006 .

[27]  P. Lions,et al.  User’s guide to viscosity solutions of second order partial differential equations , 1992, math/9207212.

[28]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[29]  Xiaobing Feng,et al.  Analysis of Galerkin Methods for the Fully Nonlinear Monge-Ampère Equation , 2007, J. Sci. Comput..

[30]  Xiaobing Feng,et al.  Mixed Finite Element Methods for the Fully Nonlinear Monge-Ampère Equation Based on the Vanishing Moment Method , 2007, SIAM J. Numer. Anal..

[31]  A. Majda.  Introduction to PDEs and Waves in Atmosphere and Ocean , 2003 .