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 991 records · Page 55

Determining critical points to control electric lighting to meet circadian lighting requirements and minimize energy use

Designing electric lighting systems to meet circadian lighting requirements may raise light levels and consequently energy use compared to existing practices. To reduce energy use, electric lighting can be controlled to be dimmed or turned off when sufficient daylight levels are available in space. This requires input from one or a few critical measurement points. However, it is unclear how critical points can be determined to ensure that all occupants receive the needed light levels while reducing electric lighting energy. This paper discusses three approaches for selecting critical points and utilizes annual daylight simulations modified to account for sky spectra, and coupled with spectral electric lighting simulations. Among the three evaluated approaches, the use of continuous daylight autonomy (modified to use EML measured at eye positions) is helpful for estimating electric lighting energy for dimmable electric lighting systems, and for identifying energy-saving strategies.

Abboushi, Belal K.↗

Computational Modeling and Experimental Characterization of Martensitic Transformations in Nicoal for Self-Sensing Materials

Fundamental changes to aero-vehicle management require the utilization of automated health monitoring of vehicle structural components. A novel method is the use of self-sensing materials, which contain embedded sensory particles (SP). SPs are micron-sized pieces of shape-memory alloy that undergo transformation when the local strain reaches a prescribed threshold. The transformation is a result of a spontaneous rearrangement of the atoms in the crystal lattice under intensified stress near damaged locations, generating acoustic waves of a specific spectrum that can be detected by a suitably placed sensor. The sensitivity of the method depends on the strength of the emitted signal and its propagation through the material. To study the transition behavior of the sensory particle inside a metal matrix under load, a simulation approach based on a coupled atomistic-continuum model is used. The simulation results indicate a strong dependence of the particle's pseudoelastic response on its crystallographic orientation with respect to the loading direction and suggest possible ways of optimizing particle sensitivity. The technology of embedded sensory particles will serve as the key element in an autonomous structural health monitoring system that will constantly monitor for damage initiation in service, which will enable quick detection of unforeseen damage initiation in real-time and during onground inspections.

Wallace, T. A.↗

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗

Multi-Physics Coupled Seismic Safety Analysis of Molten-Salt Reactors

Pool-type Molten Salt Reactors (MSRs) are a promising advanced nuclear reactor concept relying on passive physical behavior to provide increased safety characteristics. The high heat capacity, high boiling point, unpressurized liquid salt contains the nuclear fuel and as it passes through the vessel achieves criticality and produces nuclear power. The salt then carries away the heat produced along with the delayed neutron precursors towards heat exchangers and primary pumps before entering the core. This forms a highly coupled multiphysics problem between neutronics and fluid flow, which makes modeling these reactors challenging. The safety/licensing case heavily relies on modeling and simulation to prove their safety, and significant effort has been expanded by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to develop the appropriate simulation tools. Their safety/licensing will require a thorough analysis of their behavior during earthquake-based transients. Earthquakes are a transient condition commonly experienced in many regions, and are safely survived by dozens of reactors every year. Nonetheless, significant conservatisms were introduced with limited predictive modeling and simulation available. The first-of-a-kind capability developed by this research will allow for engineers to analyze ranges of design-basis and beyond-design basis accident scenarios to assess the integrity of the reactor vessel and surrounding system, possibly allowing for limiting the costly conservatisms in the design. This projects uses MASTODON, Griffin, and Pronghorn in the NEAMS ecosystem to model the coupled multiphysics problem posed by an earthquake in a molten salt reactor. It marches through a number of number of coupling schemes, from uncoupled simulations to tightly two-way coupled simulations. The appropriate level of coupling for these simulations will then be determined based on the accuracy in the quantities of interest (QoI) and computational effort involved in each scheme. The QoIs chosen for this paper include the mass flow rate across the heat exchanger (for Pronghorn simulations), power output of the reactor (for Griffin simulations), and the maximum Von Mises stress (for MASTODON simulations). This summary paper reports on the current progress of this project, with 2D tightly coupled multiphysics results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BISON fuel performance modeling optimization for experiment X447 and X447A using axial swelling and cladding strain measurements

