SEARCH · Search NASA

Results for “RANS modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Nek5000 developments in support of industry and the NRC

This year, the Nuclear Energy Advanced Modeling Simulation program (NEAMS) thermal-hydraulics verification and validation (V&V) work has focused in three areas of Nek5000 V&V-driven development. First, in a close collaborative effort with the U. S. Nuclear Regulatory Commission (NRC) staff, we have continued V&V efforts for the HYMERES-2 project using the OECD/NEA sponsored testing in the PSI PANDA facility. This year’s focus of ANL-NRC collaboration involves Nek5000 setups and validation for a range of problems relevant to and including the HYMERES-2 benchmark from PSI. The primary outcome of this year efforts is a more efficient geometry and inlet modeling simplification after a careful sensitivity study of the inlet profiles and pipe geometries. The resulting modeling choice of a short recycling/fully-developed turbulent inlet is within the experimental uncertainty estimate. This finding simplifies the next step of the cross-V&V HYMERES-2 project. In addition, the ANL team continue to provide assistance to the NRC staff in the form of Nek5000 application support in general and on the use of the HPC platforms of ALCF and INL in particular. This supports the NRC’s assessment of Nek5000 for use with the NRC Blue CRAB code suite. Second, we have implemented and tested more robust model of URANS, namely the k – τ model, a variant of the k-ω model, along with other improvements to RANS Nek5000 modeling in general. Because of its demonstrated robustness and stability, the k – τ model is the only RANS model that has been implemented in the new GPU version of the Nek5000 code, nekRS. Lastly, we report the initial implementation of Jacobian-free Newton Krylov approach to the direct Newton method for steady fluid solvers aimed at acceleration of RANS modeling and at IC improvement for LES campaigns. Also leveraging the Exascale Computing Project (ECP) ANL/CEED & SMR team’s software development effort to support NEAMS problems at large scale of the advanced computing architectures, NekRS, a GPU variant of Nek5000, built on top of kernels from libParanumal using OCCA for portability, has been successfully run on the full system of Summit (4608 nodes, 27648 GPUs).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗

DOE OSTI2021

Initial use of Nek5000/Cardinal to improve closure models in Pronghorn

Heat transfer coefficient closure models for pebble bed reactors are built using a data-driven approach by leveraging the capabilities of an Evolutionary Algorithm entitled Particle Swarm Optimization (PSO). In the present work, the Computational Fluid Dynamics code nekRS was used in order to collect the high-fidelity flow data for a core with 1,568 pebbles. To characterize the heat transfer, multiple concentric regions were considered to extract the physical quantities of interest, e.g./ the Reynolds number. The PSO algorithm is employed as part of an inverse problem targeting determine what are the coefficients for a Nusselt number correlation to match the collected data. Such correlation should follow any given format that is defined a priori. Finally, two correlations are proposed, one with an implicit dependence on the pebbles’ wall temperatures and another expressed as a fully explicit correlation depending on the flow conditions and the position within the core. Anyway, given the generic nature of the proposed approach, correlations following different formats could be tested. Preliminary results for the high-fidelity simulation of a fast MSR core are presented. The target Reynolds number is currently 20K, with the expectation that this will increase, pending the availability of further computational resources. These simulations will be used to inform lower fidelity models, including a coarse CFD turbulence model in Pronghorn. Additionally, they will serve as a reference for the RANS models in Nek5000/NekRS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Validation of URANS and STRUCT-ε turbulence models for stratified sodium flow

Simulations of a transient stratified sodium experiment are carried out using a classic unsteady RANS model and the second-generation URANS model, STRUCT-ε. Turbulence modeling challenges and their implications to stratified flow prediction are discussed in the context of other sources of error. Input errors are discussed and addressed; discretization error is calculated to be less than 5% of the inlet velocity, for 80% of the domain; and remaining errors in temperature distributions are attributed to the turbulence model. Qualitative flow features from the simulations are presented and discussed. Compared to the experiment, the STRUCT-ε turbulence model provides a more physically accurate prediction of temperature and momentum mixing in key regions of the domain. Quantitative measures such as the L 2 norm of the temperature discrepancy demonstrate the improved performance of the STRUCT-ε approach. In conclusion, the magnitude of the temperature fluctuations is very well-predicted by the STRUCT-ε, while URANS overpredicts them by approximately 50%.

42 ENGINEERING↗

NASA NTRSDate not supplied

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Adverse Pressure Gradients due to Concave Streamline Curvature

