Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “second-order methods”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

On the accuracy of Whitham's method

The steady flow of an ideal gas past a conical body is studied by the method of matched asymptotic expansions and by Whitham's method in order to assess the accuracy of the latter. It is found that while Whitham's method does not yield a correct asymptotic representation of the perturbation field to second order in regions where the flow ahead of the Mach cone of the apex is disturbed, it does correctly predict the changes of the second-order perturbation quantities across a shock (the first-order shock strength). The results of the analysis are illustrated by a special case of a flat, rectangular plate at incidence.

Zahalak, G. I.↗
NASA NTRS1977

Polynomial interpolation methods for viscous flow calculations

Higher-order collocation procedures which result in block-tridiagonal matrix systems are derived from (1) Taylor series expansions and from (2) polynomial interpolation, and the relationships between the two formulations, called respectively Hermite and spline collocation, are investigated. A Hermite block-tridiagonal system for a nonuniform mesh is derived, and the Hermite approach is extended in order to develop a variable-mesh sixth-order block-tridiagonal procedure. It is shown that all results obtained by Hermite development can be recovered by appropriate spline polynomial interpolation. The additional boundary conditions required for these higher-order procedures are also given. Comparative solutions using second-order accurate finite difference and spline and Hermite formulations are presented for the boundary layer on a flat plate, boundary layers with uniform and variable mass transfer, and the viscous incompressible Navier-Stokes equations describing flow in a driven cavity.

Rubin, S. G.↗
NASA NTRS1978

Marching methods for elliptic problems. II, III

Higher-order operators for marching methods for elliptic equations are considered. Higher-order is understood in the sense of higher-order accuracy solutions to second-order Poisson equations, and in the sense of higher-order elliptic equations such as the biharmonic equation. The use of deferred corrections for overcoming stability problems is illustrated. Direct and iterative methods of extending the mesh size are considered. Multiple marching, patching, and influence extending techniques are described.

Roache, P. J.↗
NASA NTRS1976

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.↗
NASA NTRS1982

Solution of second-order linear system by matched asymptotic expansions

Matched asymptotic expansions (MAE) are used to obtain a first order approximation to the solution of a singularly perturbed second order system. A special case is considered in which the uniform asymptotic solution obtained by MAE is shown to converge to the exact solution. Ways in which the method can be used to sole higher-order linear systems, including those which are not singularly perturbed, are also discussed.

Ardema, M. D.↗
NASA NTRS1983

A preconditioned formulation of the Cauchy-Riemann equations

A preconditioning of the Cauchy-Riemann equations which results in a second-order system is described. This system is shown to have a unique solution if the boundary conditions are chosen carefully. This choice of boundary condition enables the solution of the first-order system to be retrieved. A numerical solution of the preconditioned equations is obtained by the multigrid method.

Phillips, T. N.↗
NASA NTRS1982

Second-order sensitivity derivatives in structural analysis

Haug's (1981) expressions for the calculation of static constraint function second derivatives have resulted in an order-of-magnitude saving over conventional methods for the calculation of such derivatives. A similar but simplified procedure is presented which results in a reduction of computational cost by an additional factor of two. The results are presented in a notation that is different from that of Haug, but more common in other work on structural synthesis and analysis.

Haftka, R. T.↗
NASA NTRS1983

Embedding methods for the steady Euler equations

An approach to the numerical solution of the steady Euler equations is to embed the first-order Euler system in a second-order system and then to recapture the original solution by imposing additional boundary conditions. Initial development of this approach and computational experimentation with it were previously based on heuristic physical reasoning. This has led to the construction of a relaxation procedure for the solution of two-dimensional steady flow problems. The theoretical justification for the embedding approach is addressed. It is proven that, with the appropriate choice of embedding operator and additional boundary conditions, the solution to the embedded system is exactly the one to the original Euler equations. Hence, solving the embedded version of the Euler equations will not produce extraneous solutions.

Chang, S. H.↗
NASA NTRS1984

Empirically derived second order schemes for the advection equations with minimum dissipation and dispersion errors

Two second-order advection schemes of the Lax-Wendroff type are empirically derived which have accuracy and phase characteristics similar to that of a third-order scheme. The new schemes are compared with other currently used methods, and are shown to have superior behavior in simulating the advection of localized disturbances. The schemes are derived for constant flow and generalized to two-dimensional non-uniform flow.

Takacs, L. L.↗
NASA NTRS1986

An embedding method for the steady Euler equations

Certain difficulties arise in connection with the numerical solution of a direct finite difference representation of the steady Euler equations. Johnson (1979, 1981, 1982) has, therefore, proposed a surrogate-equation technique, in which the first-order steady Euler equations are embedded in a certain second-order system of equations. The present paper is concerned with the theoretical justification for such an embedding approach. For the numerical solution of the two-dimensional steady Euler equations, it is shown that, under a continuity restriction, it is possible to solve a second-order embedded system together with appropriate additional boundary conditions. The result indicates that a more direct and potentially more efficient approach to the steady solutions exists than the alternative of solving the unsteady equations.

