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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A method to generate small-scale, high-resolution sedimentary bedform architecture models representing realistic geologic facies

Small-scale (mm to m) sedimentary structures (e.g. ripple lamination, cross-bedding) have received a great deal of attention in sedimentary geology. The influence of depositional heterogeneity on subsurface fluid flow is now widely recognized, but incorporating these features in physically-rational bedform models at various scales remains problematic. The current investigation expands the capability of an existing set of open-source codes, allowing generation of high-resolution 3D bedform architecture models. The implemented modifications enable the generation of 3D digital models consisting of laminae and matrix (binary field) with characteristic depositional architecture. The binary model is then populated with petrophysical properties using a textural approach for additional analysis such as statistical characterization, property upscaling, and single and multiphase fluid flow simulation. One example binary model with corresponding threshold capillary pressure field and the scripts used to generate them are provided, but the approach can be used to generate dozens of previously documented common facies models and a variety of property assignments. An application using the example model is presented simulating buoyant fluid (CO 2 ) migration and resulting saturation distribution.

58 GEOSCIENCES↗
DOE OSTI2026

Direct observation of diamond formation in a shock-compressed high explosive

Understanding the formation timescale and structure of carbonaceous reaction products is critical for modeling the high-pressure equation-of-state of organic materials. We use the National Ignition Facility to shock-compress polycrystalline TATB (C 6 H 6 N 6 O 6 ) samples to ~70–130 GPa and ~4000–5500 K, employing in situ nanosecond X-ray diffraction to probe reaction products and velocimetry to measure transmitted compression wave profiles. Our diffraction data is consistent with the formation of diamond over timescales less than ~60 ns. This represents carbon condensation from a molecular explosive on timescales three times faster than previously reported and the earliest observation of diamond produced from reacting TATB. Reactive flow simulations with explicit chemistry reproduce the observed temporal structure within wave profiles to inform the distribution of P-T states. These findings provide direct evidence of ultrafast diamond formation in a reactive system at extreme conditions and provide new constraints for models of shock and detonation chemistry.

