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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Simulation of transient flow in a shock tunnel and a high Mach number nozzle

A finite volume Navier-Stokes code was used to simulate the shock reflection and nozzle starting processes in an axisymmetric shock tube and a high Mach number nozzle. The simulated nozzle starting processes were found to match the classical quasi-1-D theory and some features of the experimental measurements. The shock reflection simulation illustrated a new mechanism for the driver gas contamination of the stagnated test gas.

Jacobs, P. A.↗
NASA NTRS1986

Tip vortices of wings in subsonic and transonic flow: A numerical simulation

Thin layer Navier-Stokes and Euler equations are numerically solved using a multi-block zonal approach to simulate the formation and roll up of tip vortices of wings in subsonic and transonic flows. Several wing planforms were considered to examine the influence of tip-cap shape, planform geometry and free stream Mach number on the formation process. A good definition of the formation and qualitative roll up of tip vortices was achieved.

Srinivasan, G. R.↗
NASA NTRS1992

Evaluation of thermochemical models for particle and continuum simulations of hypersonic flow

Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered that are dominated in turn by vibrational relaxation, weak dissociation, strong dissociation and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.

Boyd, Iain D.↗
NASA NTRS1996

Simulation and modeling of the elliptic streamline flow

Direct numerical simulations are performed for the elliptic streamline flow, which is a homogeneous turbulent flow that combines the effects of solid body rotation and strain. Simulations are run over a range of parameters in order to determine the effect of changing rotation and strain separately. For early times the nonlinear cascade is suppressed, but then is re-established at later times. The growth rate of turbulent kinetic energy agrees at early times with the trends from linear theory, but at later times the flow seems to approach an asymptotic state that is independent of the ratio of mean flow rotation rate to strain rate. A comparison with standard Reynolds stress turbulence models is made. It is found that for strong rotation rates, the models predict decay of the turbulence, while the simulations show exponential growth. Close examination of the simulation results shows that they are affected by excessively low Reynolds numbers. Suggestions for reducing low Reynolds number effects in future simulations is given.

Blaisdell, G. A.↗
DOE OSTI2016

A High-Fidelity Approach for the Simulation of Flow-Induced Vibration

Flow-induced vibration (FIV) is a widespread problem in energy systems becausethey rely on fluid movement for energy conversion. Vibrating structures may be damaged as fatigue or wear occur. Given the importance of reliable components in the nuclear industry, flow-induced vibrations have long been a major concern in the safety and operation of nuclear reactors. In particular,nuclear fuel and steam generators have been known to suffer from flow-induced vibrations and related failures. Over the past five years, the Nuclear EnergyAdvanced Modeling and Simulation program has developed the integrated multiphysics code suite SHARP. The goal of developing such a tool is to perform multiphysics modeling of the components inside a reactor core, the full reactor core or portions of it, and be able to achieve that with various levelsof fidelity. This flexibility allows users to select the appropriate level of fidelity for their computational resources and design constraints.In particular SHARP contains high-fidelity single-physics codes for structural mechanics and fluid mechanics calculations: the structural mechanics implicit code Diablo and the computational fluid dynamics spectral element code Nek5000. Both codes are state-of-the-art.highly scalable (up to millions of processors in the case of Nek5000) tools that have been extensively validated. These tools form a strong basis on which to build an FIV modeling capability.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1984

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗
DOE OSTI2018

Finite-Element Simulations of Elastoplastic Flow during Compression of a Sample in a Diamond Anvil Cell under Extremely High Pressure: Effects of Geometry and Material Properties

Finite-element simulations are conducted to investigate large elastoplastic deformations of rhenium under multimegabar pressures in a diamond anvil cell (DAC), with an emphasis on the effects of geometric and material properties. The following published experimental phenomena are reproduced: (1) the pressure distribution at the sample or diamond contact surface and the final sample thickness at pressure up to 300 GPa; (2) the cupping (i.e., appearance of a cuplike concave shape of the contact diamond-sample surface) and double cupping phenomena at megabar pressures; (3) three stages at the curve of the maximum pressure versus compressive force; (4) stages of material flow with increasing load; (5) pressure drop at the periphery after cupping in that region; and (6) change in direction of material flow to the center without change in the sign of the pressure gradient. The effects of the culet geometry, bevel angle, sample thickness, and sample or gasket system are analyzed in detail. The obtained results improve the understanding of the material mechanical response under extreme pressures and large elastoplastic deformations and are beneficial for the optimum design of a DAC with the goal of reaching the high and record high pressures once or multiple times.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

Large-Eddy Simulation of Conductive Flows at Low Magnetic Reynolds Number