Chang, S.-H.↗
NASA NTRS1987

Development of upwind schemes for the Euler equations

Described are many algorithmic and computational aspects of upwind schemes and their second-order accurate formulations based on Total-Variation-Diminishing (TVD) approaches. An operational unification of the underlying first-order scheme is first presented encompassing Godunov's, Roe's, Osher's, and Split-Flux methods. For higher order versions, the preprocessing and postprocessing approaches to constructing TVD discretizations are considered. TVD formulations can be used to construct relaxation methods for unfactored implicit upwind schemes, which in turn can be exploited to construct space-marching procedures for even the unsteady Euler equations. A major part of the report describes time- and space-marching procedures for solving the Euler equations in 2-D, 3-D, Cartesian, and curvilinear coordinates. Along with many illustrative examples, several results of efficient computations on 3-D supersonic flows with subsonic pockets are presented.

Chakravarthy, Sukumar R.↗
NASA NTRS1985

Numerical solution of the unsteady Navier-Stokes equation

The construction and the analysis of nonoscillatory shock capturing methods for the approximation of hyperbolic conservation laws are discussed. These schemes share many desirable properties with total variation diminishing schemes, but TVD schemes have at most first-order accuracy, in the sense of truncation error, at extrema of the solution. In this paper a uniformly second-order approximation is constructed, which is nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time. This is achieved via a nonoscillatory piecewise linear reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell.

Osher, Stanley J.↗
NASA NTRS1987

Exposure time considerations in high temperature low cycle fatigue

The Conventional Strainrange Partitioning (CSRP) method for High-Temperature, Low-Cycle Fatigue (HTLCF) life prediction has its origins in the modeling of first-order, creep-fatigue waveform effects while treating as second-order effects, the influence of metallurgical or environmental time dependencies. Procedures are proposed to include the latter explicitly in the inelastic strainrange--life relations. For brevity, only the CP life relation will be presented in detail. The exposure-time effect within the CP inelastic strainrange (tensile creep reversed by compressive plasticity) was determined by tensile stresshold-time experiments for 316 SS at 816 C. Reductions in CP cyclic life of a factor of about two were observed with an increase in exposure time or a corresponding decrease in creep rate by a factor of about 100. The CP life relation has been modified to be expressed in terms of either Steady State Creep Rate (SSCR) or Exposure Time (ET). The applicability and accuracy of the time-dependent CP life relations is demonstrated by conducting verification experiments involving complex hysteresis loops. Metallographic examination revealed time-dependent degradation attributable to oxide formation and precipitation of carbides along grain boundaries.

Kalluri, S.↗
NASA NTRS1987

Study of instability and active control of boundary layers using Navier-Stokes equations

Spatial stability of a three-dimensional growing boundary layer is studied using Navier-Stokes equations. The finite difference method used is an explicit, predictor-corrector type which is fourth-order accurate in space and second-order in time. Active control of two-dimensional disturbances is simulated by localized surface heating. The computations show features of the flow relevant to early stages of transition such as lifting of vorticity away from the wall and localized regions of intense shear. Active control with single heating strip is only marginally effective, but may be improved by multiple strips.

Parikh, Paresh↗
NASA NTRS1987

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.↗
NASA NTRS1988

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.↗
NASA NTRS1988

Exposure time considerations in high temperature low cycle fatigue

The Conventional Strainrange Partitioning (CSRP) method for High-Temperature, Low-Cycle Fatigue (HTLCF) life prediction has its origins in the modeling of first-order, creep-fatigue waveform effects while treating as second-order effects, the influence of metallurgical or environmental time dependencies. Procedures are proposed to include the latter explicitly in the inelastic strainrange-life relations. For brevity, only the CP life relation will be presented in detail. The exposure-time effect within the CP inelastic strainrange (tensile creep reversed by compressive plasticity) was determined by tensile stresshold-time experiments for 316 SS at 816 C. Reductions in CP cyclic life of a factor of about two were observed with an increase in exposure time or a corresponding decrease in creep rate by a factor of about 100. The CP life relation has been modified to be expressed in terms of either Steady State Creep Rate (SSCR) or Exposure Time (ET). The applicability and accuracy of the time-dependent CP life relations is demonstrated by conducting verification experiments involving complex hysteresis loops. Metallographic examination revealed time-dependent degradation attributable to oxide formation and precipitation of carbides along grain boundaries.

Kalluri, S.↗
NASA NTRS1951

On the particular integrals of the Prandtl-Busemann iteration equations for the flow of a compressible fluid

The particular integrals of the second-order and third-order Prandtl-Busemann iteration equations for the flow of a compressible fluid are obtained by means of the method in which the complex conjugate variables are utilized as the independent variables of the analysis. The assumption is made that the Prandtl-Glauert solution of the linearized or first-order iteration equation for the two-dimensional flow of a compressible fluid is known. The forms of the particular integrals, derived for subsonic flow, are readily adapted to supersonic flows with only a change in sign of one of the parameters of the problem.

Kaplan, Carl↗