SEARCH · Search NASA

Results for “Multiphysics 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 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS2016

A Virtual Laboratory for the 4 Bed Molecular Sieve of the Carbon Dioxide Removal Assembly

Ongoing work to improve water and carbon dioxide separation systems to be used on crewed space vehicles combines sub-scale systems testing and multi-physics simulations. Thus, as part of NASA's Advanced Exploration Systems (AES) program and the Life Support Systems Project (LSSP), fully predictive COMSOL Multiphysics models of the Four Bed Molecular Sieve (4BMS) of the Carbon Dioxide Removal Assembly (CDRA) on the International Space Station (ISS) have been developed. This Virtual Laboratory is being used to help reduce mass, power, and volume requirements for exploration missions. In this paper we describe current and planned modeling developments in the area of carbon dioxide removal to support future missions as well as the resolution of anomalies observed in the ISS CDRA.

Coker, Robert↗

DOE OSTI2016

An Overview of Material Models in the BISON Fuel Performance Code

BISON is a modern finite element-based nuclear fuel performance code that has been under development at Idaho National Laboratory (INL) since 2009. The code is applicable to both steady and transient fuel behavior and supports the use of one-, two-, and three-dimensional meshes. BISON is designed to be a multi-fuels code. Although primarily applied to Light Water Reactor (LWR) fuel to date, the code has been used to analyze TRISO-coated particle fuel and metal fuel in rod and plate form, design and interpret fuel irradiation experiments and investigate novel fuel concepts. Nuclear fuel operates in an extreme environment that induces complex multiphysics phenomena occurring over distances ranging from inter-atomic spacing to meters, and times scales ranging from microseconds to years. To simulate this behavior requires a wide variety of material models which are often complex and nonlinear. Nearly all of the material models currently available in BISON are empirical, using correlations from existing fuel performance codes or material property libraries. Although these models perform adequately under many conditions, they have well-understood limitations. One obvious weakness is that the models are only valid in the domain of the associated experimental data, and extrapolating beyond this domain can result in substantial errors. A further weakness is that these models are often correlated in terms of imprecise quantities, such as burnup. Materials can, for example, reach a given burnup under very different circumstances, leading to very different material behavior. To overcome these weaknesses, efforts are underway with BISON to develop models based on a multi-scale approach, where material behavior is described with a series of interdependent microstructure based models. This paper provides an overview of the material models available in BISON. To keep the scope manageable, only LWR fuel will be considered. The first section focuses on UO{sub 2} fuel, identifying the models generally required and briefly describing those that have been implemented in the code. The second section repeats this approach for Zircaloy cladding. A final section provides examples of areas where multi-scale material modeling has been applied to BISON. BISON is a modern 2D and 3D nuclear fuel performance code that is applicable to multiple fuel forms. A wide variety of complex and nonlinear material models are required to simulate fuel behavior. Focusing on UO{sub 2} fuel and Zircaloy cladding (LWR fuel), the set of empirical models currently implemented in BISON have been reviewed. To overcome some of the deficiencies associated with empirical models, multi-scale material models are being developed. Two multi-scale models that have been used in BISON were reviewed. (authors)

36 MATERIALS SCIENCE↗

DOE OSTI2013

Realistic multisite lattice-gas modeling and KMC simulation of catalytic surface reactions: Kinetics and multiscale spatial behavior for CO-oxidation on metal (100) surfaces

