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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Supercritical wing sections II: A handbook

The numerical aspect of theoretical work on transonics and supercritical wing sections are compiled. A model of the trailing edge is introduced which eliminates the loss of 15 to 20 percent experienced with heavily aft-loaded models, and it is indicated how drag creep can be reduced at off-design conditions. A rotated finite difference scheme is presented which can handle supersonic as well as subsonic free stream Mach numbers and leads to an effective three-dimensional program for the computation of transonic flow past an oblique wing. In the case of two-dimensional flow, the method is extended to take into account the displacement thickness computed by a semiempirical turbulent boundary layer correction. A series of supercritical wing sections is discussed together with comparisons between experimental and theoretical data. Computer programs and a brief manual for their operation are listed. It is shown that the programs furnish a physically adequate computer simulation of the compressible flows that arise in problems of transonic aerodynamics.

Bauer, F.↗
NASA NTRS1975

Viscous shock layer solutions for turbulent flow of radiating gas mixtures in chemical equilibrium

The viscous shock layer equations for hypersonic laminar and turbulent flows of radiating or nonradiating gas mixtures in chemical equilibrium are presented for two-dimensional and axially symmetric flow fields. Solutions are obtained using an implicit finite difference scheme and results are presented for hypersonic flow over spherically blunted cone configurations at free stream conditions representative of entry into the atmosphere of Venus. These data are compared with solutions obtained using other methods of analysis.

Anderson, E. C.↗
NASA NTRS1976

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 1: Computational technique

A numerical procedure has been developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second order accurate finite difference scheme is used to integrate the three dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗
NASA NTRS1976

Development of a computer code for calculating the steady super/hypersonic inviscid flow around real configurations. Volume 2: Code description

A numerical procedure was developed to compute the inviscid super/hypersonic flow field about complex vehicle geometries accurately and efficiently. A second-order accurate finite difference scheme is used to integrate the three-dimensional Euler equations in regions of continuous flow, while all shock waves are computed as discontinuities via the Rankine-Hugoniot jump conditions. Conformal mappings are used to develop a computational grid. The effects of blunt nose entropy layers are computed in detail. Real gas effects for equilibrium air are included using curve fits of Mollier charts. Typical calculated results for shuttle orbiter, hypersonic transport, and supersonic aircraft configurations are included to demonstrate the usefulness of this tool.

Marconi, F.↗
NASA NTRS1975

Computation of three-dimensional, inviscid supersonic flows

The paper sets forth in detail a method for the finite-difference computation of three-dimensional supersonic fields in an Eulerian mesh. First-, second-, and third-order finite difference schemes are examined. Attention is given to proper treatment of the impermeable and permeable boundaries encompassing the computational plane. Numerical results are presented for certain specific configurations: a conical wing-body combination, internal corner flow, a two-dimensional blunt body, an interfering shock problem, and three-dimensional inviscid supersonic flow past a shuttle-orbiter type vehicle.

Kutler, P.↗
NASA NTRS1976

An implicit finite-difference algorithm for hyperbolic systems in conservation-law form

An implicit finite-difference scheme is developed for the efficient numerical solution of nonlinear hyperbolic systems in conservation-law form. The algorithm is second-order time-accurate, noniterative, and in a spatially factored form. Second- or fourth-order central and second-order one-sided spatial differencing are accommodated within the solution of a block tridiagonal system of equations. Significant conceptual and computational simplifications are made for systems whose flux vectors are homogeneous functions (of degree one), e.g., the Eulerian gasdynamic equations. Conservative hybrid schemes, which switch from central to one-sided spatial differencing whenever the local characteristic speeds are of the same sign, are constructed to improve the resolution of weak solutions. Numerical solutions are presented for a nonlinear scalar model equation and the two-dimensional Eulerian gasdynamic equations.

Beam, R. M.↗
NASA NTRS1976

A three-dimensional, time-dependent model of Mobile Bay

