Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow 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 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

An Adaptive Semi-Implicit Scheme for Simulations of Unsteady Viscous Compressible Flows

A numerical scheme for simulation of unsteady, viscous, compressible flows is considered. The scheme employs an explicit discretization of the inviscid terms of the Navier-Stokes equations and an implicit discretization of the viscous terms. The discretization is second order accurate in both space and time. Under appropriate assumptions, the implicit system of equations can be decoupled into two linear systems of reduced rank. These are solved efficiently using a Gauss-Seidel method with multigrid convergence acceleration. When coupled with a solution-adaptive mesh refinement technique, the hybrid explicit-implicit scheme provides an effective methodology for accurate simulations of unsteady viscous flows. The methodology is demonstrated for both body-fitted structured grids and for rectangular (Cartesian) grids.

Steinthorsson, Erlendur↗
NASA NTRS1984

Finite difference simulation of unsteady interactive flows

The development of time-dependent numerical simulations of unsteady interactive flows of an aerodynamic nature is reviewed with emphasis on compressible flows at flight Reynolds numbers and noniterative schemes based on Navier-Stokes equations. The importance of writing the equations in strong conservation-law form for a generalized body-oriented coordinate system is pointed out. The discussion covers time and length scales and numerical methods currently in use. Some computed results are presented and compared with experimental data.

Deiwert, G. S.↗
NASA NTRSDate not supplied

Flow and Acoustics Simulations of a Supersonic Engine Inlet

An inlet for a supersonic aircraft engine was previously studied in a geometrically simplified form. The present paper continues that work and specifically investigates three-dimensional effects on the combined inlet and fan system. Three-dimensional flow, inlet components and viscous effects are studied for their effect on fan response and acoustics. Operating condition was identified as the most important parameter affecting noise, which is discussed in a companion paper. The details of the inlet geometry presented here are a smaller effect. Additionally, the results indicate that an axisymmetric two-dimensional computation is sufficiently accurate for some cases while a three-dimensional computation may reveal particular details.

Supersonics↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
NASA NTRS1981

A computer program for the simulation of heat and moisture flow in soils

A computer program that simulates the flow of heat and moisture in soils is described. The space-time dependence of temperature and moisture content is described by a set of diffusion-type partial differential equations. The simulator uses a predictor/corrector to numerically integrate them, giving wetness and temperature profiles as a function of time. The simulator was used to generate solutions to diffusion-type partial differential equations for which analytical solutions are known. These equations include both constant and variable diffusivities, and both flux and constant concentration boundary conditions. In all cases, the simulated and analytic solutions agreed to within the error bounds which were imposed on the integrator. Simulations of heat and moisture flow under actual field conditions were also performed. Ground truth data were used for the boundary conditions and soil transport properties. The qualitative agreement between simulated and measured profiles is an indication that the model equations are reasonably accurate representations of the physical processes involved.

Camillo, P.↗
NASA NTRS1987

On direct numerical simulations of turbulent reacting flows

A description of the emerging field of direct numerical simulations of turbulent, chemically reacting flows is presented. The types of direct numerical simulations, physical issues related to implementing the simulations, as well as the various numerical methods used are described. Examples are presented of recent applications of direct numerical simulations to a variety of problems, displaying both the potential of the method and also some of its limitations. Finally, our view of the potential role of direct numerical simulations in future research on turbulent, chemically reacting flows is presented.

Jou, W.-H.↗
NASA NTRS1996

Simulations of Turbine Cooling Flows Using a Multiblock-Multigrid Scheme

Results from numerical simulations of air flow and heat transfer in a 'branched duct' geometry are presented. The geometry contains features, including pins and a partition, as are found in coolant passages of turbine blades. The simulations were performed using a multi-block structured grid system and a finite volume discretization of the governing equations (the compressible Navier-Stokes equations). The effects of turbulence on the mean flow and heat transfer were modeled using the Baldwin-Lomax turbulence model. The computed results are compared to experimental data. It was found that the extent of some regions of high heat transfer was somewhat under predicted. It is conjectured that the underlying reason is the local nature of the turbulence model which cannot account for upstream influence on the turbulence field. In general, however, the comparison with the experimental data is favorable.

Steinthorsson, Erlendur↗
DOE OSTI2019

MULTIDIMENSIONAL NUMERICAL SIMULATIONS OF REACTING FLOW IN A NON-PREMIXED ROTATING DETONATION ENGINE

