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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.↗
NASA NTRS1984

Numerical Simulations of Convectively Excited Gravity Waves

Magneto-convection and gravity waves are numerically simulated with a nonlinear, three-dimensional, time-dependent model of a stratified, rotating, spherical fluid shell heated from below. A Solar-like reference state is specified while global velocity, magnetic field, and thermodynamic perturbations are computed from the anelastic magnetohydrodynamic equations. Convective overshooting from the upper (superadiabatic) part of the shell excites gravity waves in the lower (subadiabatic) part. Due to differential rotation and Coriolis forces, convective cell patterns propagate eastward with a latitudinally dependent phase velocity. The structure of the excited wave motions in the stable region is more time-dependent than of the convective motions above. The magnetic field tends to be concentrated over giant-cell downdrafts in the convective zone but is affected very little by the wave motion in the stable region.

Glatzmaier, G. A.↗
NASA NTRS1987

Numerical simulations of turbulent spots in plane Poiseuille and boundary-layer flow

Direct numerical simulations of turbulent spots in plane Poiseuille and boundary-layer flows are performed. Mature, self-similar spots are obtained. The propagation velocities and spreading angles are found to compare well with corresponding experiments. The difference in shape of the two spots is also clearly discernible: the turbulent parts are contained within arrowhead regions that point in opposite directions for the two cases. The wing-tip region of the Poiseuille spot is also found to consist of a large-amplitude semiturbulent wave packet.

Henningson, Dan↗
NASA NTRS1988

Turbulence spectrum of a passive temperature field - Results of a numerical simulation

The spectrum of a passive temperature field G(k) is determined via numerical simulation using three kinds of isotropic turbulent velocity fields. For a velocity field developed concurrently with the temperature field, G(k) collapses onto the spectrum obtained using a frozen developed velocity field only for sufficiently high conductivity values. For lower conductivity values, the power law behavior of G(k) is less clear.

Chasnov, J.↗
NASA NTRS1990

Numerical simulations of the rebound shock model for solar spicules

Using time-dependent numerical simulations, the proposed rebound shock mechanism for spicules has been examined. At temperatures above a critical value, T(c), the radiation is characteristic of the conditions in the optically thin corona and near optically thin transition region. When T less than T(c), the atmosphere has a radiative cooling time, tau(rad) characteristic of chromosphere. The spicule is initiated with a quasi-impulsive force in the low chromosphere, which drives a train of upward propagating rebound shocks along the rigid magnetic flux tube. These shocks then move the transition region upward. The material below the displaced transition region has temperatures and densities similar to those of spicules when T(c) = 20,000 K or more and tau(rad) = 500 s or more, but not when T(c) = 10,000 K, and probably not when tau(rad) = 100 s. For all the cases where the cross sectional area diverges rapidly with height, the upward velocity of the transition region is less than that of spicules. Moreover, the maximum height is less than that of average spicules. Taller, higher velocity spicules result when the magnetic field cross sectional area is constant. In all cases, the rebound shock mechanism produces substantial motions and temperature and density variations in chromospheric and transition region material. It is suggested that this may be a partial explanation for the continuous dynamic state of the lower solar atmosphere.

Sterling, Alphonse C.↗
NASA NTRS1993

A numerical simulation of galaxy subcluster mergers

We present preliminary results of a 3-D numerical simulation of two merging subclusters of galaxies. By self-consistently modelling the intracluster gas and dark matter dynamics, we hope to gain insight as to how the dynamics of both relate to such observables as the cluster x-ray emission, radio source morphology, and velocity dispersions.

Roettiger, Kurt↗
NASA NTRS1993

Direct numerical simulation of turbulent H2-O2 combustion using reduced chemistry

