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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗
DOE OSTI2022

Toward control co-design of utility-scale wind turbines: Collective vs. individual blade pitch control

A large-eddy simulation framework has been coupled with controller modules to systematically investigate the impacts of collective (CPC) and individual (IPC) pitch control strategies on utility-scale wind turbine energy production and fatigue loads. Wind turbine components were parameterized using an actuator surface model to simulate the rotor blades and the turbine nacelle. The baseline CPC and IPC algorithms, consisting of single-input single-output proportional–integral controllers and two integral controllers, respectively, were incorporated into the numerical framework. A series of simulations were carried out to investigate the relative performance of the two controllers under various turbulent inflow conditions, spanning hub-height velocities of 7 to 14 m/s. The numerical simulation results of this study showed that, in comparison to the CPC, the IPC controller could successfully reduce the damage equivalent loads of utility-scale turbines at regions 2 and 3 of turbine operation by about 3% and 40%, respectively, without any penalty on the power production of the turbine. It was also shown that, despite its minor impact on the turbulence kinetic energy of the wake, the IPC controller did not influence the recovery of the turbine wake.

17 WIND ENERGY↗
DOE OSTI2023

Analysis of ultrafast magnetization switching dynamics in exchange-coupled ferromagnet–ferrimagnet heterostructures

Magnetization switching in ferromagnets has so far been limited to the current-induced spin–orbit–torque effects. Recent observation of helicity-independent all-optical magnetization switching (HI-AOS) in an exchange-coupled ferromagnet–ferrimagnet (FM-FEM) heterostructures expanded the range and applicability of such ultrafast heat-driven magnetization switching. Here we report the element-resolved HI-AOS dynamics of such an exchange-coupled system, using a modified microscopic three-temperature model. We have studied the effect of (i) the Curie temperature of the FM, (ii) FEM composition, (iii) the long-range Ruderman–Kittel–Kasuya–Yosida (RKKY) exchange-coupling strength, and (iv) the absorbed optical energy on the element-specific time-resolved magnetization dynamics. The phase-space of magnetization illustrates how the RKKY coupling strength and the absorbed optical energy influence the switching time. Our analysis demonstrates that the threshold switching energy depends on the composition of the FEM and the switching time depends on the Curie temperature of the FM as well as RKKY coupling strength. Importantly, this simulation anticipates new insights into developing faster and more energy-efficient spintronics devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗
DOE OSTI2025

Microreactor Automated Control System Test Bed Digital Architecture for Real-Time, Hardware-in-the-Loop Simulation

This work describes progress made towards the development of a real-time hardware-in-the-loop (HIL) test bed for non-nuclear testing of microreactor control schemes and failure modes. Non-nuclear testing is a crucial step in developing robust control algorithms for managing microreactor dynamics. The creation of an HIL simulation harnesses the realistic dynamics of physical analogue systems while additionally considering the challenges of variable communication delay. This collaborative effort between Oak Ridge National Laboratory and Idaho National Laboratory has resulted in a LabVIEW-based gRPC communication protocol which couples a TRANSFORM Modelica simulation of nuclear components to the ViBRANT physical hardware for realistic feedback and visual representation of control action in real time. A modular python client structure is developed to manage FMU-based Modelica simulation and real-time gRPC communication. HIL testing suggests that the modeled reactor with natural convection molten salt loop coolant configuration responds well to PID control of drum positioning for modulation of reactor core power, however, future efforts will be made to explore the added thermal inertial delay of system level control and downstream demand changes. Development of this platform with a generalized methodology provides a foundation for exploring a variety of reactor configurations and failure modes in rapid order to provide insight into the most effective avenues of study for further research and development.

McConnell, Jono [ORNL] (ORCID:0000000238984741)↗
DOE OSTI2020

Origin and effect of film sub-stoichiometry on ultraviolet, ns-laser damage resistance of hafnia single layers

Understanding the origin of laser damage-prone precursors in high index materials such as hafnia holds the key to the development of laser damage-resistant multilayer dielectric coated optics for high power and energy laser systems. In this study, we investigate the source of sub-stoichiometry, a potent laser damage precursor, in hafnia films produced by an ion beam sputtering (IBS) deposition method and the effect of such defects on the film performance upon ns ultraviolet (UV) laser (8 ns, 355 nm) exposure. Chemical analysis of data obtained via Rutherford backscattering spectroscopy (RBS) suggests that hafnia films deposited at two different planetary locations from the same deposition run exhibit anisotropic and location-dependent stoichiometries. While the oxygen-to-hafnium ratio is at the stoichiometric value of 2 for the hafnia film at the edge location, the ratio is significantly deviated and is 1.7 for that deposited at the planetary center. The sub-stoichiometric hafnia films display a much lower 1-on-1 damage onset at 1.6 ± 0.2 J/cm 2 compared to 2.3 ± 0.2 J/cm 2 in a stoichiometric film. The low damage performance films also have an over three times higher damage density at fluences above initiation. Coupled with Monte Carlo simulations, we reveal that sub-stoichiometry is primarily attributed to preferential removal of oxygen during film deposition by the bombardment of energetic reflected argon neutrals. The resulting oxygen deficiencies create the sub-bandgap states which facilitate the strong laser energy coupling and reduce the resistance to laser-induced damage in the hafnia single layer films.