A three-dimensional, time-variant mathematical model for momentum and mass transport in estuaries was developed and its solution implemented on a digital computer. The mathematical model is based on state and conservation equations applied to turbulent flow of a two-component, incompressible fluid having a free surface. Thus, bouyancy effects caused by density differences between the fresh and salt water, inertia from thare river and tidal currents, and differences in hydrostatic head are taken into account. The conservation equations, which are partial differential equations, are solved numerically by an explicit, one-step finite difference scheme and the solutions displayed numerically and graphically. To test the validity of the model, a specific estuary for which scaled model and experimental field data are available, Mobile Bay, was simulated. Comparisons of velocity, salinity and water level data show that the model is valid and a viable means of simulating the hydrodynamics and mass transport in non-idealized estuaries.

Pitts, F. H.↗
NASA NTRS1977

Numerical solutions of Navier-Stokes equations for the structure of a trailing vortex

The structure and decay of a trailing vortex were analyzed during the numerical solutions of the full Navier-Stokes equations. Unsteady forms of the governing equations were recast in terms of circulation, vorticity, and stream function as dependent variables, and a second upwind finite difference scheme was used to integrate them with prescribed initial and boundary conditions. The boundary conditions at the outer edge and at the outflow section of the trailing vortex were considered. Different models of the flow were postulated, and solutions were obtained describing the development of the flow as integration proceeds in time. A parametric study was undertaken with a view to understanding the various phenomena that may possibly occur in the trailing vortex. Using the Hoffman and Joubert law of circulation at the inflow section, the results of this investigation were compared with experimental data for a Convair 990 wind model and a rectangular wing. With an exponentially decaying law of circulation at the inflow section and an adverse pressure gradient at the outer edge of the trailing vortex, solutions depict vortex bursting through the sudden expansion of the core and/or through the stagnation and consequent reversal of the flow on the axis. It was found that this bursting takes place at lower values of the swirl ratio as the Reynolds number increases.

Jain, A. C.↗
NASA NTRS1977

Interactive calculation procedure for supersonic flows

An interactive procedure was developed for supersonic viscous flows that can be used for either two-dimensional or axisymmetric configurations. The procedure is directed to supersonic internal flows as well as those supersonic external flows that require consideration of mutual interaction between the outer flow and the boundary layer flow. The flow field is divided into two regions: an inner region which is highly viscous and mostly subsonic and an outer region where the flow is supersonic and in which viscous effects are small but not negligible. For the outer region a numerical solution is obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines. The inner region is treated by a system of equations of the boundary layer type that includes higher order effects such as longitudinal and transverse curvature and normal pressure gradients. These equations are coupled and solved simultaneously in the physical coordinates by using an implicit finite difference scheme. This system can also be used to calculate laminar and turbulent boundary layers using a scalar eddy viscosity concept.

Tassa, Y.↗
NASA NTRS1976

Steady supersonic flowfields with embedded subsonic regions

Supersonic flow past a blunt body is considered, where the flow contains an embedded subsonic region which lies between the shock wave and the body surface and is bounded by sonic lines from the body to the shock. A numerical approach is taken, which uses a basic finite difference scheme that solves the unsteady fluid dynamic equations in integral form. The unsteady equations are everywhere hyperbolic in time so no distinction need be made between subsonic and supersonic regions. Solutions to the mixed elliptic and hyperbolic steady flow equations are approached asymptotically in time. The method is illustrated for two-dimensional flows.

Maccormack, R. W.↗
NASA NTRS1976

Improved interactive calculation procedure for supersonic flows

An interactive numerical procedure has been developed for supersonic viscous flows (either two-dimensional or axisymmetric configurations). The flow field is divided into two regions: (1) an inner region which is highly viscous and mostly subsonic, and (2) an outer region where the flow is supersonic and in which viscous effects are small, but not negligible. This paper presents a detailed description of: I. Outer Region - numerical solution obtained by applying the method of characteristics to a system of equations which includes viscous and conduction transport terms only normal to the streamlines; II. Inner Region - treated by a system of equations of the boundary-layer type that includes higher order effects, such as longitudinal and transverse curvature and normal pressure gradients (equations are coupled and solved simultaneously in physical coordinates, using an implicit finite-difference scheme); III. Interactive Procedure - in the interaction mode, the two regions are coupled iteratively along a matching line, where the Mach number is of the order of 1.2.

Tassa, Y.↗
NASA NTRS1977

Prediction of outer layer mixing lengths in turbulent boundary layers

