Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor - Extended Abstract

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗
DOE OSTI2021

Coupled Multiphysics Simulations of Heat Pipe Microreactors Using DireWolf

DireWolf is a multiphysics software driver application designed to simulate heat pipe–cooled nuclear microreactors. Developed under the U.S. Department of Energy, Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the DireWolf software application’s objective is to provide the nuclear community with a design and safety analysis simulation capability. Based upon the NEAMS program Multiphysics Object-Oriented Simulation Environment (MOOSE) computational framework, DireWolf tightly couples nuclear microreactor physics, reactor physics, radiation transport, nuclear fuel performance, heat pipe thermal hydraulics, power generation, and structural mechanics to resolve the interdependent nonlinearities. DireWolf is capable of simulating both steady and transient normal reactor operation and several postulated failure scenarios. We will present the fundamental physics of heat pipe–cooled nuclear microreactors and the MOOSE-based software employed in DireWolf. Both steady and transient results for coupled reactor physics, radiation transport, and nuclear fuel performance are demonstrated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Fully Implicit Conjugate Heat Transfer Analysis of the ARC-Class Vacuum Vessel

The coupled simulation of fusion reactor blankets including neutronics, thermal-hydraulics and thermo-mechanics is expected to speed up the design cycle of fusion reactor design concepts. In this work we demonstrate tight implicit coupling of conjugate heat transfer using the open-source Computational Fluid Dynamics software OpenFOAM for thermo-fluid mechanics and Diablo for thermo-solid mechanics. The heat transfer analysis is augmented by volumetric energy deposition from neutronic calculations using the Monte Carlo N-particle code on both solid and fluid parts of the vacuum vessel. An additional heat flux is imposed on the first wall estimated from the design power of the reactor. The tight coupling is realized through the open-source coupling library, preCICE, and tested on the vacuum vessel of the affordable, robust, compact reactor design by Commonwealth Fusion Systems. The features of the coupling and the influence of different coupling parameters such as coupling schemes, acceleration techniques and convergence criterion are discussed. The coupled simulation results are compared to a thermal-hydraulics simulation which includes only the fluid domains (the liquid immersion molten salt blanket and cooling channel) to demonstrate usefulness of a coupled simulation. Further analysis is performed to identify regions of hot spots for subsequent design improvement. This introduces the outline for integrating conjugate electromagnetics and fluid/solid mechanics (e.g., allow for deformation of the cooling channel walls) with our present approach for future analysis.

Sircar, Arpan↗
DOE OSTI2025

In-Transit Data Transport Strategies for Coupled AI-Simulation Workflow Patterns

Coupled AI-Simulation workflows are becoming the major workloads for HPC facilities, and their increasing complexity necessitates new tools for performance analysis and prototyping of new in-situ workflows. We present SimAI-Bench, a tool designed to both prototype and evaluate these coupled workflows. In this paper, we use SimAI-Bench to benchmark the data transport performance of two common patterns on the Aurora supercomputer: a one-to-one workflow with co-located simulation and AI training instances, and a many-to-one workflow where a single AI model is trained from an ensemble of simulations. For the one-to-one pattern, our analysis shows that node-local and DragonHPC data staging strategies provide excellent performance compared Redis and Lustre file system. For the many-to-one pattern, we find that data transport becomes a dominant bottleneck as the ensemble size grows. Our evaluation reveals that file system is the optimal solution among the tested strategies for the many-to-one pattern.