36 MATERIALS SCIENCE↗
DOE OSTI2021

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗
DOE OSTI2021

Inference of Rock Flow and Mechanical Properties from Injection-Induced Microseismic Events During Geologic CO 2 Storage

Monitoring microseismic activities during CO 2 injection into geologic formations is important for ensuring the safety of the storage operations. The resulting data provide insight into the response of the storage formation to CO 2 injection and can be used to infer the underlying rock flow and mechanical properties. In this paper, assimilation of microseismic data is performed for dynamic characterization of the storage formation by using a stochastic simulation model to forecast the microseismic response of a geologic formation during CO 2 injection. Two modeling approaches are adopted to predict the space-time distribution of the injection-induced microseismicity. The first model is based on pore pressure relaxation assumption, while the second model uses coupled flow and geomechanics simulation to establish the complex physical relation between the storage formation properties and the corresponding microseismic responses during CO 2 injection. The stochastic predictive models in each case are used in ensemble data assimilation frameworks to estimate rock properties from the observed microseismic data. Two data assimilation methods are considered: (i) a new ensemble-based stochastic point process filter (EnPPF) that can directly integrate discrete microseismic events, and (ii) a variant of ensemble smoother, known as the ensemble smoother with multiple data assimilation (ES-MDA), which requires continuous representation of microseismic events for assimilation. The two methods are successfully applied to a geologically realistic model of the Farnsworth Field in Texas, with complex geologic flow units and interacting fault systems.

42 ENGINEERING↗
DOE OSTI2025

Aeroelastic Modelling of Large Wind Turbines: Towards a Unified OpenFAST-SEAHOWL Approach

In recent years, the scale of wind turbines has significantly increased to maximize energy capture for a given site (particularly offshore), presenting new challenges in terms of structural design and dynamics. As towers grow taller and blades grow longer, flexion and torsion of the latter have a non-negligible impact on the behavior and performance of the turbine in terms of overall loads, power production, and control. When representing large-scale wind turbines numerically to capture these important effects, particular attention must therefore be given to the level of fidelity for representing each structural component, as well as the coupling scheme used between them to keep simulations accurate, stable, and efficient. To address this issue, we combine here the two following tools: (1) OpenFAST, the reference whole-turbine simulation tool from NREL with standalone modules covering each physics and the choice between loose coupling and a new tight coupling scheme for structural dynamics, and (2) SEAHOWL, the whole-turbine simulation tool from TotalEnergies with monolithic coupling of structural dynamics through Project Chrono and partitioned coupling for multiphysics interactions.

17 WIND ENERGY↗
DOE OSTI2024

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗
DOE OSTI2020

Stability Analysis of Interface Conditions for Ocean-Atmosphere Coupling

In this paper, we analyze the stability of different coupling strategies for multidomain PDEs that arise in general circulation models used in climate simulations. We focus on fully coupled ocean–atmosphere models that are needed to represent and understand the complicated interactions of these two systems, becoming increasingly important in climate change assessment in recent years. Numerical stability issues typically arise because of different time-stepping strategies applied to the coupled PDE system. In particular, the contributing factors include using large time steps, lack of accurate interface flux, and single-iteration coupling. We investigate the stability of the coupled ocean–atmosphere models for various interface conditions such as the Dirichlet–Neumann condition and the bulk interface condition, which is unique to climate modeling. By analyzing a simplified model, we demonstrate here how the parameterization of the bulk condition and other numerical and physical parameters affect the coupling stability and establish stability conditions for different coupling strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Thermodynamic stability of a spin microemulsion in Rashba spin-orbit-coupled bosons

