Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - BOUNDARY LAYER - PRANDTL”

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 109 records · Page 6

Development and Verification of the Charring, Ablating Thermal Protection Implicit System Simulator

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver (CATPISS) is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method (FEM) with first and second order fully implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton s method, while the linear system is solved via the Generalized Minimum Residual method (GMRES). Verification results from exact solutions and Method of Manufactured Solutions (MMS) are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗

Development and Verification of the Charring Ablating Thermal Protection Implicit System Solver

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method with first and second order implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton's method, while the fully implicit linear system is solved with the Generalized Minimal Residual method. Verification results from exact solutions and the Method of Manufactured Solutions are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗

A study of the possible characteristics of a low-altitude electron layer in the Martian atmosphere

The apparent diurnal Martian surface pressure variations, as deduced from radio occultation experiments, is discussed and explained as possibly arising from the effect of a low-altitude electron layer. Possible source and loss mechanisms for the low altitude electron layer are presented and discussed. Time dependent differential equations describing the electron layer are derived, and then integrated to investigate the electron distribution resulting from several processes that might occur in the atmosphere. It is concluded that the source mechanism is the sublimation of alkali atoms from a permanent dust layer, and that the dominant loss process must involve CO2 clustering about the alkali atoms. An electron layer is developed which explains the apparent diurnal surface pressure variation.

Wallio, H. A.↗

Inertial Taylor columns on a beta plane

An investigation is conducted concerning the lowest-order effect of variable Coriolis parameters on Taylor-column formation, taking into account the flow over a bump on a beta plane. The model considered involves a two-layer fluid on a beta plane. The governing equations for the two layers are discussed along with the characteristics of the general solutions, solutions for two retrograde currents, solutions for two prograde currents, experiments related to the single-layer problem, and the geophysical implications of the results of the investigation.

Mccartney, M. S.↗

LANDSAT-D investigations in snow hydrology

Two stream methods provide rapid approximate calculations of radiative transfer in scattering and absorbing media. Although they provide information on fluxes only, and not on intensities, their speed makes them attractive to more precise methods. The methods provide a comprehensive, unified review for a homogeneous layer, and solve the equations for reflectance and transmittance for a homogeneous layer over a non reflecting surface. Any of the basic kernels for a single layer can be extended to a vertically inhomogeneous medium over a surface whose reflectance properties vary with illumination angle, as long as the medium can be subdivided into homogeneous layers.

Dozier, J.↗

Application of program LAURA to three-dimensional AOTV flowfields

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is a robust, finite volume, single-level storage, implicit upwind differencing algorithm which has been documented and tested on several three-dimensional blunt-body flows. The algorithm can run at unlimited Courant numbers (relaxing the steady-state equations) but requires the inversion of only a 5 x 5 matrix per computational cell. An alternating directional sweep Gauss-Seidel substitution strategy is used to relax the governing equations. At present, the Euler and thin-layer Navier-Stokes equations using Sutherland's law for viscosity have been modeled for a perfect gas, equilibrium air, and nonequilibrium air chemistry neglecting diffusion. The equilibrium and nonequilibrium air chemistry options have been described in a companion paper. Good comparisons with experimental data and another calculation method for pressure distributions, aerodynamic coefficients, and heat-transfer distributions have been demonstrated for three-dimensional blunt-body flows.

Gnoffo, P. A.↗

Refractive Index Effects on Radiation in an Absorbing, Emitting, and Scattering Laminated Layer

A simple set of equations is derived for predicting temperature radiative energy flow in a two-region semitransparent laminated layer in the limit of zero heat conduction. The composite is heated on its two sides by unequal amounts of incident radiation. The two layers of the composite have different refractive indices, and each material absorbs, emits, and isotropically scatters radiation. The interfaces are diffuse, and all interface reflections are included. To illustrate the thermal behavior that is readily calculated from the equations, typical results an given for various optical thicknesses and refractive indices of the layers. Internal reflections have a substantial effect on the temperature distribution and radiative heat flow.

Siegel, R.↗

Hypersonic ionizing air viscous shock-layer flows over nonanalytic blunt bodies

