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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a drop in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. The terms in the budget of the Reynolds stress were computed. In general, the individual terms do not show a significant change. However, an increase in dissipation rate rather than reduction in total production appears to be responsible for the drop in turbulent kinetic energy.

Moin, P.↗
NASA NTRS1990

Direct numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a decrease in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. Terms in the Reynolds stress transport equations were computed. The balances indicate that the decrease in turbulent kinetic energy is a result of a decrease in turbulence production, along with an increase in turbulent dissipation. Intuitive reasoning and current turbulence models would predict an increase in kinetic energy along with increases in production and dissipation rates as a result of increased mean-flow strain rate. Later in the evolution of the flow, both turbulence production and dissipation increase.

Moin, P.↗
NASA NTRS2017

Direct Numerical Simulation of Hypersonic Turbulent Boundary Layer inside an Axisymmetric Nozzle

As a first step toward a study of acoustic disturbance field within a conventional, hypersonic wind tunnel, direct numerical simulations (DNS) of a Mach 6 turbulent boundary layer on the inner wall of a straight axisymmetric nozzle are conducted and the results are compared with those for a flat plate. The DNS results for a nozzle radius to boundary-layer thickness ratio of 5:5 show that the turbulence statistics of the nozzle-wall boundary layer are nearly unaffected by the transverse curvature of the nozzle wall. Before the acoustic waves emanating from different parts of the nozzle surface can interfere with each other and undergo reflections from adjacent portions of the nozzle surface, the rms pressure fluctuation beyond the boundary layer edge increases toward the nozzle axis, apparently due to a focusing effect inside the axisymmetric configuration. Spectral analysis of pressure fluctuations at both the wall and the freestream indicates a similar distribution of energy content for both the nozzle and the flat plate, with the peak of the premultiplied frequency spectrum at a frequency of [(omega)(delta)]/U(sub infinity) approximately 6.0 inside the free stream and at [(omega)(delta)]/U(sub infinity) approximately 2.0 along the wall. The present results provide the basis for follow-on simulations involving reverberation effects inside the nozzle.

Huang, Junji↗
NASA NTRS2002

Numerical Simulation of the Combustion of Fuel Droplets: Finite Rate Kinetics and Flame Zone Grid Adaptation (CEFD)

The NASA Nebraska Space Grant (NSGC) & EPSCoR programs have continued their effort to support outstanding research endeavors by funding the Numerical Simulation of the Combustion of Fuel Droplets study at the University of Nebraska at Lincoln (UNL). This team of researchers has developed a transient numerical model to study the combustion of suspended and moving droplets. The engines that propel missiles, jets, and many other devices are dependent upon combustion. Therefore, data concerning the combustion of fuel droplets is of immediate relevance to aviation and aeronautical personnel, especially those involved in flight operations. The experiments being conducted by Dr. Gogos and Dr. Nayagam s research teams, allow investigators to gather data for comparison with theoretical predictions of burning rates, flame structures, and extinction conditions. The consequent improved hndamental understanding droplet combustion may contribute to the clean and safe utilization of fossil hels (Williams, Dryer, Haggard & Nayagam, 1997, 72). The present state of knowledge on convective extinction of he1 droplets derives fiom experiments conducted under normal gravity conditions. However, any data obtained with suspended droplets under normal gravity are grossly affected by gravity. The need to obtain experimental data under microgravity conditions is therefore well justified and addresses one of the goals of NASA s Human Exploration and Development of Space (HEDS) microgravity combustion experiment.

Gogos, George↗
NASA NTRS1995

Numerical simulation of supersonic and hypersonic inlet flow fields

This report summarizes the research performed by North Carolina State University and NASA Ames Research Center under Cooperative Agreement NCA2-719, 'Numerical Simulation of Supersonic and Hypersonic Inlet Flow Fields". Four distinct rotated upwind schemes were developed and investigated to determine accuracy and practicality. The scheme found to have the best combination of attributes, including reduction to grid alignment with no rotation, was the cell centered non-orthogonal (CCNO) scheme. In 2D, the CCNO scheme improved rotation when flux interpolation was extended to second order. In 3D, improvements were less dramatic in all cases, with second order flux interpolation showing the least improvement over grid aligned upwinding. The reduction in improvement is attributed to uncertainty in determining optimum rotation angle and difficulty in performing accurate and efficient interpolation of the angle in 3D. The CCNO rotational technique will prove very useful for increasing accuracy when second order interpolation is not appropriate and will materially improve inlet flow solutions.

Mcrae, D. Scott↗
NASA NTRS2020

