Search NASA⌕ Search

SEARCH · Search NASA

Results for “NUMERICAL ANALYSIS”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Numerical analysis of natural convection in two-dimensional square and circular containers in low gravity

A numerical study of natural convection in circular cylinder and square enclosures shows that the analytic low Rayleigh number theory of previous investigators is valid for Rayleigh numbers up to 1000. For a Rayleigh number of 5000, steady state values of maximum fluid velocity differ by 20 percent. This deviation between analytic theory and numerical results increases for higher Rayleigh numbers. In addition, the low Rayleigh number theory is shown to be valid for higher Rayleigh numbers for a portion of the transient phase before significant deviation becomes apparent. It is also shown that square shaped experimental configurations may be analytically approximated with good accuracy by circular cylinders of equal cross sectional area for the prediction of convection velocities and flow patterns at low Rayleigh number.

Robertson, S. J.↗
NASA NTRS2003

Numerical Analysis of the SCHOLAR Supersonic Combustor

The SCHOLAR scramjet experiment is the subject of an ongoing numerical investigation. The facility nozzle and combustor were solved separate and sequentially, with the exit conditions of the former used as inlet conditions for the latter. A baseline configuration for the numerical model was compared with the available experimental data. It was found that ignition-delay was underpredicted and fuel-plume penetration overpredicted, while the pressure rise was close to experimental values. In addition, grid-convergence by means of grid-sequencing could not be established. The effects of the different turbulence parameters were quantified. It was found that it was not possible to simultaneously predict the three main parameters of this flow: pressure-rise, ignition-delay, and fuel-plume penetration.

Rodriguez, Carlos G.↗
NASA NTRS2011

Numerical Analysis of a Pulse Detonation Cross Flow Heat Load Experiment

A comparison between experimentally measured and numerically simulated, time-averaged, point heat transfer rates in a pulse detonation (PDE) engine is presented. The comparison includes measurements and calculations for heat transfer to a cylinder in crossflow and to the tube wall itself using a novel spool design. Measurements are obtained at several locations and under several operating conditions. The measured and computed results are shown to be in substantial agreement, thereby validating the modeling approach. The model, which is based in computational fluid dynamics (CFD) is then used to interpret the results. A preheating of the incoming fuel charge is predicted, which results in increased volumetric flow and subsequent overfilling. The effect is validated with additional measurements.

Paxson, Daniel E.↗
NASA NTRS1984

Numerical analysis of flow and heat transfer in the VAFB LOX storage Dewar tank

The present report describes numerical simulation of three-dimensional transient distributions of velocity and temperature of liquid oxygen (LOX) in the LOX Dewar tank of Vendenberg Air Force Base (VAFB). The present analyses cover the replenish time period only. Four test cases have been considered. For all four cases, the input boundary conditions are comprised of LOX facility heat loads, drain flow rates, recirculation flow rates and dewar heating. All the quantities are prescribed as functions of time. The first two test cases considered sensitivity of results to the computational grid. In Case 3, system heat load was changed, while in Case 4, a lower LOX level was specified. Cases 1 and 2 showed that the temperatures were not sensitive to the grid refinement. This provided a basic check on the numerical model. Cases 3 and 4 showed that the thermal boundary layer motion near the tank surface becomes more significant at the late time, e.g., 5 1/2 hours from replenish start. Comparison between results of Cases 3 and 4 showed, as expected, that the smaller initial LOX volume given in Case 4, results in higher temperature level. All calculated velocity and temperature distributions were found to be plausible.

Tam, L. T.↗
NASA NTRS1991

Numerical analysis of a nuclear fuel element for nuclear thermal propulsion