Direct numerical simulations (DNS) of adverse-pressure-gradient turbulent boundary layers over a planar concave wall are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The wall geometry and flow conditions of the DNS are representative of the experimental data for a Mach 4.9 turbulent boundary layer that was tested on a two-dimensional planar concave wall model in the high-speed blow-down wind tunnel located at the National Aerothermochemistry Laboratory at Texas A&M University (TAMU). The DNS was validated against the experimental results of TAMU for the same flow conditions and wall geometry. An analysis of the DNS datasets was also conducted to provide an assessment of the validity of Morkovin’s hypothesis and the strong Reynolds analog for turbulence subject to mechanical nonequilibrium. In addition to the DNS results, RANS predictions are obtained by using the Baldwin-Lomax (BL), Spalart-Allmaras (SA), and the k - w SST turbulence models. The comparisons between RANS and DNS showed little impact of an adverse pressure gradient on the accuracy of these models, at least up to an incompressible Clauser pressure gradient parameter of beta(sub inc) 1.22. While the Boussinesq assumption provided reasonable predictions for the Reynolds shear stress, it failed to adequately predict the normal components of the Reynolds stress.

turbulent boundary layers↗

DOE OSTI2025

Towards stabilization of pressurized oxy-coal diffusion flames

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. Turbulent jet diffusion burners are adopted for the pressurized oxy-combustor in this research and a bluff-body is employed to create a recirculation zone to stabilize the flame in pressurized oxy-combustor. The objective of this numerical work is to perform a systematic analysis of the characteristics of this pressurized non-premixed flame. Specifically, a 15-bar POFC combustor of power 100 kWth is modeled employing the Ansys FLUENT commercial platform, using Reynolds-averaged Navier-Stokes (RANS) modeling. The present computational work focuses on identifying the aerothermodynamic features of the isothermal and reacting flow with a disk-shape bluff-body. Results show that fuel-oxidizer stream momentum ratio has great impacts on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. Parametric study of the width of the burner tip identifies an optimal range for this pressurized burner.

Li, Lei↗

DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗

DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗

NASA NTRS2001

Turbulence Modeling Workshop

Advances in turbulence modeling are needed in order to calculate high Reynolds number flows near the onset of separation and beyond. To this end, the participants in this workshop made the following recommendations. (1) A national/international database and standards for turbulence modeling assessment should be established. Existing experimental data sets should be reviewed and categorized. Advantage should be taken of other efforts already under-way, such as that of the European Research Community on Flow, Turbulence, and Combustion (ERCOFTAC) consortium. Carefully selected "unit" experiments will be needed, as well as advances in instrumentation, to fill the gaps in existing datasets. A high priority should be given to document existing turbulence model capabilities in a standard form, including numerical implementation issues such as grid quality and resolution. (2) NASA should support long-term research on Algebraic Stress Models and Reynolds Stress Models. The emphasis should be placed on improving the length-scale equation, since it is the least understood and is a key component of two-equation and higher models. Second priority should be given to the development of improved near-wall models. Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) would provide valuable guidance in developing and validating new Reynolds-averaged Navier-Stokes (RANS) models. Although not the focus of this workshop, DNS, LES, and hybrid methods currently represent viable approaches for analysis on a limited basis. Therefore, although computer limitations require the use of RANS methods for realistic configurations at high Reynolds number in the foreseeable future, a balanced effort in turbulence modeling development, validation, and implementation should include these approaches as well.

Rubinstein, R.↗

NASA NTRS2000

Towards LES Models of Jets and Plumes

