Search NASA⌕ Search

SEARCH · Search NASA

Results for “Effective generation time”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Implementation of Perturbation Theory and Sensitivity Capabilities in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor Multiphysics analysis application, jointly developed by Argonne and Idaho National Laboratories under the DOE-NE NEAMS program. This fiscal year, capabilities for reactivity and sensitivity evaluation using perturbation methods were implemented and verified. The First Order Perturbation Method (FOPT) was employed to compute reactivity worth resulting from small perturbations in input parameters, while the Generalized Perturbation Theory (GPT) was used to evaluate sensitivities of a range of response types, including reaction rate ratio, k-eigenvalue, neutron generation time, and effective delayed neutron fraction. These perturbation methods enable users to quantify how response quantities change due to a perturbation in a input parameter without explicitly performing an additional transport simulation for each perturbed state. In particular, the GPT formulation accounts for indirect effects arising from flux changes by solving generalized inhomogeneous equations, for which a Neumann series-based iterative solution method was developed and implemented in Griffin. The implemented reactivity and sensitivity evaluation capabilities were verified using two test problems: an infinite homogeneous system and a two-dimensional hexagonal core. The results showed excellent agreement with reference solutions obtained by a direct method based on finite difference approximation as well as GPT-based results from the PERSENT code, confirming the accuracy of both reactivity and sensitivity evaluations. Additionally, preliminary uncertainty quantification (UQ) results were obtained by combining the sensitivity values computed using GPT and external covariance data, demonstrating that the implemented sensitivity results can be reliably used for uncertainty calculations. To further demonstrate the generality and practical strength of the implementation, the sensitivity evaluation capability was successfully applied to the Empire microreactor with a geometrically complex design that poses significant modeling challenges. The results confirm that Griffin enables sensitivity evaluations even for irregular and highly heterogeneous reactor configurations, thereby establishing a foundation for UQ applications in advanced reactor designs and analyses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Preservation of kinetics parameters generated by Monte Carlo calculations in two-step deterministic calculations

The generation of accurate kinetic parameters such as mean generation time Λ and effective delayed neutron fraction β eff via Monte Carlo codes is established. Employing these in downstream deterministic codes warrants another step to ensure no additional error is introduced by the low-order transport operator when computing forward and adjoint fluxes for bilinear weighting of these parameters. Another complexity stems from applying superhomogenization (SPH) equivalence in non-fundamental mode approximations, where reference and low-order calculations rely on a 3D full core model. In these cases, SPH factors can optionally be computed for only part of the geometry while preserving reaction rates and K-effective, but the impact of such approximations on kinetics parameters has not been thoroughly studied. This paper aims at studying the preservation of bilinearly-weighted quantities in the Serpent–Griffin calculation procedure. Diffusion and transport evaluations of IPEN/MB-01, Godiva, and Flattop were carried out with the Griffin reactor physics code, testing available modeling options using Serpent-generated multigroup cross sections and equivalence data. Verifying Griffin against Serpent indicates sensitivities to multigroup energy grid selection and regional application of SPH equivalence, introducing significant errors; these were demonstrated to be reduced through the use of a transport method together with a finer energy grid.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Improvements to the Modeling of the TREAT Reactor and Experiments