Over the last two decades, detonation based propulsion has received a great deal of attention as a potential means to achieve significant improvement in the performance of air-breathing and rocket engines. Detonative combustion mode is particularly interesting due to the resulting pressure gain from reactants to products, faster heat release, decreased entropy generation, more available work and higher thrust compared to conventional deflagrative combustion. Rotating detonation engine (RDE) is one such novel combustor concept. Realistic RDE configurations utilize separate fuel and air injection schemes, hence are not perfectly premixed. Moreover, RDE performance is governed by a large number of design parameters and operating conditions. In this context, computational fluid dynamics (CFD) has the potential to enhance the understanding of RDE combustion and aid future development/optimization of this technology. In the present work, a CFD model was developed to simulate a representative non-premixed RDE combustor. Unsteady Reynolds-Averaged Navier-Stokes (RANS) simulations were performed for the full combustor geometry (including the separate fuel and air injection ports), with hydrogen as fuel and air as the oxidizer. Adaptive mesh refinement (AMR) was incorporated to achieve a trade-off between model accuracy and computational expense. A finite-rate chemistry model along with a 10-species detailed kinetic mechanism was employed to describe the H2-Air combustion chemistry. Two operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel and air mass flow rates. For both conditions, the capability of the model to capture the essential detonation wave dynamics was assessed. A validation study was performed against experimental data available on detonation wave frequency/height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model predicted the sensitivity of these wave characteristics to the operating conditions with good accuracy, both qualitatively and quantitatively. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

Pal, Pinaki↗
DOE OSTI2019

High-order two-fluid plasma solver for direct numerical simulations of plasma flows with full transport phenomena

The two-fluid plasma equations for a single ion species, with full transport terms, including temperature and magnetic field dependent ion and electron viscous stresses and heat fluxes, frictional drag force, and ohmic heating terms, have been implemented in the CFDNS code and solved by using sixth-order non-dissipative compact finite differences for plasma flows in several different regimes. In order to be able to fully resolve all the dynamically relevant time and length scales, while maintaining computational feasibility, the assumptions of infinite speed of light and negligible electron inertia have been made. Non-dimensional analysis of the two-fluid plasma equations shows that, by varying the characteristic/background number density, length scale, temperature, and magnetic strength, the corresponding Hall, resistive, and ideal magnetohydrodynamic equations can be recovered as limiting cases. The accuracy and robustness of this two-fluid plasma solver in handling plasma flows in different regimes have been validated against four canonical problems: Alfven and whistler dispersion relations, electromagnetic plasma shock, and magnetic reconnection. For all test cases, by using physical dissipation and diffusion, with negligible numerical dissipation/diffusion, fully converged Direct Numerical Simulation (DNS)-like solutions are obtained when the ion Reynolds number based on the grid size is smaller than a threshold value which is about 2.3. For the magnetic reconnection problem, the results show that the magnetic flux saturation time and value converge when the ion and magnetic Reynolds numbers are large enough. Thus, the DNS-like results become relevant to practical problems with much larger Reynolds numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Simulation of fluid flows during growth of organic crystals in microgravity

Several counter diffusion type crystal growth experiments were conducted in space. Improvements in crystal size and quality are attributed to reduced natural convection in the microgravity environment. One series of experiments called DMOS (Diffusive Mixing of Organic Solutions) was designed and conducted by researchers at the 3M Corporation and flown by NASA on the space shuttle. Since only limited information about the mixing process is available from the space experiments, a series of ground based experiments was conducted to further investigate the fluid dynamics within the DMOS crystal growth cell. Solutions with density differences in the range of 10 to the -7 to 10 to the -4 power g/cc were used to simulate microgravity conditions. The small density differences were obtained by mixing D2O and H2O. Methylene blue dye was used to enhance flow visualization. The extent of mixing was measured photometrically using the 662 nm absorbance peak of the dye. Results indicate that extensive mixing by natural convection can occur even under microgravity conditions. This is qualitatively consistent with results of a simple scaling analysis. Quantitave results are in close agreement with ongoing computational modeling analysis.

Roberts, Gary D.↗
NASA NTRS1987

Numerical simulation of the flow field and fuel sprays in an IC engine

A two-dimensional model for axisymmetric piston-cylinder configurations is developed to study the flow field in two-stroke direct-injection Diesel engines under motored conditions. The model accounts for turbulence by a two-equation model for the turbulence kinetic energy and its rate of dissipation. A discrete droplet model is used to simulate the fuel spray, and the effects of the gas phase turbulence on the droplets is considered. It is shown that a fluctuating velocity can be added to the mean droplet velocity every time step if the step is small enough. Good agreement with experimental data is found for a range of ambient pressures in Diesel engine-type microenvironments. The effects of the intake swirl angle in the spray penetration, vaporization, and mixing in a uniflow-scavenged two-stroke Diesel engine are analyzed. It is found that the swirl increases the gas phase turbulence levels and the rates of vaporization.

