A multidomain method for the solution of elliptic CFD problems with an ADI scheme is described. Two methods of treatment of internal boundary conditions for ADI functions are discussed, namely an explicit and a semiimplicit method. Stability conditions for the proposed methods are derived theoretica
A factored implicit scheme for numerical weather prediction
โ Scribed by J. M. Augenbaum; S. E. Cohn; E. Isaacson; D. P. Dee; D. Marchesin
- Publisher
- John Wiley and Sons
- Year
- 1985
- Tongue
- English
- Weight
- 662 KB
- Volume
- 38
- Category
- Article
- ISSN
- 0010-3640
No coin nor oath required. For personal study only.
โฆ Synopsis
The nonlinear partial differential equations of atmospheric dynamics govern motion on two time scales, a fast one and a slow one. Only the slow-scale motions are relevant in predicting the evolution of large weather patterns. Implicit numerical methods are therefore attractive for weather prediction, since they permit a large time step chosen to resolve only the slow motions, To develop an implicit method which is efficient for problems in more than one spatial dimension, one must approximate the problem by smaller, usually one-dimensional problems. A popular way to do so is to approximately factor the multidimensional implicit operator into one-dimensional operators. The factorization error incurred in such methods, however, is often unacceptably large for problems with multiple time scales.
We propose a new factorization method for numerical weather prediction which is based on factoring separately the fast and slow parts of the implicit operator. We show analytically that the new method has small factorization error, which is comparable to other discretization errors of the overall scheme. The analysis is based on properties of the shallow water equations, a simple two-dimensional version of the fully three-dimensional equations of atmospheric dynamics.
๐ SIMILAR VOLUMES
Optimal approximations are derived for locally linearized partial differential equations with matrices of the same form as Galerkin B-spline approximations, but not necessarily having Galerkin equivalents. The matrix coefficients are chosen for correct frequency and poup velocity of long waves and a
Possible extensions of the present scheme to further improve efficiency are also discussed.
A fully coupled, implicit, numerical scheme has been developed for solving highly stiff systems of parabolic conservation equations. The finite-domain equations are formed by integration of the governing conservation equations, expressed in vector notation, over control volumes. The central idea is