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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Numerical Simulation of Dual-Mode Scramjet Combustors

Results of a numerical investigation of a three-dimensional dual-mode scramjet isolator-combustor flow-field are presented. Specifically, the effect of wall cooling on upstream interaction and flow-structure is examined for a case assuming jet-to-jet symmetry within the combustor. Comparisons are made with available experimental wall pressures. The full half-duct for the isolator-combustor is then modeled in order to study the influence of side-walls. Large scale three-dimensionality is observed in the flow with massive separation forward on the side-walls of the duct. A brief review of convergence-acceleration techniques useful in dual-mode simulations is presented, followed by recommendations regarding the development of a reliable and unambiguous experimental data base for guiding CFD code assessments in this area.

Rodriguez, C. G.↗
NASA NTRS1992

Numerical simulation of compressor endwall and casing treatment flow phenomena

A numerical study is presented of the flow in the endwall region of a compressor blade row, in conditions of operation with both smooth and grooved endwalls. The computations are first compared to velocity field measurements in a cantilevered stator/rotating hub configuration to confirm that the salient features are captured. Computations are then interrogated to examine the tip leakage flow structure since this is a dominant feature of the endwall region. In particular, the high blockage that can exist near the endwalls at the rear of a compressor blade passage appears to be directly linked to low total pressure fluid associated with the leakage flow. The fluid dynamic action of the grooved endwall, representative of the casing treatments that have been most successful in suppressing stall, is then simulated computationally and two principal effects are identified. One is suction of the low total pressure, high blockage fluid at the rear of the passage. The second is energizing of the tip leakage flow, most notably in the core of the leakage vortex, thereby suppressing the blockage at its source.

Crook, A. J.↗
NASA NTRS1987

Numerical simulation of polar lows and comma clouds using simple dry models

Linear and nonlinear numerical channel models are used to simulate polar low/comma cloud evolution. The basic state of the model is idealized and can be specified analytically. It has strong horizontal as well as vertical wind shear and implicitly includes the effect of sensible heat fluxes from the lower boundary. This is done by having a layer of reduced stability present in the lower troposphere. The results of the linear model study are presented and used to provide input for the nonlinear model integration. It is shown that, given sufficient horizontal resolution, many features of polar low/comma clouds can be simulated by the nonlinear primitive equations model and the idealized basic state, in particular a rapidly growing disturbance with a wavelength on the order of 1000 km.

Mudrick, Stephen E.↗
NASA NTRS1989

Numerical simulation and comparison of symmetrical/supercritical airfoils for the near tip region of a helicopter in forward flight

Aerodynamic loads on a multi-bladed helicopter rotor in forward flight at transonic tip conditions are calculated. The unsteady, three-dimensional, time-accurate compressible Reynolds-averaged thin layer Navier-Stokes equations are solved in a rotating coordinate system on a body-conformed, curvilinear grid of C-H topology. Detailed boundary layer and global numerical comparisons of NACA-0012 symmetrical and CAST7-158 supercritical airfoils are made under identical forward flight conditions. The rotor wake effects are modeled by applying a correction to the geometric angle of attack of the blade. This correction is obtained by computing the local induced downwash velocity with a free wake analysis program. The calculations are performed on the Numerical Aerodynamic Simulation Cray 2 and the VPS32 (a derivative of a Cyber 205 at the Langley Research Center) for a model helicopter rotor in forward flight.

Badavi, F. F.↗
NASA NTRS2004

The QuakeSim Project: Numerical Simulations for Active Tectonic Processes

In order to develop a solid earth science framework for understanding and studying of active tectonic and earthquake processes, this task develops simulation and analysis tools to study the physics of earthquakes using state-of-the art modeling, data manipulation, and pattern recognition technologies. We develop clearly defined accessible data formats and code protocols as inputs to the simulations. these are adapted to high-performance computers because the solid earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes.

tectonic processes↗
NASA NTRS1985

