Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Numerical simulation and optimization of red mud separation thickener with self-dilute feed

In order to acquire the flow pattern and investigate the settling behavior of the red mud in the separation thickener, computational fluid dynamics (CFD), custom subroutines and agglomerates settling theory were employed to simulate the three-dimensional flow field in an industrial scale thickener with the introduction of a self-dilute feed system. Our simulation results show good agreement with the measurement onsite and the flow patterns of the thickener are presented and discussed on both velocity and concentration field. Optimization experiments on feed well and self-dilute system were also carried out, and indicate that the optimal thickener system can dilute the solid concentration in feed well from 110 g/L to 86 g/L which would help the agglomerates’ formation and improve the red mud settling speed. The additional power of recirculation pump can be saved and flocculants dosage was reduced from 105g/t to 85g/t in the operation.

42 ENGINEERING↗
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↗
DOE OSTI2019

Numerical Simulations of the Supercritical Carbon Dioxide Round Turbulent Jet

We explore the fundamental connection between properties of a supercritical fluid and observed behavior of the flow by comparing simulations of a supercritical carbon dioxide (sCO2) round jet to canonical simulations using an ideal gas model. sCO2 has desirable behavior that can improve power density compared to traditional fluids for applications involving closed-cycle gas turbines, heat transmission, and hydraulic fracturing while facilitating carbon sequestration. While literature exists focusing on the prospects of sCO2 in various industries, the quantities of interest studied tend to be application specific whereas our investigation generalizes across multiple applications. The Soave-Redlich-Kwong equation of state is utilized to close our system of equations; we use a second order finite volume method in conjunction with adaptive mesh refinement as implemented in PeleC. The jet is at p=10 MPa and T=600 K in order to maintain a single-phase fluid. Quantities of interest for this study include the mean axial velocity and Reynolds stresses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Numerical simulation of seawater intrusion to coastal aquifers and brine water/freshwater interaction in south coast of Laizhou Bay, China

Seawater intrusion and brine water/freshwater interaction have significantly affected agriculture, industry and public water supply at Laizhou Bay, Shandong Province, China. A two-dimensional SEAWAT model is developed to simulate the seawater intrusion to coastal aquifers and brine water/fresh water interaction in the south of Laizhou Bay. This model is applied to predict the seawater intrusion and brine water/freshwater interface development in the coming years. The model profile is perpendicular to the coastal line with two interfaces, freshwater-saline water interface near the shore and inland brine water-saline water-seawater interface. The hydrogeological parameters in the SEAWAT-2000 model are calibrated by the head and salinity measurements. The precipitation infiltration coefficient, boundary conditions and thicknesses of aquifers are studied in a sensitivity analysis. The predicted results indicate that equivalent freshwater head in shallow freshwater-saline water area will decline 2.0 m by the end of the forecasting period, caused by groundwater over-pumping for farmland irrigation. The groundwater head in the brine-saline water area will also decrease about 1.8 m by the end of forecasting period, caused by excessive brine mining. Salinity decreases below 105 g/L in the brine area, but increases in other areas and contaminates fresh groundwater resources.

54 ENVIRONMENTAL SCIENCES↗
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.↗
DOE OSTI2015

Numerical simulations of vertical growth of hydraulic fractures and brine migration in geological formations above the Marcellus shale

One of the critical environmental questions about hydraulic fracturing in shales is the potential for contamination of ground and surface water. There are two specific concerns arising from hydraulic treatments: 1) whether hydraulic fractures extend upward through overlying strata to reach overlying aquifers containing drinking water, and 2) whether injected fluids push native fluids upward into these overlying aquifers. In this work, the extent of likely fracture growth through overlying layers during hydraulic treatment of the Marcellus shale was estimated using a hydraulic fracture model. A wide range of material and fluid flow properties in a multi-layered geologic model was considered. The model was based on conditions and characteristics applicable to the Marcellus shale in that part of the Appalachian basin within southwestern Pennsylvania. Predictions of vertical termination frequencies for hydraulic fractures were used in a multi-layer model of the strata and natural fractures for studying brine migration through the natural and induced fracture network. NFFLOW, the software for explicitly modeling flow within networks of fractures, was utilized to compute transient flow rates according to the schedule of injected fluid during hydraulic fracturing. To aid our analysis, the modeled sequence of geologic strata was capped with a fictitious unfractured, but moderately-permeable layer, which serves as a monitoring zone. The analysis assumes one well lateral was placed in the middle of the Marcellus shale with hydraulic fractures penetrating layers in the model. The newly-developed geomechanical module within NFFLOW was used to represent stress-sensitivity of the fractures. This allows the opening and closing of fracture apertures with changes in fluid pressures within fracture segments. Pressure increases in the formations overlying the Tully limestone, indicating fluid flow, was observed due to the hydraulic stimulation; and the impact of these increased pressures on brine migration towards the surface was considered.

03 NATURAL GAS↗