The equations which govern the viscous shock-layer flow are presented and the method by which the equations are solved is discussed. The predictions of the present finite-difference method are compared with other numerical predictions as well as with experimental data. The principal emphasis is placed on predictions of the viscous flowfield for the windward plane of symmetry of the space shuttle orbiter and other axisymmetric bodies which approximate the shuttle orbiter geometry. Experimental data on two slender sphere-cones at hypersonic conditions are also considered. The present predictions agreed well with experimental data and with the past predictions. Substantial differences were found between present predictions and more approximate methods.

Miner, E. W.↗

Recent advances in Runge-Kutta schemes for solving 3-D Navier-Stokes equations

A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.

Vatsa, Veer N.↗

The Proteus Navier-Stokes code

An effort is currently underway at NASA Lewis to develop two and three dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. Proteus solves the Reynolds-averaged, unsteady, compressible Navier-Stokes equations in strong conservation law form. Turbulence is modeled using a Baldwin-Lomax based algebraic eddy viscosity model. In addition, options are available to solve thin layer or Euler equations, and to eliminate the energy equation by assuming constant stagnation enthalpy. An extensive series of validation cases have been run, primarily using the two dimensional planar/axisymmetric version of the code. Several flows were computed that have exact solution such as: fully developed channel and pipe flow; Couette flow with and without pressure gradients; unsteady Couette flow formation; flow near a suddenly accelerated flat plate; flow between concentric rotating cylinders; and flow near a rotating disk. The two dimensional version of the Proteus code has been released, and the three dimensional code is scheduled for release in late 1991.

Towne, Charles E.↗

Numerical studies of major and minor stratospheric warmings caused by orographic forcing

A primitive equation spectral model using spherical harmonics is formulated to study dynamic interactions between the troposphere and stratosphere in association with sudden stratospheric warmings. Using sigma coordinates for five tropospheric layers and log-pressure coordinates for 26 stratospheric and mesospheric layers, separate model equations for each system are combined to form single matrix governing equations. The gradual introduction of large scale topography to balanced initial states representative of observed mean winter conditions in the Northern Hemisphere is used for the generation of planetary waves during 40-day time integrations. Results of these integrations indicate that stratospheric warmings can be simulated by this orographic forcing and that mean momentum flux divergence due to zonal mean motion appears to be an essential mechanism of these simulated sudden warmings. It was found that the strength of the polar night jet can be a determining factor whether a warming becomes 'major' or 'minor'.

Koermer, J. P.↗

Solutions for Reacting and Nonreacting Viscous Shock Layers with Multicomponent Diffusion and Mass Injection

Numerical solutions are presented for the viscous shocklayer equations where the chemistry is treated as being either frozen, equilibrium, or nonequilibrium. Also the effects of the diffusion model, surface catalyticity, and mass injection on surface transport and flow parameters are considered. The equilibrium calculations for air species using multicomponent: diffusion provide solutions previously unavailable. The viscous shock-layer equations are solved by using an implicit finite-difference scheme. The flow is treated as a mixture of inert and thermally perfect species. Also the flow is assumed to be in vibrational equilibrium. All calculations are for a 45 deg hyperboloid. The flight conditions are those for various altitudes and velocities in the earth's atmosphere. Data are presented showing the effects of the chemical models; diffusion models; surface catalyticity; and mass injection of air, water, and ablation products on heat transfer; skin friction; shock stand-off distance; wall pressure distribution; and tangential velocity, temperature, and species profiles.

Moss, J. N.↗

Sensitivity of the Tropical Pacific Ocean to Precipitation Induced Freshwater Flux

We have performed a series of experiments using an ocean model to study the sensitivity of tropical Pacific Ocean to variations in precipitation induced freshwater fluxes. Variations in these fluxes arise from natural causes on all time scales. In addition, estimates of these fluxes are uncertain because of differences among measurement techniques. The model used is a quasi-isopycnal model, covering the Pacific from 40 S to 40 N. The surface forcing is constructed from observed wind stress, evaporation, precipitation, and surface temperature (SST) fields. The heat flux is produced with an iterative technique so as to maintain the model close to the observed climatology, but with only a weak damping to that climatology. Climatological estimates of evaporation are combined with various estimates of precipitation to determine the net surface freshwater flux. Results indicate that increased freshwater input decreases salinity as expected, but increases temperatures in the upper ocean. Using the freshwater flux estimated from the Microwave Sounding Unit leads to a warming of up to 0.6 C in the western Pacific over a case with zero net freshwater flux. SST is sensitive to the discrepancies among different precipitation observations, with root-mean-square differences in SST on the order of 0.2-0.3 C. The change in SST is more pronounced in the eastern Pacific, with differences of over 1 C found among the various precipitation products. Interannual variation in precipitation during El Nino events leads to increased warming. During the winter of 1982-83, freshwater flux accounts for about 0.4 C (approximately 10-15% of the maximum warming) of the surface warming in the central-eastern Pacific. Thus, the error of SST caused by the discrepancies in precipitation products is more than half of the SST anomaly produced by the interannual variability of observed precipitation. Further experiments, in which freshwater flux anomalies are imposed in the western, central, and eastern Pacific, show that the influence of net freshwater flux is also spatially dependent. The imposition of freshwater flux in the far western Pacific leads to a trapping of salinity anomaly to the surface layers near the equator. An identical flux imposed in the central Pacific produces deeper and off-equatorial salinity anomalies. The contrast between these two simulations is consistent with other simulations of the western Pacific barrier layer information.

