Search NASA⌕ Search

SEARCH · Search NASA

Results for “second-order accuracy”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

Second-order non-iterative ADI solution of non-linear partial differential equations

A new method for the solution of non-linear partial differential equations by an ADI procedure is described. Although the method is second order accurate in time, it does not require either iterations or predictor corrector methods to overcome the nonlinearity of the equations. Thus the computational effort required for the solution of the non-linear problem becomes similar to that required for the linear case. The method is applied to a two-dimensional 'extended Burgers equation'. Linear stability is studied, and some numerical solutions obtained. The improved accuracy obtained by the 2nd order truncation error is clearly manifested.

Wolfshtein, M.↗

Multidimensional difference schemes with fourth-order accuracy

An explicit finite-difference algorithm is presented for the solution of quasilinear divergence free multidimensional hyperbolic systems. The method consists of four steps per time level. The resulting scheme is fourth-order accurate in both space and time, though the intermediate steps are only first-order accurate. The family of schemes introduced is dissipative, and hence, suitable for both smooth flows and flows containing shocks. This method is compared, in several numerical examples, with both second-order schemes and others that are fourth order in space, but second order in time.

Turkel, E.↗

A critical study of higher-order numerical methods for solving the boundary-layer equations

A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Principles of equivalence, Eotvos experiments, and gravitational red-shift experiments - The free fall of electromagnetic systems to post-post-Coulombian order

Free fall in an external static gravitational field of a composite test body consisting of electromagnetically interacting charged particles is analyzed using the TH-epsilon-mu formalism developed by Lightman and Lee (1973). The Lightman and Lee calculation that was carried out to 'post-Coulombian' order is extended to 'post-post-Coulombian' order, and violations of the weak equivalence principle are found which result from the magnetostatic self-energy of the test body. The current generation of Eotvos experiments is shown to test the cited formalism at first order in the external gravitational potential to an accuracy of 4 by 10 to the -10th power for electrostatic energies and 6 by 10 to the -6th power for magnetostatic energies. It is concluded that although the formalism can be tested to an accuracy of 0.04 for electrostatic energies at second order in the gravitational potential, a significant second-order test for magnetostatic energies would require Eotvos experiments accurate to better than one part in 10 to the 15th power.

Haugan, M. P.↗

Application of higher-order numerical methods to the boundary-layer equations

A fourth-order method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations for both attached and separated flows. The efficiency of the present method is compared with other higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, the three-point spline methods, and a modified finite-element method. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗

Second-order estimates in singular perturbation of output regulator problems

A method of order reduction in singular perturbation problems due to Chow and Kokotovic is studied. The method is used for an output regulator problem in control theory. The quadratic term in the asymptotic expansion of Chow and Kokotovic is calculated and it is shown by some examples how this term can give evidence of the accuracy of the procedure.

Rublein, G.↗

Polynominal Interpolation Methods for Viscous Flow Calculations

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and by Hermitian (Taylor series) finite-difference discretization. The similarities and special features of these different developments are discussed. The governing systems apply for both uniform and variable meshes. Hybrid schemes resulting from two different polynomial approximations for the first and second derivatives lead to a nonuniform mesh extension of the so-called compact or Pad? difference technique (Hermite 4). A variety of fourth-order methods are described and the Hermitian approach is extended to sixth-order (Hermite 6). The appropriate spline boundary conditions are derived for all procedures. For central finite differences, this leads to a two-point, second-order accurate generalization of the commonly used three-point end-difference formula. Solutions with several spline and Hermite procedures are presented for the boundary layer equations, with and without mass transfer, and for the incompressible viscous flow in a driven cavity. Divergence and nondivergence equations are considered for the cavity. Among the fourth-order techniques, it is shown that spline 4 has the smallest truncation error. The spline 4 procedure generally requires one-quarter the number of mesh points in a given coordinate direction as a central finite-difference calculation of equal accuracy. The Hermite 6 procedure leads to remarkably accurate boundary layer solutions.

Rubin, S. G.↗

The effect of accuracy, conservation and filtering on numerical weather forecasting