Clarke, Samantha M. [Lawrence Livermore National L↗
DOE OSTI2020

Reaction heterogeneity in practical high-energy lithium–sulfur pouch cells

The lithium–sulfur (Li–S) battery is a promising next-generation energy storage technology because of its high theoretical energy and low cost. Extensive research efforts have been made on new materials and advanced characterization techniques for mechanistic studies. However, it is uncertain how discoveries made on the material level apply to realistic batteries due to limited analysis and characterization of real high-energy cells, such as pouch cells. Evaluation of pouch cells (>1 A h) (instead of coin cells) that are scalable to practical cells provides a critical understanding of current limitations which enables the proposal of strategies and solutions for further performance improvement. Herein, we design and fabricate pouch cells over 300 W h kg -1 , compare the cell parameters required for high-energy pouch cells, and investigate the reaction processes and their correlation to cell cycling behavior and failure mechanisms. Spatially resolved characterization techniques and fluid-flow simulation reveal the impacts of the liquid electrolyte diffusion within the pouch cells. We found that catastrophic failure of high-energy Li–S pouch cells results from uneven sulfur/polysulfide reactions and electrolyte depletion for the first tens of cycles, rather than sulfur dissolution as commonly reported in the literature. The uneven reaction stems from limited electrolyte diffusion through the porous channels into the central part of thick cathodes during cycling, which is amplified both across the sulfur electrodes and within the same electrode plane. A combination of strategies is suggested to increase sulfur utilization, improve nanoarchitectures for electrolyte diffusion and reduce consumption of the electrolytes and additives.

25 ENERGY STORAGE↗
DOE OSTI2014

Vibrational and electronic collisional-radiative model in air for Earth entry problems

The two-temperature collisional-radiative model CoRaM-AIR, working over a wide range for pressure and temperatures, has been developed for the flow conditions around a space vehicle entering the Earth's atmosphere. The species N{sub 2}, O{sub 2}, NO, N, O, Ar, N{sub 2}{sup +}, O{sub 2}{sup +}, NO{sup +}, N{sup +}, O{sup +}, Ar{sup +}, and free electrons are taken into account. The model is vibrationally specific on the ground electronic state of N{sub 2}, O{sub 2}, and NO, and electronically specific for all species, with a total of 169 vibrational states and 829 electronic states, respectively. A wide set of elementary processes is considered under electron and heavy particle impact given the temperatures involved (up to 30 000 K). This set corresponds to almost 700 000 forward and backward elementary processes. The relaxation from initial thermal or chemical nonequilibrium is studied for dissociation-ionization situations in conditions related to the FIRE II flight experiment. Boltzmann plots clearly prove that the vibrational and electronic excitation distributions are far from being Boltzmanian. In particular, high-lying vibrational levels remain underpopulated for most of the duration of the relaxation. This relaxation can be separated in a first phase characterized by the dissociation and the excitation of the molecular species, and a second phase leading to the excitation and the ionization of the dissociation products. Owing to the vibrational relaxation, the time scales are slightly higher than the ones predicted by former kinetic mechanisms usually used in flow simulations. In the present FIRE II conditions, radiation does not play a significant role.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Radial–axial transport coordination enhances sugar translocation in the phloem vasculature of plants

Abstract Understanding mass transport of photosynthates in the phloem of plants is necessary for predicting plant carbon allocation, productivity, and responses to water and thermal stress. Several hypotheses about optimization of phloem structure and function and limitations of phloem transport under drought have been proposed and tested with models and anatomical data. However, the true impact of radial water exchange of phloem conduits with their surroundings on mass transport of photosynthates has not been addressed. Here, the physics of the Munch mechanism of sugar transport is re-evaluated to include local variations in viscosity resulting from the radial water exchange in two dimensions (axial and radial) using transient flow simulations. Model results show an increase in radial water exchange due to a decrease in sap viscosity leading to increased sugar front speed and axial mass transport across a wide range of phloem conduit lengths. This increase is around 40% for active loaders (e.g. crops) and around 20% for passive loaders (e.g. trees). Thus, sugar transport operates more efficiently than predicted by previous models that ignore these two effects. A faster front speed leads to higher phloem resiliency under drought because more sugar can be transported with a smaller pressure gradient.

58 GEOSCIENCES↗
DOE OSTI2024

A Novel Framework to Quantify Power Grid Resilience

The quantification of an operating power grid’s resilience is highly significant today, given its criticality as an enabler of other infrastructures, complexity, and the threat it faces due to a wide range of detrimental events, from extreme climate to cyber attacks. Currently, there exist no standardized definitions and metrics for measuring the resilience of an operating grid. In this paper, we introduce a novel resilience quantification framework and demonstrate a method to measure the flexibility towards topological/structural changes due to potential failures in the power grid to assess operational resilience. We start with the state estimation data from a large utility and use the graph analysis methods and power flow simulation tools to compute the identified resilience parameters.

Yoginath, Srikanth↗
DOE OSTI2023

Constraints on OPF Surrogates for Learning Stable Local Volt/Var Controllers

We consider the problem of learning local Volt/Var controllers in distribution grids (DGs). Our approach starts from learning separable surrogates that take both local voltages and reactive powers as arguments and predict the reactive power setpoints that approximate optimal power flow (OPF) solutions. We propose an incremental control algorithm and identify two different sets of slope conditions on the local surrogates such that the network is collectively steered toward desired configurations asymptotically. Our results reveal the trade-offs between each set of conditions, with coupled voltage-power slope constraints allowing an arbitrary shape of surrogate functions but risking limitations on exploiting generation capabilities, and reactive power slope constraints taking full advantage of generation capabilities but constraining the shape of surrogate functions. AC power flow simulations on the IEEE 37-bus feeder illustrate their guaranteed stability properties and respective advantages in two DG scenarios.

asymptotic stability↗
DOE OSTI2020

Analytical Voltage Sensitivity Analysis for Unbalanced Power Distribution System

Large scale integration of distributed energy resources and electric vehicles in a transactive energy environment present new challenges in terms of voltage stability and fluctuations in a power distribution system. The impact of different level of DER/EV penetration on the voltages across the network is typically quantified through voltage sensitivity analyses. Existing methods of voltage sensitivity analysis are computationally expensive and prior efforts to develop analytical approximation lacks generality and have not been effectively validated. The objective of this work is to provide a new analytical method of voltage sensitivity analysis that has low computational cost and also allows for stochastic analysis of voltage change. This paper first derives an analytical approximation of change in voltage at a particular bus due to change in power consumption at other bus in a radial three phase unbalanced power distribution system. Then, the proposed method is shown to be valid for different load configurations, which demonstrates its generality. The results from our analytical approach is validated via classical load flow simulation of the test system based on IEEE 37 bus network. The proposed method is shown to have good accuracy, and computation complexity is of order O(1), compared to O(n3) in classical sensitivity analysis approaches.

Munikoti, Sai↗
DOE OSTI2018

Chance-constrained Service Restoration for Distribution Networks with Renewables

Power systems operation is facing great challenges from natural disasters and cyber-attacks. It is critical but also difficult to enhance the reliability and resilience against extreme events. To better response to inevitable outages or blackouts, service restoration in distribution networks is important to minimize the disastrous impacts of catastrophic events. The increasing penetration of distributed energy resources (DERs) provides new opportunities to expedite the restoration process. However, the coordination with conventional distribution system control devices and the uncertainty and variability of intermittent renewable energy resources requires new operation and control strategies for distribution service restoration (DSR). This paper develops an optimal bottom-up DSR strategy by coordinating DERs with voltage regulators and capacitor banks. The chance-constrained (CC) programming approach is used to model the probabilistic output limit of solar radiation and PV generation. The Markov's inequality and Latin hypercube sampling techniques are applied to convert and incorporate the chance constraints into the DSR optimization problem. The CC-DSR problem is formulated as a mixed integer convex programming problem, considering various operational cost functions and bidirectional three-phase unbalanced load flow. Simulation results on the modified IEEE 13-node test feeder system demonstrate the effectiveness and flexibility of the bottom-up DSR strategy.

Chance constrained optimization, Distribution serv↗
DOE OSTI2026

Closed-Loop Wind Turbine Controllers for Active Wake Mixing Strategies

Wind turbine active wake mixing (AWM) is an exciting new field of research where dynamic actuation, usually on the blade pitch angles, is used to increase wind farm-wide power production. From a controls perspective, the current state-of-the-art AWM strategies are very simple: a dynamic, usually periodic, reference signal is prescribed to actuators in an open-loop (OL) fashion. The actuation is then presumed to have a certain desired effect on the system, i.e., the wind turbine and the flow it affects. However, this system is highly nonlinear and experiences disturbances in the form of wind variations that are not known a priori. As a result, the OL approach might not yield optimal results. In this article, a novel approach is presented, which closes the loop on AWM controllers. A feedback loop is implemented, which uses measurements of the individual blade bending moments that are widely available on modern wind turbines to determine the individual blade pitch angles. A proof of concept of this implementation is presented in this article, and a thorough comparison with the OL method is executed using high-fidelity flow simulations. These simulations show that the novel closed-loop controller does not substantially increase wake mixing but achieves similar performance as the OL controller while reducing fatigue loads on the controlled turbine.

17 WIND ENERGY↗
DOE OSTI2024

Study of eVTOL Charging Impact on Airport Electrical Grids

Growth in commercial air traffic has increased the aviation industry's contribution to global emissions. Aircraft electrification can reduce greenhouse gas emissions along with being efficient and reliable. Advanced air mobility introduces the idea of an advanced and more efficient air transportation system. Electric and hybrid aircraft designs, such as electric vertical takeoff and landing (eVTOL), are being considered for a number of applications, such as urban air mobility, cargo delivery, and medical emergencies. Although the electrification of aircraft provides environmental benefits, their charging needs will place significant demand on the electrical grid. This work presents a grid impact study and infrastructure upgrade requirements analysis of eVTOL charging on a realistic distribution grid through power flow simulation analysis. The case study results show that the eVTOL charging can lead to power quality issues, such as undervoltage and equipment capacity overloading. Grid infrastructure, such as lines and transformers, would need to be upgraded, and voltage management strategies would need to be deployed to accommodate eVTOL charging at an existing airport.

advanced air mobility (AAM)↗
DOE OSTI2024

Fracture length data for geothermal applications

Fracture lengths govern permeability and are unknowns in geothermal assessment. Along their lengths, fracture widths vary due to growth by linkage. Under the influence of diagenesis, narrow widths seal, breaking porosity continuity and reducing open length. The largest range of widths and thus susceptibility to fill occurs where fractures are linked by narrow segments. Outcrops of a geothermal target, Cambrian Potsdam quartz arenite, contain opening-mode fractures having lengths spanning five orders of magnitude from 0.082 mm to 17.9 m. Combined lengths measured at a range of scales can be described by power laws, but at a given image resolution, lengths are best fit by exponential functions. Owing to preferential sealing of small fractures, open fractures follow exponential functions, but values depend on rules for designating fractures as continuous. En échelon segments, offset 10 mm, are connected by narrow fractures or microfractures (hard linked) not evident on outcrop 1 m-elevation LiDAR or 30 m-height drone images. A rule that identifies where narrow and likely connected segments are located can yield lengths meaningful for flow simulation. Depending on diagenesis, continuity rules can halve or double average and maximum length values. Length values from outcrop for geothermal applications should be adjusted based on wellsite-specific diagenesis information.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

PeleLMeX [SWR-22-48]

PeleLMeX is a solver for high fidelity reactive flow simulations, namely direct numerical simulation (DNS) and large eddy simulation (LES). The solver combines a low Mach number approach, adaptive mesh refinement (AMR), embedded boundary (EB) geometry treatment and high performance computing (HPC) to provide a flexible tool to address research questions on platforms ranging from small workstations to the world's largest GPU-accelerated supercomputers. PeleLMeX has been used to study complex flame/turbulence interactions in RCCI engines and hydrogen combustion or the effect of sustainable aviation fuel on gas turbine combustion. PeleLMeX is part of the Pele combustion Suite (https://amrex-combustion.github.io/)

Day, Marcus↗
DOE OSTI2023

MARBLES (Multi-scale Adaptively Refined Boltzmann LatticE Solver) [SWR-23-37]

MARBLES (Multi-scale Adaptively Refined Boltzmann LatticE Solver) is an open-source computational fluid dynamics package powered by the lattice Boltzmann equations and built on AMReX. In the lattice Boltzmann method, local collisions between meso-scale fictitious particles drive the governing equations which enables MARBLES to easily simulate flow around complex and/or moving geometry without the generation of a body-conforming mesh. Using AMReX data structures and operations ensures a high level of computational performance and parallel scaling on heterogenous architectures while also naturally supporting locally enhanced grid resolution and fidelity through automatic mesh refinement. New domains and problem definitions are easily specified through an input file with examples and guidance on all options and variables provided in the MARBLES documentation.

Henry de Frahan, Marc↗
DOE OSTI2025

pdas-experiments

SAND2025-04589O pdas-experiments automates computational experiments of fluid flow simulations. It uses the pressio-demoapps-schwarz package as a basis to break down complex simulations into smaller, manageable parts. This application is an extension of the Sandia Pressio software which uses domain decomposition to work with complex simulations more efficiently. Users can test different simulation setups, while keeping a detailed record of their experiments so they can be reproduced later. The software includes a C++ program that runs individual experiments based on user-defined settings in a YAML file, as well as a Python script that can manage multiple simulations at once. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Tezaur, Irina [Sandia National Lab. (SNL-CA), Live↗