With the recent need to qualify new reactor designs such as the Versatile Test Reactor (VTR), fuel performance calculations need to be performed to determine safety criteria of the proposed designs. In order to validate the fuel performance results obtained by a fuel performance code, BISON, for new reactor designs, legacy fuel from EBR-II and FFTF MFF with Post -Irradiation Examination (PIE) data need to be used as validation cases to benchmark models. Here in this work, BISON has been paired with the Fuels Irradiation & Physics Database (FIPD) and IFR Materials Information System (IMIS) to supply PIE data for comparison with simulations of EBR-II experiments X447/X447A. X447/X447A were assessed by implementing models for Fuel Cladding Chemical Interaction (FCCI) within BISON and optimizing the friction coefficient between the fuel surface and the cladding, the anisotropic swelling factor, and the HT9 first thermal creep scalar (which scales the first term in the HT9 creep equation) to best match the PIE axial fuel swelling height and cladding profilometry for all pins in X447/X447A. The optimal values were found using a generic algorithm developed to select different values for the three parameters until end criteria was met and error couldn’t be reduced further. The BISON-simulated cladding profilometry was evaluated using Standard Error of the Estimate (SEE) to account for the profile shape of the cladding profilometry. Optimal values for the friction coefficient, anisotropic fuel swelling factor, and HT9 first thermal creep scalar were found to best fit the BISON simulation results to the PIE measurements found in IMIS and FIPD. Improvements to current models are suggested to account for the underprediction of fuel swelling at low burnups and the overprediction of fuel swelling at higher burnups observed for the axial fuel swelling height. Although two pins in EBR-II X447/X447A (DP70 and DP75) were known to fail due to FCCI, none of the pins simulated in BISON reached a cumulative damage fraction (CDF) above 0.008 with FCCI correlations coupled in the BISON simulations. The error estimate generated for all pins in X447/X447A using optimal values was 209 µm, which is deemed acceptable.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Zero-Field NMR and Millitesla-SLIC Spectra for >200 Molecules from Density Functional Theory and Spin Dynamics

NMR is usually performed at magnetic fields of 1 T and above to obtain sufficient sensitivity and spectral dispersion to identify chemicals based on chemical shifts and J couplings. At lower fields, the advent of hyperpolarization technologies and sensitive detectors can address sensitivity concerns. However, it remains disputed whether spectral signatures at zero and ultra-low fields are sufficient for chemical identification. Here, we report an all–electron DFT-based batch calculation of J-coupling constants, which are used to generate J coupling NMR spectra at zero field and 6.5 mT for over 200 small molecules. In the developed computational tool chain, we first used the all-electron FHI-aims code to calculate the molecular J couplings and chemical shifts. We then fed the calculated NMR parameters into the NMR simulation package SPINACH to simulate both heteronuclear J coupling spectra at zero-field, and homonuclear J coupling spectra as spin-lock induced crossing (SLIC) spectra at ultra-low field (6.5 mT). The resulting spectra demonstrate that zero and ultra-low field NMR spectra can represent unique identifiers of chemical structure for small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulated Space Radiation and Weightlessness: Vascular-Bone Coupling Mechanisms to Preserve Skeletal Health

