Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite difference scheme”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Second-order upwind differencing in a recirculating flow

The accuracy and stability of the second order upwind differencing scheme was investigated. The solution algorithm employed is based on a coupled solution of the nonlinear finite difference equations by the multigrid technique. Calculations have been made of the driven cavity flow for several Reynolds numbers and finite difference grids. In comparison with the hybrid differencing, the second order upwind differencing is somewhat more accurate but it is not monotonically accurate with mesh refinement. Also, the convergence of the solution algorithm deteriorates with the use of the second order upwind differencing.

Vanka, S. P.↗
NASA NTRS1993

Higher-order accurate Osher schemes with application to compressible boundary layer stability

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.

Van Der Vegt, J. J. W.↗
NASA NTRS1993

Higher-order accurate Osher schemes with application to compressible boundary layer stability

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.

Vandervegt, J. J. W.↗
NASA NTRS1992

Numerical simulation of entry flow over blunt swept-wing planes

Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Li, C. P.↗
NASA NTRS1994

On the Dynamics of Computing a Chemically Relaxed Nonequilibrium Flow

In order to gain insights into the strong dependence of numerical solutions on initial data for finite time steps, a set of nonlinear test problems rich enough to capture the behavior of difference schemes were recently identified and the numerical basins of attraction for these problems were computed using commonly used time discretizations in CFD. Our study revealed a wealth of surprisingly nonlinear behavior of numerical schemes that were not observed before, in particular for the implicit time discretizations that are commonly used in CFD. The goal of this work is to apply these tools to study a practical model from non-equilibrium flowfield relaxation. This type of problem arises in chemically nonequilibrium hypersonic flows such as in a shock tube experiment or an expanding nozzle. Here we consider a reacting mixture of (N2, N) for an inviscid one-dimensional steady model. Preliminary numerical results indicate that, aside from the possibility of spurious numerical solutions being introduced by the time discretizations, limitations on the model for physical or accurate solutions may also play a part in the dynamics observed.

Sweby, P. K.↗
NASA NTRS1995

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗
NASA NTRS1979

Fast, conservative schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) are investigated for the solution of matrix equations resulting from finite-difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, is introduced by an upwind bias of the density. Two implicit AF procedures are presented, and their convergence performance is compared with that of the standard transonic solution procedure: successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. Results indicate a substantial improvement in convergence rate for AF schemes relative to SLOR.

Holst, T. L.↗
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 NTRS1983

An improved finite-difference analysis of uncoupled vibrations of tapered cantilever beams

An improved finite difference procedure for determining the natural frequencies and mode shapes of tapered cantilever beams undergoing uncoupled vibrations is presented. Boundary conditions are derived in the form of simple recursive relations involving the second order central differences. Results obtained by using the conventional first order central differences and the present second order central differences are compared, and it is observed that the present second order scheme is more efficient than the conventional approach. An important advantage offered by the present approach is that the results converge to exact values rapidly, and thus the extrapolation of the results is not necessary. Consequently, the basic handicap with the classical finite difference method of solution that requires the Richardson's extrapolation procedure is eliminated. Furthermore, for the cases considered herein, the present approach produces consistent lower bound solutions.

Subrahmanyam, K. B.↗
NASA NTRS1987

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗
NASA NTRS1987

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗
NASA NTRS1988

An improved flux-split algorithm applied to hypersonic flows in chemical equilibrium

An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.

Palmer, Grant↗
NASA NTRS1975

Inverse solutions for laminar boundary-layer flows with separation and reattachment

Numerical solutions of the laminar, incompressible boundary layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified; the pressure is deduced from the solution. A vorticity-stream-function formulation of the boundary layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite difference equations, in which the flow direction is taken into account, are obtained by several global iteration schemes which are stable and have unconditional diagonal dominance. Results are compared with Klineberg and Steger's separated boundary layer calculations, and with Briley's solution of Navier-Stokes equations for a separated region. In addition, an approximate technique is presented in which the streamwise convection of vorticity is set equal to zero in the reversed flow region; such a technique results in a quick forward marching procedure for separated flows.

