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 NTRS1972

Further analytical study of hybrid rocket combustion

Analytical studies of the transient and steady-state combustion processes in a hybrid rocket system are discussed. The particular system chosen consists of a gaseous oxidizer flowing within a tube of solid fuel, resulting in a heterogeneous combustion. Finite rate chemical kinetics with appropriate reaction mechanisms were incorporated in the model. A temperature dependent Arrhenius type fuel surface regression rate equation was chosen for the current study. The governing mathematical equations employed for the reacting gas phase and for the solid phase are the general, two-dimensional, time-dependent conservation equations in a cylindrical coordinate system. Keeping the simplifying assumptions to a minimum, these basic equations were programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and, the Crank-Nicolson scheme for the solid phase.

Hung, W. S. Y.↗
NASA NTRS1973

Three-dimensional compressible laminar boundary layers on sharp and blunt circular cones at angle of attack

A method for solving the three-dimensional compressible laminar boundary layer equations for the case of a circular cone and a sphere-cone body at an angle of attack is presented. The governing equations are modified by a similarity type transformation and then transformed into a Crocco-type form. The resulting set of equations is solved simultaneously by an iterative method using an implicit finite difference scheme by means of an efficient algorithm for equations of tridiagonal form. The effects of streamline swallowing on a sharp cone are included by introducing the true inviscid edge conditions at the distance from the wall equal to the boundary layer thickness. The validity of the approach was established by comparison of the computational results with similar results by other methods and with experimental data. It was concluded that at sufficiently high Mach number and moderate to large angles of attack, the streamline swallowing effects on a sharp cone result in higher values of skin friction and heat transfer as compared with the classical results for constant entropy.

Popinski, Z.↗
NASA NTRS1973

The use of Levy-Lees variables in three-dimensional boundary-layer flows

A method for solving a general class of three-dimensional boundary layer flows is developed. In the development, Levy-Lees variables are extended to three dimensions and equations are placed in these similarity variables. An implicit finite difference scheme which is stable for negative transverse velocities is used to solve these equations. The method developed is applied to obtain solutions for sharp and spherically blunted circular cones at angle of attack. Longitudinal and transverse distributions are presented for these cases. Good agreement is found with the results obtained by other numerical schemes and the experimental data of Tracy, for sharp circular cones at angle of attack. For spherically blunted cones at angle of attack, the results are in good agreement with axisymmetric sphere results up to the region where spherical symmetry holds.

Vatsa, V. N.↗
NASA NTRS1973

An analysis of a charring ablator with thermal nonequilibrium, chemical kinetics, and mass transfer

The differential equations governing the transient response of a one-dimensional ablative thermal protection system are presented for thermal nonequilibrium between the pyrolysis gases and the char layer and with finite rate chemical reactions occurring. The system consists of three layers (the char layer, the uncharred layer, and an optical insulation layer) with concentrated heat sinks at the back surface and between the second and third layers. The equations are solved numerically by using a modified implicit finite difference scheme to obtain solutions for the thickness of the charred and uncharred layers, surface recession and pyrolysis rates, solid temperatures, porosity profiles, and profiles of pyrolysis-gas temperature, pressure, composition, and flow rate. Good agreement is obtained between numerical results and exact solutions for a number of simplified cases. The complete numerical analysis is used to obtain solutions for an ablative system subjected to a constant heating environment. Effects of thermal, chemical, and mass transfer processes are shown.

Clark, R. K.↗
NASA NTRS1973

Impact of space shuttle orbiter reentry on mesospheric NOx.

Shuttle orbiter reentry will produce large amounts of nitric oxide (NO) through shock heating at mesospheric altitudes. The effects of this reentry are modeled by considering an impulsive source of NO in time-dependent equations including photochemistry and transport. The model uses an odd-nitrogen approach similar to Strobel's. Parameterized flux boundary conditions are imposed at the mesopause and stratopause and concentrations of the odd nitrogen species (NO, NO2, NO+, N(4S) and N(2D) are found numerically by a finite difference scheme. A localized disturbance can last many hours, depending on the effectiveness of eddy transport. Although the model is one-dimensional, effects of horizontal transport are also included approximately.

Stolarski, R. S.↗
NASA NTRS1974

The computation of steady nozzle flow by a time-dependent method

The steady flow in two-dimensional and axisymmetric nozzles was computed using a time-dependent method. In this method the interior mesh points were computed using the MacCormack finite-difference scheme, while a characteristic scheme was used to calculate the boundary mesh points. No explicit artificial viscosity term was included. The fluid was assumed to be a perfect gas. This method was used to compute the flow in a 45 deg - 15 deg conical, converging-diverging nozzle, a 15 deg conical, converging nozzle, and a 10 deg conical, plug nozzle. Good agreement between the numerical solution and experimental data was found. In contrast to previous time-dependent methods, the computational times were less than one minute on a CDC 6600 computer.

Cline, M. C.↗
NASA NTRS1974

Numerical calculations of two dimensional, unsteady transonic flows with circulation

The feasibility of obtaining two-dimensional, unsteady transonic aerodynamic data by numerically integrating the Euler equations is investigated. An explicit, third-order-accurate, noncentered, finite-difference scheme is used to compute unsteady flows about airfoils. Solutions for lifting and nonlifting airfoils are presented and compared with subsonic linear theory. The applicability and efficiency of the numerical indicial function method are outlined. Numerically computed subsonic and transonic oscillatory aerodynamic coefficients are presented and compared with those obtained from subsonic linear theory and transonic wind-tunnel data.