Weightlessness causes a cephalad fluid shift and reduction in mechanical stimulation, adversely affecting both cortical and trabecular bone tissue in astronauts. In rodent models of weightlessness, the onset of bone loss correlates with reduced skeletal perfusion, reduced and rarified vasculature and lessened vasodilation, which resembles blood-bone symbiotic events that can occur with fracture repair and aging. These are especially serious risks for long term, exploration class missions when astronauts will face the challenge of increased exposure to space radiation and abrupt transitions between different gravity environments upon arrival and return. Previously, we found using the mouse hindlimb unloading model and exposure to heavy ion radiation, both disuse and irradiation cause an acute bone loss that was associated with a reduced capacity to produce bone-forming osteoblasts from the bone marrow. Together, these findings led us to hypothesize that exposure to space radiation exacerbates weightlessness-induced bone loss and impairs recovery upon return, and that treatment with anti-oxidants may mitigate these effects. The specific aims of this recently awarded grant are to: AIM 1 Determine the functional and structural consequences of prolonged weightlessness and space radiation (simulated spaceflight) for bone and skeletal vasculature in the context of bone cell function and oxidative stress. AIM 2 Determine the extent to which an anti-oxidant protects against weightlessness and space radiation-induced bone loss and vascular dysfunction. AIM 3 Determine how space radiation influences later skeletal and vasculature recovery from prolonged weightlessness and the potential of anti-oxidants to preserve adaptive remodeling.

Alwood, J. S.↗

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Approach to Integrate Global-Sun Models of Magnetic Flux Emergence and Transport for Space Weather Studies

The Sun lies at the center of space weather and is the source of its variability. The primary input to coronal and solar wind models is the activity of the magnetic field in the solar photosphere. Recent advancements in solar observations and numerical simulations provide a basis for developing physics-based models for the dynamics of the magnetic field from the deep convection zone of the Sun to the corona with the goal of providing robust near real-time boundary conditions at the base of space weather forecast models. The goal is to develop new strategic capabilities that enable characterization and prediction of the magnetic field structure and flow dynamics of the Sun by assimilating data from helioseismology and magnetic field observations into physics-based realistic magnetohydrodynamics (MHD) simulations. The integration of first-principle modeling of solar magnetism and flow dynamics with real-time observational data via advanced data assimilation methods is a new, transformative step in space weather research and prediction. This approach will substantially enhance an existing model of magnetic flux distribution and transport developed by the Air Force Research Lab. The development plan is to use the Space Weather Modeling Framework (SWMF) to develop Coupled Models for Emerging flux Simulations (CMES) that couples three existing models: (1) an MHD formulation with the anelastic approximation to simulate the deep convection zone (FSAM code), (2) an MHD formulation with full compressible Navier-Stokes equations and a detailed description of radiative transfer and thermodynamics to simulate near-surface convection and the photosphere (Stagger code), and (3) an MHD formulation with full, compressible Navier-Stokes equations and an approximate description of radiative transfer and heating to simulate the corona (Module in BATS-R-US). CMES will enable simulations of the emergence of magnetic structures from the deep convection zone to the corona. Finally, a plan will be summarized on the development of a Flux Emergence Prediction Tool (FEPT) in which helioseismology-derived data and vector magnetic maps are assimilated into CMES that couples the dynamics of magnetic flux from the deep interior to the corona.

Mansour, Nagi N.↗

Monte Carlo solution of Boltzmann equation for a simple model of highly nonequilibrium diatomic gases - Translational rotational energy relaxation

Theoretical results pertaining to internally excited translational-rotational energy relaxation in a spatially uniform diatomic gas far removed from solid boundaries are obtained by solving the Boltzmann equation by means of the Monte Carlo direct simulation method. The analysis is based on calculations involving three different types of initial conditions: equilibrium, nonequilibrium-equipartition (i.e., equipartition is satisfied, but the distributions are perturbed), and nonequilibrium-nonequipartition (i.e., both equipartition and the distributions are perturbed). Results of monatomic-gas simulations are also included to facilitate comparisons with the coupled translational-rotational relaxation simulations, and some simulations for a normal shock-wave structure are briefly examined. The results show that: (1) single-step transitions are the significant mechanisms of intermodal energy transfer; (2) translational-rotational transitions are coupled most efficiently for low-lying states of rotationally excited molecules and least efficiently for highly rotationally excited molecules; and (3) relaxation occurs via a successive set of distributions that are not Maxwell-Boltzmann (nonlocal Maxwellian).

Yoshikawa, K. K.↗

Design Power Control Strategies of Grid-Forming Inverters for Microgrid Application