Direct Numerical Simulation of Fog Formation Within the Surface Layer

We investigate the effects of shear-driven turbulence on fog formation in the stably stratified surface layer, using recent direct numerical simulations (DNS) of dry and moist open-channel flows. A constant cooling rate applied at the ground mimics longwave radiative cooling, leading to either turbulent (weakly stable) or laminar (very stable) flows. The transition of the Obukhov, Ozmidov and buoyancy length scales are investigated for different cooling rates and Reynolds numbers. Compared to the dry case, the condensation of liquid water in the moist case supports slightly higher cooling rates before leading to turbulence collapse. In the very stable moist cases, runaway cooling leads to substantial condensation of liquid water close to the ground and fog (visibility less than 1 km) results over much of the domain. In the weakly stable cases, turbulent mixing narrowly yields visibilities of 1 km close to the ground over a similar time period. Despite the idealised nature of the system, the present results suggest that turbulence impedes, although will not necessarily inhibit, fog formation.

Teixeira, Joao↗
NASA NTRS1977

Dynamical response of the solar corona. III - Numerical simulation of the 1973 June 10 coronal transient

A spherically symmetric adiabatic single-fluid model is outlined for simulating the nonlinear time-dependent response of the corona to solar events that are simulated by perturbations in the appropriate physical variables from their steady-state values at the coronal base. Several observed features of the coronal transient that occurred on June 10, 1973, are simulated by using a particular steady-state solar wind, a specific combination of density and temperature perturbations, and a particular time dependence of the perturbations. A different steady-state solar wind, a perturbation of shorter duration, and other perturbation combinations are also employed to simulate the same transient so that the effect of each quantity can be determined. It is found that the model cannot adequately simulate all the observational results for the investigated transient, that the steady-state solar wind is relatively unimportant in such numerical simulations, and that studies which attempt only to reproduce observed shock trajectories may lead to erroneous conclusions regarding the physics of the solar event that produced the transient.

Steinolfson, R. S.↗
NASA NTRS1991

Spatial Direct Numerical Simulation of Boundary-Layer Transition Mechanisms: Validation of PSE Theory

A study of instabilities in incompressible boundary-layer flow on a flat plate is conducted by spatial direct numerical simulation (DNS) of the Navier-Stokes equations. Here, the DNS results are used to critically evaluate the results obtained using parabolized stability equations (PSE) theory and to study mechanisms associated with breakdown from laminar to turbulent flow. Three test cases are considered: two-dimensional Tollmien-Schlichting wave propagation, subharmonic instability breakdown, and oblique-wave break-down. The instability modes predicted by PSE theory are in good quantitative agreement with the DNS results, except a small discrepancy is evident in the mean-flow distortion component of the 2-D test problem. This discrepancy is attributed to far-field boundary- condition differences. Both DNS and PSE theory results show several modal discrepancies when compared with the experiments of subharmonic breakdown. Computations that allow for a small adverse pressure gradient in the basic flow and a variation of the disturbance frequency result in better agreement with the experiments.

Joslin, R. D.↗
NASA NTRS1986

Numerical simulation of tip vortices of wings in subsonic and transonic flows

A multi block zonal algorithm which solves the thin-layer Navier-Stokes and the Euler equations is used to numerically simulate the formation and roll-up of the tip vortex in both subsonic and transonic flows. Four test cases which used small and large aspect ratio wings have been considered to examine the influence of the tip-cap shape, the tip planform and the free-stream Mach number. It appears that both the tip-planform and the tip-cap shape have some influence on the formation of the tip vortex, but its subsequent roll-up seems to be more influenced by the tip-planform shape. In general, a good definition of the formation and the roll-up of the tip vortex has been observed for all the cases considered here. Comparions of the numerical results with the limited, available experimental data show good agreement with both the surface pressures and the tip-vortex strength.

Srinivasan, G. R.↗
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
DOE OSTI2024

Numerical Simulations in Support of a Long-Term Test of Gas Production From Hydrate Accumulations on the Alaska North Slope: Water Production and Associated Design and Management Issues

Here, we investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.

02 PETROLEUM↗
NASA NTRS1980

Numerical simulation of sudden stratospheric warmings

The paper presents a mechanistic, quasi-geostrophic, semi-spectral model with a self-consistent calculation of the mean zonal flow fields used to numerically simulate sudden stratospheric warmings generated by a single zonal harmonic (m) planetary wave. The development of a warming depends on the strength of the westerly winds in the lower stratosphere, and the growth rate of a m=2 warming is accelerated during the westerly zonal wind phase of the quasi-biennial oscillation. An m=1 warming is characterized by oscillation of wave amplitude and mean flow from resonantly trapped, westward, propagating planetary waves moving in and out of phase with the tropospherically forced stationary planetary wave. An m=2 critical level develops in the equatorial region and advances poleward, and its warming is characterized by a sudden intensification after an initially slow growth in contrast to slowly developing m=1 warmings.

