SEARCH · Engineering Papers

Results for “Flow simulations”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

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 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↗

DOE OSTI2024

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

CO2 Plume Imaging with Accelerated Deep Learning-based Data Assimilation Considering Multiple Realizations: Application to the Illinois Basin-Decatur Carbon Sequestration Project

We propose a fast and efficient deep learning workflow for near real-time data assimilation, forecasting and visualization of CO2 plume evolution in saline aquifer and demonstrate its application at a field site. Unlike the previous work, this study incorporates the impact of spatial heterogeneity using multiple realizations. In the proposed workflow, a neural network model utilizes available monitoring data such as downhole pressure measurements as input and predicts the propagating pressure ‘front’ using the diffusive time of flight (DTOF) map which is considered as representative reservoir image of the flow field. The DTOF is the arrival time of pressure front propagation, which can be computed by the Fast Marching Method rapidly without flow simulations. Reservoir model calibration can be implemented by selecting the training data samples that describe the predicted DTOF map based on observed data. The power and efficacy of our workflow is demonstrated by application to the Illinois Basin-Decatur Project.

CO2 plume imaging↗

DOE OSTI2026

Ultra-High Vacuum Outgassing Characterization of Thermally Processed Low-Carbon Steel for Advanced Particle Accelerator and Gravitational Wave Detector Applications

This dissertation investigated AISI 1020 low-carbon steel as an alternative vacuum chamber material to conventional stainless steel for ultra-high vacuum (UHV) and extreme-high vacuum (XHV) applications. After a 400 °C/48 h bake, AISI 1020 tube chambers achieved a hydrogen outgassing rate of 2.4 × 10¿¹6 Torr·L·s¿¹·cm¿², approximately 2,300 times lower than the prebaked 316L stainless-steel comparator, among the lowest hydrogen outgassing rates ever reported for an uncoated metallic vacuum chamber. Bare and magnetite-coated AISI 1020 chambers were then compared using throughput and rate-of-rise methods. The magnetite coating yielded 5× lower water outgassing at room temperature, but this advantage disappeared after 80 °C baking. After full thermal conditioning (400 °C/48 h prebake followed by 150 °C/96 h and 200 °C/110 h), bare steel achieved 25× lower hydrogen outgassing than the magnetite-coated chamber (9.6 × 10¿¹6 Torr·L·s¿¹·cm¿²) and >99% H2 purity with carbon species below RGA detection. Monte Carlo molecular flow simulations of a CEBAF photogun beamline (96 scenarios) showed that replacing 304L stainless steel with AISI 1020 reduces equilibrium H2 pressure by a factor of 833; a single 304L electrode contributes 98.8% of the gas load despite occupying only 9.1% of the internal surface area. A 500-m Einstein Telescope beampipe screening showed that corrugated bellows contribute 18% of the gas load from only 0.7% of the surface area. A five-model adsorption isotherm framework applied to 22 pumpdown datasets (164 fits with AR(1)-GLS correction) established that the experimental protocol, not the material, controls isotherm identifiability: Dubinin–Radushkevich wins isothermal pipe pumpdowns; Langmuir wins thermally dominated chamber bakes. Cross-dataset joint fitting of the AISI 1020 pipe pumpdowns yielded an H2 diffusion activation energy Ed = 7.24 ± 1.28 kcal·mol¿¹, consistent with trap dominated diffusion in commercial low-carbon steels. Two companion innovations were developed: a Variable Conductance Device (VCD, patent pending IDF-00723) for XHV outgassing measurement, and VacuumDesignerPro (VDP), a MATLAB-based design tool validated against LIGO benchmarks.

Al-Allaq, Aiman H [Old Dominion University]↗

DOE OSTI2024

StOKeDMD: Streaming Occupation kernel dynamic mode decomposition