Tummalapalli, Harikrishna [Argonne National Labora↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2022

Quantum Simulations of Vibrational Strong Coupling via Path Integrals

A quantum simulation of vibrational strong coupling (VSC) in the collective regime via thermostated ring-polymer molecular dynamics (TRPMD) is reported. For a collection of liquid-phase water molecules resonantly coupled to a single lossless cavity mode, the simulation shows that as compared with a fully classical calculation, the inclusion of nuclear and photonic quantum effects does not lead to a change in the Rabi splitting but does broaden polaritonic line widths roughly by a factor of 2. Moreover, under thermal equilibrium, both quantum and classical simulations predict that the static dielectric constant of liquid water is largely unchanged inside vs outside the cavity. Furthermore, this result disagrees with a recent experiment demonstrating that the static dielectric constant of liquid water can be resonantly enhanced under VSC, suggesting either limitations of our approach or perhaps other experimental factors that have not yet been explored.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

NPP Simulators for Coupled Thermal and Electric Power Dispatch

The Light Water Reactor Sustainability (LWRS) program within the United States Department of Energy supports extending the operation of the U.S. commercial nuclear power plant (NPP) fleet. Within the LWRS program, the Flexible Plant Operation and Generation (FPOG) Pathway works to diversify the revenue streams of light water reactors (LWRs) by opening opportunities for the co-generation of non-electric products in addition to supplying electrical power to the grid. Recent events have added greater motivation to these efforts. For example, the recent Inflation Reduction Act (IRA) passed by the U.S. federal government offers substantial tax incentives for producing clean hydrogen, the technology readiness level of dispatchable and high-efficiency hydrogen production has dramatically increased in a short time, and societal response to world climate change is driving a transition away from fossil fuels. Producing hydrogen with maximum efficiency using nuclear power requires dispatching both electrical and thermal power from the nuclear plant to the hydrogen plant, so testing concepts of operations for combined electrical and thermal power dispatch (TPD) from an NNP to a hydrogen plant is of interest. This report documents achievement of the Light Water Reactor Sustainability (LWRS) program milestone “Install and demonstrate a vendor-developed simulator on the Human Systems Simulation Laboratory (HSS) for dispatch of LWR electrical power to a close-coupled electrolysis plant” with a due date of Dec. 22, 2022. Several factors provide motivation for this effort. Coupling the power generation deck of a nuclear power plant to a hydrogen production facility introduces new possibilities for operational transients that must be addressed. In particular, the performance of the integrated system during startup and shutdown of the hydrogen production facility, as well as offnormal conditions, need to be evaluated to ensure there are no adverse effects on the operation of the existing NPP. The concept of operations involving the NPP, the hydrogen plant, and the electric power grid must be tested using NPP simulators and operating procedures that have been modified for TPD operations. These tests must also include dynamic simulations of the coupled tertiary thermal and electric loads as well as coordinated activities with NPP operators, tertiary load operators and grid power coordinators. The report summarizes progress in developing and testing full-scope NPP simulators at the HSSL, including a generic BWR simulator from GSE Systems, Inc. and generic PWR simulator from Westinghouse. In the case of the TPD-GBWR Simulator from GSE Systems, Inc., a BWR is thermally coupled to a high temperature electrolysis (HTE) plant that produces hydrogen and oxygen from de-ionized water. The hydrogen plant is not explicitly simulated but only included as a transient heat sink. A thermal power dispatch (TPD) system transfers heat between the steam systems at the BWR and the hydrogen plant. Operational results from two versions of the modified simulator are presented. The first version uses synthetic oil as a heat transfer fluid in a closed delivery heat loop (DHL) that generates steam at the hydrogen plant. The second version uses steam as the heat transfer fluid in a delivery steam line (DSL) to provide steam to the hydrogen plant. For both versions, the estimated thermal power delivery distance is approximately one kilometer. The amount of thermal power dispatched in the simulators is 15% of the total reactor thermal power such that the simulators provide a tool to study the feasibility of coupling a BWR to industrial processes that benefit from a combination electrical and thermal power dispatch. Ongoing work within a CRADA is also developing a full-scope PWR simulator provided by Westinghouse for both thermal and electric power coupling. This simulator is based on a PWR plant with two three-loop Westinghouse reactors. Westinghouse PWRs are sufficiently similar that a simulator of a three-loop reactor is an appropriate representation for two-loop and four-loop PWR reactors. The three-loop simulator will initially be modified for close-coupling to a 100 MW HTE hydrogen production plant that will require approximately 25 MW of thermal power while operating at its maximum rated capacity. The simulator testing will include full coupling to dynamic simulations of a hydrogen production plant and a representative bulk electric grid. The simulator provided by Westinghouse is similar to the GPWR simulator that INL has already obtained from GSE Systems but has a few important added benefits. First, the Westinghouse simulator is based on digital controls and has additional screens that can be called up to show parameter trends to assist operators in decision-making. The Westinghouse simulator also has upgrades to the controls and hardware representations, such as valve actuators, that make it more realistic and flexible in terms of accurately sim

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

CFD Simulation of Floating Body Motion with Mooring Dynamics: Coupling MoorDyn with OpenFOAM

It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.

floating body↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Mutually guided light and particle beam propagation

The polarizability of atoms and molecules gives rise to optical forces that trap particles and a refractive index that guides light beams, potentially leading to a self-guided laser and particle beam propagation. In this paper, the mutual interactions between an expanding particle beam and a diffracting light beam are investigated using an axisymmetric particle-light coupled simulation. The nonlinear coupling between particles and photons is dependent on the particle beam radius, particle density, particle velocity and temperature, polarizability, light beam waist, light frequency (with respect to the resonance frequency), and light intensity. The computational results show that the maximum propagation distance is achieved when the waveguiding effect is optimized to single-mode operation. The application of the coupled beam propagation as a space propulsion system is discussed.

42 ENGINEERING↗
DOE OSTI2025

Incorporating Elevation in Traffic-Vehicle CO-Simulation: Issues, Impacts, and Solutions

Traffic-vehicle co-simulation couples microscopic traffic simulation with full-body vehicle dynamics to assess system-level impacts on mobility, energy, and safety with greater realism. Incorporating elevation is critical for accurately modeling vehicle behavior and energy use, especially for gradient-sensitive vehicles such as electric and heavy-duty trucks. However, raw elevation data often contain noise, discontinuities, and inconsistencies. While such issues may be negligible in traditional traffic simulations, they significantly affect traffic-vehicle co-simulations where vehicle dynamics are sensitive to road grade variations. This paper investigates the impact of unprocessed elevation data on vehicle behavior and energy consumption using a 42-mile simulation along Interstate 81. We propose an elevation processing workflow that can mitigate the effects stem from elevation data issues, improving the realism and stability of traffic-vehicle co-simulation. Results show that the method effectively removes noise and abrupt elevation transitions while preserving roadway geometry.

Xu, Guanhao [ORNL] (ORCID:0000000214326357)↗
DOE OSTI2020

Electron beam irradiated foil expansion simulated by a coupled HYDRA/MCNP6 application

Simulating the hydrodynamic evolution of a temporally pulsed electron beam with multiple metal foils requires coupling the nearly monoenergetic, relativistic electron beam to the foil target. This deposition may be effectively calculated by Monte Carlo transport. Once this energy is deposited on a computational mesh representing the foil target, highly developed radiation-hydrodynamic codes, such as HYDRA, can describe the subsequent evolution. Currently, HYDRA does not have an internal Monte Carlo electron transport capability. To partially rectify this deficiency, a PYTHON-based scheme to link MCNP6 to HYDRA has been developed for quadrilateral meshes. Results from a seven-foil tantalum target configuration are presented to demonstrate this capability.

42 ENGINEERING↗
DOE OSTI2025

Simulating nationwide coupled disease and fear spread in an agent-based model

Human cognitive responses, behavioral responses, and disease dynamics co-evolve over the course of any disease outbreak, and can result in complex feedbacks. We present a dynamic agent-based model that explicitly couples the spread of disease with the spread of fear surrounding the disease, implemented within the EpiCast simulation framework. EpiCast models transmission within a realistic synthetic population, capturing individual-level interactions. In our model, fear propagates through both in-person contact and broadcast media, prompting individuals to adopt protective behaviors that reduce disease spread. In order to better understand these coupled dynamics, we create and compare a range of compartmental models to ensure that introducing additional disease states does not prevent the emergence of multiple waves in these simpler models. Additionally, we compare a range of behavioral scenarios within EpiCast, varying the level and intensity of fear and behavior change. Our results show that the addition of asymptomatic, exposed, and pre-symptomatic disease states can impact both the rate at which an outbreak progresses and its overall trajectory in compartmental models. In EpiCast, the combination of non-local fear spread via broadcasters and strong behavioral responses by fearful individuals generally leads to multiple epidemic waves, an outcome that occurs only within a narrow parameter range when fear spreads purely through local contact. Accounting for the coupled spread of fear and disease is critical for understanding disease dynamics and designing timely, targeted responses to emerging infectious threats.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2020

Multi-Physics Simulations of Molten Chloride Fast Reactor using Nek5000 and PROTEUS-NODAL

This report documents the FY20 work to demonstrate a specific multi-physics capability within NEAMS tools to model molten chloride fast reactors (MCFR). Specifically, the CFD code Nek5000 and the neutronics code PROTEUS-NODAL were coupled off-line to calculate the delayed neutron precursor distribution within the flowing fuel salts in an MCFR core. An MCFR benchmark problem and a 3D test problem based on a conceptual MSR design were successfully solved through the coupled simulation of Nek5000 and PROTEUS-NODAL. In the coupled simulations for these test problems, Nek5000 provided high-quality flow field information, and PROTEUS-NODAL performed the neutron diffusion calculation along with the precursor drift model. The impact of precursor drift on the core status was also investigated. Additionally, the precursor drift model implemented in PROTEUS-NODAL was verified by making use of the reference precursor distribution generated from Nek5000. Results of these test problems suggest that the multi-physics capability of Nek5000 and PROTEUS-NODAL was properly developed for tracking delayed neutron precursor drift.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Bridging particle and reactor scales in the simulation of biomass fast pyrolysis by coupling particle resolved simulation and coarse grained CFD-DEM

Efficient utilization of biomass is challenging due to the complexity of the feedstock at molecular, particle, and reactor scales. Direct coupling of particle scale model and reactor scale simulation is computationally infeasible. In this work, the results from a particle resolved simulation were used to calibrate the heat transfer and reaction kinetics of particles. A hybrid drag model was used to consider the different fluidization properties of sands and biomass particles. The computation cost was reduced using a coarse-grained Discrete Element Method (DEM). The simulated conversion rate, char yield, tar yield, and pyrolysis gas yield compare well with experiment. The residence time and axial distributions of biomass were also analyzed. This novel multi-scale method provides an efficient and accurate tool for the modeling of biomass pyrolysis reactors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Time-dependent equation-of-motion coupled-cluster simulations with a defective Hamiltonian

Here simulations of laser-induced electron dynamics in a molecular system are performed using time-dependent (TD) equation-of-motion (EOM) coupled-cluster (CC) theory. The target system has been chosen to highlight potential shortcomings of truncated TD-EOM-CC methods [represented in this work by TD-EOM-CC with single and double excitations (TD-EOM-CCSD)], where unphysical spectroscopic features can emerge. Specifically, we explore driven resonant electronic excitations in magnesium fluoride in the proximity of an avoided crossing. Near the avoided crossing, the CCSD similarity-transformed Hamiltonian is defective, meaning that it has complex eigenvalues, and oscillator strengths may take on negative values. When an external field is applied to drive transitions to states exhibiting these traits, unphysical dynamics are observed. For example, the stationary states that make up the time-dependent state acquire populations that can be negative, exceed one, or even complex-valued.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