This paper summarizes the latest improvements and lessons learned from the modeling and simulation of the transient test reactor at Idaho National Laboratory using the MAM- MOTH reactor physics application. MAMMOTH is a MOOSE-based, Finite Element Method application that specializes in the analysis of the spatial dynamics behavior of nuclear reactors. Since early 2018 several transient tests have been conducted at TREAT, thus providing the opportunity to apply and benchmark modern modeling and simulation tools. MAMMOTH was used to provide predictions of the power coupling factor between the core and the experiment for various experiments. Even though the power coupling factor predictions agree very well with the experimental data, within the bounds of the experimental uncertainty, one shortcoming was the underprediction of the total energy deposited in the core and in the experiment. Determination of the sources for this discrepancy are ongoing, but several key problems have been identified and resolved, thus providing valuable insights for future research. This paper discusses several of these lessons learned. First, the heat capacity data for the TREAT fuel has some significant problems due to limitations of the measurement techniques used circa 1960s. The sensitivity of the peak power and the total energy deposition to various representations of the heat capacity is approximately 5%. Second, the effects of the biological shield and thermal column on the modeling of the core are non-negligible, since they affect the mean generation time and the effective reflection of neutrons back into the core, which is suspected to be important during the core heat up. Matching the reactor period resolves the fact that the reduced spatial domain used in the MAMMOTH model underpredicts the mean generation time. The neutron reflection from these regions is marginally improved with the use of an albedo boundary condition. Third, modeling of the control rod movement with a multi-scheme method is introduced and its current limitations are exposed. Fourth, we explore the effects of using a homogenized model with Superhomogenization equivalence and how that differs from fully heterogeneous simulations. Finally, the energy condensation effects for this graphite core are significant. Solutions with 10 and 26 energy groups show the benefits of using a finer coarse group structure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Seismic Recovery of Acoustic Source Height During the OSIRIS-REx Sample Return Capsule Re-Entry

Controlled hypersonic reentries with well‐constrained trajectories are exceedingly rare, yet they provide valuable reference cases for validating geophysical source characterization methods applicable to bolides and space debris. The 24 September 2023 re‐entry of the OSIRIS‐REx Sample Return Capsule, the first interplanetary capsule to return over the contiguous United States in nearly two decades, presented such an opportunity. Here, we present seismic observations of this event recorded at 16 stations across three sites at ground distances of 7–20 km from the trajectory. Using eigenray tracing through ground‐to‐space atmospheric specifications, we invert picked acoustic‐to‐seismic travel times for effective shock‐generated acoustic source altitudes along the known trajectory. The method recovers source heights of 57.9–59.5 km with travel‐time residuals <1 s across all stations. The seismic‐derived source heights, when compared at their along‐track positions, fall within the range of previously published infrasound‐derived heights for the corresponding trajectory segment, with intermethod differences comparable with the timing‐based uncertainty. These results suggest that air‐to‐ground coupling does not introduce a resolvable systematic timing bias in the picked near‐field seismic arrivals, supporting seismic observations as a complementary tool for altitude estimation of bolides, re‐entry vehicles, and other hypersonic atmospheric sources.

Silber, Elizabeth Allaryce [Sandia National Labora↗
DOE OSTI2021

Fuel properties effects on wall-wetting and fuel-in-lube dilution under stochastic pre ignition prone conditions

Fuel-engine interactions have been shown to influence stochastic pre-ignition (SPI), however there is uncertainty on the root causes, participially if effects are primarily physical or chemical in nature. To better delineate the physical and chemical dependencies of fuel properties on SPI, this work studies interactions between aromatics, distillation, and particle matter index on SPI with a laser induced fluorescence (LIF) diagnostic to observe wall wetting. Fuels with high aromatic components or higher distillation temperatures tend to increase SPI events, but it is not well understood if the increased fuel dilution or a chemical aspect of this process is the primary mechanism for SPI. By employing the novel LIF in this work we measure fuel dilution in real-time in a running single cylinder research engine operated at SPI conditions. Three fuels were used in the present work, with the LIF diagnostic. Each fuel was operated at two injection timings, generating fuel specific effects on wall wetting. The LIF measurements are coupled with engine combustion measurements and SPI event counts. Results enable a deconvolution of fuel physical and chemical properties on SPI propensity, providing unique and novel insights into the dominant fuel effects affecting SPI.

Splitter, Derek↗
DOE OSTI2025

Uncertainty quantification and sensitivity analysis of a nuclear thermal propulsion reactor startup sequence

The research presented in this article describes progress in applying stochastic methods, uncertainty quantification, parametric studies, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated via the radiation transport code Griffin to simulate neutronics. Our goal is to develop a reduced-order (surrogate) model that can be rapidly sampled with perturbations to multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums (CDs), which is itself controlled by a hybrid proportional-integral-derivative (PID) controller actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time [Λ] and effective delayed neutron fraction [ β eff ] from a transient Griffin simulation executed via Griffin’s improved quasi-static solver to provide the kinetic parameters) as inputs to functions that control the CD rotation angle. By investigating numerous stochastic approaches, we developed a dual-purpose surrogate model of the NTP system, using polynomial regression in the Multiphysics Object-Oriented Simulation Environment (MOOSE) Stochastic Tools Module (STM). The trained model can be rapidly sampled while simultaneously perturbing various input parameters, such as coefficients on the PID control or temperature (directly affecting the neutron cross section). The surrogate model delivers accurate (within 5%) results at speeds orders of magnitude faster (minutes, not days of computational time) than the base model. Once the surrogate model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs as well as the effects of these inputs on the model output. For example, coefficients used in the PID control system may vary due to some type of physical interference, or uncertainty may exist in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters, and operational boundaries can be determined. The goal of this work is to support development of an advanced control system for operating CDs in a functioning NTP system. This work is a scoping study of the MOOSE STM.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Hierarchical Control of Utility-Scale Solar PV Plants for Mitigation of Generation Variability and Ancillary Service Provision