Results of direct numerical simulations of hydrogen-oxygen combustion using a partial-equilibrium chemistry scheme in constant density, decaying, isotropic turbulence are reported. The simulations qualitatively reproduce many features of experimental results, such as superequilibrium radical species mole fractions, with temperature and major species mole fractions closer to chemical equilibrium. It was also observed that the peak reaction rates occur in narrow zones where the stoichiometric surface intersects regions of high scalar dissipation, as might be expected for combustion conditions close to chemical equilibrium. Another finding was that high OH mole fraction correspond more closely to the stoichiometric surface than to areas of high reaction rate for conditions of the simulations. Simulation results were compared to predictions of the Conditional Moment Closure model. This model was found to give good results for all quantities of interest when the conditionally averaged scalar dissipation was used in the prediction. When the nonconditioned average dissipation was used, the predictions compared well to the simulations for most of the species and temperature, but not for the reaction rate. The comparison would be expected to improve for higher Reynolds number flows, however.

Montgomery, Christopher J.↗
NASA NTRS1999

Numerical Simulation of Low-Density Shock-Wave Interactions

Computational Fluid Dynamics (CFD) numerical simulations of low-density shock-wave interactions for an incident shock impinging on a cylinder have been performed. Flow-field density gradient and surface pressure and heating define the type of interference pattern and corresponding perturbations. The maximum pressure and heat transfer level and location for various interaction types (i.e., shock-wave incidence with respect to the cylinder) are presented. A time-accurate solution of the Type IV interference is employed to demonstrate the establishment and the steadiness of the low-density flow interaction.

Glass, Christopher E.↗
NASA NTRS2002

Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer

The results from direct numerical simulations of a spatially evolving, supersonic, flat-plate turbulent boundary-layer flow, with free-stream Mach number of 2.25 are presented. The simulated flow field extends from a transition region, initiated by wall suction and blowing near the inflow boundary, into the fully turbulent regime. Distributions of mean and turbulent flow quantities are obtained and an analysis of these quantities is performed at a downstream station corresponding to Re(sub x)= 5.548 x10(exp 6) based on distance from the leading edge.

Gatski, T. B.↗
NASA NTRS2003

Numerical Simulation of the Water Cycle Change Over the 20th Century

We have used numerical models to test the impact of the change in Sea Surface Temperatures (SSTs) and carbon dioxide (CO2) concentration on the global circulation, particularly focusing on the hydrologic cycle, namely the global cycling of water and continental recycling of water. We have run four numerical simulations using mean annual SST from the early part of the 20th century (1900-1920) and the later part (1980-2000). In addition, we vary the CO2 concentrations for these periods as well. The duration of the simulations is 15 years, and the spatial resolution is 2 degrees. We use passive tracers to study the geographical sources of water. Surface evaporation from predetermined continental and oceanic regions provides the source of water for each passive tracer. In this way, we compute the percent of precipitation of each region over the globe. This can also be used to estimate precipitation recycling. In addition, we are using the passive tracers to independently compute the global cycling of water (compared to the traditional, Q/P calculation).

Bosilovich, Michael G.↗
DOE OSTI2024

Direct numerical simulation of a high-pressure hydrogen micromix combustor: Flame structure and stabilisation mechanism

A high-pressure hydrogen micromix combustor has been investigated using direct numerical simulation with detailed chemistry to examine the flame structure and stabilisation mechanism. The configuration of the combustor was based on the design by Schefer et al., using numerical periodicity to mimic a large square array. A precursor simulation of an opposed jet-in-crossflow was first conducted to generate appropriate partially-premixed inflow boundary conditions for the subsequent reacting simulation. The resulting flame can be described as an predominantly-lean inhomogeneously-premixed lifted jet flame. Five main zones were identified: a jet mixing region, a core flame, a peripheral flame, a recirculation zone, and combustion products. The core flame, situated over the jet mixing region, was found to burn as a thin reaction front, responsible for over 85% of the total fuel consumption. The peripheral flame shrouded the core flame, had low mean flow with high turbulence, and burned at very lean conditions (in the distributed burning regime). It was shown that turbulent premixed flame propagation was an order-of-magnitude too slow to stabilise the flame at these conditions. Stabilisation was identified to be due to ignition events resulting from turbulent mixing of fuel from the jet into mean recirculation of very lean hot products. Ignition events were found to correlate with shear-driven Kelvin-Helmholtz vortices, and increased in likelihood with streamwise distance. At the flame base, isolated events were observed, which developed into rapidly burning flame kernels that were blown downstream. Further downstream, near-simultaneous spatially-distributed ignition events were observed, which appeared more like ignition sheets. The paper concludes with a broader discussion that considers generalising from the conditions considered here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Direct numerical simulation of incompressible axisymmetric flows