Dynamic mode decomposition (DMD) has become a common technique for constructing surrogate models for dynamical systems from observed system states. The Occupation Kernel DMD (OKDMD) method proposed in (Rosenfeld et al., 2022) and (Rosenfeld et al., 2024) is a Liouville operator based method that builds surrogate models from system state trajectories. Here, this paper proposes an extension of OKDMD to the case when the system states are observed in a streaming fashion, i.e., only a small fraction of the state trajectory is available at a given time. The developed method, Streaming Occupation Kernel DMD (StOKeDMD), accommodates the streaming data input by leveraging properties of specific choices of kernel functions and occupation kernels. We apply the StoKeDMD method as a compression method for streaming data, analyze the memory complexity, and demonstrate the performance of StoKeDMD in the compression of streaming data generated from a Lorenz system and a fluid flow simulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Mitigating Impact Through Community-Engaged Flood Modeling

Urban pluvial flooding poses a growing threat to the city of Baltimore, driven by heavy rainfall, increased impervious area, and aging infrastructure. Adapting to the risks posed by pluvial flooding is critical for building greater climate resiliency in Baltimore's Inner Harbor Watershed. This study addresses these challenges through community-informed decision analysis, which uses hydrologic modeling and optimization tools to identify robust flooding adaptation pathways. We will collaborate with community partners to identify key concerns and objectives regarding flooding. These concerns have been purposefully built in to a combined surface-subsurface dynamic flow simulation model. Model outputs are used to identify flooding locations within the Inner Harbor, and to test adaptation methods. Machine learning will be used search for solutions which meet diverse environmental, financial, and social goals, and solution performance will be examined under a wide range of potential future climatic conditions and integrated with an adaptive planning approach. This novel set of adaptation pathways will enhance the City's capacity to respond to evolving pluvial flood risk.

climate resilience↗

DOE OSTI2026

Three-Dimensional Grid Visualization for Planning Activities: A Dubai Case Study

National Laboratory of the Rockies (NLR), in collaboration with the Dubai Electricity and Water Authority (DEWA) and Infra-X, has undertaken the Energy Visualization Analysis Project. The aim of this project is to enhance analytical and 3D visualization capabilities for distribution network planning and renewable energy integration. As modern grid continues to evolve with large-scale solar PV deployment and emerging distributed energy resources (DERs), the ability to effectively analyze, visualize, and communicate complex grid behaviors has become increasingly critical. The project focuses on developing empirical use cases based on real distribution feeder data and engineering workflows, ensuring the outcomes are directly aligned with operational environment. Through time-series power flow simulations and nodal hosting capacity analysis, the study quantifies the impacts of high PV penetration on voltage and thermal limits within representative 11 kV feeders. These analyses identify specific nodes and conditions where DER integration challenges arise. Furthermore, a Battery Energy Storage System (BESS) optimization algorithm was applied to determine the optimal size and placement of storage systems that can mitigate network constraints and enhance hosting capacity. The comparative results between base-case and BESS-augmented scenarios clearly demonstrate improvements in network stability and load management efficiency. In parallel, the NLR team developed an immersive 3D visualization framework, enabling interactive exploration of grid simulations using commodity head-mounted display (HMD) systems. This framework transforms conventional 2D simulation data into spatially intuitive visual environments - allowing engineers to analyze feeder conditions, PV hosting potential, and BESS effects in real time. This report represents the first foundational phase in establishing a visualization-driven analytical ecosystem. It provides a methodological foundation for data integration, visualization architecture, and simulation-based decision support, paving the way for large-scale adoption of immersive visualization across DEWA's Smart Grid Initiative, R&D activities, and future network resilience studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Unsupervised Learning for Equitable DER Control: Preprint

