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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Numerical simulation of low gravity draining

A boundary value problem was solved numerically for a liquid that is assumed to be inviscid and incompressible, having a motion that is irrotational and axisymmetric, and having a constant (5 degrees) solid-liquid contact angle. The avoidance of excessive mesh distortion, encountered with strictly Lagrangian or Eulerian kinematics, was achieved by introducing an auxiliary kinematic velocity field along the free surface in order to vary the trajectories used in integrating the ordinary differential equations simulating the moving boundary. The computation of the velocity potential was based upon a nonuniform triangular mesh which was automatically revised to varying depths to accommodate the motion of the free surface. These methods permitted calculation of draining induced axisymmetric slosh through the many (or fractional) finite amplitude oscillations that can occur depending upon the balance of draining, gravitational, and surface tension forces. Velocity fields, evolution of the free surface with time, and liquid residual volumes were computed for three and one half decades of Weber number and for two Bond numbers, tank fill levels, and drain radii. Comparisons with experimental data are very satisfactory.

Bizzell, G. D.↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗
NASA NTRS1990

Numerical simulation of the compressible Orszag-Tang vortex 2. Supersonic flow

The numerical investigation of the Orszag-Tang vortex system in compressible magnetofluids will consider initial conditions with embedded supersonic regions. The simulations have initial average Mach numbers 1.0 and 1.5 and beta 10/3 with Lundquist numbers 50, 100, or 200. The behavior of the system differs significantly from that found previously for the incompressible and subsonic analogs. Shocks form at the downstream boundaries of the embedded supersonic regions outside the central magnetic X-point and produce strong local current sheets which dissipate appreciable magnetic energy. Reconnection at the central X-point, which dominates the incompressible and subsonic systems, peaks later and has a smaller impact as M increases from 0.6 to 1.5. Similarly, correlation between the momentum and magnetic field begins significant growth later than in subsonic and incompressible flows. The shocks bound large compression regions, which dominate the wavenumber spectra of autocorrelations in mass density, velocity, and magnetic field.

Picone, J. M.↗
DOE OSTI2024

Bulk cloud microphysical properties as seen from numerical simulation and remote sensing products: case study of a hailstorm event over the La Plata Basin

Hailstorms develop over the La Plata Basin, in south-eastern South America, more often during later winter and early austral spring, between September and October. These systems have significant socioeconomic impacts over the region. Thus, a better understanding of how atmospheric drivers modulate the formation of hailstorms is important to improve the forecast of such phenomena. In this study, we selected a hailstorm event observed over the eastern La Plata Basin during 14–15 July 2016 to evaluate the performance of the Brazilian developments on the Regional Atmospheric Modelling System (BRAMS) model. The ability of the model in simulating cloud microphysical properties was evaluated by comparing simulations driven by different global forcings against in situ and remote sensing observations. The model results showed good skill in capturing the basic characteristics of the thunderstorm, particularly in terms of the spatial distribution of hydrometeors. The simulated spatial distribution of hail covers locations where hail fall was reported. The BRAMS simulations suggest that, despite relatively low values of the convective available potential energy (CAPE) (700–1000 J kg -1 ), environments with strong 0–8-km bulk shear (60–70 kt, ~30.9–36.0 m s –1 ) can promote the formation of ice clouds and hail fall over the eastern La Plata Basin. To be more conclusive, however, further research is needed to understand how different combinations of CAPE and shear affect hail formation over the region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1978

A numerical simulation of seasonal stratospheric climate. I - Zonal temperatures and winds

Two numerical experiments were designed to isolate the effects of seasonal and latitudinal variations in longwave radiative transfer processes. The first experiment aimed at computing the rate of temperature change in the stratosphere due to longwave radiative transfer by CO2, H2O and O3 from Ramanathan's (1976) radiative-convective model applied to the troposphere-stratosphere circulation. The second experiment employs the Newtonian cooling approximation in which the rate of temperature change in the stratosphere by longwave radiative transfer is set equal to the product of the Newtonian cooling coefficient 'h' and the departure of the local temperature from a reference temperature. It is shown that the latitudinal temperature distribution of the lower stratosphere is maintained by the combined effects of dynamics, O3 solar heating and the longwave radiative coupling between troposphere and stratosphere. The latitudinal gradient in the troposphere-stratosphere longwave radiative coupling is maximum during winter and spring. The radiative response time (1/h) of the middle and upper stratosphere undergoes significant latitudinal and seasonal variations, largely due to the temperature dependence of h.

Ramanathan, V.↗
NASA NTRS1984

A numerical simulation of barotropic instability. I Wave-mean flow interaction