This paper develops and compares two control schemes in the application control layer of a non-phase-locked loop (non-PLL) grid-forming (GFM) inverter to gain insight and understanding into how the two schemes affect the dynamic responses of GFM inverters and the transition operation of microgrids. The first scheme adopts power tracking based on an outer current loop in grid-connected mode and droop control in islanded mode, and the second uses droop control in both grid-connected and islanded modes. Analytical study is developed to compare the performance of these two strategies from various aspects, including fundamental differences, transition operation, power tracking, and P-Q capability with low point-of-common-coupling (PCC) voltages. Extensive simulation testing is carried out to validate the control performance of these two control strategies, and the simulation testing confirms the findings of the analytic study. Finally, the second strategy is recommended because of its superior control performance and ease of implementation. The analysis and results are useful in developing reliable control schemes for non-PLL GFM inverters because increasing numbers of inverters will work as non-PLL grid-forming sources in future grids because of their improved stability and reliability.

grid-forming↗

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗

A new dynamic zOnal model with air-diffuser (DOMA) - Application to thermal comfort prediction

A new Dynamic zOnal Model with Air-diffuser (DOMA) was developed. Several case studies were investigated and tested to evaluate and validate this program using measurement data. This new model was integrated into a TRaNsient SYstems Simulation program library and coupled with the multi-zone thermal model. The DOMA/TRNSYS coupled model was then used to predict room temperature distribution over an entire day of a single-zone building. The results show that increasing the heating outputs of the electric floor system, for example, from 75 to 200 W/m 2 , would not effectively improve the indoor thermal comfort, since the thermostat will reach the set point first and then turn off the system before the room gets enough heat and reach a comfortable level. This indicates the importance of selecting an appropriate location and set point for the thermostat when using a floor heating system. This potential thermal comfort issue can only be identified through the two-node model with a dynamic zonal model rather than the conventional PMV model, which thus suggests that for optimizing indoor thermal comfort of a building equipped with a time-sensitive control strategy and/or HVAC system, the TSENS results obtained from the two-node model integrated with DOMA are more appropriate than PMVs.

Construction & Building Technology↗

A Hybrid CFD/Engineering Model Tool for Lunar Lander Surface Erosion Prediction

Plume-Surface Interaction (PSI) between lander engine plumes and landing area regolith poses risks to the landers and space exploration missions through view obscuration and high-energy ejecta impacts from eroded particle clouds and vehicle tilting from plume-induced craters. To address these risks, NASA MSFC Fluid Dynamics Branch (ER42) has developed a cascade of predictive simulation capabilities for PSI, including the recent development of a hybrid CFD/engineering model tool that allows for rapid PSI simulations with a descending/ascending vehicle and an eroding surface. The development of this fully-coupled, moving body/eroding surface simulation capability is detailed along with a demonstration of the tool’s capabilities on the Apollo 12 Lunar Module landing. Predictions of eroded mass flow rates as a function of time are verified against Apollo 12 flight data, from which the developed, viscous erosion model has been calibrated.

Plume Surface Interaction↗

A Hybrid CFD/Engineering Model Tool for Lunar Lander Surface Erosion Prediction

Plume-Surface Interaction (PSI) between lander engine plumes and landing area regolith poses risks to the landers and space exploration missions through view obscuration and high-energy ejecta impacts from eroded particle clouds and vehicle tilting from plume-induced craters. To address these risks, NASA MSFC Fluid Dynamics Branch (ER42) has developed a cascade of predictive simulation capabilities for PSI, including the recent development of a hybrid CFD/engineering model tool that allows for rapid PSI simulations with a descending/ascending vehicle and an eroding surface. The development of this fully-coupled, moving body/eroding surface simulation capability is detailed along with a demonstration of the tool’s capabilities on the Apollo 12 Lunar Module landing. Predictions of eroded mass flow rates as a function of time are verified against Apollo 12 flight data, from which the developed, viscous erosion model has been calibrated.

Plume Surface Interaction↗

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