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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Material Characterization and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

TPS↗

NASA NTRSDate not supplied

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

PATO↗

DOE OSTI2022

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Numerical Simulations of Noisy Variational Quantum Eigensolver Ansatz Circuits

This is a case study of the variational quantum eigensolver (VQE) method using numerical simulations to test the influence of noise on the accuracy of the underlying circuit ansatz. We investigate a computational chemistry application of VQE to calculate the electronic ground state and its energy for Sodium Hydride (NaH), a prototypical two-electron problem. Using a one-parameter ansatz derived from unitary coupled cluster (UCC) theory, we simulate the effects of noise on the energy expectation value and variance with respect to the ansatz parameter. These numerical simulations provide insights into the accuracy of the prepared quantum state and the efficiency of the classical optimizer that iteratively refines the ansatz. We conduct a comparative study between analytical results derived for the UCC ansatz in the absence of noise and the noisy numerical simulation results obtained using an isotropic depolarizing noise model for each gate. We also compare the relative increase in noise on logically equivalent UCC ansatz circuits generated by randomized compiling. Notably, we observe that the intrinsic variance in the energy due to the simplicity of the ansatz itself compares with the noise induced by the bare circuit.

Gowrishankar, Meenambika↗

NASA NTRS1983

Simulation of transport and removal processes of the Saharan dust

A planetary boundary layer (PBL) model has been modified to include a Saharan air layer containing the bulk of Saharan dust. The Saharan air layer is recognized as a deep mixed layer, which extends up to 4-6 km during hot summer months and is characterized by high potential temperature and high dust concentration. Microphysical processes of particle coagulation and sedimentation have been coupled with dynamic processes to simulate the evolution of dust particles. In agreement with observations, simulations indicate that dust particles greater than 10 micron radius can be produced and transported in the Sahara air layer across the Atlantic Ocean, while dust concentrations in the PBL over the Atlantic Ocean are smaller by a factor of two or more than in the Saharan air layer.

Lee, I.-Y.↗

DOE OSTI2026

Differentiable multiphase flow model for physics-informed machine learning in reservoir pressure management

Accurate subsurface reservoir pressure control is extremely challenging due to geological heterogeneity and multiphase fluid-flow dynamics. Predicting behavior in this setting relies on high-fidelity physics-based simulations that are computationally expensive. Yet, the uncertain, heterogeneous properties that control these flows make it necessary to perform many of these expensive simulations, which is often prohibitive. To address these challenges, we introduce a physics-informed machine learning workflow that couples a fully differentiable multiphase flow simulator, which is implemented in the DPFEHM framework with a convolutional neural network (CNN). The CNN learns to predict fluid extraction rates from heterogeneous permeability fields to enforce pressure limits at critical reservoir locations. By incorporating transient multiphase flow physics into the training process, our method enables more practical and accurate predictions for realistic injection-extraction scenarios compared to previous works. To speed up training, we pretrain the model on single-phase, steady-state simulations and then finetune it on full multiphase scenarios, which dramatically reduces the computational cost. We demonstrate that high-accuracy training can be achieved with fewer than three thousand full-physics multiphase flow simulations – compared to previous estimates requiring up to ten million. This drastic reduction in the number of simulations is achieved by leveraging transfer learning from much less expensive single phase simulations.

25 ENERGY STORAGE↗

NASA NTRS1999

Dynamical Signatures of Living Systems

One of the main challenges in modeling living systems is to distinguish a random walk of physical origin (for instance, Brownian motions) from those of biological origin and that will constitute the starting point of the proposed approach. As conjectured, the biological random walk must be nonlinear. Indeed, any stochastic Markov process can be described by linear Fokker-Planck equation (or its discretized version), only that type of process has been observed in the inanimate world. However, all such processes always converge to a stable (ergodic or periodic) state, i.e., to the states of a lower complexity and high entropy. At the same time, the evolution of living systems directed toward a higher level of complexity if complexity is associated with a number of structural variations. The simplest way to mimic such a tendency is to incorporate a nonlinearity into the random walk; then the probability evolution will attain the features of diffusion equation: the formation and dissipation of shock waves initiated by small shallow wave disturbances. As a result, the evolution never "dies:" it produces new different configurations which are accompanied by an increase or decrease of entropy (the decrease takes place during formation of shock waves, the increase-during their dissipation). In other words, the evolution can be directed "against the second law of thermodynamics" by forming patterns outside of equilibrium in the probability space. Due to that, a specie is not locked up in a certain pattern of behavior: it still can perform a variety of motions, and only the statistics of these motions is constrained by this pattern. It should be emphasized that such a "twist" is based upon the concept of reflection, i.e., the existence of the self-image (adopted from psychology). The model consists of a generator of stochastic processes which represents the motor dynamics in the form of nonlinear random walks, and a simulator of the nonlinear version of the diffusion equation which represents the mental dynamics. It has been demonstrated that coupled mental-motor dynamics can simulate emerging self-organization, prey-predator games, collaboration and competition, "collective brain," etc.

