Search NASASearch

Engineering topics

Marchesin, D.

Publications and source records attributed to Marchesin, D..

A factored implicit scheme for numerical weather prediction

An implicit method is proposed to factor the nonlinear partial differential equations governing fast and slow modes of dynamic motion in numerical weather prediction schemes. The method permits separate factorization of the slow and fast modes of the implicit operator. A simple two-dimensional version of the system of three-dimensional equations governing atmospheric dynamics over shallow water was analyzed to assess the accuracy of the proposed method. It is shown that the method has a small error which is comparable to other discretization errors in the overall scheme.

Augenbaum, J. M.

A fully implicit scheme for the barotropic primitive equations

An efficient implicit finite-difference method is developed and tested for a global barotropic model. The scheme requires, at each time step, the solution of only one-dimensional block-tridiagonal linear systems. This additional computation is offset by the use of a time step chosen independently of the mesh spacing. The method is second-order accurate in time and fourth-order accurate in space. Present experience indicates that this implicit method is practical for numerical simulation on fine meshes.

Cohn, S. E.

A Factored Implicit Scheme for Numerical Weather Prediction with Small Factorization Error

Numerical results show that, for large time steps, the factorization error can be significant, even for the slowly propagating Rossby modes. A new scheme is formulated based on a more accurate factorization of the equations. By grouping separately the terms of the equations which give rise to the fast and slow motion, the equations are factored more accurately. The fast-slow factorization eliminated the factorization error. If each of the fast and slow factors are factored again according to spatial components, the resulting scheme only involves the solution of one dimensional linear systems, and computational efficient. It is shown that the factorization error for the slow made component is negligible for this new scheme.

Augenbaum, J. M.

A fully implicit scheme for global numerical weather prediction

A fast-slow factored scheme is presented for use with shallow-water primitive equation numerical weather prediction models. The technique was developed to reduce the rotational mode errors which arise when the fast and slow terms of the governing differential equations are treated simultaneously. The method factors out the fast and slow terms along the coordinate directions by means of a modified Crank-Nicolson scheme. A finite-difference spatial discretization is carried out in the zonal and meridional directions to reduce the factorization error to near-zero, and that time steps of 60-90 min can be used to obtain acceptably accurate results, even in the presence of fine spatial structures in the flow.

Augenbaum, J. M.

Using exact solutions to develop an implicit scheme for the baroclinic primitive equations

The exact solutions presently obtained by means of a novel method for nonlinear initial value problems are used in the development of numerical schemes for the computer solution of these problems. The method is applied to a new, fully implicit scheme on a vertical slice of the isentropic baroclinic equations. It was not possible to find a global scale phenomenon that could be simulated by the baroclinic primitive equations on a vertical slice.

Marchesin, D.

The effect of compact implicit differencing in a baroclinic primitive equations model

In an effort to improve forecast accuracy, the horizontal accuracy of gridpoint forecast models at a number of numerical weather prediction (NWP) centers was increased from second order to fourth order. The current GLAS model uses second-order explicit finite-difference formulas in the vertical and fourth-order explicit finite-difference formulas in the horizontal. Experiments are described which indicate the increased forecast accuracy gained by use of compact fourth-order finite differences in a simple baroclinic model. Tables show an increase in forecast accuracy by a factor of 100-150 by use of fourth-order instead of second-order vertical discretization.

Augenbaum, J. M.

Numerical analysis of spectral properties of coupled oscillator Schroedinger operators. I - Single and double well anharmonic oscillators

Several methods for computing many eigenvalues and eigenfunctions of a single anharmonic oscillator Schroedinger operator whose potential may have one or two minima are described. One of the methods requires the solution of an ill-conditioned generalized eigenvalue problem. This method has the virtue of using a bounded amount of work to achieve a given accuracy in both the single and double well regions. Rigorous bounds are given, and it is proved that the approximations converge faster than any inverse power of the size of the matrices needed to compute them. The results of computations for the g:phi(4):1 theory are presented. These results indicate that the methods actually converge exponentially fast.

Isaacson, D.

Numerical methods for studying anharmonic oscillator approximations to the phi super 4 sub 2 quantum field theory

This paper is an expanded version of a talk given at the 1979 T.I.C.O.M. conference. It is a self-contained introduction, for applied mathematicians and numerical analysts, to quantum mechanics and quantum field theory. It also contains a brief description of the authors' numerical approach to the problems of quantum field theory, which may best be summarized by the question; Can we compute the eigenvalues and eigenfunctions of Schrodinger operators in infinitely many variables.

Isaacson, D.