A numerical model is used to study the evolution of the barotropic point jet instability as it interacts with the mean flow. The linearized instability solution agrees well with the recent analytical solutions of Lindzen. Stabilization of the point jet instability occurs as the mean flow is modified by wave vorticity transport. Assuming stabilization occurs when the meridional gradient of the zonal mean vorticity is no longer negative, the maximum integrated wave enstrophy can be predicted. In addition, an estimate of the integrated wave enstrophy at steady state can be made by balancing the generation of vorticity against dissipation. These limits are found to be in good agreement with the numerical results.

Schoeberl, M. R.↗
NASA NTRS1986

Numerical simulation of synthetic aperture radar image spectra for ocean waves

A numerical model for predicting the synthetic aperture radar (SAR) image of a moving ocean surface is described, and results are presented for two SIR-B data sets collected off the coast of Chile. Wave height spectra measured by the NASA radar ocean wave spectrometer (ROWS) were used as inputs to this model, and results are compared with actual SIR-B image spectra from orbits 91 and 106. Additional parametric variations are presented to illustrate the effects of nonlinearities in the imaging process.

Lyzenga, D. R.↗
NASA NTRS1987

Investigation of turbulent shear flows by numerical simulation

A data base (Kim et al., 1986) generated by direct numerical integration of the Navier-Stokes equations for unsteady three-dimensional turbulent channel flow at low Reynolds numbers is applied to investigate the physical processes acting in wall-bounded shear flows. The techniques used to construct the data base are outlined; the modeling of turbulence, streaks and streamwise vortices, horseshoe and inverted-horseshoe vortices, and vortex rings is described; the conditional sampling techniques employed are discussed; and typical results are presented graphically. An organized vorticity-field structure is revealed and found to be associated with the bursting process and a local peak in the wall pressure, indicating a possible approach to reducing viscous drag.

Kim, J.↗
NASA NTRS1988

Numerical simulation of unsteady incompressible viscous flows in generalized coordinate systems

Several numerical solutions of the three dimensional unsteady incompressible Navier-Stokes equations in generalized coordinate systems are presented. The governing equations are discretized by finite volumes with special care to the accurate approximation of the geometric quantities. The unknowns are the pressure and the volume fluxes over the computational cell faces. This formulation results in a robust fractional step solution method for solving discrete equations. Although this method is formulated for the three dimensional case, only two dimensional unsteady results are given. Results are presented for a lid driven two dimensional cavity flow at Reynolds number of 10,000, and for the flow over a circular cylinder with vortex shedding for several Reynolds numbers in the range 100 less than Re less than 1000.

Rosenfeld, Moshe↗
NASA NTRS1989

Numerical simulation of unsteady incompressible viscous flows in generalized coordinate systems

Several numerical solutions of the three-dimensional unsteady incompressible Navier-Stokes equations in generalized coordinate systems are presented. The governing equations are discretized by finite volumes with special care to the accurate approximation of the geometric quantities. The unknowns are the pressure and the volume fluxes over the computational cell faces. This formulation results in a robust fractional step solution method for solving discrete equations. Although this method is formulated for the three-dimensional case, only two-dimensional unsteady results are given. Results are presented for a lid driven two-dimensional cavity flow at Reynolds number of 10,000 and for the flow over a circular cylinder with vortex shedding for several Reynolds numbers in the range 100 less than Re less than 1000.

Rossenfeld, Moshe↗
NASA NTRS1990

Numerical simulation of 3-D Benard convection with gravitational modulation

In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.

Biringen, S.↗
NASA NTRS1990

Numerical simulations of three-dimensional rotor blade-vortex interactions

A numerical procedure which has been developed (Hassan and Charles, 1989) for the prediction of helicopter blade loads during self-generated blade-vortex interactions (BVI) is described. These BVIs result in significant and frequent impulsive changes in the rotor blade aerodynamic loads and moments that are known to lead to vibrations and severe noise problems. Potential BVIs are identified and tracked in time at even increments of rotor azimuth employing the lifting line helicopter/rotor trim code CAMRAD. The accuracy of the predicted leading edge pressures for self-generated interactions is shown to rely heavily on the user-specified vortex core radius, and more significantly on the CAMRAD-predicted geometry of the interaction vortex elements and their relative orientation with respect to the blade. Thus, the robustness of the surface transpiration approach for modeling subcritical-type BVIs is proven.

Hassan, A. A.↗
NASA NTRS1992

Numerical simulation of vortex unsteadiness on a slender body at high incidence

