Algorithms for Vector Field Generation in Mass Consistent Models

作者：Flores Ciro; Juarez Hector^*; Nunez Marco A; Luisa Sandoval Maria

来源：Numerical Methods for Partial Differential Equations, 2010, 26(4): 826-842.

DOI：10.1002/num.20458

• mass consistent model
• variational formulation
• finite element method
• elliptic problem
• saddle-point problem
• conjugate gradient

Diagnostic models in meteorology are based on the fulfillment of some time independent physical constraints as, for instance, mass conservation. A successful method to generate an adjusted wind field, based on mass conservation equation, was proposed by Sasaki and leads to the solution of an elliptic problem for the multiplier. Here we study the problem of generating an adjusted wind field from given horizontal initial velocity data, by two ways. The first one is based on orthogonal projection in Hilbert spaces and leads to the same elliptic problem but with natural boundary conditions for the multiplier. We derive from this approach the so called E algorithm. An innovative alternative proposal is obtained from a second approach where we consider the saddle point formulation of the problem avoiding boundary conditions for the multiplier and solving this problem by iterative conjugate gradient methods. This leads to an algorithm that we call the CG-algorithm, which is inspired from Glowinski's approach to solve Stokes like problems in computational fluid dynamics. Finally, the introduction of new boundary conditions for the multiplier in the elliptic problem generates better adjusted fields than those obtained with the original boundary conditions.

• 出版日期2010-7