In this paper we study the LES method with dynamic procedure in the context of conductive flows subject to an applied external magnetic field at low magnetic Reynolds number R(sub m). These kind of flows are encountered in many industrial applications. For example, in the steel industry, applied magnetic fields can be used to damp turbulence in the casting process. In nuclear fusion devices (Tokamaks), liquid-lithium flows are used as coolant blankets and interact with the surrounding magnetic field that drives and confines the fusion plasma. Also, in experimental facilities investigating the dynamo effect, the flow consists of liquid-sodium for which the Prandtl number and, as a consequence, the magnetic Reynolds number is low. Our attention is focused here on the case of homogeneous (initially isotropic) decaying turbulence. The numerical simulations performed mimic the thought experiment described in Moffatt in which an initially homogeneous isotropic conductive flow is suddenly subjected to an applied magnetic field and freely decays without any forcing. Note that this flow was first studied numerically by Schumann. It is well known that in that case, extra damping of turbulence occurs due to the Joule effect and that the flow tends to become progressively independent of the coordinate along the direction of the magnetic field. Our comparison of filtered direct numerical simulation (DNS) predictions and LES predictions show that the dynamic Smagorinsky model enables one to capture successfully the flow with LES, and that it automatically incorporates the effect of the magnetic field on the turbulence. Our paper is organized as follows. In the next section we summarize the LES approach in the case of MHD turbulence at low R(sub m) and recall the definition of the dynamic Smagorinsky model. In Sec. 3 we describe the parameters of the numerical experiments performed and the code used. Section 4 is devoted to the comparison of filtered DNS results and LES results. Conclusions are presented in Sec. 5.

Knaepen, B.↗
DOE OSTI2017

Final Report for ALCC Allocation: Predictive Simulation of Complex Flow in Wind Farms

This report documents work performed using ALCC computing resources granted under a proposal submitted in February 2016, with the resource allocation period spanning the period July 2016 through June 2017. The award allocation was 10.7 million processor-hours at the National Energy Research Scientific Computing Center. The simulations performed were in support of two projects: the Atmosphere to Electrons (A2e) project, supported by the DOE EERE office; and the Exascale Computing Project (ECP), supported by the DOE Office of Science. The project team for both efforts consists of staff scientists and postdocs from Sandia National Laboratories and the National Renewable Energy Laboratory. At the heart of these projects is the open-source computational-fluid-dynamics (CFD) code, Nalu. Nalu solves the low-Mach-number Navier-Stokes equations using an unstructured- grid discretization. Nalu leverages the open-source Trilinos solver library and the Sierra Toolkit (STK) for parallelization and I/O. This report documents baseline computational performance of the Nalu code on problems of direct relevance to the wind plant physics application - namely, Large Eddy Simulation (LES) of an atmospheric boundary layer (ABL) flow and wall-modeled LES of a flow past a static wind turbine rotor blade. Parallel performance of Nalu and its constituent solver routines residing in the Trilinos library has been assessed previously under various campaigns. However, both Nalu and Trilinos have been, and remain, in active development and resources have not been available previously to rigorously track code performance over time. With the initiation of the ECP, it is important to establish and document baseline code performance on the problems of interest. This will allow the project team to identify and target any deficiencies in performance, as well as highlight any performance bottlenecks as we exercise the code on a greater variety of platforms and at larger scales. The current study is rather modest in scale, examining performance on problem sizes of O(100 million) elements and core counts up to 8k cores. This will be expanded as more computational resources become available to the projects.

17 WIND ENERGY↗
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.↗
DOE OSTI2022

Multiscale Experiments and Multiphysics Simulation of Multiphase Flow for Transportable Small Modular Reactors

A new type of safe, small, transportable nuclear reactor would address the intense and ever-growing global demand for energy produced via a resilient, carbon-free energy source. In this regard, transportable small modular reactors (SMRs) are being designed and developed for electricity generation within small/micro-grid/off-grid isolated systems, as well as for heat generation in industrial/residential applications. These reactors feature the capability to be fully factory fabricated and then directly transported to utilities’ sites as “plug-and-play” systems. Research and development (R&D) programs are underway at Idaho National Laboratory (INL) to successfully design, develop, and demonstrate such safe-by-design mobile reactor technologies, in collaboration with partner organizations. Multiscale experimental facilities and multiphysics simulation tools are required for reactor design verification and validation (V&V), and licensing. These advanced reactors are intended to feature passive safety systems such as passive containment cooling systems (PCCS), which consist of multiphase flows and multispecies distributions. This seminar talk will focus on designing and analyzing transportable SMR PCCS by using multiscale experiments and multiphysics computational fluid dynamics (CFD) simulations to support reactor licensing and safety. The corresponding research challenges are addressed via supportive verification and validation results generated by the models and simulation tools in combination with selective parametric and uncertainty analysis. This solution approach could blaze the trail for commercial adoption of such technologies. The facilities, simulation capabilities, and research opportunities available at INL in regard to such reactors and the integrated energy systems with which they go hand in hand are also discussed briefly, and may spark interest in deeper research as well as new collaborative projects.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2018

An Approach for Utilizing Power Flow Modeling for Simulations of Hybrid Electric Propulsion Systems

This paper describes an approach to creating simulations of the electric components for a hybrid electric propulsion system. The proposed modeling technique is based on power/load flow modeling and is designed to provide a modular framework that includes buses, lines, and other electrical components that can be connected together to form the electrical distribution system. The purpose of this paper is to detail an electric distribution system modeling technique and to demonstrate how these models may be integrated with turbomachinery simulations. These general modeling techniques were created to be utilized for system and control design studies. Additionally, steady-state and dynamic performance for a proposed model example is compared with data from a hardware in the loop simulation.