This paper presents a hierarchical control system to mitigate the variability of solar photovoltaic (PV) power plant and provide ancillary services to the electric grid without the need for additional non-solar resources. With coordinated management of each inverter in the system, the control system commands the power plant to proactively curtail a small fraction of its instantaneous maximum power potential, which gives the plant enough headroom to ramp up production from the overall power plant, for a service such as regulation reserve. This control system is practical for continuously changing cloud cover conditions in partially cloudy days. A case study from a site in Hawaii with one-second resolution solar irradiance data is used to verify the efficacy of the proposed control system. The proposed control algorithm is subsequently compared with the alternative control technology from the literature, the grouping control algorithm; the results show that the proposed hierarchical control system is over 10 times more effective in reducing generator mileage to support power fluctuations from solar PV power plants.

14 SOLAR ENERGY↗
DOE OSTI2020

Modification of terahertz emission spectrum using microfabricated spintronic emitters

Terahertz (THz) radiation with sub-millimeter wavelength falls in the gap between optical and radio frequencies. Conventional THz emitters do not intertwine with spin degrees of freedom. However, it was recently shown that broadband THz radiation can be efficiently created also by exploiting spin-based effects on ultrafast time scales. Here, we demonstrate the generation and control of THz radiation from microstructured spintronic THz emitters based on the inverse spin-Hall effect. Using time-domain THz spectroscopy, we compare the THz spectra of different stripe patterns made of Fe/Pt bilayers with a spectrum obtained from an extended Fe/Pt bilayer film. It is found that the THz spectrum can be altered by a proper choice of the microstructure dimensions. The experimentally observed spectra are interpreted in terms of a simplified multi-slit interference model, which captures the main experimental features. Our results pave the way for an efficient control of THz light emitted from magnetic heterostructures. This is a crucial step forward for the design and realization of directional THz sources.

Wu, Weipeng↗
DOE OSTI2022

Dynamic inflow and unsteady aerodynamics models for modal and stability analyses in OpenFAST

This work presents two aerodynamic models that are used to perform modal and stability analyses of wind turbines. Including aerodynamic effects in stability analyses is important to capture aerodynamic inertia, stiffness, and damping. The two models consist of a dynamic inflow model based on Øye’s model and an unsteady aerodynamics model based on a simplified Beddoes-Leishman dynamic stall model. Both models are written in a continuous-time state-space form, which is the form needed to perform linearizations within OpenFAST. In this article, we present the new models, perform some simple verifications against HAWC2 results, and then perform stability analyses. We investigate the contribution of the aerodynamic effects on the turbine modes, particularly the flapwise modes, which are the most affected. Because of the development presented in this work, Campbell diagrams accounting for dynamic stall and dynamic wake effects were generated for the first time using OpenFAST.

17 WIND ENERGY↗
DOE OSTI2020

Rupture Dynamics Simulations of Shallow Crustal Earthquakes on Reverse Slip Faults