Nguyen, H. L.↗
NASA NTRS1989

Multifractal spectra in shear flows

Numerical simulations of three-dimensional homogeneous shear flow and fully developed channel flow, are used to calculate the associated multifractal spectra of the energy dissipation field. Only weak parameterization of the results with the nondimensional shear is found, and this only if the flow has reached its asymptotic development state. Multifractal spectra of these flows coincide with those from experiments only at the range alpha less than 1.

Keefe, L. R.↗
DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Simulation and Modeling of the Elliptic Streamline Flow

Direct numerical simulation (DNS) was used to study homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. There are many important reasons for studying the elliptic streamline flow. This flow contains the effects of both rotation and strain and is therefore similar to the mean flow in a vortex strained in the plane perpendicular to its axis. Such flows provide insight into the fundamental vortical interactions within turbulence. A strained vortex also occurs in airplane wakes, in which each wingtip vortex induces a strain field on the other. The strain field can affect the stability of these vortices and thereby their turbulent structure downstream. The ability to understand and predict the turbulent structure of the vortices is important to the wake hazard problem which is of major concern for the safety of commercial aircraft.

Blaisdell, Gregory A.↗
DOE OSTI2016

Numerical simulation of break flows occurring in the reactor cavity cooling system of a very high temperature reactor using a system CFD approach - 18437

To avoid accidents due to human error or the failure of components, the reactor cavity cooling system (RCCS) of a high temperature reactor must be able to operate naturally or passively safe. The system under consideration in this investigation is an air-cooled system which consists of an annular downcomer which is connected via a bottom-header to riser tubes. The downcomer and riser tubes are connected to the atmosphere via inlet and outlet manifold systems respectively. In this study the entire system is simulated using the 1D network codes GAMMA+ and Flownex. The study investigates the effect that breaks, causing a short circuit between one half of the inlet manifold and the corresponding half of the outlet manifold, will have on the performance of the cooling system. The agreement between the results obtained with the two codes is good. When the break size exceeds a certain value the flow may reverse in the relevant half of the outlet manifold. The critical break size is dependent on the elevation at which the break occurs. The reversal in the flow results in a drop of 10% in the heat removal capacity of the RCCS. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Modeling and simulation of multiphase flows

This presentation provides an overview of the National Energy Technology Laboratory’s (NETL) multiphase computational fluid dynamics codes. The highly successful Multiphase Flows with Interphase eXchanges (MFIX) suite has been used to model a wide range of applications including post-combustion carbon capture, bioreactor optimization, and bio-FCC regeneration. MFIX-Exa, a state-of-the-art CFD code, developed under DOE’s Exascale Computing Project, is built on the AMReX software framework (https://amrex-codes.github.io/) and is designed to leverage modern accelerator-based compute architectures. This presentation further reviews the underlying physical models of both MFIX and MFIX-Exa and contrasts their similarities and differences. Examples of past and present CFD simulations will illustrate how scientific computing at NETL is being used not only for scientific exploration but also for design, optimization and scale-up of multiphase flow devices.

Musser, Jordan [NETL]↗
NASA NTRS2015

Detached Eddy Simulation for the F-16XL Aircraft Configuration

Numerical simulations for the flow around the F-16XL configuration as a contribution to the Cranked Arrow Wing Aerodynamic Project International 2 (CAWAPI-2) have been performed. The NASA Langley Tetrahedral Unstructured Software System (TetrUSS) with its USM3D solver was used to perform the unsteady flow field simulations for the subsonic high angle-of-attack case corresponding to flight condition (FC) 25. Two approaches were utilized to capture the unsteady vortex flow over the wing of the F-16XL. The first approach was to use Unsteady Reynolds-Averaged Navier-Stokes (URANS) coupled with standard turbulence closure models. The second approach was to use Detached Eddy Simulation (DES), which creates a hybrid model that attempts to combine the most favorable elements of URANS models and Large Eddy Simulation (LES). Computed surface static pressure profiles are presented and compared with flight data. Time-averaged and instantaneous results obtained on coarse, medium and fine grids are compared with the flight data. The intent of this study is to demonstrate that the DES module within the USM3D solver can be used to provide valuable data in predicting vortex-flow physics on a complex configuration.

Elmiligui, Alaa↗