Direct numerical simulation of equilibrium turbulent boundary layers

This paper describes the simulation of turbulent boundary layers by direct numerical solution of the three-dimensional, time-dependent Navier-Stokes equations, using a spectral method. The flow is incompressible, with Re sub delta = 1000. The equations are written in the self-similar coordinate system and periodic streamwise and spanwise boundary conditions are imposed. A family of nine 'equilibrium' boundary layers, from the strongly accelerated 'sink' flow to Stratford's separating boundary layer is treated. Good general agreement with experiments is observed. The effects of pressure gradients on the structures and statistics, both in the wall and wake regions are discussed.

Spalart, P. R.↗
NASA NTRS1993

Numerical simulation of upstream disturbance on flows around a slender body

Numerical solutions of the thin-layer approximation of the compressible Navier-Stokes equations have been obtained for flows around an ogive-cylinder body with and without a small fixed disturbance placed upstream of the body tip. Locating the disturbance at positions in the flow field upstream of the tip provokes the same range of behavior of the asymmetric flow that was numerically produced earlier by use of a geometrical disturbance on the body tip. Results remain consistent with the presence of a convective instability mechanism, and demonstrate the potential for a precise mapping of the body's receptivity to fixed disturbances in the flow field. Numerical solutions were also obtained for the flow-field responses to impulsive upstream disturbances to determine whether there is a growing response to an asymmetric impulsive disturbance that is consistent with presence of a convective instability mechanism. Results for surface pressure are interpreted with the aid of a mathematical model. The model suggests that the observed growth of surface pressure gradient with time and distance along a ray in response to an asymmetnc impulsive disturbance is in accord with the solution of a Ginzberg-Landau equation, with distinguishing features of the solution being consistent with the convective instability mode of behavior.

Degani, David↗
NASA NTRS1994

Fundamental mechanisms in premixed flame propagation via vortex-flame interactions: Numerical simulations

The goal of the present study is to assess numerically the ability of single-step and two-step chemical models to describe the main features encountered during the interaction between a two-dimensional vortex pair and a premixed laminar flame. In the two-step mechanism, the reaction kinetics are represented by a first chain branching reaction A + X yields 2X and a second chain termination reaction X + X yields P. This paper presents the fundamental mechanisms occurring during vortex-flame interactions and the relative impact of the major parameters encountered in turbulent premixed flames and suspected of playing a role in quenching mechanism: (1) Influence of stretch is investigated by analyzing the contribution of curvature and tangential strain on the local structure of the flame. The effect of Lewis number on the flame response to a strained field is analyzed. (2) Radiative heat losses which are suspected to be partially or totally responsible for quenching are also investigated. (3) The effect of the diffusion of the radicals is studied using a two-step mechanism in which an intermediate species is present. The parameters of the two-step mechanism are entirely determined from physical arguments. (4) Precise quantitative comparisons between the DNS and the experimental results of Samaniego et al are performed. These comparisons concern the evolution of the minimum heat release rate found along the flame front during the interaction and the distribution of the heat release rate along the flame front.

Mantel, Thierry↗
NASA NTRS2010

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Shih, Tsan-Hsing↗
NASA NTRS1988

Numerical simulations of an oblique detonation wave engine

An account is given of the numerical methods employed in a code for the simulation of supersonic combustion, which is then applied to the simulation of attached detonations and flames associated with the oblique-detonation wave supersonic combustor concept. The addition of heat by a detonation wave results in a shorter combustor than can be obtained in more conventional scramjet designs. Pure oblique detonations have been produced in a stoichiometric, uniformly mixed hydrogen/air stream; the wave rotates upstream with energy release, according to simple analytical arguments. Flow visualization maps for Mach number and temperature are presented.

Cambier, Jean-Luc↗
NASA NTRS1990

Full numerical simulation of coflowing, axisymmetric jet diffusion flames

