Study of radiative heat transfer and flow physics from medium-scale methanol pool fire simulations
Not Available
Kirsch, Jared (ORCID:0000000268852942)↗
SEARCH · Search NASA
Results for “flow physics”
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.
ExaWind: Exascale Predictive Wind Plant Flow Physics Modeling
The scientific goal of the ExaWind project is to advance our fundamental understanding of the flow physics governing whole wind plant performance, including wake formation, complex terrain impacts, and turbine-turbine-interaction effects. The primary application codes in the ExaWind environment are Nalu-Wind, an unstructured-grid computational fluid dynamics (CFD) code, AMR-Wind, a structured-grid CFD code, and OpenFAST, a whole-turbine simulation code. In this poster we present our progress towards simulating the ExaWind challenge problem, which is a predictive simulation of a wind farm with tens of megawatt-scale wind turbines dispersed over an area of 50 square kilometers.
49 EE - Wind and Water Power Program - Wind (EE-4W↗
Mixed precision efficiency improvements for multi-physics flow simulations
This report summarizes research activities at Cascade Technologies during February 2020 - November 2020 to support DOE SBIR phase I contract DE-SC0020548. The primary technical objectives of the work are two-fold: a) to establish a theoretical framework to understand the propagation and impact of finite-precision-related errors and b) to quantify the trade-offs between accelerator throughput and solution accuracy through a series of validation tests involving multi-physics flow phenomena.
97 MATHEMATICS AND COMPUTING↗
Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer
This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Egg yolk as a model for gelation: From rheometry to flow physics
Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.
Marsh, Maxwell C. [Department of Mechanical Scienc↗
Wind Plant Flow Physics and Power Performance in Complex Environments: Cooperative Research and Development (Final Report)
Cornell University will partner with NLR on the topic of wind farm wake effects to improve understanding of interactions between complex atmospheric flows, terrain, and wind turbine wakes and plant efficiency. Wind plant flow simulation tools will also be validated. The work performed will help improve wind farm modeling by analyzing data, applying models, designing and performing experiments to acquire additional wind farm data, and develop better models.
17 WIND ENERGY↗
MRT 7365: Power flow physics and key physics phenomena
The Z accelerator at Sandia National Laboratories conducts z-pinch experiments at 26 MA in support of DOE missions in stockpile stewardship, dynamic materials, fusion, and other basic sciences. Increasing the current delivered to the z-pinch would extend our reach in each of these disciplines. To achieve increases in current and accelerator efficiency, a fraction of Z’s shots are set aside for research into transmission-line power flow. These shots, with supporting simulations and theory, are incorporated into this Advanced Diagnostics milestone report. The efficiency of Z is reduced as some portion of the total current is shunted across the transmission-line gaps prior to the load. This is referred to as “current loss”. Electrode plasmas have long been implicated in this process, so the bulk of dedicated power-flow experiments are designed to measure the plasma environment. The experimental analyses are enhanced by simulations conducted using realistic hardware and Z voltage pulses. In the same way that diagnostics are continually being improved for sensitivity and resolution, the modeling capability is continually being improved to provide faster and more realistic simulations. The specifics of the experimental hardware, diagnostics, simulations, and algorithm developments are provided in this report. The combined analysis of simulation and data confirms that electrode plasmas have the most detrimental impact on current delivery. Experiments over the last three years have tested the theoretical current-loss mechanisms of enhanced ion current, plasma gap closure, and Hall-related current. These mechanisms are not mutually exclusive and may be coincident in the final feed as well as in upstream transmission lines. The final-feed geometries tested here, however, observe lower-density plasmas without dominant ion currents which is consistent with a Hall-related current. The picture of plasma formation and transport formed from experiment and simulation is informing hardware designs being fielded on Z now and being proposed for the Next-Generation Pulsed Power (NGPP) facility. In this picture, the strong magnetic fields that heat the electrodes above particle emission thresholds also confine the charged particles near the surface. Some portion of the plasmas thus formed is transported into the transmission-line gap under the force of the electric field, with aid from plasma instabilities. The gap plasmas are then transported towards the load by a cross-field drift, where they accumulate and contribute to a likely Hall-related cross-gap current. The achievements in experimental execution, model validation, and physical analysis presented in this report set the stage for continued progress in power flow and load diagnostics on Z. The planned shot schedule for Z and Mykonos will provide data for extrapolation to higher current to ensure the predicted performance and efficiency of a NGPP facility.
43 PARTICLE ACCELERATORS↗
MFC 5.0: An exascale many-physics flow solver
Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions
Computational fluid dynamics↗
Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics
n this paper, we introduce a multigrid block-based preconditioner for solving linear systems arising from a Discontinuous Galerkin discretization of the all-speed Navier-Stokes equations with phase change. The equations are discretized in conservative form with a reconstructed Discontinuous Galerkin (rDG) method and integrated with fully-implicit time discretization schemes. To robustly converge the numerically stiff systems, we use the Newton-Krylov framework with a primitive-variable formulation (pressure, velocity, and temperature), which is better conditioned than the conservative-variable form at low-Mach number. In the limit of large acoustic CFL number and viscous Fourier number, there is a strong coupling between the velocity-pressure system and the linear systems become non-diagonally dominant. To effectively solve these ill-conditioned systems, an approximate block factorization preconditioner is developed, which uses the Schur complement to reduce a 3 x 3 block system into a sequence of two 2 x 2 block systems: velocity-pressure,v P , and velocity-temperature, v T . We compare the performance of the v P -v T Schur complement preconditioner to classic preconditioning strategies: monolithic algebraic multigrid (AMG), element-block SOR, and primitive variable block Gauss-Seidel. The performance of the preconditioned solver is investigated in the limit of large CFL and Fourier numbers for low-Mach lid-driven cavity flow, Rayleigh-Bénard melt convection, compressible internally heated convection, and 3D laser-induced melt pool flow. Here, numerical results demonstrate that the v P -v T Schur complement preconditioned solver scales well both algorithmically and in parallel, and is robust for highly ill-conditioned systems, for all tested rDG discretization schemes (up to 4th-order).
97 MATHEMATICS AND COMPUTING↗
Complementary Study of Radiative Heat Transfer and Flow Physics from Moderate-scale Hydrocarbon Pool Fire Simulations
As part of the Advanced Simulation and Computing Verification and Validation (ASCVV) program, a 0.3-m diameter hydrocarbon pool fire with multiple fuels was modeled and simulated. In the study described in this report, systematic examination was performed on the radiation model used in a series of coupled Fuego/Nalu simulations. A calibration study was done with a medium-scale methanol pool fire and the effect of calibration traced throughout the radiation model. This analysis provided a more detailed understanding of the effect of radiation model parameters on each other and on other quantities in the simulations. Heptane simulation results were also examined using this approach and possible areas for further improvement of the models were identified. The effect of soot on radiative losses was examined by comparing heptane and methanol results.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
ExaWind: Exascale Predictive Wind Plant Flow Physics Modeling.
Abstract not provided.
Sprague, Michael↗
ExaWind: Exascale Predictive Wind Plant Flow Physics Modeling.
Abstract not provided.
Sprague, Michael A.↗
Study of Radiative Heat Transfer and Flow Physics from Medium-scale Methanol Pool Fire Simulations [Slides]
Analysis of methanol pool fire conducted as part of validation study for SIERRA/Fuego. Radiation model was effectively calibrated by modifying radiation model parameters for methanol. Computing integrated buoyancy flux, entrainment rate, and turbulent kinetic energy allowed for evaluation of less typical quantities in this validation study. Quantities were compared with experimental data or correlations and generally showed agreement. TKE needed fine mesh to be computed accurately. Predicted flame height less sensitive to variations in mixture fraction than temperature. Mixture fraction is a preferable threshold variable for this application.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
MRT 7365 Power flow physics and key physics phenomena: EMPIRE verification suite
This milestone work baselines electromagnetic particle-in-cell capability of the EMPIRE plasma simulation code to model key processes germane to the physics of electrode plasmas arising in magnetically-insulated transmission lines operating at or near 20 MA. This evaluation is done so through the provision of benchmark verification problems designed to exercise the individual and combined physics models on a small-scale surrogate geometry for the final-feed-to-load region of the Z accelerator under representative operating conditions. In this report, we overview our test designs, and present a portfolio of simulation results along with performance assessments which altogether establish state-of-the-art. In particular, two main verification categories are covered this report: (1) Z-relevant desorption physics (Temkin isotherm), and (2) two approaches to simulate electrode plasma creation and dynamics (automatic creation versus self-consistent creation through direct simulation Monte Carlo collisions).
43 PARTICLE ACCELERATORS↗
Mixed precision efficiency improvements for multi-physics flow simulations
Explore the source record for details and available documents.
97 MATHEMATICS AND COMPUTING↗
Cylinders and jets in crossflow: Wake formations as a result of varying geometric proximities
The combined flow physics of several canonical flow configurations is experimentally studied. Here, we analyze an array of jets issuing into a crossflow, then immediately navigating past an array of cylinders. This is achieved with a 2 × 3 triangular pattern of jets and symmetric cylinders at three jets to crossflow velocity ratios, enabling near-complete optical access of each jet, with velocities measured by time-resolved particle image velocimetry. Jet trajectories reveal that each configuration adheres to a power-law trend and that greater penetration is achieved by the downstream and confined jets compared to the more conventional upstream one. Recirculation regions of the upstream and downstream jets are nearly independent, with the confined jet encompassing regions of overlap with both. Turbulent statistics reveal the influence of geometric placement and velocity ratio on the time-averaged velocity, anisotropy, and Reynolds stresses incurred by each jet. Galilean decomposition utilizes a supplemental crossflow-only velocity field to delineate the influence of each jet's low- and high-pressure regions on the otherwise uniform stream. Proper orthogonal decomposition suggests that increased jet penetration decreases the number of modes required for truncation in the investigated spanwise plane. Vortex identification algorithms are applied to the reconstructed flow fields, reaffirming that with increasing velocity ratio, the jets generate vortices of their own in similar statistical formations as the cylinders. This investigation provides a foundation to aid future modeling efforts toward characterizing flow physics of importance in designing and passively controlling transverse jets with varying blockage proximities in a crossflow.
Kristo, Paul J. (ORCID:0000000225605175)↗
Physical and Flow Properties of Glass Forming Chemicals (V2O5, SnO, SnO2, Cr2O3, FeCr2O4, and ZrSiO4) and Mixtures
For a sustainable nuclear waste vitrification process at the Hanford Tank Waste Treatment and Immobilization Plant (WTP), proper selection and consistent supply of glass-forming chemicals (GFCs) are crucial. Establishing rigorous acceptance criteria for the characterization of GFCs will be required to operate the vitrification facility and to mitigate any processing issues or failures. Low-activity wastes (LAW) are blended with GFCs to form slurry melter feeds and vitrified in a melter. To enhance properties of waste glasses, new chemicals are being introduced to the current GFC mixture (Vienna et al. 2016; Muller et al. 2017, 2019). In this study, three new GFCs were evaluated for enhanced LAW glass formulations: chromium oxide (Cr 2 O 3 ), vanadium oxide (V 2 O 5 ), and stannic oxide (SnO 2 ). These three oxide components are included in enhanced waste glass (EWG) formulations and GFCs with the appropriate physical and flow properties are needed. As a starting point, single metal oxide GFCs: Cr 2 O 3 , V 2 O 5 , and SnO 2 were sourced and tested. To characterize these new individual GFCs and mixtures of GFCs, the industrial bulk characterization consultant, Jenike and Johanson, was employed to measure physical and flow properties of individual GFCs and their mixtures. Pacific Northwest National Laboratory (PNNL) also measured several selected physical properties for data evaluation as a quality assurance step. In addition, PNNL measured physical and rheological properties of slurry melter feeds containing those GFCs. Subsequent data analyses and verification were conducted. The purpose of this report is to assess the applicability of these GFCs for LAW vitrification based on their properties. This report will help understand measured data and evaluate new GFCs for use. Moreover, this report may give useful insights to help troubleshoot any GFC and melter feed transport and mixing issues that arise during processing, leading to a successful cleanup mission at WTP.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