Use of the tensor product for numerical weather prediction by the finite element method
Abstract
This is Part 2 of a report-pair concerning application of the
tensor product in solving large sets of simultaneous linear
equations arising in finite element formulations of Numerical
Weather Prediction problems. A rectangular region having a graded mesh with Dirichlet boundary conditions on all four edges
is considered. Coefficient matrices are the "mass" matrix and
the "stiffness" matrix of the finite element method. For the stiffness matrix, which appears in Poisson's equation, operation counts and storage requirements are compared with corresponding numbers for solutions by successive over-relaxation and Gaussian elimination. FORTRAN programs for implementation of the tensor product formulations are given.
Rights
This publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.NPS Report Number
NPS69-84-005Related items
Showing items related by title, author, creator and subject.
-
Robust Approaches to Handling Complex Geometries with Galerkin Difference Methods
Kozdon, Jeremy E.; Wilcox, Lucas C.; Hagstromb, Thomas; Banks, Jeffrey W. (Cornell University Library, 2018-06-14);The Galerkin difference (GD) basis is a set of continuous, piecewise polynomials defined using a finite difference like grid of degrees of freedom. The one dimensional GD basis functions are naturally extended to multiple ... -
An LES-like stabilization of the spectral element solution of the euler equations for atmospheric flows
Marras, Simone; Giraldo, Francis X.; Müller, Andreas (2014);The solution of the Euler equations by the spectral element method (SEM) is subject to oscillatory behavior if the high-frequency modes are not damped in some way. In this analysis, we extend to high order spectral elements ... -
Some consequences of a covariant differential coordinate system
Comstock, Craig (Monterey, California: U.S. Naval Postgraduate School, 1961);The assumption is made that the concept of a covariant coordinate system is meaningful and that such a system must satisfy the differential equations a) dx(I) = g(ij)dx (j) where the g(ij) are the components of the ...