Search NASA⌕ Search

Engineering topics

Wornom, S. F.

Publications and source records attributed to Wornom, S. F..

17 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Relaxation algorithms for the Euler equations

The alternating direction implicit central-difference scheme algorithms may be replaced by successive-line-relaxation (SLR) procedures that are stable in three dimensions. Several Beamand Warming-type codes are presently modified on the basis of SLR solution procedures in order to compute the flows over a cylinder, a NACA-0012 airfoil, and a shock wave reflecting off a plate. Since the codes are written in a 'delta form', the implementation of relaxation procedures that rely only on the corrections is relatively easy; the steady state residual computation and any artificial viscosity terms may be left unchanged.

Wornom, S. F.↗
NASA NTRS1984

Calculation of quasi-one-dimensional flows with shocks

The two-point subsonic method of Wornom (1983) is modified to permit calculation of transonic flows with shocks using artificial density and/or artificial pressure. This method, which requires no numerical boundary conditions, is then used to calculate transonic flow for large and small CFL (Courant-Friedrichs-Lewy) numbers. Convergence acceleration is achieved for small CFL numbers by using multi-grid.

Wornom, S. F.↗
NASA NTRS1984

A rule for selecting analytical boundary conditions for the conservative quasi-one-dimensional nozzle flow equations

For the one-dimensional nozzle flow equations, the number of analytical boundary conditions that may be applied corresponds to the number of exterior characteristics reaching the boundary points. Since the slopes of the characteristics are the characteristic speeds u, u+c, and u-c, one concludes that two analytical boundary conditions are required at a subsonic inflow boundary and one at a subsonic outflow boundary. However, using this guideline alone can lead to nonunique solutions in subsonic cases. This paper presents a rule for selecting analytical boundary conditions for the dependent variables which will yield unique solutions if they exist.

Wornom, S. F.↗
NASA NTRS1984

A two-point difference scheme for computing steady-state solutions to the conservative one-dimensional Euler equations

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central-difference scheme for the flux derivatives with dissipation added at supersonic points via the retarded density concept. Application of the method to 1-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. For 1-dimensional Euler calculations, it is shown that the scheme offers two advantages over the more widely-used three-point schemes. The first is in regard to application of boundary conditions, and the second relates to the fact that the two-point algorithm is well-conditioned for large time steps.

Wornom, S. F.↗
NASA NTRS1983

Implicit conservative characteristic modeling schemes for the Euler equations - A new approach

An implicit characteristic-modeling solution scheme for the Euler equations is presented. The scheme does not require the governing equations to be written in characteristic variables or the flux terms to be split into positive and negative contributions. For the two-dimensional problem of a shock wave reflecting from a flat plate, this feature and the simple solution algorithm combine to reduce the computational work per mesh point by 40 percent from that required by a standard, central-difference, implicit, solution algorithm. Application of the method to the quasi-one-dimensional nozzle flow equations for subsonic and supersonic flows without shocks shows the method to be well-conditioned for large time steps.

Wornom, S. F.↗
NASA NTRS1982

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite-difference method is presented for obtaining steady-state solutions to the time-dependent, conservative Euler equations for flows containing shocks. The method uses a two-point central difference scheme with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions show the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely-used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗
NASA NTRS1982

Application of two-point difference schemes to the conservative Euler equations for one-dimensional flows

An implicit finite difference method is presented for obtaining steady state solutions to the time dependent, conservative Euler equations for flows containing shocks. The method used the two-point differencing approach of Keller with dissipation added at supersonic points via the retarded density concept. Application of the method to the one-dimensional nozzle flow equations for various combinations of subsonic and supersonic boundary conditions shows the method to be very efficient. Residuals are typically reduced to machine zero in approximately 35 time steps for 50 mesh points. It is shown that the scheme offers certain advantages over the more widely used three-point schemes, especially in regard to application of boundary conditions.