Predictions of mixing length by Pletcher's (1976) method (using a two-layer eddy viscosity model for a turbulent boundary layer at low Reynolds numbers) are compared to some values derived from turbulent boundary-layer profiles by other authors. The model is incorporated into a finite-difference scheme, to accurately predict low Reynolds number skin friction in supersonic flow.

Watson, R. D.↗
NASA NTRS1976

Numerical problems in the advection of pollutants

Prediction of pollutant transport involves solution of the advection equation numerically. Truncation errors arising from integrating the advection equation by using finite difference schemes are discussed. Numerical errors manifest themselves in artificial dissipation and dispersion of pollutants. Criteria for measuring the dissipation and dispersion are presented, and minimization of those errors is discussed. An alternative procedure by integration in the wavenumber space is examined. Examples illustrating quantitative discrepancy in actual forecasting are also shown.

Liu, C. Y.↗
NASA NTRS1976

Numerical solution of the steady-state Navier-Stokes equations for hypersonic flow about blunt axisymmetric bodies

The steady-state Navier-Stokes equations are solved for hypersonic flow about blunt axisymmetric bodies. The equations of motion are solved by successive approximations using an implicit finite-difference scheme. The results are compared with viscous shock-layer theory, experimental data, and time-dependent solutions of the Navier-Stokes equations. It is demonstrated that viscous shock-layer theory is sufficiently accurate for the range of flight conditions normally encountered by entry vehicles.

Anderson, E. C.↗
NASA NTRS1977

Calculation of unsteady two-dimensional laminar and turbulent boundary layers with fluctuations in external velocity

A numerical method is proposed for calculating unsteady two-dimensional laminar and turbulent boundary layers with allowance for fluctuations in external flow velocity. An eddy-viscosity formulation is used to model the Reynolds shear stress term appropriate to turbulent flow, and the governing boundary layer equations are solved by a two-point finite-difference scheme. For illustration, the method is applied to the calculation of the phase angles between the shear stress at the wall and an oscillating external laminar boundary layer on a flat plate. The effectiveness of the method (for engineering purposes) is demonstrated by comparing calculations for turbulent flows with experimental data.

Cebeci, T.↗
NASA NTRS1977

A solution to the Navier-Stokes equations based upon the Newton Kantorovich method

An implicit finite difference scheme based on the Newton-Kantorovich technique was developed for the numerical solution of the nonsteady, incompressible, two-dimensional Navier-Stokes equations in conservation-law form. The algorithm was second-order-time accurate, noniterative with regard to the nonlinear terms in the vorticity transport equation except at the earliest few time steps, and spatially factored. Numerical results were obtained with the technique for a circular cylinder at Reynolds number 15. Results indicate that the technique is in excellent agreement with other numerical techniques for all geometries and Reynolds numbers investigated, and indicates a potential for significant reduction in computation time over current iterative techniques.

Davis, J. E.↗
NASA NTRS1978

Solution of transonic flows by an integro-differential equation method

Solutions of steady transonic flow past a two-dimensional airfoil are obtained from a singular integro-differential equation which involves a tangential derivative of the perturbation velocity potential. Subcritical flows are solved by taking central differences everywhere. For supercritical flows with shocks, central differences are taken in subsonic flow regions and backward differences in supersonic flow regions. The method is applied to a nonlifting parabolic-arc airfoil and to a lifting NACA 0012 airfoil. Results compare favorably with those of finite-difference schemes.

Ogana, W.↗
NASA NTRS1977

Computational design of the basic dynamical processes of the UCLA general circulation model

The 12-layer UCLA general circulation model encompassing troposphere and stratosphere (and superjacent 'sponge layer') is described. Prognostic variables are: surface pressure, horizontal velocity, temperature, water vapor and ozone in each layer, planetary boundary layer (PBL) depth, temperature, moisture and momentum discontinuities at PBL top, ground temperature and water storage, and mass of snow on ground. Selection of space finite-difference schemes for homogeneous incompressible flow, with/without a free surface, nonlinear two-dimensional nondivergent flow, enstrophy conserving schemes, momentum advection schemes, vertical and horizontal difference schemes, and time differencing schemes are discussed.

Arakawa, A.↗