In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Loulou, Patrick↗
NASA NTRS1991

Numerical simulations of the nonlinear kink modes in linearly stable supersonic slip surfaces

High-resolution numerical simulations of supersonic slip surfaces were performed using the Cray 2 at the Minnesota Supercomputer Center to confirm earlier analytic nonlinear stability calculations of such structures. The present study presents a numerical solution for the problem solved analytically by Artola and Majda (1987) and follows growing modes well into the large-amplitude regime. The response of an equal-density Mach-4 slip surface to a variety of impinging nonlinear sound wave trains is computed using the piecewise-parabolic method (PPM). The nature of the nonlinear kink modes observed in the simulations was found to be similar to that discussed by Artola and Majda. Most of the speeds moved in either direction, with speeds near the predicted value. It is suggested that the stationary modes will eventually dominate the flow at much later times, and it is argued that the stationary modes are more disruptive than the propagating kink modes.

Pedelty, Jeffrey A.↗
NASA NTRS1990

Efficient numerical simulation of electron states in quantum wires

A new algorithm is presented for the numerical simulation of electrons in a quantum wire as described by a two-dimensional eigenvalue problem for Schroedinger's equation coupled with Poisson's equation. Initially, the algorithm employs an underrelaxed fixed point iteration to generate an approximation which is reasonably close to the solution. Subsequently, this approximate solution is employed as an initial guess for a Jacobian-free implementation of an approximate Newton method. In this manner the nonlinearity in the model is dealt with effectively. The effectiveness of this approach is demonstrated in a set of numerical experiments which study the electron states on the cross section of a quantum wire structure based on III-V semiconductors at 4.2 and 77 K.

Kerkhoven, Thomas↗
NASA NTRS1977

Determination of thermospheric quantities from simple ionospheric observations using numerical simulation

Measured ionospheric electron content and peak electron concentration data are introduced into a numerical simulation of the ionosphere to yield values of induced plasma drifts and exospheric neutral temperatures consistent with the observations. Data collected on 23-24 March 1970 on the East Coast of the U.S.A. are analyzed and the results are in agreement with incoherent radar measurements at Millstone Hill, Massachusetts. Neutral winds and meridional exospheric temperature gradients that give rise to the computed plasma drifts are calculated through the use of a dynamic model of the thermosphere.

Antoniadis, D. A.↗
NASA NTRS2014

Analysis of Numerical Simulation Database for Pressure Fluctuations Induced by High-Speed Turbulent Boundary Layers

Direct numerical simulations (DNS) of Mach 6 turbulent boundary layer with nominal freestream Mach number of 6 and Reynolds number of Re(sub T) approximately 460 are conducted at two wall temperatures (Tw/Tr = 0.25, 0.76) to investigate the generated pressure fluctuations and their dependence on wall temperature. Simulations indicate that the influence of wall temperature on pressure fluctuations is largely limited to the near-wall region, with the characteristics of wall-pressure fluctuations showing a strong temperature dependence. Wall temperature has little influence on the propagation speed of the freestream pressure signal. The freestream radiation intensity compares well between wall-temperature cases when normalized by the local wall shear; the propagation speed of the freestream pressure signal and the orientation of the radiation wave front show little dependence on the wall temperature.

Duan, Lian↗