Zak, M.↗

DOE OSTI2022

Newton trust-region methods with primary variable switching for simulating high temperature multiphase porous media flow

Coupling multiphase flow with energy transport due to high temperature heat sources introduces significant new challenges since boiling and condensation processes can lead to dry-out conditions with subsequent re-wetting. The transition between two-phase and single-phase behavior can require changes to the primary dependent variables adding discontinuities as well as extending constitutive nonlinear relations to extreme physical conditions. Practical simulations of large-scale engineered domains lead to Jacobian systems with a very large number of unknowns that must be solved efficiently using iterative methods in parallel on high-performance computers. Performance assessment of potential nuclear repositories, carbon sequestration sites and geothermal reservoirs can require numerous Monte-Carlo simulations to explore uncertainty in material properties, boundary conditions, and failure scenarios. Due to the numerical challenges, standard NR iteration may not converge over the range of required simulations and require more sophisticated optimization method like trust-region. In this study, we use the open-source simulator PFLOTRAN for the important practical problem of the safety assessment of future nuclear waste repositories in the U.S. DOE geologic disposal safety assessment Framework. The simulator applies the PETSc parallel framework and a backward Euler, finite volume discretization. We demonstrate failure of the conventional NR method and the success of trust-region modifications to Newton’s method for a series of test problems of increasing complexity. Trust-region methods essentially modify the Newton step size and direction under some circumstances where the standard NR iteration can cause the solution to diverge or oscillate. Furthermore, we show how the Newton Trust-Region method can be adapted for Primary Variable Switching (PVS) when the multiphase state changes due to boiling or condensation. The simulations with high-temperature heat sources which led to extreme nonlinear processes with many state changes in the domain did not converge with NR, but they do complete successfully with the trust-region methods modified for PVS. This implementation effectively decreased weeks of simulation time needing manual adjustments to complete a simulation down to a day. Finally, we show the strong scalability of the methods on a single node and multiple nodes in an HPC cluster.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗

NASA NTRS2007

Future Climate Impacts of Direct Radiative Forcing Anthropogenic Aerosols, Tropospheric Ozone, and Long-lived Greenhouse Gases

Long-lived greenhouse gases (GHGs) are the most important driver of climate change over the next century. Aerosols and tropospheric ozone (O3) are expected to induce significant perturbations to the GHG-forced climate. To distinguish the equilibrium climate responses to changes in direct radiative forcing of anthropogenic aerosols, tropospheric ozone, and GHG between present day and year 2100, four 80-year equilibrium climates are simulated using a unified tropospheric chemistry-aerosol model within the Goddard Institute for Space Studies (GISS) general circulation model (GCM) 110. Concentrations of sulfate, nitrate, primary organic (POA) carbon, secondary organic (SOA) carbon, black carbon (BC) aerosols, and tropospheric ozone for present day and year 2100 are obtained a priori by coupled chemistry-aerosol GCM simulations, with emissions of aerosols, ozone, and precursors based on the Intergovernmental Panel on Climate Change (IPCC) Special Report on Emissions Scenario (SRES) A2. Changing anthropogenic aerosols, tropospheric ozone, and GHG from present day to year 2100 is predicted to perturb the global annual mean radiative forcing by +0.18 (considering aerosol direct effects only), +0.65, and +6.54 W m(sup -2) at the tropopause, and to induce an equilibrium global annual mean surface temperature change of +0.14, +0.32, and +5.31 K, respectively, with the largest temperature response occurring at northern high latitudes. Anthropogenic aerosols, through their direct effect, are predicted to alter the Hadley circulation owing to an increasing interhemispheric temperature gradient, leading to changes in tropical precipitation. When changes in both aerosols and tropospheric ozone are considered, the predicted patterns of change in global circulation and the hydrological cycle are similar to those induced by aerosols alone. GHG-induced climate changes, such as amplified warming over high latitudes, weakened Hadley circulation, and increasing precipitation over the Tropics and high latitudes, are consistent with predictions of a number of previous GCM studies. Finally, direct radiative forcing of anthropogenic aerosols is predicted to induce strong regional cooling over East and South Asia. Wintertime rainfall over southeastern China and the Indian subcontinent is predicted to decrease because of the increased atmospheric stability and decreased surface evaporation, while the geographic distribution of precipitation is also predicted to be altered as a result of aerosol-induced changes in wind flow.

