Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulation”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Numerical simulation of shock-induced combustion/detonation

A numerical study is conducted to simulate the shock-induced combustion in premixed H2-air mixtures. Two types of bodies, blunt (spherical projectile) and sharp (wedge), are considered in the study. A nine-species, 18-step finite-rate H2-air chemical reaction mechanism coupled with the Navier-Stokes equations is solved. The flow field over the blunt body is found to be unsteady, when the projectile velocity is same as the Chapman-Jouget velocity of the mixture. The unsteadiness is caused by the periodic instabilities originating in the stagnation zone. Numerical results show good qualitative agreement with the ballistic range shadowgraph. In addition, the frequency of oscillations, determined by using the Fourier power spectrum, is found to be in good agreement with the experiment. The flow field over the wedge is found to be stable for the conditions considered in this study. The oblique detonation wave structure is investigated and the important flow features are discussed.

Singh, D. J.↗
NASA NTRS1992

Three-dimensional numerical simulation of the 20 June 1991, Orlando microburst

On 20 June 1991, NASA's Boeing 737, equipped with in-situ and look-ahead wind-shear detection systems, made direct low-level penetrations (300-350 m AGL) through a microburst during several stages of its evolution. This microburst was located roughly 20 km northeast of Orlando International Airport and was monitored by a Terminal Doppler Weather Radar (TDWR) located about 10 km south of the airport. The first NASA encounter with this microburst (Event 142), at approximately 2041 UTC, was during its intensification phase. At flight level, in-situ measurements indicated a peak 1-km (averaged) F-factor of approximately 0.1. The second NASA encounter (Event 143) occurred at approximately 2046 UTC, about the time of microburst peak intensity. It was during this penetration that a peak 1-km F-factor of approximately 17 was encountered, which was the largest in-situ measurement of the 1991 summer deployment. By the third encounter (Event 144), at approximately 2051 UTC, the microburst had expanded into a macroburst. During this phase of evolution, an in-situ 1-km F-factor of 0.08 was measured. The focus of this paper is to examine this microburst via numerical simulation from an unsteady, three-dimensional meteorological cloud model. The simulated high-resolution data fields of wind, temperature, radar reflectivity factor, and precipitation are closely examined so as to derive information not readily available from 'observations' and to enhance our understanding of the actual event. Characteristics of the simulated microburst evolution are compared with TDWR and in-situ measurements.

Proctor, Fred H.↗
NASA NTRS1993

Numerical simulation of non-Newtonian free shear flows

Free shear flows, like those of mixing layers, are encountered in aerodynamics, in the atmosphere, and in the ocean as well as in many industrial applications such as flow reactors or combustion chambers. It is, therefore, crucial to understand the mechanisms governing the process of transition to turbulence in order to predict and control the evolution of the flow. Delaying transition to turbulence as far downstream as possible allows a gain in energy expenditure while accelerating the transition can be of interest in processes where high mixing is desired. Various methods, including the use of polymer additives, can be effective in controlling fluid flows. The drag reduction obtained by the addition of small amounts of high polymers has been an active area of research for the last three decades. It is now widely believed that polymer additives can affect the stability of a large variety of flows and that dilute solutions of these polymers have been shown to produce drag reductions of over 80 percent in internal flows and over 60 percent in external flows under a wide range of conditions. The major thrust of this work is to study the effects of polymer additives on the stability of the incompressible mixing layer through large scale numerical simulations. In particular, we focus on the two dimensional flow and examine how the presence of viscoelasticity may affect the typical structures of the flow, namely roll-up and pairing of vortices.

Homsy, G. M.↗
NASA NTRS2013

Numerical Simulations of Instabilities in Single-Hole Office Elements

