Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

An analysis of the relaxation of laminar boundary layer on a flat plate after passage of an interface with application to expansion-tube flows

The relaxation of the accelerating-gas boundary layer to the test-gas boundary layer over a flat plate in an expansion tube is analyzed. Several combinations of test gas and acceleration gas are considered. The problem is treated in two conically similar limits: (1) when the time lag between the arrival of the shock and the interface at the leading edge of the plate is very large, and (2) when this lag is negligible. The time-dependent laminar-boundary-layer equations of a binary mixture of perfect gases are taken as the flow-governing equations. This coupled set of differential equations, written in terms of the Lam-Crocco variables, has been solved by a line-relaxation finite-difference techniques. The results presented include the Stanton number and the local skin-friction coefficient as functions of shock Mach number and the nondimensional distance-time variable. The results indicate that more than 95 percent of the test-gas boundary layer exists over a length, measured from the leading edge of the plate, equal to about three-tenths of the distance traversed by the interface in the free stream.

Gupta, R. N.↗
NASA NTRS1972

The effect of gravity induced free convection upon the melting phenomena of a finite paraffin slab for thermal control

Two theoretical models were developed to predict the thermal response of the phase change material to a given hot plate temperature. A two-dimensional pure conduction model was developed to predict the melting of the phase change material when heat transfer was a function of conduction. A combined conduction-convection model, also two-dimensional, was developed to predict the phase change phenomena when heat transfer was a function of conduction and gravity-induced free convection. Both models were solved using explicit finite difference approximations on a digital computer. The experimental equipment consisted of a rectangular cell utilizing a heat chamber, an expansion chamber, and a test chamber; a sixteen channel multipoint recorder, and a fluid flow system. The recorder monitered hot and cold plate temperatures and interior node temperatures at two second intervals. A comparison of theoretical temperature profiles and experimental temperature profiles is presented for six runs at various angles of inclination of the test cell with respect to the horizontal direction. A detailed discussion of results is presented.

Bain, R. L.↗
NASA NTRS1972

Thermal modeling of phase change solidification in thermal control devices including natural convection effects

Natural convection effects in phase change thermal control devices were studied. A mathematical model was developed to evaluate natural convection effects in a phase change test cell undergoing solidification. Although natural convection effects are minimized in flight spacecraft, all phase change devices are ground tested. The mathematical approach to the problem was to first develop a transient two-dimensional conduction heat transfer model for the solidification of a normal paraffin of finite geometry. Next, a transient two-dimensional model was developed for the solidification of the same paraffin by a combined conduction-natural-convection heat transfer model. Throughout the study, n-hexadecane (n-C16H34) was used as the phase-change material in both the theoretical and the experimental work. The models were based on the transient two-dimensional finite difference solutions of the energy, continuity, and momentum equations.

Ukanwa, A. O.↗
NASA NTRS1973

Steady inviscid transonic flows over planar airfoils: A search for a simplified procedure

A finite difference procedure based upon a system of unsteady equations in proper conservation form with either exact or small disturbance steady terms is used to calculate the steady flows over several classes of airfoils. The airfoil condition is fulfilled on a slab whose upstream extremity is a semi-circle overlaying the airfoil leading edge circle. The limitations of the small disturbance equations are demonstrated in an extreme example of a blunt-nosed, aft-cambered airfoil. The necessity of using the equations in proper conservation form to capture the shock properly is stressed. Ability of the steady relaxation procedures to capture the shock is briefly examined.

Magnus, R.↗
NASA NTRS1973

Thermal correction of deformations in a telescope mirror

Orbiting astronomical observatories have the potential for making observations far superior to those from earth-based mirrors. In order for this performance to be realized, the contour of the primary mirror must be very accurately controlled. A preliminary investigation of the use of thermally induced elastic strains for correcting axisymmetric deformations in space telescope mirrors has been presented. The relation between axial deformation and thermal inputs was determined by a finite difference solution of the equations for thin elastic shells. The use of this technique was demonstrated analytically on a beryllium paraboloid. This mirror had 10 equally spaced thermal inputs and results are presented which show the nature of the temperature distribution required to correct deformations due to an acceleration-type loading.

Rhodes, M. D.↗
NASA NTRS1973

Preliminary study of a wall stabilized constructed arc

An iterative, implicit, finite-difference numerical technique is described which is suitable for obtaining solutions to the governing equations for a gas flowing in an axially symmetric constricted-arc heater. The method is shown to provide adequate solutions for three cases of simple pipe flows found in the literature, and for flow in a constricted-arc heater by direct comparison with experimental data. The comparison with arc-heater data includes static pressure, arc voltage, and wall heat flux, all as a function of axial location, and a radial temperature profile at one axial station. The arc-heater data were taken with air as the test gas at a heater inlet pressure approximately 0.40 atm and two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 584 amperes.

Graves, R. A.↗
NASA NTRS1973

Vertical distribution of vibrational energy of molecular nitrogen in a stable auroral red arc and its effect on ionospheric electron densities

Previous solutions of the problem of the distribution of vibrationally excited molecular nitrogen in the thermosphere have either assumed a Boltzmann distribution and considered diffusion as one of the loss processes or solved for the energy level populations and neglected diffusion. Both of the previous approaches are combined by solving the time dependent continuity equations, including the diffusion process, for the first six energy levels of molecular nitrogen for conditions in the thermosphere corresponding to a stable auroral red arc. The primary source of molecular nitrogen excitation was subexcitation, and inelastic collisions between thermal electrons and molecular nitrogen. The reaction rates for this process were calculated from published cross section calculations. The loss processes for vibrational energy were electron and atomic oxygen quenching and vibrational energy exchange. The coupled sets of nonlinear, partial differential equations were solved numerically by employing finite difference equations.