Beam, R. M.↗
NASA NTRS1974

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 NTRS1974

Analysis of dynamic stall using unsteady boundary-layer theory

The unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself. The predicted delay in flow-reversal onset, and its variation with pitch rate, are shown to be in reasonable agreement with experimental data relating to the delay in dynamic stall. From the comparisons it can be concluded (a) that the effects of time-dependence are sufficient to explain the failure of the boundary layer to separate during the dynamic overshoot, and (b) that there may be some link between forward movement of the reversal point and dynamic stall.

Scruggs, R. M.↗
NASA NTRS1974

Analysis of flow-reversal delay for a pitching airfoil

The coupled unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself.

Scruggs, R. M.↗
NASA NTRS1973

Computer simulation of shock wave boundary layer interaction

Development of a boundary layer correction for an analysis program previously used by Bauer et al. (1972) to study two-dimensional flow past transonic airfoils at off-design conditions where shocks which interact with the boundary layer appear. Using a method which combines conformal mapping with a finite-difference scheme due to Murman and Cole (1971), an analysis is made of the interaction between shock waves on a supercritical wing section and a turbulent boundary layer, assuming that the shocks are weak and that separation is insignificant. It is shown that relatively simple shock wave/boundary layer interactions of the kind envisioned in this study can be treated by determining the displacement thickness iteratively together with the flow and altering the profile accordingly. It is concluded that the proposed boundary layer correction appears to furnish an essential improvement in cases where both pronounced aft loading and a large enclosed supersonic region of flow are present.

Bauer, F.↗
NASA NTRS1974

Development of inverse inviscid transonic solution methods

A numerical method suitable for the analysis and/or design of supercritical transonic airfoils is reported. In order to achieve accuracy, the method utilizes the full inviscid potential flow equations; and in order to remain simple it solves the problem in a stretched Cartesian grid system. The resulting computer program has several advantages over others of its type -- its use in either the direct analysis mode in which the airfoil shape is prescribed and the flow field and surface pressures are determined, or in the inverse mode in which the surface pressures are given and the airfoil shape and flow field are computed. Other advantages of the program include its use in a design program, the rotated finite difference scheme and its determination of the airfoil shape simultaneously with the flow field relaxation solution.

Carlson, L. A.↗
NASA NTRS1975

Numerical solution of the hypersonic viscous-shock-layer equations for laminar, transitional, and turbulent flows of a perfect gas over blunt axially symmetric bodies

The viscous shock layer equations applicable to hypersonic laminar, transitional, and turbulent flows of a perfect gas over two-dimensional plane or axially symmetric blunt bodies are presented. The equations are solved by means of an implicit finite difference scheme, and the results are compared with a turbulent boundary layer analysis. The agreement between the two solution procedures is satisfactory for the region of flow where streamline swallowing effects are negligible. For the downstream regions, where streamline swallowing effects are present, the expected differences in the two solution procedures are evident.

Anderson, E. C.↗
NASA NTRS1975

Supercritical wing sections 2, volume 108

A mathematical theory for the design and analysis of supercritical wing sections was previously presented. Examples and computer programs showing how this method works were included. The work on transonics is presented in a more definitive form. For design, a better model of the trailing edge is introduced which should eliminate a loss of fifteen or twenty percent in lift experienced with previous heavily aft loaded models, which is attributed to boundary layer separation. How drag creep can be reduced at off-design conditions is indicated. A rotated finite difference scheme is presented that enables the application of Murman's method of analysis in more or less arbitrary curvilinear coordinate systems. This allows the use of supersonic as well as subsonic free stream Mach numbers and to capture shock waves as far back on an airfoil as desired. Moreover, it 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 semi-empirical turbulent boundary layer correction.

Bauer, F.↗
NASA NTRS1971

A computer program for two-dimensional and axisymmetric nonreacting perfect gas and equilibrium chemically reacting laminar, transitional and-or turbulent boundary layer flows

A computer program is described in detail for laminar, transitional, and/or turbulent boundary-layer flows of non-reacting (perfect gas) and reacting gas mixtures in chemical equilibrium. An implicit finite difference scheme was developed for both two dimensional and axisymmetric flows over bodies, and in rocket nozzles and hypervelocity wind tunnel nozzles. The program, program subroutines, variables, and input and output data are described. Also included is the output from a sample calculation of fully developed turbulent, perfect gas flow over a flat plate. Input data coding forms and a FORTRAN source listing of the program are included. A method is discussed for obtaining thermodynamic and transport property data which are required to perform boundary-layer calculations for reacting gases in chemical equilibrium.

Miner, E. W.↗
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 were obtained using an implicit finite-difference scheme and results are presented for hypersonic flow over spherically-blunted cone configurations at freestream 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 NTRS1975

Numerical simulation of turbulence in the presence of shear

The numerical calculations are presented of the large eddy structure of turbulent flows, by use of the averaged Navier-Stokes equations, where averages are taken over spatial regions small compared to the size of the computational grid. The subgrid components of motion are modeled by a local eddy-viscosity model. A new finite-difference scheme is proposed to represent the nonlinear average advective term which has fourth-order accuracy. This scheme exhibits several advantages over existing schemes with regard to the following: (1) the scheme is compact as it extends only one point away in each direction from the point to which it is applied; (2) it gives better resolution for high wave-number waves in the solution of Poisson equation, and (3) it reduces programming complexity and computation time. Examples worked out in detail are the decay of isotropic turbulence, homogeneous turbulent shear flow, and homogeneous turbulent shear flow with system rotation.

Shaanan, S.↗