A realistic molecular-level description of catalytic reactions on single-crystal metal surfaces can be provided by stochastic multisite lattice-gas (msLG) models. This approach has general applicability, although in this report, we will focus on the example of CO-oxidation on the unreconstructed fcc metal (100) or M(100) surfaces of common catalyst metals M = Pd, Rh, Pt and Ir (i.e., avoiding regimes where Pt and Ir reconstruct). These models can capture the thermodynamics and kinetics of adsorbed layers for the individual reactants species, such as CO/M(100) and O/M(100), as well as the interaction and reaction between different reactant species in mixed adlayers, such as (CO + O)/M(100). The msLG models allow population of any of hollow, bridge, and top sites. This enables a more flexible and realistic description of adsorption and adlayer ordering, as well as of reaction configurations and configuration-dependent barriers. Adspecies adsorption and interaction energies, as well as barriers for various processes, constitute key model input. The choice of these energies is guided by experimental observations, as well as by extensive Density Functional Theory analysis. Model behavior is assessed via Kinetic Monte Carlo (KMC) simulation. We also address the simulation challenges and theoretical ramifications associated with very rapid diffusion and local equilibration of reactant adspecies such as CO. These msLG models are applied to describe adsorption, ordering, and temperature programmed desorption (TPD) for individual CO/M(100) and O/M(100) reactant adlayers. In addition, they are also applied to predict mixed (CO + O)/M(100) adlayer structure on the nanoscale, the complete bifurcation diagram for reactive steady-states under continuous flow conditions, temperature programmed reaction (TPR) spectra, and titration reactions for the CO-oxidation reaction. Extensive and reasonably successful comparison of model predictions is made with experimental data. Furthermore, we discuss the possible transition from traditional mean-field-type bistability and reaction kinetics for lower-pressure to multistability and enhanced fluctuation effects for moderate- or higher-pressure. Behavior in the latter regime reflects a stronger influence of adspecies interactions and also lower diffusivity in the higher-coverage mixed adlayer. We also analyze mesoscale spatiotemporal behavior including the propagation of reaction diffusion fronts between bistable reactive and inactive states, and associated nucleation-mediated transitions between these states. This behavior is controlled by complex surface mass transport processes, specifically chemical diffusion in mixed reactant adlayers for which we provide a precise theoretical formulation. The msLG models together with an appropriate treatment of chemical diffusivity enable equation-free heterogeneous coupled lattice-gas (HCLG) simulations of spatiotemporal behavior. In addition, msLG + HCLG modeling can describe coverage variations across polycrystalline catalysts surfaces, pressure variations across catalyst surfaces in microreactors, and could be incorporated into a multiphysics framework to describe mass and heat transfer limitations for high-pressure catalysis. (C) 2013 Elsevier Ltd. All rights reserved.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Development of a BlueCRAB/MELCOR Framework for Supporting Realistic Mechanistic Source Term Calculations in Microreactors

Efforts are currently underway to deploy microreactor modeling and simulation tools to better support vendors and regulatory authorities in submitting and reviewing licensing applications. In particular, the Nuclear Regulatory Commission is expected to rely on the Comprehensive Reactor Analysis Bundle (BlueCRAB) multiphysics toolset in performing design- and beyond-design-basis accident analyses. In addition, the Nuclear Regulatory Commission has been using the MELCOR code to estimate mechanistic source terms during accidents. As MELCOR relies on isotopic inventory and reactor temperature/power evolution profiles during accident conditions—all of which can theoretically be obtained from BlueCRAB—the ultimate goal of this activity is to establish a common BlueCRAB-MELCOR framework. However, prior to the present research, BlueCRAB had never been used to calculate such quantities of interest at the full-core level. While there are many Monte Carlo (MC) codes capable of computing such quantities of interest, they are unable to readily account for multiphysics feedback. BlueCRAB allows for the coupling of different physics codes together to perform multiphysics-informed calculations. Therefore, the purpose of this fiscal year 2023 work is to investigate the feasibility and challenges of performing such calculations within BlueCRAB so as to generate the data that MELCOR relies on. To demonstrate the methodology, the proposed workflow was applied to a prototypical heat pipe-cooled microreactor model. To predict isotopic concentrations (taking into account the ac- cumulation of fission products during operation), the necessary microscopic cross sections were generated via OpenMC and tabulated with respect to temperature and burnup. Next, a recently developed capability in Griffin (the reactor physics application in BlueCRAB) was used to convert the OpenMC output format into the ISOXML format used by Griffin. A multiphysics microscopic depletion calculation that involved performing a coupled full-core, heterogeneous neutron trans- port and thermal calculation at each depletion step was conducted to deplete the core to end of life (EOL) conditions so as to provide both isotopics and the initial condition for the transient calculation. Following a brief null-transient to verify that the initial condition had been properly restarted and was indeed in thermal equilibrium, a heat pipe failure transient was simulated. Thus, the entire workflow of using BlueCRAB to generate MELCOR inputs, from cross-section generation to producing isotopic inventory and power/temperature evolution profiles during transients, is demonstrated. This report also details the identified gaps in the workflow and how they were (for the most part) addressed. Future work should focus on directly including MEL- COR into the workflow by performing a MELCOR calculation using the BlueCRAB-generated input data. In addition, the heat pipe reactor design should be improved so as to reflect more prototypical burnup characteristics at EOL.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2016