Wornom, S. F.↗
NASA NTRS1978

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, which is the natural extension of the second order box scheme to fourth order, was demonstrated with application to the incompressible, laminar and turbulent, boundary layer equations. The efficiency of the present method is compared with two point and three point higher order methods, namely, the Keller box scheme with Richardson extrapolation, the method of deferred corrections, a three point spline method, and a modified finite element method. For equivalent accuracy, numerical results show the present method to be more efficient than higher order methods for both laminar and turbulent flows.

Wornom, S. F.↗
NASA NTRS1978

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

A fourth-order box method for solving the boundary layer equations

A fourth order box method for calculating high accuracy numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations is presented. 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. Numerical results for high accuracy test cases show the method to be significantly faster than other higher order and second order methods.

Wornom, S. F.↗
NASA NTRS1977

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

Forward marching procedure for separated boundary-layer flows

A forward-marching procedure for separated boundary-layer flows which permits the rapid and accurate solution of flows of limited extent is presented. The streamwise convection of vorticity in the reversed flow region is neglected, and this approximation is incorporated into a previously developed (Carter, 1974) inverse boundary-layer procedure. The equations are solved by the Crank-Nicolson finite-difference scheme in which column iteration is carried out at each streamwise station. Instabilities encountered in the column iterations are removed by introducing timelike terms in the finite-difference equations. This provides both unconditional diagonal dominance and a column iterative scheme, found to be stable using the von Neumann stability analysis.

Carter, J. E.↗
NASA NTRS1975

Solutions for incompressible separated boundary layers including viscous-inviscid interaction

Numerical solutions are presented for the laminar and turbulent boundary-layer equations for incompressible flows with separation and reattachment. The separation angularity is avoided by using an inverse technique in which the displacement thickness is prescribed and the pressure is deduced from the resulting solution. The turbulent results appear qualitatively correct despite the use of a two-layer eddy-viscosity model which is generally assumed appropriate only for mild-pressure-gradient flows. A new viscous-inviscid interaction technique is presented in which the inviscid flow is solved inversely by prescribing the pressure from the boundary-layer solution and deducing the new displacement thickness from the solution of a Cauchy integral. Calculations are presented using this interaction procedure for a laminar flow in which separation and reattachment occur on a solid surface.

Carter, J. E.↗
NASA NTRS1972

Longitudinal curvature and displacement speed effects on incompressible laminar boundary layers.

The title problem is considered for the case of flow past a circular cylinder placed normal to a uniform mainstream with Reynolds numbers from 40 to 200. Implicit finite difference numerical solutions are obtained for a set of boundary-layer equations that account for the second order effects associated with surface curvature and displacement speed. It was found that both of these contributors have a significant influence on the internal structure of the viscous region and that an accurate estimate of the surface pressure distribution is essential for estimating the surface shear stress.

Werle, M. J.↗
NASA NTRS1972

Displacement interaction and surface curvature effects on hypersonic boundary layers.

The title problem was studied employing implicit finite-difference methods to obtain numerical solutions to a composite set of boundary-layer equations valid to second order. Results are given for flow up a two-dimensional cubic compression ramp for free-stream Mach numbers of 6, 8, and 12.25 and for free-stream Re/inch equal to 85,800 and 25,800 at a wall-to-stagnation temperature ratio of 0.223. Comparisons with independent theories and experimental results are given. Nonsingular separation was produced at a free-stream Mach number of 12.25. For all cases considered, displacement and curvature effects canceled one another when a consistent treatment of inviscid and viscous curvature corrections was employed - the second-order theory virtually reproducing the first-order results.

Wornom, S. F.↗
NASA NTRS1971

A Numerical Study of Displacement Body and Curvature Effects on Incompressible and Compressible Laminar Boundary Layers

This technique has been applied to study such effects on incompressible flow around cylinders at moderate to low Reynolds numbers and for compression ramps at hypersonic Mach numbers by employing a finite difference method to obtain numerical solutions. The results indicate that the technique can be applied successfully in both regimes and does predict the correct trend in regions of large curvature and displacement body effects. It was concluded that curvature corrections should only be attempted in cases where all displacement effects can be fully accounted for.

Wornom, S. F.↗