Considerations leading to the numerical design of the GLAS fourth-order global atmospheric model are discussed, including changes recently introduced into the model. The computation time and memory requirements for the fourth-order model are similar to those of the present second-order GLAS model with the same 4 deg latitude, 5 deg longitude, and 9 vertical-level resolution. However, the fourth-order model forecast skill is significantly better than that of the current GLAS model, and after three days it is comparable to the 2.5 by 3 deg version of the GLAS model in the sea level pressure maps, and has less phase errors in the 500 mb maps.

Kalnay-Rivas, E.↗

Implicit Numerical Methods in Meteorology

The development of a fully implicit finite-difference model, whose time step is chosen solely to resolve accurately the physical flow of interest is discussed. The method is based on an operator factorization which reduces the dimensionality of the implicit approach: at each time step only (spatially) one-dimensional block-tridiagonal linear systems must be solved. The scheme uses two time levels and is second-order accurate in time. Compact implicit spatial differences are used, yielding fourth-order accuracy both vertically and horizontally. In addition, the development of a fully interactive computer code is discussed. With this code the user will have a choice of models, with various levels of accuracy and sophistication, which are imbedded, as subsets of the fully implicit 3D code.

Augenbaum, J.↗

Vibration analysis of rotating turbomachinery blades by an improved finite difference method

The problem of calculating the natural frequencies and mode shapes of rotating blades is solved by an improved finite difference procedure based on second-order central differences. Lead-lag, flapping and coupled bending-torsional vibration cases of untwisted blades are considered. Results obtained by using the present improved theory have been observed to be close lower bound solutions. The convergence has been found to be rapid in comparison with the classical first-order finite difference method. While the computational space and time required by the present approach is observed to be almost the same as that required by the first-order theory for a given mesh size, accuracies of practical interest can be obtained by using the improved finite difference procedure with a relatively smaller matrix size, in contrast to the classical finite difference procedure which requires either a larger matrix or an extrapolation procedure for improvement in accuracy.

Subrahmanyam, K. B.↗

An efficient method for solving the steady Euler equations

An efficient numerical procedure for solving a set of nonlinear partial differential equations is given, specifically for the steady Euler equations. Solutions of the equations were obtained by Newton's linearization procedure, commonly used to solve the roots of nonlinear algebraic equations. In application of the same procedure for solving a set of differential equations we give a theorem showing that a quadratic convergence rate can be achieved. While the domain of quadratic convergence depends on the problems studied and is unknown a priori, we show that firstand second-order derivatives of flux vectors determine whether the condition for quadratic convergence is satisfied. The first derivatives enter as an implicit operator for yielding new iterates and the second derivatives indicates smoothness of the flows considered. Consequently flows involving shocks are expected to require larger number of iterations. First-order upwind discretization in conjunction with the Steger-Warming flux-vector splitting is employed on the implicit operator and a diagonal dominant matrix results. However the explicit operator is represented by first- and seond-order upwind differencings, using both Steger-Warming's and van Leer's splittings. We discuss treatment of boundary conditions and solution procedures for solving the resulting block matrix system. With a set of test problems for one- and two-dimensional flows, we show detailed study as to the efficiency, accuracy, and convergence of the present method.

Liou, M. S.↗

Finite difference operators on unstructured triangular meshes

A new form of the Laplace operator is derived which is shown to be second-order accurate when calculated at mesh triangle nodes if the cells corresponding to the nodes have a high degree of symmetry. It is demonstrated that the Laplace operator directly derived from Stokes' integral theorem is only zeroth-order accurate on these same cells. The construction of gradient and Laplace operators on unstructured triangular grids are discussed from the variational point of view. Two discrete representations of the Laplacian are compared to each other, and conclusions are drawn regarding their accuracy on a pointwise basis. Finally, geometric relations valid on arbitrary polygons are given for completeness in an appendix.

Erlebacher, G.↗

High resolution upwind schemes for the three-dimensional incompressible Navier-Stokes equations