A computational fluid dynamics model with porosity and permeability formulations in the transport equations has been developed to study the concept of nuclear thermal propulsion through the analysis of a pulsed irradiation of a particle bed element (PIPE). The numerical model is a time-accurate pressure-based formulation. An adaptive upwind scheme is employed for spatial discretization. The upwind scheme is based on second- and fourth-order central differencing with adaptive artificial dissipation. Multiblocked porosity regions have been formulated to model the cold frit, particle bed, and hot frit. Multiblocked permeability regions have been formulated to describe the flow shaping effect from the thickness-varying cold frit. Computational results for several zero-power density PIPEs and an elevated-particle-temperature PIPE are presented. The implications of the computational results are discussed.

Wang, Ten-See↗
NASA NTRS1986

Numerical analysis of flow in the hot gas manifold of the Space Shuttle main engine

This paper presents a numerical model and results of analyses carried out to characterize the flow through the two duct configuration of the Hot Gas Manifold of the Space Shuttle Main Engine. Three dimensional computations have been performed for a half-scale air test model using a nonorthogonal body-fitted coordinate system. The calculation domain is extended from the inlet of the turnaround duct to the exit of the transfer duct. Three test cases, one for laminar flow and two for turbulent flow, have been considered. For turbulent flows, constant eddy viscosity and the k-epsilon model of turbulence were employed. As expected, laminar flow calculation predicts much larger pressure drop than turbulent flow cases. The turbulent flow results are in good agreement with the available flow-visualization data. This study and experimental data indicate that the two-transfer duct design will significantly improve the flow distribution in the Hot Gas Manifold and thereby enhance the overall performance of the SSME.

Owens, S. F.↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗
NASA NTRS1986

Numerical analysis of confined recirculating gas-solid turbulent flows

Particle dispersion in confined recirculating turbulent flows has been investigated numerically. The present computational scheme utilizes Eulerian forms of the governing equations and allows two-way coupling between gas and solid phases. A recently developed two-phase closure model based on the multiple-scale turbulence model was used for the estimation of turbulent viscosities and diffusivities. For the particle size and loading considered in this study, the effect of particulate phase on the conveying gas is small, the nonequilibrium between the two phases is significant. Effects of recirculation, expanded chamber size and secondary annular jet momentum on the particle mixing rate are also investigated. In general, the present numerical results are in reasonably good agreement with the available experimental data.

Chen, C. P.↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.

Bune, Andris V.↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the perioidic boundary conditions. Approximate nonreflecting boundary conditons were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distribution for certain flow conditions.

Huff, Dennis L.↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and non-zero inter-blade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions have been coded for the inlet and exit boundary conditions. Sample unsteady solutions have been performed for an oscillating cascade and compared to experimental data. Also, test cases were fun for a flat plate cascade to compare with an unsteady, small-perturbation, subsonic analysis. The predictions for oscillating cascades with non-zero inter-blade phase angles are in good agreement with experimental data and small-perturbation theory. The zero degree inter-blade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus non-reflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions. This code is expected to be used as a tool for reviewing simpler models that do not include the full nonlinear aerodynamics or as a final check for designs against flutter in turbomachinery.

Huff, Dennis L.↗
NASA NTRS1989

Numerical analysis of flow through oscillating cascade sections

The design of turbomachinery blades requires the prevention of flutter for all operating conditions. However, flow field predictions used for aeroelastic analysis are not well understood for all flow regimes. The present research focuses on numerical solutions of the Euler and Navier-Stokes equations using an ADI procedure to model two-dimensional, transonic flow through oscillating cascades. The model prescribes harmonic pitching motions for the blade sections for both zero and nonzero interblade phase angles. The code introduces the use of a deforming grid technique for convenient specification of the periodic boundary conditions. Approximate nonreflecting boundary conditions were coded for the inlet and exit boundary conditions. Sample unsteady solutions were performed for an oscillating cascade and compared to experimental data. Also, test cases were run for a flat plate cascade to compare with the unsteady, small-perturbation, subsonic analyis. The predictions for oscillating cascades with nonzero interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. The zero degree interblade phase angle cases, which were near a resonant condition, differ from the experiment and theory. Studies on reflecting versus nonreflecting inlet and exit boundary conditions show that the treatment of the boundary can have a significant effect on the first harmonic, unsteady pressure distributions for certain flow conditions.