Recent finite-temperature numerical simulations have unveiled a quantum “spin” microemulsion analog, found by raising the temperature of a stripe supersolid phase in a Rashba spin-orbit coupled Bose gas. This microemulsion state is a highly correlated, isotropic normal fluid where atoms self-arrange based on their internal pseudospin into patterns that resemble bicontinuous microemulsions. This finding leaves several open questions regarding the broader accessibility of this phase in experiments. Here, we use equilibrium finite-temperature numerical simulations based on a coherent-state path integral representation to perform a computational investigation into the thermodynamic stability of the spin microemulsion state. Numerical simulations emphasize the requirement of a nearly, but not perfectly, isotropic spin-orbit coupling in order to achieve the microemulsion phase in cold-atom experiments. Moreover, the microemulsion state exists independent of miscibility of the pseudospin components and for a wide range of pseudospin population imbalance, suggesting a high degree of flexibility in choosing the atom and hyperfine states in an experimental realization. Lastly, we demonstrate this feasibility by mimicking a Rashba spin-orbit-coupled 87 Rb experiment in an isotropic harmonic trap, where we confirm the microemulsion's existence via its density profile and equilibrium quasimomentum distribution.

Complex Langevin dynamics↗
DOE OSTI2020

Identifying Characteristics of Wildfire Towers and Troughs

Wildfire behavior is dictated by the complex interaction of numerous physical phenomena including dynamic ambient and fire-induced winds, heat transfer, aerodynamic drag on the wind by the fuel and combustion. These phenomena create complex feedback effects between the fire and its surroundings. In this study, we aim to study the mechanisms by which buoyant flame dynamics along with vortical motions and instabilities control wildfire propagation. Specifically, this study employs a suite of simulations conducted with the physics-based coupled fire-atmosphere behavior model (FIRETEC). The simulations are initialized with a fire line and the fires are allowed to propagate on a grass bed, where the fuel heights and wind conditions are varied systematically. Flow variables are extracted to identify the characteristics of the alternating counter-rotational vortices, called towers and troughs, that drive convective heat transfer and fire spread. These vortices have previously been observed in wildfires and laboratory fires, and have also been observed to arise spontaneously in FIRETEC due to the fundamental physics incorporated in the model. However, these past observations have been qualitative in nature and no quantitative studies can be found in the literature which connected these coherent structures fundamental to fire behavior with the constitutive flow variables. To that end, a variety of state variables are examined in the context of these coherent structures under various wind profile and grass height conditions. Identification of various correlated signatures and fire-atmosphere feedbacks in simulations provides a hypothesis that can be tested in future observational or experimental efforts, potentially assisting experimental design, and can aid in the interpretation of data from in situ detectors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Reionization time of the Local Group and Local-Group-like halo pairs

ABSTRACT Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present-day observations. Several numerical simulation-based studies show correlations between the reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired haloes within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation, and isolation criteria but in other environments, alongside isolated haloes within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out, given our reionization model. Within this model, haloes in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: Pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ∼625 Myr against 660 ± 4 Myr (z ∼ 8.25 against 7.87 ± 0.02) on average.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Impact of solvent forces and broken symmetry on the assembly of designed proteins at a liquid-solid interface

The era of protein design has enabled the creation of hybrid protein-inorganic interfaces, leading to both surface-directed self-assembly of de novo protein architectures and protein-directed formation of inorganic materials. However, the resulting patterns of protein assembly are often unexpected, implying that essential interactions are not accounted for in current design platforms. Here, we use high-speed atomic force microscopy (AFM) analyzed through machine learning to follow the assembly of protein nanorods in aqueous electrolytes on two types of mica exhibiting disparate symmetry elements, which are imprinted on the overlying hydration structure. Using Monte Carlo simulations, we reproduce the observed phases and show that an observed smectic phase, previously thought to be unstable for non-interacting rods in two dimensions, emerges when crystal symmetry introduces a directional bias. The findings demonstrate the importance of incorporating solvent forces as modulated by the hydration structure inherent to interfacial systems when designing protein assemblies at liquid-crystal interfaces. Coupling physics-based simulations that can account for these factors to de novo protein design algorithms can lead to improved design platforms for bio-inspired, hybrid materials.

bioinspired materials↗
DOE OSTI2021

A Co-design Framework for Online Data Analysis and Reduction

Science applications preparing for the exascale era are increasingly exploring in situ computations comprising of simulation-analysis-reduction pipelines coupled in-memory. Efficient composition and execution of such complex pipelines for a target platform is a codesign process that evaluates the impact and tradeoffs of various application- and system-specific parameters. In this article, we describe a toolset for automating performance studies of composed HPC applications that perform online data reduction and analysis. We describe Cheetah, a new framework for composing parametric studies on coupled applications, and Savanna, a runtime engine for orchestrating and executing campaigns of codesign experiments. Furthermore, this toolset facilitates understanding the impact of various factors such as process placement, synchronicity of algorithms, and storage versus compute requirements for online analysis of large data. Ultimately, we aim to create a catalog of performance results that can help scientists understand tradeoffs when designing next-generation simulations that make use of online processing techniques. We illustrate the design of Cheetah and Savanna, and present application examples that use this framework to conduct codesign studies on small clusters as well as leadership class supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A numerical study of process complexity in permafrost dominated regions

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Calibration↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