Numerical results obtained with the present3D thin-layer Navier-Stokes solver are shown to qualitatively follow experimental results. They also indicate the existence of at least two main frequencies on an ogive cylinder at 40 and 60 deg angles-of-attack: (1) a low frequency associated with fluctuations of the primary vortices, and leads to vortex shedding, and (2) a high frequency associated with shear-layer fluctuations that are concentrated in several small areas above the leeward side of the body's cylindrical portion.

Degani, David↗
DOE OSTI2025

Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems

Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Numerical simulation of hypersonic inlet flows with equilibrium or finite rate chemistry

An efficient numerical program incorporated with comprehensive high temperature gas property models has been developed to simulate hypersonic inlet flows. The computer program employs an implicit lower-upper time marching scheme to solve the two-dimensional Navier-Stokes equations with variable thermodynamic and transport properties. Both finite-rate and local-equilibrium approaches are adopted in the chemical reaction model for dissociation and ionization of the inlet air. In the finite rate approach, eleven species equations coupled with fluid dynamic equations are solved simultaneously. In the local-equilibrium approach, instead of solving species equations, an efficient chemical equilibrium package has been developed and incorporated into the flow code to obtain chemical compositions directly. Gas properties for the reaction products species are calculated by methods of statistical mechanics and fit to a polynomial form for C(p). In the present study, since the chemical reaction time is comparable to the flow residence time, the local-equilibrium model underpredicts the temperature in the shock layer. Significant differences of predicted chemical compositions in shock layer between finite rate and local-equilibrium approaches have been observed.

Yu, Sheng-Tao↗
NASA NTRS1980

Full-coverage film cooling. II - Heat transfer data and numerical simulation

Experimental research into heat transfer from full-coverage film-cooled surfaces with three injection geometries was described in Part I. This part has two objectives. The first is to present a simple numerical procedure for simulation of heat transfer with full-coverage film cooling. The second objective is to present some of the Stanton number data that was used in Part I of the paper. The data chosen for presentation are the low-Reynolds number, heated-starting-length data for the three injection geometries with five-diameter hole spacing. Sample data sets with high blowing ratio and with ten-diameter hole spacing are also presented. The numerical procedure has been successfully applied to the Stanton number data sets.

Crawford, M. E.↗
NASA NTRS2002

Numerical simulation of the influence of gravity and posture on cardiac performance

A numerical model of the cardiovascular system was used to quantify the influences on cardiac function of intrathoracic pressure and intravascular and intraventricular hydrostatic pressure, which are fundamental biomechanical stimuli for orthostatic response. The model included a detailed arterial circulation with lumped parameter models of the atria, ventricles, pulmonary circulation, and venous circulation. The venous circulation was divided into cranial, central, and caudal regions with nonlinear compliance. Changes in intrathoracic pressure and the effects of hydrostatic pressure were simulated in supine, launch, sitting, and standing postures for 0, 1, and 1.8 G. Increasing intrathoracic pressure experienced with increasing gravity caused 12% and 14% decreases in cardiac output for 1 and 1.8 G supine, respectively, compared to 0 G. Similar results were obtained for launch posture, in which the effects of changing intrathoracic pressure dominated those of hydrostatic pressure. Compared to 0 G, cardiac output decreased 0.9% for 1 G launch and 15% for 1.8 G launch. In sitting and standing, the position of the heart above the hydrostatic indifference level caused the effects of changing hydrostatic pressure to dominate those of intrathoracic pressure. Compared to 0 G, cardiac output decreased 13% for 1 G sitting and 23% for 1.8 G sitting, and decreased 17% for 1 G standing and 31% for 1.8 G standing. For a posture change from supine to standing in 1 G, cardiac output decreased, consistent with the trend necessary to explain orthostatic intolerance in some astronauts during postflight stand tests. Simulated lower body negative pressure (LBNP) in 0 G reduced cardiac output and mean aortic pressure similar to I G standing, suggesting that LBNP provides at least some cardiovascular stimuli that may be useful in preventing postflight orthostatic intolerance. A unifying concept, consistent with the Frank-Starling mechanism of the heart, was that cardiac output was proportional to cardiac diastolic transmural pressure for all postures and gravitational accelerations.

Non-NASA Center↗
NASA NTRS1984

Numerical simulation of reconnection in an emerging magnetic flux region

The resistive MHD equations are numerically solved in two dimensions for an initial-boundary-value problem which simulates reconnection between an emerging magnetic flux region and an overlying coronal magnetic field. The solution involves both ideal-MHD and resistive-MHD processes, and the solution shows an evolution which is remarkably suggestive of the preflare, impulsive, and main phases of the flare-cycle.

Forbes, T. C.↗