In the context of managing distributed energy resources (DERs) within distribution networks (DNs), this work focuses on the task of developing local controllers. We propose an unsupervised learning framework to train functions that can closely approximate optimal power flow (OPF) solutions. The primary aim is to establish specific conditions under which these learned functions can collectively guide the network towards desired configurations asymptotically, leveraging an incremental control approach. The flexibility of the proposed methodology allows to integrate fairness-driven components into the cost function associated with the OPF problem. This addition seeks to mitigate power curtailment disparities among DERs, thereby promoting equitable power injections across the network. To demonstrate the effectiveness of the proposed approach, power flow simulations are conducted using the IEEE 37-bus feeder. The findings not only showcase the guaranteed system stability but also underscore its improved overall performance.

asymptotic stability↗

DOE OSTI2024

Integrating a Water Tracer Model Into WRF‐Hydro for Characterizing the Effect of Lateral Flow in Hydrologic Simulations

Abstract Most current land models approximate terrestrial hydrological processes as one‐dimensional vertical flow, neglecting lateral water movement from ridges to valleys. Such lateral flow is fundamental at catchment scales and becomes crucial for finer‐scale land models. To test the effect of incorporating lateral flow toward three‐dimensional representations of hydrological processes in the next generation land models, we integrate a water tracer model into the WRF‐Hydro framework to track water movement from precipitation to discharge and evapotranspiration. This hydrologic‐tracer integrated system allows us to identify the key mechanisms by which lateral flow affects the flow paths and transit times in WRF‐Hydro. By comparing modeling experiments with and without lateral routing in two contrasting catchments, we determine the impacts of lateral flow on the transit times of precipitation event‐water. Results show that with limited hydrologic connectivity, lateral flow extends the transit times by reducing (increasing) event‐water drainage loss (accumulation) in ridges (valleys) and allowing reinfiltration of infiltration‐excess flow, which is missing in most land models. On the contrary with high hydrologic connectivity, lateral flow can effectively accelerate the water release to streams and reduce the transit time. However, the transit times are substantially underestimated by the model compared with isotope‐derived estimates, indicating model limitations in representing flow paths and transit times. This study provides some insights on the fundamental differences in terrestrial hydrology simulated by land models with and without lateral flow representation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Counter-Current Flow Limitation Studies in Complex Geometries Utilizing Interface Capturing Simulations Coupled with PID Flow Rate Controller

In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.

Counter-current flow limitation↗

DOE OSTI2026

An Approach for Simulating Transport of Nonspherical, Sublimating Particles in Rarefied Gas Flows via the Direct Simulation Monte Carlo Method

This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Anomaly detection with flow-based fast calorimeter simulators

Recently, several normalizing flow-based deep generative models have been proposed to accelerate the simulation of calorimeter showers. Using caloflow as an example, we show that these models can simultaneously perform unsupervised anomaly detection with no additional training cost. As a demonstration, we consider electromagnetic showers initiated by one (background) or multiple (signal) photons. The caloflow model is designed to generate single-photon showers, but it also provides access to the shower likelihood. We use this likelihood as an anomaly score and study the showers tagged as being unlikely. As expected, the tagger struggles when the signal photons are nearly collinear but is otherwise effective. This approach is complementary to a supervised classifier trained on only specific signal models using the same low-level calorimeter inputs. While the supervised classifier is also highly effective at unseen signal models, the unsupervised method is more sensitive in certain regions, and, thus, we expect that the ultimate performance will require a combination of these approaches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Large-Eddy Simulation Study of Flow and Combustion Dynamics in a Full-Scale Hydrogen–Air Rotating Detonation Combustor-Stator Integrated System

In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.

Combustion↗

DOE OSTI2025

Modeling and simulation of multiphase flows

Multiphase flow can be found in a wide range of natural and industrial systems, such as sediment flow in rivers and coasts to pipeline transport, and reactors of all types. Computational fluid dynamics (CFD) modeling provides a rigorous physics-based method to visualize and quantify these flows and to assess their influence on their surroundings. This presentation will highlight computational research on two varied technologies in the energy sector: carbon capture and wellbore acoustics.

acoustic analysis↗

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)↗

DOE OSTI2025

Nonlinear simulation of under-resolved flows with shocks

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

97 MATHEMATICS AND COMPUTING↗