The main goal of this study was to perform dynamic rupture simulations in order constrain the characterization of kinematic rupture models implemented in techniques for simulating strong ground motion for crustal earthquakes. First, we investigated the rupture process of the 2008 Iwate-Miyagi Nairiku earthquake by developing a dynamic rupture model on a reverse fault using a trial-and-error technique that produced a slip distribution and near-fault ground motion that matched the recorded ones. The simulations were performed in the frequency range 0-2 Hz, using a 3D staggered grid finite-difference method and a linear slip weakening friction law. Constrained by the observed slip distribution and consistent with dynamic rupture models, the derived kinematic rupture model of the Iwate-Miyagi Nairiku earthquake contains areas with relatively large slip rate, representing strong motion generation areas, set against lower amplitude heterogeneous background slip, and a relatively low slip rate in the weak zone of the top-most crust (upper 3km). Second, we performed rupture dynamics modeling to constrain shallow slip characterization in rupture models for strike-slip crustal earthquakes. The objective was to establish general rules about the characterization of slip rate function and slip at shallow depths and in the strong motion generation areas (SMGAs). The simulations of spontaneous rupture were performed in the frequency range 0-2.0Hz, using a 3D staggered-grid finite-difference method and a layer over half-space 1D crustal velocity model with a minimum shear-wave velocity of 2.8 km/s. In order to account for changes in material ductility and reduction of stress drop, observed in the shallow crust (upper 3-5 km), and the transition from ductile state to brittle state in the upper seismogenic zone, in our stress models we included a shallow weak zone (<4km). In this zone the stress drop was set to zero at the free surface and gradually increased with depth, while the slip weakening distance was set to 75 cm at the free surface and decreased to 50 cm at the base of the weak zone. As in the case of reverse faulting during the Iwate-Miyagi Nairiku earthquake, from these computations we found a systematic change in the shape of the slip-rate function, from Kostrov-type in the deeper part of the fault to a more symmetric cosine-type in the upper few kilometers, near the free-surface. Moreover, the average slip duration in the weak zone, with respect to slip duration in the deeper parts of the fault, increases by at most a factor of 1.5. We found a systematic and gradual change in the shape of the slip-rate function from Kostrov-type in the asperity areas (SMGAs), to more symmetric cosine-type in the upper few km near the free-surface, on the long period motion generation area (LMGA). The effective rise time in the LMGA, with respect to that in the SMGAs, increases by at most a factor of 2. Effective rise time is the time difference between the time at which the slip rate drops to a level that is equal to 25% of its peak and the onset time of the signal. In addition, the slip in the LMGA located above the SMGAs is almost the same as the one in the SMGAs, and about 1.5 times larger than the average slip.

58 GEOSCIENCES↗
DOE OSTI2021

Pulsar Based Timing for Grid Synchronization

Existing synchronization systems in the power grid, such as the global positioning system, are susceptible to temporary or permanent failures due to various unpredictable and uncontrollable factors such as cyber-attack and electromagnetic interferences, thus affecting the accuracy and reliability of generated timing signal. In this article, a pulsar astronomy-based timing system is proposed to provide an alternative synchronization signal. Further, this clock will offer significant security improvements to power grid applications, such as a wide-area monitoring system, which depends on a precise timing signal. The hardware and software frameworks are described in detail. First, a high-speed sampling hardware platform is designed to collect signals from radio telescopes. Then a periodic pulse extraction method with three steps is proposed to process the pulsar signal, including polyphase filterbanks, incoherent de-dispersion, and sliding window folding. Lastly, three experiments are conducted to verify the effectiveness of the frameworks. The generated pulsar timing pulse is presented, and the factors affecting its accuracy are also discussed. The analysis results demonstrate that the pulsar signals can provide high-accurate timing pulses for grid synchronization.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Enhancement of reactive oxygen species production by ultra-short electron pulses

The development of laser-driven accelerators-on-chip has provided an opportunity to miniaturize devices for electron radiotherapy delivery. Laser-driven accelerators produce highly time-compressed electron pulses, on the 100 fs to 1 ps scale. This delivers electrons at high peak power yet low average beam current compared with conventional delivery devices, which generate pulses of approximately 3 µs. The biophysical effects of this time structure, however, are unclear. Here, we use a Monte Carlo simulation approach to explore the effects of the electron beam time structure on the production of reactive oxygen species (ROS) in water. Our results show a power law increase in the generation of hydroxyl ions per deposited electron with decreasing pulse length over the pulse length range of 10 µs to 100 fs. Similar trends were observed for hydrogen peroxide, superoxide, hydroperoxyl, hydronium and solvated electrons. In practical terms, this indicates a fourfold increase in the efficiency of free radical production for sub-picosecond pulses, relative to that of conventional microsecond pulses, for the same number of deposited electrons.

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

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2024