Chen, Wei-Ting↗

NASA NTRS2017

Smectite Formation in the Presence of Sulfuric Acid: Implications for Acidic Smectite Formation on Early Mars

The excess of orbital detection of smectite deposits compared to carbonate deposits on the martian surface presents an enigma because smectite and carbonate formations are both favored alteration products of basalt under neutral to alkaline conditions. We propose that Mars experienced acidic events caused by sulfuric acid (H2SO4) that permitted phyllosilicate, but inhibited carbonate, formation. To experimentally verify this hypothesis, we report the first synthesis of smectite from Mars-analogue glass-rich basalt simulant (66 wt% glass, 32 wt% olivine, 2 wt% chromite) in the presence of H2SO4 under hydrothermal conditions (approximately 200 degC). Smectites were analyzed by X-ray diffraction, Mossbauer spectroscopy, visible and near-infrared reflectance spectroscopy and electron microprobe to characterize mineralogy and chemical composition. Solution chemistry was determined by Inductively Coupled Plasma Mass Spectrometry. Basalt simulant suspensions in 11-42 mM H2SO4 were acidic with pH less than or equal to 2 at the beginning of incubation and varied from acidic (pH 1.8) to mildly alkaline (pH 8.4) at the end of incubation. Alteration of glass phase during reaction of the basalt simulant with H2SO4 led to formation of the dioctahedral smectite at final pH approximately 3 and trioctahedral smectite saponite at final pH approximately 4 and higher. Anhydrite and hematite formed in the final pH range from 1.8 to 8.4 while natroalunite was detected at pH 1.8. Hematite was precipitated as a result of oxidative dissolution of olivine present in Adirondack basalt simulant. Formation of secondary phases, including smectite, resulted in release of variable amounts of Si, Mg, Na and Ca while solubilization of Al and Fe was low. Comparison of mineralogical and solution chemistry data indicated that the type of smectite (i.e., dioctahedral vs trioctahedral) was likely controlled by Mg leaching from altering basalt and substantial Mg loss created favorable conditions for formation of dioctahedral smectite. We present a model for global-scale smectite formation on Mars via acid-sulfate conditions created by the volcanic outgassing of SO2 in the Noachian and early Hesperian.

Smectite↗

DOE OSTI2026

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

DOE OSTI2021

Alpha particle driven Alfvénic instabilities in ITER post-disruption plasmas

Fusion-born alpha particles in ITER disruption simulations are investigated as a possible drive of Alfvénic instabilities. The ability of these waves to expel runaway electron (RE) seed particles is explored in the pursuit of a passive, inherent RE mitigation scenario. The spatiotemporal evolution of the alpha particle distribution during the disruption is calculated using the linearized Fokker–Planck solver CODION coupled to a fluid disruption simulation. These simulations are done in the limit of no alpha particle transport during the thermal quench, which can be seen as a most pessimistic situation where there is also no RE seed transport. Under these assumptions, the radial anisotropy of the resulting alpha population provides free energy to drive Alfvénic modes during the quench phase of the disruption. We use the linear gyrokinetic magnetohydrodynamic code LIGKA to calculate the Alfvén spectrum and find that the equilibrium is capable of sustaining a wide range of modes. The self-consistent evolution of the mode amplitudes and the alpha distribution is calculated utilizing the wave-particle interaction tool HAGIS. Intermediate mode number (n = 7–15, 22–26) toroidal Alfvén eigenmodes are shown to saturate at an amplitude of up to δB/B ≈ 0.1% in the spatial regimes crucial for RE seed formation. We find that the mode amplitudes are predicted to be sufficiently large to permit the possibility of significant radial transport of REs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

A New Efficient Quantitative Multi-component Phase Field: Lattice Boltzmann Model for Simulating Ti6Al4V Solidified Dendrite Under Forced Flow