An orifice element is commonly used in liquid rocket engine test facilities either as a flow metering device, a damper for acoustic resonance or to provide a large reduction in pressure over a very small distance in the piping system. While the orifice as a device is largely effective in stepping down pressure, it is also susceptible to a wake-vortex type instability that generates pressure fluctuations that propagate downstream and interact with other elements of the test facility resulting in structural vibrations. Furthermore in piping systems an unstable feedback loop can exist between the vortex shedding and acoustic perturbations from upstream components resulting in an amplification of the modes convecting downstream. Such was the case in several tests conducted at NASA as well as in the Ariane 5 strap-on P230 engine in a static firing test where pressure oscillations of 0.5% resulted in 5% thrust oscillations. Exacerbating the situation in cryogenic test facilities, is the possibility of the formation of vapor clouds when the pressure in the wake falls below the vapor pressure leading to a cavitation instability that has a lower frequency than the primary wake-vortex instability. The cavitation instability has the potential for high amplitude fluctuations that can cause catastrophic damage in the facility. In this paper high-fidelity multi-phase numerical simulations of an orifice element are used to characterize the different instabilities, understand the dominant instability mechanisms and identify the tonal content of the instabilities.

Ahuja, Vineet↗
NASA NTRS1986

Numerical simulations of a siphon mechanism for quiescent prominence formation

Quiescent prominences represent a significant challenge to our understanding of the flow of mass and energy in the outer layers of the solar atmosphere. A small number of quiescent prominences contain as much mass as the entire corona (Athay, 1976). The problem then is how to get that much material into the relatively small volume of a prominence and maintain it at a temperature of 10,000 K in close proximity to material at one million K. The thermal insulation to conduction provided by the magnetic field explains the disparate temperatures. The mass source problem is less well understood. One method for supplying mass to the prominence is to siphon it from the chromosphere. The siphon mechanism begins with a magnetic loop that evolves into a configuration with a gravitational well, such as that described by Kippenhahn and Schluter (1957). This could be formed, for example, by a twist in the magnetic field. A gravitational well could also be formed by a condensation induced sag in the field. This could further enhance the condensation process. Once this well has formed, or as it is forming, the material in the well area of the loop must cool and condense to the point where radiative losses exceed any heat input. Additional material must also flow into the well from the underlying chromosphere to supply the mass required to form the prominence. One example from a series of numerical simulations that were performed to study the formation of quiescent prominences is presented.

Poland, A. I.↗
NASA NTRS1991

Procedures for Including Secondary Electron Emission in Numerical Simulations of Plasma-Insulator Interactions

Previous Monte Carlo simulations provide a data base for properties of secondary electron emission (SEE) from insulators and metals. Incident primary electrons are considered at energies up to 1200 eV. The behavior of secondary electrons is characterized by (1) yield vs. primary energy E(sub p), (2) distribution vs. secondary energy E(sub s), and (3) distribution vs. angle of emission theta. Special attention is paid to the low energy range E(sub p) up to 50 eV, where the number and energy of secondary electrons is limited by the finite band gap of the insulator. For primary energies above 50 eV the SEE yield curve can be conveniently parameterized by a Haffner formula. The energy distribution of secondary electrons is described by an empirical formula with average energy about 8.0 eV. The angular distribution of secondaries is slightly more peaked in the forward direction than the customary cos theta distribution. Empirical formulas and parameters are given for all yield and distribution curves. Procedures and algorithms are described for using these results to find the SEE yield, and then to choose the energy and angle of emergence of each secondary electron. These procedures can readily be incorporated into numerical simulations of plasma-solid surface interactions in low earth orbit.

Beyst, Brian↗
NASA NTRS1995

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Lin, Yuh-Lang↗
NASA NTRS1993

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The first section is on 3-D numerical modeling of terrain-induced circulations and covers the following: (1) additional insights into gravity wave generation mechanisms based on the control simulation; (2) ongoing nested-grid numerical simulations; (3) work to be completed during the remainder of FY-93; and (4) work objectives for FY-94. The second section is on linear theory and theoretical modeling and covers the following: (1) the free response of a uniform barotropic flow to an initially stationary unbalanced (ageostrophic) zonal wind anomaly; and (2) the free response of a uniform barotropic flow to an initially stationary balanced zonal wind anomaly.