Chapman, Jeffryes W.↗
NASA NTRS1994

Direct simulation Monte Carlo and Navier-Stokes simulations of blunt body wake flows

Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

WAKE STRUCTURE↗
DOE OSTI2026

Simulation of Channel Flow with Square Ribs for Blanket First-Wall Cooling: Geometry-Specific Tuning of k-ω Model Using Adjoint Method

Cooling of the plasma-facing first wall is challenging in the design of blanket components because of the high heat flux (on the order of 𝑀𝑊/𝑚2) from the plasma, especially when a low thermal mass medium like helium is chosen as the coolant. Therefore, heat transfer enhancement in which the convective heat transfer rate is augmented by the addition of turbulence-promoting structures becomes a key initiative for providing sufficient cooling capability with helium. Previously, computational fluid dynamics simulations had been performed on pipe flows with different transverse and longitudinal ribbed geometries at Oak Ridge National Laboratory to compare the enhancement performance among different ribbed geometries. Rib shape morphing had been conducted to obtain an optimized rib profile. In the work presented here, the adjoint method is adopted in the ANSYS Fluent solver for turbulence model augmentation, and the Generalized k-ω (GEKO) turbulence model is employed because of its ability of tuning the turbulence model. The Nusselt number and pressure drop obtained from the channel flow with bottom ribbed wall experiments are used as the targets. Sensitivity analysis provides information as guidance to improve the turbulence model accuracy. The augmented GEKO model is tuned for the studied ribbed channel geometry and flow conditions, providing improved predictive accuracy within this context. Extension to other configurations offers potential but may require additional tuning and validation.

Xu, Tracy [ORNL] (ORCID:0009000193700887)↗
NASA NTRS1984

Direct simulation of transitional flow for hypersonic reentry conditions

This paper presents results of flowfield calculations for typical hypersonic reentry conditions encountered by the nose region of the Space Shuttle Orbiter. Most of the transitional flow regime is covered by the altitude range of 150 to 92 km. Calculations were made with the Direct Simulation Monte Carlo (DSMC) method that accounts for translational, rotational, vibrational, and chemical nonequilibrium effects. Comparison of the DSMC heating results with both Shuttle flight data and continuum predictions showed good agreement at the lowest altitude considered. However, as the altitude increased, the continuum predictions, which did not include slip effects, departed rapidly from the DSMC results by overpredicting both heating and drag. The results demonstrate the effects of rarefaction on the shock and the shock layer, along with the extent of the slip and temperature jump at the surface. Also, the sensitivity of the flow structure to the gas-surface interaction model, thermal accommodation, and surface catalysis are studied.

Moss, J. N.↗
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 NTRS1993

Kinetic simulation of plasma flows in the inner magnetosphere

A one-dimensional hybrid particle code is used to study the interactions between upflowing thermal ions from conjugate ionospheres. The simulation model allows for multiple species, convection of plasmaspheric flux tubes, and Coulomb self-collisions which conserve momentum and energy locally. The model incorporates a variable-flux boundary condition where the flux, at the boundaries, approaches zero as the plasmasphere fills and equilibrium conditions are reached. The effects of two important processes on plasmaspheric refilling have been considered. The first includes convection of the plasmaspheric flux tube. The second is the interaction of ionospheric thermal plasma and particle injection from an external source. Particle injection seems to play an important role in the evolution of the total particle distribution on the early timescales (t less than 1 hour); however, for late timescales (t larger than 8 days) the thermal plasma from the ionosphere dominates the particle distribution.

Miller, Ronald H.↗
DOE OSTI2016

High Order Numerical Methods for the Dynamic SGS Model of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. This paper addresses the accuracy improvement of LES of turbulent flows in two ways: (a) from the SGS model standpoint and (b) from the numerical method improvement standpoint. In an internal report, Kotov et al. ( “High Order Numerical Methods for large eddy simulation (LES) of Turbulent Flows with Shocks”, CTR Tech Brief, Oct. 2014, Stanford University), we performed a preliminary comparative study of different approaches to reduce the loss of accuracy within the framework of the dynamic Germano SGS model. The high order low dissipative method of Yee & Sjögreen (2009) using local flow sensors to control the amount of numerical dissipation where needed is used for the LES simulation. The considered improved dynamics model approaches include applying the one-sided SGS test filter of Sagaut & Germano (2005) and/or disabling the SGS terms at the shock location. For Mach 1.5 and 3 canonical shock-turbulence interaction problems, both of these approaches show a similar accuracy improvement to that of the full use of the SGS terms. The present study focuses on a five levels of grid refinement study to obtain the reference direct numerical simulation (DNS) solution for additional LES SGS comparison and approaches. As a result, one of the numerical accuracy improvements included here applies Harten's subcell resolution procedure to locate and sharpen the shock, and uses a one-sided test filter at the grid points adjacent to the exact shock location.

97 MATHEMATICS AND COMPUTING↗