As pointed out by Rodi standard integral solutions for jets and plumes developed for discharge into infinite, quiescent ambient are difficult to extend to complex situations, particularly in the presence of boundaries such as the sea floor or ocean surface. In such cases the assumption of similarity breaks down and it is impossible to find a suitable entrainment coefficient. The models are also incapable of describing any but the most slowly varying unsteady motions. There is therefore a need for full time-dependent modeling of the flow field for which there are three main approaches: (1) Reynolds averaged numerical simulation (RANS), (2) large eddy simulation (LES), and (3) direct numerical simulation (DNS). Rodi applied RANS modeling to both jets and plumes with considerable success, the test being a match with experimental data for time-averaged velocity and temperature profiles as well as turbulent kinetic energy and rms axial turbulent velocity fluctuations. This model still relies on empirical constants, some eleven in the case of the buoyant jet, and so would not be applicable to a partly laminar plume, may have limited use in the presence of boundaries, and would also be unsuitable if one is after details of the unsteady component of the flow (the turbulent eddies). At the other end of the scale DNS modeling includes all motions down to the viscous scales. Boersma et al. have built such a model for the non-buoyant case which also compares well with measured data for mean and turbulent velocity components. The model demonstrates its versatility by application to a laminar flow case. As its name implies, DNS directly models the Navier-Stokes equations without recourse to subgrid modeling so for flows with a broad spectrum of motions (high Re) the cost can be prohibitive - the number of required grid points scaling with Re(exp 9/4) and the number of time steps with Re(exp 3/4). The middle road is provided by LES whereby the Navier-Stokes equations are formally filtered with the filter chosen to only exclude the smallest turbulent motions. If successful, LES should provide much of the detail available to DNS but at more bearable cost. Fatica et al. in comparing LES with DNS for a low Reynolds number jet showed that the LES could simulate the temporally evolving behavior including growth of the jet thickness. It is the intention of this report to explore the application of an LES model to jets and plumes. As always, before tackling complex situations, the model must be tested for the simplest of cases and so we address only two, a non-buoyant axisymmetric jet issuing steadily from an orifice into a semi-infinite stationary environment and a buoyant jet in the same environment. The work is a continuation of Basu and Mansour.

Webb, A. T.↗

NASA NTRS2012

Investigation of Transonic Wake Dynamics for Mechanically Deployable Entry Systems

A numerical investigation of transonic flow around a mechanically deployable entry system being considered for a robotic mission to Venus has been performed, and preliminary results are reported. The flow around a conceptual representation of the vehicle geometry was simulated at discrete points along a ballistic trajectory using Detached Eddy Simulation (DES). The trajectory points selected span the low supersonic to transonic regimes with freestream Mach numbers from 1:5 to 0:8, and freestream Reynolds numbers (based on diameter) between 2:09 x 10(exp 6) and 2:93 x 10(exp 6). Additionally, the Mach 0:8 case was simulated at angles of attack between 0 and 5 . Static aerodynamic coefficients obtained from the data show qualitative agreement with data from 70deg sphere-cone wind tunnel tests performed for the Viking program. Finally, the effect of choices of models and numerical algorithms is addressed by comparing the DES results to those using a Reynolds Averaged Navier-Stokes (RANS) model, as well as to results using a more dissipative numerical scheme.

Stern, Eric↗

NASA NTRS2017

Design of an Axisymmetric Afterbody Test Case for CFD Validation

As identified in the CFD Vision 2030 Study commissioned by NASA, validation of advanced RANS models and scale-resolving methods for computing turbulent flow fields must be supported by continuous improvements in fundamental, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, the underpinnings of a new test platform referred to herein as the NASA Axisymmetric Afterbody are presented. The devised body-of-revolution is a modular platform consisting of a forebody section and afterbody section, allowing for a range of flow behaviors to be studied on interchangeable afterbody geometries. A body-of-revolution offers advantages in shape definition and fabrication, in avoiding direct contact with wind tunnel sidewalls, and in tail-sting integration to facilitate access to higher Reynolds number tunnels. The current work is focused on validation of smooth-body turbulent flow separation, for which a six-parameter body has been developed. A priori RANS computations are reported for a risk-reduction test configuration in order to demonstrate critical variation among turbulence model results for a given afterbody, ranging from barely-attached to mild separated flow. RANS studies of the effects of forebody nose (with/without) and wind tunnel boundary (slip/no-slip) on the selected afterbody are presented. Representative modeling issues that can be explored with this configuration are the effect of higher Reynolds number on separation behavior, flow physics of the progression from attached to increasingly-separated afterbody flows, and the effect of embedded longitudinal vortices on turbulence structure.

Disotell, Kevin J.↗

DOE OSTI2023

Simulations for Planning of Liquid Hydrogen Spill Test