Huff, Dennis L.↗
NASA NTRS2019

Numerical analysis of unsteady flow in a converging-diverging nozzle

This study investigates the physics of jet noise by the numerical solution of the three dimensional unsteady Navier-Stokes equations for a low-supersonic jet (Mach numbers on the order of 1.4) from a converging-diverging (C-D) nozzle and mixing with a subsonic surrounding flow. The numerical method used for the solution of these equations is the second order accurate MacCormack explicit predictor-corrector scheme. The numerical results are compared to the experimental results of Yamamoto for the same C-D nozzle geometry and similar operating conditions. The goal of the study is to use numerical simulation techniques to resolve accurately physical processes which contribute to jet noise, and to assess whether performing these simulations in a three dimensional geometry is worth the additional cost in computing resources as compared to a comparable simulation on a two dimensional planar or axisymmetric geometry. Scientific visualization is used to observe the spreading characteristics of the jet, as well as to examine the structure of the shocks within the jet core. The time average velocity and turbulence intensity level at various locations are computed and compared with experimental results. A two point space-time cross correlation is used to compute the mean convection speed through the shear layer, and these predictions are compared with results from theory.

James N Scott↗
NASA NTRS1991

Leading-edge receptivity to a vortical freestream disturbance: A numerical analysis

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is investigated numerically. The flow is simulated by solving the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. A finite-difference scheme which is second-order accurate in both space and time is used. As a first step, the steady basic-state solution is computed. Then a small amplitude vortical disturbance is introduced at the upstream boundary and the governing equations are solved time-accurately to evaluate the spatial and temporal growth of the perturbations leading to instability waves (Tollmien-Schlichting waves) inside the boundary layer. Preliminary results for a symmetric, 2-D disturbance reveal the presence of Tollmien-Schlichting waves aft of the flat-plate/ellipse juncture.

Buter, Thomas A.↗
NASA NTRS2017

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This presentation describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL↗
NASA NTRS2017

Numerical Analysis of Mixed-Phase Icing Cloud Simulations in the NASA Propulsion Systems Laboratory

This paper describes the development of a numerical model that couples the thermal interaction between ice particles, water droplets, and the flowing gas of an icing wind tunnel for simulation of NASA Glenn Research Centers Propulsion Systems Laboratory (PSL). The ultimate goal of the model is to better understand the complex interactions between the test parameters and have greater confidence in the conditions at the test section of the PSL tunnel. The model attempts to explain the observed changes in test conditions by coupling the conservation of mass and energy equations for both the cloud particles and flowing gas mass. Model predictions were compared to measurements taken during May 2015 testing at PSL, where test conditions varied gas temperature, pressure, velocity and humidity levels, as well as the cloud total water content, particle initial temperature, and particle size distribution.

NASA PSL↗
NASA NTRSDate not supplied

Numerical Analysis of Cavitation-Induced Two-Phase Fluid Transients in Cryogenic Feedlines Using Finite Volume Procedure

This article presents a numerical investigation of the cavitation induced fluid transient in cryogenic feedline due to sudden closure of the valve A one dimensional numerical model based on finite volume procedure consisting of nodes and branches using Generalized Fluid System Simulation Program ( has been employed for the prediction of pressure variations The unsteady friction formulation accurately determines wall shear stress, and a homogeneous mixture model is used for two phase flow modeling Numerical results obtained with the present formulation are compared with experimental results and Method of Characteristic ( based results It is observed that low residual in flow variable calculation improves the solution providing numerically stable results GFSSP is able to reproduce the fluid transient peak due to vapor collapse Subsequently, it can retain the shape of pressure wave for two phase fluid transients, but it is dampening faster compared to experiments Further investigation is required with set of experiments to validate the present numerical approach for different experimental data.

Fluid Transient↗
NASA NTRS1966

Air bearing analysis

Numerical analysis of air bearing system to reduce high frequency noise from vertical vibrations

NUMERICAL ANALYSIS↗