NREL Multiphysics Modeling Tools and ISC Device for Designing Safer Li-Ion Batteries

The National Renewable Energy Laboratory has developed a portfolio of multiphysics modeling tools to aid battery designers better understand the response of lithium ion batteries to abusive conditions. We will discuss this portfolio, which includes coupled electrical, thermal, chemical, electrochemical, and mechanical modeling. These models can simulate the response of a cell to overheating, overcharge, mechanical deformation, nail penetration, and internal short circuit. Cell-to-cell thermal propagation modeling will be discussed.

25 ENERGY STORAGE↗

NASA NTRS2011

Microwave Processing of Planetary Surfaces for the Extraction of Volatiles

In-Situ Resource Utilization will be necessary for sustained exploration of space. Volatiles are present in planetary soils, but water by far has the most potential for effective utilization. The presence of water at the lunar poles, Mars, and possibly on Phobos opens the possibility of producing LOX for propellant. Water is also a useful radiation shielding material , and valuable to replenish expendables (water and oxygen) required for habitation in space. Because of the strong function of water vapor pressure with temperature, heating soil effectively liberates water vapor by sublimation. Microwave energy will penetrate soil and heat from within much more efficiently than heating from the surface with radiant heat. This is especially true under vacuum conditions since the heat transfer rate is very low. The depth of microwave penetration is a strong function of the microwave frequency and to a lesser extent on soil dielectric properties. Methods for complex electric permittivity and magnetic permeability measurement are being developed and used for measurements of lunar soil simulants. A new method for delivery of microwaves deep into a planetary surface is being prototyped with laboratory experiments and modeled with COMSOL MultiPhysics. We are planning to set up a planetary testbed in a large vacuum chamber in the coming year. Recent results are discussed.

Ethridge, Edwin C.↗

DOE OSTI2021

Conduction Heat Transfer through Solid in Porous Materials: A Comparative Study by Finite-Element Simulations and Effective Medium Approximations

The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Exploring utilization options of heterogeneous architectures for multi-physics simulations

Top supercomputers in the TOP500 list have transitioned from homogeneous node architectures toward heterogeneous manycore nodes with accelerators and CPUs. These new architectures present significant challenges to developers of large-scale multiphysics applications, especially at DOE laboratories that have invested heavily in scalable MPI codes over decades. Much of these scientific application porting efforts for the new heterogeneous architectures are focused on running the computation on the accelerators, which usually comprise > 90% of the FLOPS of the system, typically using one core per GPU to offload work. Furthermore, we describe an approach to utilizing the remaining FLOPS on a heterogeneous machine by running a portion of the computation on the CPUs cooperatively with the GPU computation. Orthogonally, we evaluate using the additional cores to also offload work to GPUs, effectively oversubscribing the GPUs. We present a proof-of-concept implementation in ARES, a multiphysics ALE-AMR code at LLNL. ARES uses a portability layer, RAJA, which enables us to utilize the same source code for both the CPU and the GPU. We develop an approach to utilize both types of processors cooperatively in a heterogeneous system. Our implementation divides the work between the computing resources via domain decomposition, and utilizes all cores of the CPU and all of the GPUs on the node for computation. Load balancing is necessary to use the heterogeneous resources effectively. We present preliminary results on early delivery pre-Sierra machines at LLNL, showing up to 40% performance benefit of oversubscribing the GPUs (with caveats), and up to an 18% performance benefit of using the CPUs on the heterogeneous nodes for computing in addition to using the GPUs.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2014

Full wave simulation of waves in ECRIS plasmas based on the finite element method