Schoeberl, M. R.↗
NASA NTRS1993

Numerical simulation of stability and stability control of high speed compressible rotating couette flow

The nonlinear temporal evolution of disturbances in compressible flow between infinitely long, concentric cylinders is investigated through direct numerical simulations of the full, three-dimensional Navier-Stokes and energy equations. Counter-rotating cylinders separated by wide gaps are considered with supersonic velocities of the inner cylinder. Initially, the primary disturbance grows exponentially in accordance with linear stability theory. As the disturbances evolve, higher harmonics and subharmonics are generated in a cascading order eventually reaching a saturation state. Subsequent highly nonlinear stages of the evolution are governed by the interaction of the disturbance modes, particularly the axial subharmonics. Nonlinear evolution of the disturbance field is characterized by the formation of high-shear layers extending from the inner cylinder towards the center of the gap in the form of jets similar to the ejection events in transitional and turbulent wall-bounded shear flows.

Biringen, Sedat↗
NASA NTRS1995

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗
NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility demonstrated. Application of the surrogate equation technique allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations or their various approximations. Computational results are presented for the full Euler equations and for the transonic disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗
NASA NTRS1992

Numerical simulations of inductive-heated float-zone growth

The present work provides an improved fluid flow and heat-transfer modeling of float-zone growth by introducing a RF heating model so that an ad hoc heating temperature profile is not necessary. Numerical simulations were carried out to study the high-temperature float-zone growth of titanium carbide single crystal. The numerical results showed that the thermocapillary convection occurring inside the molten zone tends to increase the convexity of the melt-crystal interface and decrease the maximum temperature of the molten zone, while the natural convection tends to reduce the stability of the molten zone by increasing its height. It was found that the increase of induced heating due to the increase of applied RF voltage is reduced by the decrease of zone diameter. Surface tension plays an important role in controlling the amount of induced heating. Finally, a comparison of the computed shape of the free surface with a digital image obtained during a growth run showed adequate agreement.

Chan, Y. T.↗
NASA NTRS1994

Numerical simulation of incidence and sweep effects on delta wing vortex breakdown

The structure of the vortical flowfield over delta wings at high angles of attack was investigated. Three-dimensional Navier-Stokes numerical simulations were carried out to predict the complex leeward-side flowfield characteristics, including leading-edge separation, secondary separation, and vortex breakdown. Flows over a 75- and a 63-deg sweep delta wing with sharp leading edges were investigated and compared with available experimental data. The effect of variation of circumferential grid resolution grid resolution in the vicinity of the wing leading edge on the accuracy of the solutions was addressed. Furthermore, the effect of turbulence modeling on the solutions was investigated. The effects of variation of angle of attack on the computed vortical flow structure for the 75-deg sweep delta wing were examined. At moderate angles of attack no vortex breakdown was observed. When a critical angle of attack was reached, bubble-type vortex breakdown was found. With further increase in angle of attack, a change from bubble-type breakdown to spiral-type vortex breakdown was predicted by the numerical solution. The effects of variation of sweep angle and freestream Mach number were addressed with the solutions on a 63-deg sweep delta wing.

INCIDENCE↗
NASA NTRS2011

Numerical Simulations on Origin of Galilean Moons' Magnetic Anomalies

Galileo mission detected the magnetic anomalies originated from Galilean moons. These anomalies are likely generated in the moons interiors, under the influence of a strong ambient Jovian field. Among various possible generation mechanisms of the anomalies, we focus on magneto-convection and dynamos in the interiors via numerical simulation. To mimic the electromagnetic environment of the moons, we introduce in our numerical model an external uniform magnetic field B(sub 0) with a fixed orientation but varying field strength. Our results show that a finite B(sub 0) can substantially alter the dynamo processes inside the core. When the ambient field strength B(sub 0) increases to approximately 40% of the field generated by the pure dynamo action, the convective state in the core changes significantly: the convective flow decreases by 80% in magnitude, but the differential rotation becomes stronger in much of the fluid layer, leading to a stronger field generated in the core. The field morphologies inside the core tend to align with the ambient field, while the flow patterns show the symmetry-breaking effect under the influence of B(sub 0). Furthermore, the generated field tends to be temporally more stable.

Jiao, LiQuo↗