Based on flux-difference splitting, implicit high resolution schemes are constructed for efficient computations of steady-state solutions to the three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. These schemes use first-order accurate Euler backward-time differencing and second-order central differencing for the viscous shear fluxes. Up to third-order accurate upwind differencing is achieved through a reconstruction of the solution from its cell averages. The reconstruction is accomplished by linear interpolation, where the node stencils are selected such that in regions of smooth solution the flow is highly resolved while spurious oscillations in regions of rapid changes in gradient are still suppressed. Fairly rapid convergence to steady-state solutions is attained with a completely vectorizable hybrid time-marching method. Flows around a sharp-edged delta wing are computed with the maximum accuracy of the upwind-differencing restricted to first-, second-, and third-order, to illustrate the effect of accuracy on the global and on the local vortical flow fields. The results are validated with experimental data.

Hartwich, PETER-M.↗

Third-order multi-dimensional Euler/Navier-Stokes solver

A nonoscillatory multidimensional advection code is presented which is applicable to the study of very high Reynolds number steady flow, and results are obtained for the oblique advection of a step profile. The scheme is based on an exponential upwinding or linear extrapolation refinement of a previous QUICK scheme, allowing sharp monotonic resolution of the step while retaining third-order accuracy. The resulting simple flux-limiting criteria have been applied to a higher order method and a formally third-order supercompressive algorithm. Results are presented for first-order, second-order, and third-order upwinding.

Leonard, B. P.↗

Three-dimensional computations of cross-flow injection and combustion in a supersonic flow

A low-storage version of the SPARK3D code which is based on the temporally second-order accurate MacCormack (1969) explicit scheme is used to solve the governing equations for three-dimensional chemically reacting flows with finite-rate chemistry. The code includes a fourth-order compact spatial scheme capable of providing higher order spatial accuracy, and it is used to study two-dimensional linear advection, two-dimensional Euler flow, and three-dimensional viscous flow. Also considered are the injection, mixing, and combustion of hydrogen in a supersonic cross stream.

Carpenter, M. H.↗

Solution algorithms for the two-dimensional Euler equations on unstructured meshes

The objective of the study was to analyze implicit techniques employed in structured grid algorithms for solving two-dimensional Euler equations and extend them to unstructured solvers in order to accelerate convergence rates. A comparison is made between nine different algorithms for both first-order and second-order accurate solutions. Higher-order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The discussion is illustrated by results for flow over a transonic circular arc.

Whitaker, D. L.↗

On the length of the drift region in the Ramsey cavity

The interaction of atoms in a beam with the microwave field in a separated field geometry such as a Ramsey cavity is generally described in terms of the three regions traversed successively by the atoms, namely two interaction regions of length (l) separated by a drift, or free precession, region of length L. For a monokinetic beam of velocity v, the linewidth of the central fringe in the Ramsey resonance pattern is usually expressed as delta omega equals pi v/L. A more detailed calculation shows, however, that the linewidth is equal to pi v/L asterisk, where the equivalent drift L asterisk is larger than L by an amount of the order of l/L. The correction depends on the field distribution in the interaction regions. Its origin lies in the fact that atomic precession is not limited to the field-free regions but also occurs in the interaction regions, where atomic coherence builds up or decreases continuously. Although the correction to the equivalent length of the drift region is small, it may be relevant to the evaluation of the second-order Doppler effect bias in primary cesium-beam standards to the extent that the atomic velocity is deduced from the lineshape and from the geometrical parameters of the cavity. It is shown that in current and projected standards with atoms of average thermal velocity, use of corrected dimensions may lead to a change of the calculated bias of the order of 10(exp -14), which is significant at the levels of accuracy considered nowadays.

Thomann, Pierre↗

A high-order compact numerical algorithm for supersonic flows

A dissipative compact two-four scheme (second-order time, fourth-order space) similar to the original MacCormack scheme has been developed, that exhibits greater accuracy than conventional fourth-order schemes. The dissipative nature of the scheme allows it to resolve weak discontinuities without artificial damping. A derivation of the scheme is presented, as well as the theoretical stability characteristics. The temporal scheme is then generalized into a steady-state formulation which achieves fourth-order spatial accuracy at steady-state. Several test problems are used to show that the scheme is more accurate than the traditional MacCormack scheme, and is nearly as efficient.

Carpenter, M. H.↗