Monopole Catalyzed Baryogenesis with a $\theta$ angle

Monopoles are generally expected in Grand Unified Theories (GUTs) where they can catalyze baryon decay at an unsuppressed rate by the Callan-Rubakov effect. For the first time, we show this catalysis effect can generate the observed baryon asymmetry at GeV scale temperatures. We study the minimal SU(5) GUT model and demonstrate that monopoles-fermion scattering with a $CP$-violating $\theta$-term leads to realistic baryogenesis even when $\theta\lesssim 10^{-10}$ is below the neutron EDM bound, potentially detectable in the future measurements. Our calculation also shows that to generate the observed baryon asymmetry, the abundance of the monopoles is below the current experiential bounds.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ponderomotive barriers in rotating mirror devices using static fields

Particularly for aneutronic fusion schemes, it is advantageous to manipulate the fuel species differently from one another and to expel ash promptly. The ponderomotive effect can be used to selectively manipulate particles. It is commonly a result of particle–wave interactions and has a complex dependence on the particle charge and mass, enabling species selectivity. If the plasma is rotating, e.g., due to E x B motion, the ponderomotive effect can be generated using static (i.e., time-independent) perturbations to the electric and magnetic fields, which can be significantly cheaper to produce than time-dependent waves. We propose that this feature can be particularly useful in rotating mirror machines where mirror confinement can be enhanced by rotation, both through centrifugal confinement and additionally through a ponderomotive interaction with a static azimuthal perturbation. We identify specific static perturbations that generate a ponderomotive barrier and other perturbations that can generate either a repulsive barrier or an attractive ponderomotive well, which can be used to attract particles of a certain species while repelling another. We identify the regimes in which the ponderomotive potential can enhance net plasma confinement and the regime in which plasma confinement is not enhanced. The viability of each of these effects is found to be dependent on the specifics of the rotation profile and the resultant dispersion relation in the rotating plasma.

Aneutronic fusion↗
DOE OSTI2023

