Premium
Algorithms for vector field generation in mass consistent models
Author(s) -
Flores Ciro,
Juárez Héctor,
Núñez Marco A.,
Sandoval María Luisa
Publication year - 2010
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.20458
Subject(s) - mathematics , lagrange multiplier , conjugate gradient method , vector field , multiplier (economics) , boundary value problem , partial differential equation , saddle point , boundary (topology) , mathematical analysis , mathematical optimization , geometry , economics , macroeconomics
Abstract 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 Glowinsk'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. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010