Newton, G. P.↗
NASA NTRS1973

Numerical calculations of flow fields

Numerical calculations were made of flow fields generated by various aerodynamic configurations. Data cover flow fields generated by a finitely thick lifting three dimensional wing with subsonic tips moving at supersonic speeds, cross flow instability associated with lifting delta wing configurations such as space shuttles, and flow fields produced by a lifting elliptic cone. Finite difference techniques were used to determine elliptic cone flow.

Anderson, D.↗
NASA NTRS1973

Unsteady turbulent boundary layer analysis

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of turbulent energy equations which include upstream history efforts, are solved numerically. An explicit finite difference analog to the partial differential equations is formulated and developed into a computer program. Calculations were made for a variety of unsteady flows in both two and three dimensions but primarily for two dimensional flow fields in order to first understand some of the fundamental physical aspects of unsteady turbulent boundary layers. Oscillating free stream flows without pressure gradient, oscillating retarded free stream flows and monotonically time-varying flows have all been studied for a wide frequency range. It was found that to the lowest frequency considered, the lower frequency bound being determined by economic considerations (machine time), there were significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗
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 NTRS1972

Design of the UCLA general circulation model

An edited version is reported of notes distributed at the Summer Workshop on the UCLA General Circulation Model in June 1971. It presents the computational schemes of the UCLA model, along with the mathematical and physical principles on which these schemes are based. Included are the finite difference schemes for the governing fluid-dynamical equations, designed to maintain the important integral constraints and dispersion characteristics of the motion. Also given are the principles of parameterization of cumulus convection by an ensemble of identical clouds. A model of the ground hydrology, involving the liquid, ice and snow states of water, is included. A short summary is given of the scheme for computing solar and infrared radiation transfers through clear and cloudy air.

Arakawa, A.↗
NASA NTRS1973

Analysis and testing of two-dimensional slot nozzle ejectors with variable area mixing sections

Finite difference computer techniques have been used to calculate the detailed performance of air to air two-dimensional ejectors with symmetric variable area mixing sections and coaxial coverging primary nozzles. The analysis of the primary nozzle assumed correct expansion of the flow and is suitable for subsonic and slightly supersonic velocity levels. The variation of the mixing section channel walls is assumed to be gradual so that the static pressure can be assumed uniform on planes perpendicular to the axis. A test program was run to provide two-dimensional ejector test data for verification of the computer analysis. A primary converging nozzle with a discharge geometry of 0.125 inch x 8.0 inch was supplied with 600 SCFM of air at about 35 psia and 180 F. This nozzle was combined with two mixing section geometries with throat sizes of 1.25 inch x 8.0 inch and 1.875 inch x 8.0 inch and was tested at a total of 11 operating points.

Gilbert, G. B.↗
NASA NTRS1972

Large deflection elastic-plastic dynamic response of stiffened shells of revolution

The formulation and check out porblems for a computer code DYNAPLAS, which analyzes the large deflection elastic-plastic dynamic response of stiffened shells of revolution, are presented. The formulation for special discretization is by the finite element method with finite differences being used for the evaluation of the pseudo forces due to material and geometric nonlinearities. Time integration is by the Houbolt method. The stiffeners may be due to concentrated or distributed eccentric rings and spring supports at arbitrary angles around the circumference of the elements. Check out porblems include the comparison of solutions from DYNAPLAS with experimental and other computer solutions for rings, conical and cylindrical shells and a curved panel. A hypothetical submarine including stiffeners and missile tube is studied under a combination of hydrostatic and dynamically applied asymmetrical pressure loadings.

Stricklin, J. A.↗
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

A numerical solution for thermoacoustic convection of fluids in low gravity

A finite difference numerical technique for solving the differential equations which describe thermal convection of compressible fluids in low gravity are reported. Results of one-dimensional calculations are presented, and comparisons are made to previous solutions. The primary result presented is a one-dimensional radial model of the Apollo 14 heat flow and convection demonstration flight experiment. The numerical calculations show that thermally induced convective motion in a confined fluid can have significant effects on heat transfer in a low gravity environment.

Spradley, L. W.↗
NASA NTRS1973

CONMIN: A FORTRAN program for constrained function minimization: User's manual

CONMIN is a FORTRAN program, in subroutine form, for the solution of linear or nonlinear constrained optimization problems. The basic optimization algorithm is the Method of Feasible Directions. The user must provide a main calling program and an external routine to evaluate the objective and constraint functions and to provide gradient information. If analytic gradients of the objective or constraint functions are not available, this information is calculated by finite difference. While the program is intended primarily for efficient solution of constrained problems, unconstrained function minimization problems may also be solved, and the conjugate direction method of Fletcher and Reeves is used for this purpose. This manual describes the use of CONMIN and defines all necessary parameters. Sufficient information is provided so that the program can be used without special knowledge of optimization techniques. Sample problems are included to help the user become familiar with CONMIN and to make the program operational.

Vanderplaats, G. N.↗
NASA NTRS1973

Some problems of the calculation of three-dimensional boundary layer flows on general configurations

An accurate solution of the three-dimensional boundary layer equations over general configurations such as those encountered in aircraft and space shuttle design requires a very efficient, fast, and accurate numerical method with suitable turbulence models for the Reynolds stresses. The efficiency, speed, and accuracy of a three-dimensional numerical method together with the turbulence models for the Reynolds stresses are examined. The numerical method is the implicit two-point finite difference approach (Box Method) developed by Keller and applied to the boundary layer equations by Keller and Cebeci. In addition, a study of some of the problems that may arise in the solution of these equations for three-dimensional boundary layer flows over general configurations.

Cebeci, T.↗