Determining the Effects of Neutron Irradiation on the Structural Integrity of Additively Manufactured Heat Exchangers for Very Small Modular Reactor Applications, DOE Final Report (Project # 19-16980)

Auburn University (AU) teamed with the University of Missouri Research Reactor (MURR) and Kansas State University (KSU) to determine how to best use laser-powder bed fusion (L-PBF) additive manufacturing (AM) methods for generating radiation resistant nickel-based superalloys, Inconel alloy 625 and 718, for special purpose reactor (SPR) or very small modular reactor (vSMR) heat-exchangers (HeXs). Compact, conformal, and durable HeXs that are tolerant of extreme environments are needed for supporting the technical maturity of next-generation, portable compact reactors. AM is an enabler for realizing this new wave of HeXs – providing a means to make customizable hot and cold stream architectures with novel flow path geometries (e.g., tortuous channels with non-uniform, asymmetric cross-sections) and reduced layer-to-layer contact resistance (i.e., no separate bonding procedure required). AM further enables a more time/cost efficient means for fabricating SPRs by reducing the number of suppliers required for HeX assembly and allowing for on-site HeX fabrication. The project aim has been to better understand how neutron irradiation affects the microstructure and properties of additively manufactured nickel-based superalloys, to accelerate their safe, reliable use in the modular reactor industry. The major objective was to qualify/quantify the microstructure and microhardness of nickel-based superalloys (including Inconel 718 and 625) additively manufactured via the L-PBF process in the neutron-dosed (irradiated) and non-irradiated states over a course of 3 years. Effects of build orientation during L-PBF and post-AM heat treatments on neutron resistance, microstructure and mechanical properties were also investigated. Neutron damage mechanisms via hardening were characterized. This project combined subject-matter experts in AM, mechanical/microstructure metallic part characterization, and neutron irradiation, as well as unique assets and capabilities at AU and MURR at MU, to ensure project results translated to effectively addressing known gaps in nuclear science and engineering. Parts were fabricated using L-PBF systems readily available at AU. Specimens were then irradiated using MURR facilities; a manipulator equipped hot cell was also used to measure material hardness after dosing. MURR, a 10 MW, light-water nuclear reactor, is the largest, highest powered, highest-flux university owned research reactor in the U.S. The major findings in this project provide evidence that AM can serve as an alternative way to build structural components for future advanced small modular reactors using advanced materials like Inconel 625 (IN625) or Inconel 718 (IN718). After full spectrum neutron irradiation, vertically as-built AM IN625 samples were observed to display better resistance towards radiation-induced-hardening defects relative to traditionally machined metals. A Vickers microhardness tester, using settings of 1 kgf and dwell time of 15 seconds per indentation, was used to measure hardness in this study. The as-built, vertically printed samples experienced 1.2% of radiation hardening vs. 5.25% radiation hardening observed in wrought IN625. Another set of IN625 and IN718 samples were exposed to fast neutron irradiation. It was observed that IN718 showed more resistance towards radiation hardening compared to IN625 samples indicating IN718 had a better performance. Results showed that the IN718 samples experienced less change (-2.5 to 3.24 %) in microhardness. On the other hand, IN625 samples underwent more (0.9 to 7.21%) change in microhardness after fast neutron irradiation. AM IN625 samples were irradiated using an ion (proton) beam in cyclotron. The mechanical properties of AM samples post irradiation were compared with wrought samples. The irradiated region on the samples were tested using nano-hardness indention. It was observed that the beam current and time used in this study generated an annealing effect and thus reduced the hardness of the samples. The sum of the project results provide precious insight into how one may minimize radiation hardening in AM materials while maintaining material property constraints. Results should assist engineers in selecting an appropriate heat treatment for AM nickel-based superalloys for increased radiation resistance. Results should increase confidence levels for adopting AM for building nuclear reactor components which perform the same or better than conventionally manufactured components. Fast neutron irradiation testing provided an accelerated means of obtaining radiation effects without making materials radioactive and difficult to handle.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Incorporating Full Elastodynamic Effects and Dipping Fault Geometries in Community Code Verification Exercises for Simulations of Earthquake Sequences and Aseismic Slip (SEAS)

ABSTRACT Numerical modeling of earthquake dynamics and derived insight for seismic hazard relies on credible, reproducible model results. The sequences of earthquakes and aseismic slip (SEAS) initiative has set out to facilitate community code comparisons, and verify and advance the next generation of physics-based earthquake models that reproduce all phases of the seismic cycle. With the goal of advancing SEAS models to robustly incorporate physical and geometrical complexities, here we present code comparison results from two new benchmark problems: BP1-FD considers full elastodynamic effects, and BP3-QD considers dipping fault geometries. Seven and eight modeling groups participated in BP1-FD and BP3-QD, respectively, allowing us to explore these physical ingredients across multiple codes and better understand associated numerical considerations. With new comparison metrics, we find that numerical resolution and computational domain size are critical parameters to obtain matching results. Codes for BP1-FD implement different criteria for switching between quasi-static and dynamic solvers, which require tuning to obtain matching results. In BP3-QD, proper remote boundary conditions consistent with specified rigid body translation are required to obtain matching surface displacements. With these numerical and mathematical issues resolved, we obtain excellent quantitative agreements among codes in earthquake interevent times, event moments, and coseismic slip, with reasonable agreements made in peak slip rates and rupture arrival time. We find that including full inertial effects generates events with larger slip rates and rupture speeds compared to the quasi-dynamic counterpart. For BP3-QD, both dip angle and sense of motion (thrust versus normal faulting) alter ground motion on the hanging and foot walls, and influence event patterns, with some sequences exhibiting similar-size characteristic earthquakes, and others exhibiting different-size events. These findings underscore the importance of considering full elastodynamics and nonvertical dip angles in SEAS models, as both influence short- and long-term earthquake behavior and are relevant to seismic hazard.

Geochemistry & Geophysics↗