Carter, J. E.↗
NASA NTRS1977

Fast methods incorporating direct elliptic solvers for nonlinear applications in fluid dynamics

Semidirect methods are discussed, their present role, as well as some developments for their application in computational fluid dynamics. A semidirect method is a computational scheme that uses a fast, direct, elliptic solver as the driving algorithm for the iterative solution of finite difference equations. Specific subtopics include: (1) direct Cauchy Riemann solvers for first order elliptic equations; (2) application of the semidirect method to the mixed elliptic hyperbolic problem of steady, inviscid transonic flow; and (3) the treatment of interior conditions, such as those on an airfoil or wing, in semidirect methods.

Martin, E. D.↗
NASA NTRS1977

Investigations related to the inviscid-viscous interaction in transonic flows about finite 3-D wings

Inviscid small-disturbance theory has been shown to predict three-dimensional transonic flows about finite wings reasonably well as long as viscous effects are negligible. In order to include these effects, the inviscid small-disturbance solution of Bailey and Ballhaus (1975) has been combined with a finite-difference solution for Prandtl's boundary-layer equations. This solution employs the conditionally stable Krause (1968) scheme, implicit in the direction normal to the wall, to cope with the domain-of-dependence problem that arises for reverse cross flow. To be consistent with the inviscid-flow solution, the boundary layer is computed in the representative wing planform plane which is transformed into rectangular shape in the computational domain. The flow has been assumed turbulent, and a scalar eddy-viscosity model is adopted. The interaction between inviscid and viscous flow is modeled with the help of the displacement surface which is added to the geometric wing shape. Sample distributions of displacement thickness for swept wings are presented for weak and strong interaction cases.

Kordulla, W.↗
NASA NTRS1977

The separated turbulent boundary layer over a wavy wall

A study and application of the fourth order spline collocation procedure, numerical solution of boundary layer like differential equations, is presented. A simple inversion algorithm for the simultaneous solution of the resulting difference equations is given. Particular attention is focused on the boundary condition representation for the spline second derivative approximations. Solutions using the spline procedure, as well as the three point finite difference method, are presented for several model problems in order to assess and improve the spline numerical scheme. Application of the resulting algorithm to the incompressible laminar self similar boundary layer equations is presented.

Polak, A.↗
NASA NTRS1978

Numerical study of transonic flow over oscillating airfoils using the full potential equation

The behavior of unsteady aerodynamic loadings on airfoils oscillating in transonic flow has been investigated numerically with particular attention given to supercritical airfoil sections. A previously developed finite difference method, which is based on the full potential equation and which uses a quasi-conservative scheme for proper capture of a shock wave motion, was employed for the present study. The unsteady aerodynamic pressure and load distributions on several different airfoil sections are presented with particular emphasis on the effects of free-stream Mach number, reduced frequency, and mean angle of attack. These parameters are demonstrated to have a significant effect on the behavior of the unsteady aerodynamic loadings. Comparisons of the present calculations with the exact inviscid solution and with the experimental results are also presented.

Isogai, K.↗
NASA NTRS1984

An implict LU scheme for the Euler equations applied to arbitrary cascades

An implicit scheme for solving the Euler equations is derived and demonstrated. The alternating-direction implicit (ADI) technique is modified, using two implicit-operator factors corresponding to lower-block-diagonal (L) or upper-block-diagonal (U) algebraic systems which can be easily inverted. The resulting LU scheme is implemented in finite-volume mode and applied to 2D subsonic and transonic cascade flows with differing degrees of geometric complexity. The results are presented graphically and found to be in good agreement with those of other numerical and analytical approaches. The LU method is also 2.0-3.4 times faster than ADI, suggesting its value in calculating 3D problems.

Buratynski, E. K.↗