In order to better understand the complex pooling and vaporization of a liquid hydrogen spill, Sandia National Laboratories is conducting a highly instrumented, controlled experiment inside their Shock Tube Facility. Simulations were run before the experiment to help with the planning of experimental conditions, including sensor placement and cross wind velocity. This paper describes the modeling used in this planning process and its main conclusions. Sierra Suite’s Fuego, an in-house computational fluid dynamics code, was used to simulate a RANS model of a liquid hydrogen spill with five crosswind velocities: 0.45, 0.89, 1.34, 1.79, and 2.24 m/s. Two pool sizes were considered: a diameter of 0.85 m and a diameter of 1.7. A grid resolution study was completed on the smaller pool size with a 1.34 m/s crosswind. A comparison of the length and height of the plume of flammable hydrogen vaporizing from the pool shows that the plume becomes longer and remains closer to the ground with increasing wind speed. The plume reaches the top of the facility only in the 0.45 m/s case. From these results, we concluded that it will be best for the spacing and location of the concentration sensors to be reconfigured for each wind speed during the experiment.

08 HYDROGEN↗

NASA NTRS2017

Development of an Axisymmetric Afterbody Test Case for Turbulent Flow Separation Validation

As identified in the CFD Vision 2030 Study commissioned by NASA, validation of advanced RANS models and scale-resolving methods for computing turbulent flows must be supported by improvements in high-quality experiments designed specifically for CFD implementation. A new test platform referred to as the Axisymmetric Afterbody allows for a range of flow behaviors to be studied on interchangeable afterbodies while facilitating access to higher Reynolds number facilities. A priori RANS computations are reported for a risk-reduction configuration to demonstrate critical variation among turbulence model results for a given afterbody, ranging from barely-attached to mild separated flow. The effects of body nose geometry and tunnel-wall boundary condition on the computed afterbody flow are explored to inform the design of an experimental test program.

Disotell, Kevin J.↗

NASA NTRSDate not supplied

Exploring the Accuracy of RANS Simulations for Mars Entry Vehicles

Accurate yet inexpensive predictions of aerodynamic coefficients for Mars entry vehicles have remained a consistent challenge over the past five decades. Below Mach 6, drag on the backshell becomes significant and must be accurately predicted. Steady Reynolds-averaged Navier-Stokes (RANS) models are commonly used, despite their poor predictions of backshell pressure. While scale-resolving simulations have shown promise in the past decade, there is still a need for cheap, accurate RANS predictions for large aerodynamic databases. The Mars Science Laboratory (MSL) is used as a case study to examine predictive accuracy and known shortcomings for RANS predictions of Mars entry vehicles. Several different grid generation techniques are compared, including a comparison between prismatic boundary layer grids and fully unstructured, tetrahedral grids. Comparisons are made to experimental data for Mach 2.5, 3.5, and 4.5. The accuracy of predicted aerodynamic coefficients is examined. Overpredictions in axial force and drag are explained by a closer examination of the surface pressure. These findings document sensitivities and best practices for future RANS database development of Mars entry vehicles.

RANS↗

NASA NTRSDate not supplied

Exploring the Accuracy of RANS Simulations for Mars Entry Vehicles

CFD↗

DOE OSTI2020

Gas Generation Models for Thermal Runaway in Lithium Ion Batteries

A numerical model is implemented to study cell venting, internal pressure, and gas-phase dynamics behavior of Li-ion cells undergoing thermal runaway. A k-ɛ Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer’s equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. A series of computational fluid dynamics (CFD) simulations are conducted on a single cell at various states-of-charge (100%, 50%, 25%) to demonstrate detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is shown that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

Kim, Jinyong↗

DOE OSTI2021

Modeling cell venting and gas-phase reactions in 18650 lithium ion batteries during thermal runaway

A numerical model is developed to study cell venting, internal pressure, and gas-phase dynamics behavior of 18650 Li-ion cells undergoing thermal runaway. A k-e Reynolds-Averaged Navier-Stokes (RANS) model is adopted to describe the turbulent flow out of the cells, while the fluid dynamics inside the cells is described by Darcy-Forchheimer's equation. Thermal abuse reactions and gas generation kinetics are described by a single-step lumped reaction model. Then, a series of computational fluid dynamics (CFD) simulations are conducted on a single 18650 cell at various states-of-charge (100%, 50%, 25%) to study detailed flow and thermal behavior as a function of quantity of gas generated during cell venting. Venting events are categorized into two stages: i) breaching of the cell container and ii) thermal runaway reactions. It is found that the cell response is dominated by the second stage since most of the gases are generated during thermal runaway. Also, the propensity for propagation is highly affected by state-of-charge (SOC). Cells at higher SOCs produce more heat and gas during the venting event, owing to higher mass and concentrations of reacting gases, and consequently reach higher internal cell pressures which increase the risk of side-wall breaching.

25 ENERGY STORAGE↗