Ti6Al4V is a widely used, multi-component alloy in additive manufacturing, during which the fluid flow in the molten pool significantly affects the solidified dendrites. To predict and further control the microstructure, modeling and simulating the microstructure evolution play a critical role. In this study, a newly developed, efficient, quantitative multi-component phase-field (PF) model is coupled with a lattice Boltzmann (LB) model to simulate Ti6Al4V solidified dendrite evolution under fluid flow. The accuracy and convergence behavior of the model is validated by the Gibbs–Thomson relation at the dendrite tip. Single and multiple two-dimensional (2D) equiaxed dendrite evolution cases under forced flow were simulated. Results show that the dendrite pattern is influenced remarkably by the fluid flow. Underlying mechanisms of the asymmetrical evolution are revealed by discussing the interaction among the flow, composition distribution and dendrite morphology, quantitatively. The dendrite kinetics are also derived, which ascertains the relationship between tip velocity and undercooling and inlet velocity and is the foundation for larger-scale simulation. We believe that the coupled quantitative multi-component PF–LB framework employed in this study helps in investigating the solidified dendrite morphology evolution in a deep and quantitate manner.

36 MATERIALS SCIENCE↗

NASA NTRS2020

Investigation of Atmospheric Boundary-Layer Effects on Launch-Vehicle Ground Wind Loads

A launch vehicle ground-wind-loads program is underway at the NASA Langley Transonic Dynamics Tunnel. The objectives are to quantify key aerodynamic and structural characteristics that impact the occurrence of large wind-induced oscillations of a launch vehicle when exposed to ground winds prior to launch. Of particular interest is the dynamic response of a launch vehicle when a von Kármán vortex street forms in the wake of the vehicle resulting in quasiperiodic lift and drag forces. Vehicle response to these quasiperiodic forces can become quite large when the frequency of vortex shedding nears that of a lowly-damped structural mode thereby exciting a resonant response. Wind approaching the vehicle can be characterized by a varying speed with height and turbulence content. The combination of both the varying speed and turbulence content is referred to herein as the atmospheric boundary-layer. The importance of the atmospheric boundary-layer upon launch vehicle wind-induced oscillation response has long been questioned, and its effects are not well understood. Although there are several facilities around the world dedicated to replicating atmospheric boundary layers, the development of such a boundary layer in a wind tunnel capable of producing flight-representative Reynolds numbers for aeroelastically-scaled launch vehicle models has only recently been accomplished. The NASA Langley Transonic Dynamics Tunnel is capable of simulating flight-representative Reynolds numbers of launch vehicles on the pad and is uniquely capable of replicating many fluid-structure scaling parameters typical of aeroelastic tests. Recent test efforts successfully developed representative atmospheric boundary-layers for three launch sites in the Transonic Dynamics Tunnel, thereby allowing all known aerodynamic and fluid-structure coupling parameters to be simultaneously simulated for those sites. Dynamic aeroelastically-scaled models representative of typical large launch vehicles were constructed for testing. Aeroelastic scaling includes matching geometry, mode shapes, reduced frequencies, damping, running mass ratios, and running stiffness ratios. The models were tested in smooth uniform flow and then immersed in the atmospheric boundary-layer for comparison of these effects. Dynamic data were acquired measuring unsteady pressure, acceleration, and base bending moment. It was discovered that peak dynamic loads resulting from resonant wind-induced oscillation response are similar when acquired in either smooth uniform flow or an atmospheric boundary-layer. This indicates that resonant lock-in events are minimally impacted by representative turbulence and/or wind profile. Alternately, nonresonant wind-induced oscillation response events are stronger when acquired in an atmospheric boundary-layer. This indicates that a lowly-damped structural response will increase when exposed to an increased magnitude of random excitation, which is consistent with historical comparisons. Loads created by the resonant response events were substantially stronger than those from the nonresonant response events. Therefore, if testing is done to simply identify worst-case conditions and load magnitude, then smooth uniform flow is likely an adequate test technique. However, if nonresonant response loads are of primary interest, then atmospheric boundary-layer simulation is required.

Thomas G Ivanco↗

DOE OSTI2025

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Dynamical attractors in contracting spacetimes dominated by kinetically coupled scalar fields

We present non-perturbative numerical relativity simulations of slowly contracting spacetimes in which the scalar field driving slow contraction is coupled to a second scalar field through an exponential non-linear σ model-type kinetic interaction. These models are important because they can generate a nearly scale-invariant spectrum of super-Hubble density fluctuations fully consistent with cosmic microwave background observations. We show that the non-linear evolution rapidly approaches a homogeneous, isotropic and flat Friedmann-Robertson-Walker (FRW) geometry for a wide range of inhomogeneous and anisotropic initial conditions. Ultimately, we find, the kinetic coupling causes the evolution to deflect away from flat FRW and towards a novel Kasner-like stationary point, but in general this occurs on time scales that are too long to be observationally relevant.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