This paper describes the modeling and the full wave numerical simulation of electromagnetic waves propagation and absorption in an anisotropic magnetized plasma filling the resonant cavity of an electron cyclotron resonance ion source (ECRIS). The model assumes inhomogeneous, dispersive and tensorial constitutive relations. Maxwell's equations are solved by the finite element method (FEM), using the COMSOL Multiphysics{sup ®} suite. All the relevant details have been considered in the model, including the non uniform external magnetostatic field used for plasma confinement, the local electron density profile resulting in the full-3D non uniform magnetized plasma complex dielectric tensor. The more accurate plasma simulations clearly show the importance of cavity effect on wave propagation and the effects of a resonant surface. These studies are the pillars for an improved ECRIS plasma modeling, that is mandatory to optimize the ion source output (beam intensity distribution and charge state, especially). Any new project concerning the advanced ECRIS design will take benefit by an adequate modeling of self-consistent wave absorption simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Dimming the Lights: 2D Simulations of Deflagrations of Hybrid C/O/Ne White Dwarfs Using FLASH

The dimmest and most numerous outlier of the Type Ia supernova population, the Type Iax event, is increasingly being found in the results of observational campaigns. There is currently no single accepted model to describe these events. This 2D study explores the viability of modeling Type Iax events as a hybrid C/O/Ne white dwarf progenitor undergoing a deflagration using the multiphysics software FLASH. This hybrid was created using the stellar evolution code MESA, and its C-depleted core and mixed structure have demonstrated lower yields than traditional C/O progenitors in previous deflagration-to-detonation studies. To generate a sample, 30 “realizations” of this simulation were performed, the only difference being the shape of the initial match head used to start the deflagration. Consistent with earlier work, these realizations produce the familiar hot dense bound remnant surrounded by sparse ejecta. Our results indicate that the majority of the star remains unburned (∼70%) and bound (>90%). Our realizations produce total ejecta yields on the order of 10 −2 –10 −1 M ☉ , ejected 56 Ni yields on the order of 10 −4 –10 −2 M ☉ , and ejecta kinetic energies on the order of 10 48 –10 49 erg. Compared to yields inferred from recent observations of the dimmest Type Iax events—SN 2007qd, SN 2008ha, SN 2010ae, SN 2019gsc, SN 2019muj, SN 2020kyg, and SN 2021fcg—our simulation produces comparable 56 Ni yields but too-small total yields and kinetic energies. Reignition of the remnant is also seen in some realizations.

Hydrodynamical simulations↗

DOE OSTI2015

A high-fidelity multiphysics model for the new solid oxide iron-air redox battery part I: Bridging mass transport and charge transfer with redox cycle kinetics

A high-fidelity two-dimensional axial symmetrical multi-physics model is described in this paper as an effort to simulate the cycle performance of a recently discovered solid oxide metal-air redox battery (SOMARB). The model collectively considers mass transport, charge transfer and chemical redox cycle kinetics occurring across the components of the battery, and is validated by experimental data obtained from independent research. In particular, the redox kinetics at the energy storage unit is well represented by Johnson-Mehl-Avrami-Kolmogorov (JIVIAK) and Shrinking Core models. The results explicitly show that the reduction of Fe3O4 during the charging cycle limits the overall performance. Distributions of electrode potential, overpotential, Nernst potential, and H-2/H2O-concentration across various components of the battery are also systematically investigated. (C) 2015 Elsevier B.V. All rights reserved.

Jin, XF↗

DOE OSTI2024

Coupling thermal energy storage with a thermally anisotropic building envelope for building demand-side management across various US climate conditions