The near field of a non-premixed flame in a low speed, coflowing axisymmetric jet is investigated numerically using full simulation. The time-dependent governing equations are solved by a second-order, explicit finite difference scheme and a single-step, finite rate model is used to represent the chemistry. Steady laminar flame results show the correct dependence of flame height on Peclet number and reaction zone thickness on Damkoehler number. Forced simulations reveal a large difference in the instantaneous structure of scalar dissipation fields between nonbuoyant and buoyant cases. In the former, the scalar dissipation marks intense reaction zones, supporting the flamelet concept; however, results suggest that flamelet modeling assumptions need to be reexamined. In the latter, this correspondence breaks down, suggesting that modifications to the flamelet modeling approach are needed in buoyant turbulent diffusion flames.

Mahalingam, S.↗
NASA NTRS2005

Numerical Simulation of Liquid Jet Atomization Including Turbulence Effects

This paper describes numerical implementation of a newly developed hybrid model, T-blob/T-TAB, into an existing computational fluid dynamics (CFD) program for primary and secondary breakup simulation of liquid jet atomization. This model extend two widely used models, the Kelvin-Helmholtz (KH) instability of Reitz (blob model) and the Taylor-Analogy-Breakup (TAB) secondary droplet breakup by O'Rourke and Amsden to include turbulence effects. In the primary breakup model, the level of the turbulence effect on the liquid breakup depends on the characteristic scales and the initial flow conditions. For the secondary breakup, an additional turbulence force acted on parent drops is modeled and integrated into the TAB governing equation. Several assessment studies are presented and the results indicate that the existing KH and TAB models tend to under-predict the product drop size and spray angle, while the current model provides superior results when compared with the measured data.

Trinh, Huu P.↗
NASA NTRS1984

Numerical simulations of the Ross Sea tides

Tidal currents below the floating Ross Ice shelf are reconstructed by using a numerical tidal model. They are predominantly diurnal, achieve maximum strength in regions near where the ice shelf runs aground, and are significantly enhanced by topographic Rossby wave propagation along the ice front. A comparison with observations of the vertical motion of the ice shelf surface indicates that the model reproduces the diurnal tidal characteristics within 20 percent. Similar agreement for the relatively weak semidiurnal tides was not obtained, and this calls attention to possible errors of the open boundary forcing obtained from global-ocean tidal simulations and to possible errors in mapping zones of ice shelf grounding. Air-sea contact below the ice shelf is eliminated by the thick ice cover. The dominant sub-ice-shelf circulation may thus be tidally induced. A preliminary assessment of sub-ice-shelf conditions based on the numerical tidal simulations suggests that (1) strong barotropic circulation is driven along the ice front and (2) tidal fronts may form in the sub-ice-shelf cavity where the water column is thin and where the buoyancy input is weak.

Macayeal, D. R.↗
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗
NASA NTRS1991

Numerical simulations of a nonequilibrium argon plasma in a shock-tube experiment

A code developed for the numerical modeling of nonequilibrium radiative plasmas is applied to the simulation of the propagation of strong ionizing shock waves in argon gas. The simulations attempt to reproduce a series of shock-tube experiments which will be used to validate the numerical models and procedures. The ability to perform unsteady simulations makes it possible to observe some fluctuations in the shock propagation, coupled to the kinetic processes. A coupling mechanism by pressure waves, reminiscent of oscillation mechanisms observed in detonation waves, is described. The effect of upper atomic levels is also briefly discussed.

Cambier, Jean-Luc↗
NASA NTRS1994

Oblique magnetohydrodynamic cosmic-ray-modified shocks: Two-fluid numerical simulations

We present the first results of time-dependent two-fluid cosmic-ray (CR) modified magnetohydrodynamic (MHD) shock simulations. The calculations were carried out with a new numerical code for one-dimensional ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks, including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they apppear to have converged to dynamical steady states. We also present results from an initial time-dependent simulation which extends the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescales than parallel shocks.

Frank, Adam↗