Lin, Yuh-Lang↗
NASA NTRS1987

Numerical simulation of transition

The paper presents numerical algorithms for studying the physics of transition and turbulence in simple geometries. The highly nonlinear stages of transition prior to turbulent spot formation are studied in detail. The use of simulations to study the sensitivity of laminar flow control techniques in the nonlinear regime is demonstrated. A new instability mechanism associated with the center modes in channel flows is revealed.

Hussaini, M. Y.↗
NASA NTRS1995

Software For Numerical Simulation Of Geophysical Flows

Program verified by comparisons with both experimental and numerical studies. GEOSIM implements numerical model simulating geophysical fluid flow for wide range of problems. Allows for more accurate control over experimental conditions and provides complete data source for performing diagnostic studies. Used by experienced and/or professional fluid dynamicists. Written in FORTRAN 77.

Miller, Timothy L.↗
NASA NTRS1998

Numerical Simulations of a Lobed Fuel Injector

Numerical modeling of the nonreactive mixing processes associated with a lobed fuel injector in a coflowing air stream is presented. The lobed fuel injector is a device which generates strong streamwise vorticity, producing locally high strain rates which can enhance the molecular mixing of reactants while delaying ignition in a controlled manner. Vortex element modeling is used to simulate flow field evolution and fuel element mixing characteristics for this lobed fuel injector. Quantitative predictions for vorticity generation and qualitative results for streamwise rollup compare well qualitatively with recent experimental investigations of this flow field [Smith et al, Phys. Fluids 9, 667 (1997)]. Parametric studies of the effects of lobe amplitude-to-wavelength ratio, lobe angle, and lobe shape for given flow conditions suggest that geometrical features may be optimized to enhance mixing and control reaction processes.

Strickland, J. H.↗
NASA NTRS1985

Numerically-simulated formation and propagation of interplanetary shocks

The present numerical method for simulating the formation and propagation of interplanetary shocks is based on the shock-capturing finite difference scheme of Lax (1950) and Lax and Wendroff (1960), as well as the recent method of NEAR characteristics of Nakagawa (1980, 1981). Attention is given to examples which strongly suggest that all the shocked solar wind plasma parameters due to given physical perturbations, such as flare-generated shocks, can be predicted through the use of this method; the method is, however, limited to the supersonic and super-Alfvenic flow.

Wu, S. T.↗
NASA NTRS1992

Numerical simulation of conservation laws

A new numerical framework for solving conservation laws is being developed. This new approach differs substantially from the well established methods, i.e., finite difference, finite volume, finite element and spectral methods, in both concept and methodology. The key features of the current scheme include: (1) direct discretization of the integral forms of conservation laws, (2) treating space and time on the same footing, (3) flux conservation in space and time, and (4) unified treatment of the convection and diffusion fluxes. The model equation considered in the initial study is the standard one dimensional unsteady constant-coefficient convection-diffusion equation. In a stability study, it is shown that the principal and spurious amplification factors of the current scheme, respectively, are structurally similar to those of the leapfrog/DuFort-Frankel scheme. As a result, the current scheme has no numerical diffusion in the special case of pure convection and is unconditionally stable in the special case of pure diffusion. Assuming smooth initial data, it will be shown theoretically and numerically that, by using an easily determined optimal time step, the accuracy of the current scheme may reach a level which is several orders of magnitude higher than that of the MacCormack scheme, with virtually identical operation count.

Chang, Sin-Chung↗
DOE OSTI2017

Supercritical CO 2 Water Displacements and CO 2 Capillary Trapping: Micromodel Experiment and Numerical Simulation

