We discuss a conjugate gradient type method -- the conjugate residual -- suitable for solving linear elliptic equations that result from discretization of complex atmospheric dynamical problems. Rotation and irregular boundaries typically lead to nonself-adjoint elliptic operators whose matrix representation on the grid is definite but not symmetric. On the other hand, most established methods for solving large sparse matrix equations depend on the symmetry and definiteness of the matrix. Furthermore, the explicit construction of the matrix can be both difficult and computationally expensive. An attractive feature of conjugate gradient methods in general is that they do not require any …
continued below
Publisher Info:
Los Alamos National Lab., NM (USA)
Place of Publication:
New Mexico
Provided By
UNT Libraries Government Documents Department
Serving as both a federal and a state depository library, the UNT Libraries Government Documents Department maintains millions of items in a variety of formats. The department is a member of the FDLP Content Partnerships Program and an Affiliated Archive of the National Archives.
Descriptive information to help identify this report.
Follow the links below to find similar items on the Digital Library.
Description
We discuss a conjugate gradient type method -- the conjugate residual -- suitable for solving linear elliptic equations that result from discretization of complex atmospheric dynamical problems. Rotation and irregular boundaries typically lead to nonself-adjoint elliptic operators whose matrix representation on the grid is definite but not symmetric. On the other hand, most established methods for solving large sparse matrix equations depend on the symmetry and definiteness of the matrix. Furthermore, the explicit construction of the matrix can be both difficult and computationally expensive. An attractive feature of conjugate gradient methods in general is that they do not require any knowledge of the matrix; and in particular, convergence of conjugate residual algorithms do not rely on symmetry for definite operators. We begin by reviewing some basic concepts of variational algorithms from the perspective of a physical analogy to the damped wave equation, which is a simple alternative to the traditional abstract framework of the Krylov subspace methods. We derive two conjugate residual schemes from variational principles, and prove that either definiteness or symmetry ensures their convergence. We discuss issues related to computational efficiency and illustrate our theoretical considerations with a test problem of the potential flow of a Boussinesq fluid flow past a steep, three-dimensional obstacle.
Physical Description
39 p.
Notes
OSTI as DE94007725; Paper copy available at OSTI: phone, 865-576-8401, or email, reports@adonis.osti.gov
This report is part of the following collection of related materials.
Office of Scientific & Technical Information Technical Reports
Reports, articles and other documents harvested from the Office of Scientific and Technical Information.
Office of Scientific and Technical Information (OSTI) is the Department of Energy (DOE) office that collects, preserves, and disseminates DOE-sponsored research and development (R&D) results that are the outcomes of R&D projects or other funded activities at DOE labs and facilities nationwide and grantees at universities and other institutions.
Smolarkiewicz, P. K. & Margolin, L. G.Variational elliptic solver for atmospheric applications,
report,
March 1, 1994;
New Mexico.
(https://digital.library.unt.edu/ark:/67531/metadc1317039/:
accessed April 18, 2024),
University of North Texas Libraries, UNT Digital Library, https://digital.library.unt.edu;
crediting UNT Libraries Government Documents Department.