Search NASA⌕ Search

SEARCH · Search NASA

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical simulation of fluid flow and heat transfer in a thin liquid film over a stationary and rotating disk and comparison with experimental data

In the first section, improvements in the theoretical model and computational procedure for the prediction of film height and heat-transfer coefficient of the free surface flow of a radially-spreading thin liquid film adjacent to a flat horizontal surface of finite extent are presented. Flows in the presence and absence of gravity are considered. Theoretical results are compared to available experimental data with good agreement. In the presence of gravity, a hydraulic jump is present, isolating the flow into two regimes: supercritical upstream from the jump and subcritical downstream of it. In this situation, the effects of surface tension are important near the outer edge of the disk where the fluid experiences a free fall. A region of flow separation is present just downstream of the jump. In the absence of gravity, no hydraulic jump or separated flow region is present. The variation of the heat-transfer coefficient for flows in the presence and absence of gravity are also presented. In the second section, the results of a numerical simulation of the flow field and associated heat transfer coefficients are presented for the free surface flow of a thin liquid film adjacent to a horizontal rotating disk. The computation was performed for different flow rates and rotational velocities using a 3-D boundary-fitted coordinate system. Since the geometry of the free surface is unknown and dependent on flow rate, rate of rotation, and other parameters, an iterative procedure had to be used to ascertain its location. The computed film height agreed well with existing experimental measurements. The flow is found to be dominated by inertia near the entrance and close to the free surface and dominated by centrifugal force at larger radii and adjacent to the disk. The rotation enhances the heat transfer coefficient by a significant amount.

Faghri, Amir↗
NASA NTRS1995

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗
NASA NTRS1992

Numerical simulations of cometary dust

Most observations of comets are done photometrically or spectrophotometrically. The interpretation of the aperture-averaged flux is relatively simple for an isotropic, radially expanding coma of infinite extent - the canonical model. However, the interpretation of the observations is not so clear when the motion of the dust is affected by radiation pressure, or when the emission is time-varying and anisotropic. For example, in a sample of CCD images of 10 comets, Jewitt and Meech (1987, Ap.J. 317, 992) found that the photometric profiles of only three comets were consistent, within the observational errors, with the profiles predicted form the canonical model. Photometric observations with large apertures, however, seem to suggest that the canonical model may be quite adequate (c.f. Osip, Schleicher, and Millis, 1992, Icarus 98, 115). The dust itself is characterized by a size distribution, with size dependencies on the expansion velocity, the scattered and thermal radiation, the response to radiation pressure, and probably the density. How good then are the approximations normally used in determining the production rates of the dust when these effects are present? As part of a program to better understand the dynamics of cometary dust and gas, a computer program has been developed which numerically simulates the emission of both dust and gas from a tilted rotating nucleus.

Lien, David J.↗
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↗
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.↗
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.↗