We report that the CO 2 capillary trapping is an important scientific issue in geological carbon sequestration, but few researches focus on the trapping mechanism at pore scale under supercritical CO 2 condition. In this study, based on the high-pressure fluids-microscopy-micromodel experimental system, we performed drainage experiment, i.e. supercritical CO 2 displacing water, and imbibition experiment, i.e. water displacing CO 2 , under the conditions of 45°C and 8.5 MPa. The DSLR camera was used to capture pictures of CO 2 -water two-phase immiscible flow and the microscopy was used to capture the capillary trapping behavior for the supercritical CO 2 at the pore scale. The computational fluid dynamic method was adopted to simulate the two-phase fluid flow processes. The numerical results are generally in agree ment with the experimental observations, and further provide three-dimensional geometries on the interface during the drainage-imbibition processes and the trapped supercritical CO 2 droplet/cluster. Lastly, the capillary trapping curve, i.e. the relationship between the initial CO 2 saturation and the residual saturation, was obtained from the numerical results, and we made an assessment of the three capillary trapping models, i.e. Land's, Jurauld's and Spiteri's trapping models. A comparison of the models performance indicates that Jurauld's model behaves slightly better than Land's model, whereas Spiteri's model behaves poorly. However, given that Land's model only contains one parameter of clear physical meaning, it is recommended for practical use.

58 GEOSCIENCES↗
NASA NTRS1997

A Numerical Simulation of a Normal Sonic Jet into a Hypersonic Cross-Flow

This study involves numerical modeling of a normal sonic jet injection into a hypersonic cross-flow. The numerical code used for simulation is GASP (General Aerodynamic Simulation Program.) First the numerical predictions are compared with well established solutions for compressible laminar flow. Then comparisons are made with non-injection test case measurements of surface pressure distributions. Good agreement with the measurements is observed. Currently comparisons are underway with the injection case. All the experimental data were generated at the Southampton University Light Piston Isentropic Compression Tube.

Jeffries, Damon K.↗
DOE OSTI2014

Numerical Simulation and Optimization of the Thermodynamic Process of the Molten Salt Furnace

In this paper, a numerical model of the molten salt furnace process was developed, by using computational fluid dynamics (CFD) technique and considering the gas flow, the combustion and the radiation heat transfer. The results demonstrate that the performances of the salt furnace could be improved by optimization using the numerical model. The temperatures along the circumference of the furnace coil and outside shell are quite even, and the deviant combustion phenomenon is not observed. A back-flow formed in the upper part of the furnace chamber enhances the circulation and the mixing of the gas and effectively stabilizes the combustion in the furnace. The behaviors of CO, CO 2 , NO x , and H 2 O are presented in terms of the gas flow, temperature distribution and volumetric concentration distribution. It is concluded that the furnace with the constant air flow rate of 15,500 Nm 3 /h and the guiding vane angle at 48–50 deg is optimized for the combustion effectiveness.

42 ENGINEERING↗
NASA NTRS1989

Numerical simulations of an isolated microbursts. II - Sensitivity experiments

A time-dependent, high-resolution, axisymmetric numerical model simulating isolated and stationary microbursts is used to evaluate the sensitivity of microbursts to the environment and other factors. It is found that microbursts intensity is sensitive to ther vertical distirbution of ambient temperature and humidity, the horiontal width of the precipitation shaft or downdraft, the magnitutde of precipitation loading, and the type and duration of precipitation. Scenarios for the generation of an intense microburst and the possible applications of the results are examined.

Proctor, Fred H.↗
NASA NTRS1993

Numerical simulation of X-ray-heated winds in Seyfert galaxies. I - The case of zero angular momentum

Numerical hydrodynamics simulations are presented of X-ray heated winds believed to form the 'mirrors' which reflect light from obscured Seyfert galaxy nuclei into our line of sight. Above a likely range of parameters, it is shown that steady state flows are created when a cool torus surrounds an AGN. Some of the mass evaporated off the inner edge of the torus is captured by the central gravitational well, while the rest is driven off in a warm wind. When the central gravity is strong relative to the radiative heating rate, i.e., L/L(E) is small, nearly all the injected matter is captured. When the gravity is weaker, a given combination of heating rate and injected pressure leads to a maximum mass loss rate; when the rate at which mass is injected exceeds this amount, the excess is captured. The dividing line between these two cases comes at L/L(E) = about 0.08.

Balsara, Dinshaw S.↗