Here, the thermally anisotropic building envelope (TABE) is a novel active building envelope that enhances energy efficiency and thermal comfort in buildings by transferring heat and cold between building envelopes and hydronic loops. When coupled with thermal energy storage (TES) units, the TABE + TES enables the storage of both heat and cold energy captured by the TABE roof or exterior walls. This stored energy can be later released by the TABE floor for indoor heating and cooling, benefiting both the grid and the end user. This paper evaluates the merits of TABE + TES for building demand-side management across various US climate conditions, focusing on peak load shaving, annual energy savings, and cost savings under time-of-use (TOU) electric rate schedules. Simulations were conducted by integrating time-of-day–informed, rule-based control strategies in MATLAB, TABE components and TES units in COMSOL Multiphysics, and whole-building energy analysis in EnergyPlus. A case study using the US Department of Energy’s prototype single-family detached house model in Birmingham, Alabama; Los Angeles, California; Oak Ridge, Tennessee; and Denver, Colorado, showed that the TABE + TES system achieved (1) 70 % peak load shaving in Los Angeles and Denver and 20 % in Birmingham and Oak Ridge; (2) significant peak electricity savings of 351–497 kWh, reducing peak energy consumption by 38 %–78 %; and (3) annual heating cost savings of 0.79 $\$$/m2–1.17 $\$$/m 2 and cooling cost savings of 0.60 $\$$/m 2 –1.17 $\$$/m 2 using a normal utility rate or low-TOU rate. The benefits of employing the TABE + TES system are even more significant under high TOU rates.

25 ENERGY STORAGE↗

DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗

DOE OSTI2015

Computational fluid dynamics modeling of high level waste plug formation - 15228

At Hanford, an extensive network of pipelines traversing several miles is used to transfer the high level radioactive nuclear waste from tanks to the treatment facilities. During transfer operations, however, there is a potential risk for the transport lines to plug, causing significant delays, increasing operation costs and creating hazardous conditions for personnel and the environment. In this study we investigate the phenomena of plug formation in the pipelines due to solids depositing during the waste transfer process. A three-dimensional (3-D) Computational Fluid Dynamics (CFD) model has been developed using Comsol Multiphysics 4.3b by Florida International University (FIU) to evaluate the factors that can contribute to pipeline plug formation. Several virtual scenarios are developed simulating the settling mechanics as a function of particle size, solids volume fraction and solids density. The numerical results are validated with the available experimental results and the established critical velocity correlations. The 3-D numerical results were also compared with the previous 2-D numerical studies to understand the trade-off between the two studies in terms of computing speed and numerical accuracy. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Magnetostrictive Ultrasonic Waveguide Transducer for In-Pile Thermometry

Real-time and reliable temperature measurement on nuclear fuels is crucial for the safe operation of existing pressurized-water reactors and future advanced nuclear reactors. Magnetostrictive materials deform when subjected to a magnetic field or exhibit magnetization variation when stressed. Here, based on these properties, this study prototyped an ultrasonic thermometer (UT) consisting of a magnetostrictive waveguide, a dc coil providing appropriate magnetic biasing, and an ac coil generating an acoustic impulse and detecting the resulting acoustic echoes. By tracking the time of flight between the excitation pulse and the echoes, the UT can potentially detect nuclear fuel cladding temperature from a long distance. In this study, magnetostrictive iron–gallium alloys, or Galfenol, were selected as the waveguide due to their large magnetostriction, superior temperature survivability, excellent radiation resilience, and high mechanical robustness. A multiphysics finite-element model considering electrical, magnetic, and mechanical dynamics in the magnetostrictive UT was then developed. The model exhibited an error of 0.24% in time-of-flight simulation and, therefore, enabled computer-aided design and guided signal processing. Between room temperature and 120 °C, the new Galfenol-based UT exhibits a linear sensitivity of 162.8×10 ₋6 °C ₋1 , which is 51.7% higher than a previous magnetostrictive UT based on iron–cobalt–vanadium alloys.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗

DOE OSTI2017

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: Coordinate Multiscale U 3 Si 2 Modeling

Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced accident tolerance when compared to traditional \uo~fuel zircaloy clad fuel rods. One of the potential replacement fuels is uranium silicide (\usi) for its higher thermal conductivity and uranium density. The lower melting temperature is of potential concern during postulated accident conditions. Another disadvantage for \usi~ is the lack of experimental data under power reactor conditions. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022~\cite{bragg-sitton_2014} multiscale multiphysics modeling approaches have been used to provide insight into these materials. \\ \noindent The purpose of this letter report is to highlight the multiscale modeling effort for \usi~fuel as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuel (ATF) High Impact Problem (HIP). This fiscal year two new models for \usi~fuel have been incorporated into the BISON fuel performance code~\cite{Williamson2012} based upon lower length scale simulations using the GRASS-SST rate theory code. The models are briefly described in the following sections.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