Yang, Song↗

The characteristics of a possible low altitude electron layer in the Martian atmosphere

The apparent diurnal Martian surface pressure variation, as deduced from radio occultation experiments, is discussed and explained as possibly arising from the effect of a low altitude electron layer. Possible source and loss mechanisms for the low altitude electron layer are presented and discussed. Time-dependent differential equations describing the electron layer are derived and then integrated to investigate the electron distribution resulting from the several processes that might occur in the atmosphere. It is concluded that the source mechanism is the sublimation of alkali atoms from a permanent dust layer (a dust layer of 0.2 micron particles of density 9/cu cm is sufficient), and that the dominant loss process must involve CO2 clustering to the alkali atoms. Using these processes, an electron layer is developed which would explain the apparent diurnal surface pressure.

Wallio, H. A.↗

Non-Boltzmann Modeling for Air Shock-Layer Radiation at Lunar-Return Conditions

This paper investigates the non-Boltzmann modeling of the radiating atomic and molecular electronic states present in lunar-return shock-layers. The Master Equation is derived for a general atom or molecule while accounting for a variety of excitation and de-excitation mechanisms. A new set of electronic-impact excitation rates is compiled for N, O, and N2+, which are the main radiating species for most lunar-return shock-layers. Based on these new rates, a novel approach of curve-fitting the non-Boltzmann populations of the radiating atomic and molecular states is developed. This new approach provides a simple and accurate method for calculating the atomic and molecular non-Boltzmann populations while avoiding the matrix inversion procedure required for the detailed solution of the Master Equation. The radiative flux values predicted by the present detailed non-Boltzmann model and the approximate curve-fitting approach are shown to agree within 5% for the Fire 1634 s case.

Johnston, Christopher O.↗

Onset of finger convection in a horizontal porous layer underlying a fluid layer

The problem of the onset of finger convection in a porous layer underlying a fluid layer is considered using linear stability analysis. The linear stability equations for the porous layer are formulated for temperature and salinity gradients existing in both layers. The eigenvalue problem is solved by a shooting method. The solution method and associated computer program are validated by comparison with the results of Sun (1973) for the thermal convection case. Results are also presented for the onset of salt-finger convection.

Chen, F.↗

Diagonal implicit multigrid solution of compressible turbulent flows

The Caughey (1988) multigrid diagonal implicit algorithm for Euler equations is extended to solve the two-dimensional thin layer Navier-Stokes equations for turbulent transonic flows. To further accelerate convergence to a steady state, the implicit scheme is used within the multigrid method. Results are presented for transonic flows past airfoils, and the flow-field results are compared with other computational data and experiments demonstrating the accuracy of the method.

Varma, R. R.↗

Computational fluid dynamics application to hypersonic flow over a Martian entry vehicle (MEV) - A correlation with experiment

A class of vehicles for a mission to Mars are analyzed for aerodynamic characteristics using advanced Computational Fluid Dynamics (CFD). The general configuration is a modified cone-conical-frustum geometry where the nose radius has a large influence on the flowfield. Inviscid-compressible flow using the Euler equations and viscous-compressible flow using the thin-layer Navier-Stokes equations is applied to the configuration. The surface modeling, grid generation and application of state-of-the-art CFD software are described. The effects of nose radius, angle of attack, and hypersonic velocity on the flight characteristics of the vehicle are discussed. The numerical simulations demonstrate the merits of the inviscid and viscous software. Results are compared with wind tunnel experiments.

Stewart, John E.↗