Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Full Matrix Capture”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Elastic Model Transitions: A Hybrid Approach Utilizing Quadratic Inequality Constrained Least Squares (LSQI) and Direct Shape Mapping (DSM)

A method for transitioning linear time invariant (LTI) models in time varying simulation is proposed that utilizes a hybrid approach for determining physical displacements by augmenting the original quadratically constrained least squares (LSQI) algorithm with Direct Shape Mapping (DSM) and modifying the energy constraints. The approach presented is applicable to simulation of the elastic behavior of launch vehicles and other structures that utilize discrete LTI finite element model (FEM) derived mode sets (eigenvalues and eigenvectors) that are propagated throughout time. The time invariant nature of the elastic data presents a problem of how to properly transition elastic states from the prior to the new model while preserving motion across the transition and ensuring there is no truncation or excitation of the system. A previous approach utilizes a LSQI algorithm with an energy constraint to effect smooth transitions between eigenvector sets with no requirement that the models be of similar dimension or have any correlation. This approach assumes energy is conserved across the transition, which results in significant non-physical transients due to changing quasi-steady state energy between mode sets, a phenomenon seen when utilizing a truncated mode set. The computational burden of simulating a full mode set is significant so a subset of modes is often selected to reduce run time. As a result of this truncation, energy between mode sets may not be constant and solutions across transitions could produce non-physical transients. In an effort to abate these transients an improved methodology was developed based on the aforementioned approach, but this new approach can handle significant changes in energy across mode set transitions. It is proposed that physical velocities due to elastic behavior be solved for using the LSQI algorithm, but solve for displacements using a two-step process that independently addresses the quasi-steady-state and non-steady-state contributions to the elastic displacement. For structures subject to large external forces, such as thrust or atmospheric drag, it is imperative to capture these forces when solving for elastic displacement. To simplify the mathematical formulation, assumptions are made regarding mass matrix normalization, constant external forcing, and constant viscous damping. These simplifications allow for direct solutions to the quasi-steady-state displacements through a process titled Direct Shape Mapping. DSM solves for the displacements using the eigenvalues of the elastic modes and the external forcing and returns a set of elastic displacements dictated by the eigenvectors of the post-transition mode set. For the non-steady-state contributions to displacement we formulate a LSQI problem that is constrained by energy of the non-steady state terms. The contributions from the quasi-steady-state and non-steady state solutions are then combined to obtain the physical displacements associated with the new set of eigenvectors. Results for the LSQI-DSM approach show significant reduction/complete removal of transients across mode set transitions while maintaining elastic motion from the prior state. For time propagation applications employing discrete elastic models that need to be transitioned in time and where running with full a full mode set is not feasible, the method developed offers a practical solution to simulating vehicle elasticity.

Hannan, Mike R.↗
DOE OSTI2022

Spin-phonon dispersion in magnetic materials

Microscopic coupling between the electron spin and the lattice vibration is responsible for an array of exotic properties from morphic effects in simple non magnets to magnetodielectric coupling in multiferroic spinels and hematites. Traditionally, a single spin–phonon coupling constant is used to characterize how effectively the lattice can affect the spin, but it is hardly enough to capture novel electromagnetic behaviors to the full extent. Here, we introduce a concept of spin–phonon dispersion to project the spin moment change along the phonon crystal momentum direction, so the entire spin change can be mapped out. Different from the phonon dispersion, the spin–phonon dispersion has both positive and negative frequency branches even in the equilibrium ground state, which correspond to the spin enhancement and spin reduction, respectively. Our study of bcc Fe and hcp Co reveals that the spin force matrix, that is, the second-order spatial derivative of spin moment, is similar to the vibrational force matrix, but its diagonal elements are smaller than the off-diagonal ones. This leads to the distinctive spin–phonon dispersion. Furthermore, the concept of spin–phonon dispersion expands the traditional Elliott-Yafet theory in nonmagnetic materials to the entire Brillouin zone in magnetic materials, thus opening the door to excited states in systems such as CoF 2 and NiO, where a strong spin-lattice coupling is detected in the THz regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

Delaminations Investigated With Ultrasonic Spectroscopy

A previous study suggested that the ultrasonic spectroscopy technique identified possible disbonds or delaminations in polymer matrix composite (PMC) rings sectioned from flywheel rotors (ref. 1). These results went unsubstantiated by other nondestructive evaluation (NDE) methods. To explain the results, PMC rings were further investigated with ultrasonic spectroscopy (ref. 2). The ultrasonic spectroscopy system utilizes a continuous-swept sine waveform as the input. After the swept sine wave traverses the material, the captured waveform is subjected to two fast Fourier transforms. The second fast Fourier transform along with equalization of the frequency spectrum, allows for evaluation of the fundamental resonant frequency. The full-thickness resonance, the resonance corresponding to the location of the intentional disbond, and the frequency spectrum were examined in an effort to characterize the sensitivity of the NDE method to various delamination conditions.

Cosgriff, Laura M.↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2023

Scalar Breit interaction for molecular calculations

Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Au n (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. In conclusion, the scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory.

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

A matrix completion algorithm for efficient calculation of quantum and variational effects in chemical reactions

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study (https://doi.org/10.1063/5.0018326), exploits the low-rank character of the polynomial expansion of potential energy to recover, using a small sample, vibrational frequencies (square roots of nuclear Hessian eigenvalues) constituting the reaction path. Furthermore, these frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage models at the highest levels of assembly. The effort to generate/verify these models against their detailed counterparts is a challenge not yet conclusively solved. Higher level organizations often place a limit on the node count for delivered models, assuming smaller models result in faster computation. The burden of producing and verifying the accuracy of reduced models is primarily on the lower-level organizations, which consumes resources to produce these models. Limiting the number of allowable nodes may prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which require more nodes than centroid based models. A methodology using Thermal Desktop was described in 2010 using: (1) finite elements and edge nodes for the conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. While the methodology was clear, the implementation would have had to be done manually; with the inclusion of the OpenTD Application Programming Interface, this methodology can now be implemented programmatically and be a viable approach for the generation of reduced models. The approach was implemented at NASA-Goddard for the Capture-Containment-and-Return-System (CCRS) payload, resulting in the TCYEE tool. CCRS requires delivery with node limitations through the European Space Agency to its contractors. TCYEE generated the reduced model for delivery and the predictions compared favorably to the detailed model. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements. This paper describes the methodology, implementation, and the performance of reduced models compared to their detailed counterparts.

Computer Programming and Software↗
NASA NTRS2022

Continuum Covariance Propagation for Understanding Variance Loss in Advective Systems

Motivated by the spurious variance loss encountered during covariance propagation in atmospheric and other large-scale data assimilation systems, we consider the problem for state dynamics governed by the continuity and related hyperbolic partial differential equations. This loss of variance has been attributed to reduced-rank representations of the covariance matrix, as in ensemble methods for example, or else to the use of dissipative numerical methods. Through a combination of analytical work and numerical experiments, we demonstrate that significant variance loss, as well as gain, typically occurs during covariance propagation, even at full rank. The cause of this unusual behavior is a discontinuous change in the continuum covariance dynamics as correlation lengths become small, for instance in the vicinity of sharp gradients in the velocity field. This discontinuity in the covariance dynamics arises from hyperbolicity: the diagonal of the kernel of the covariance operator is a characteristic surface for advective dynamics. Our numerical experiments demonstrate that standard numerical methods for evolving the state are not adequate for propagating the covariance, because16they do not capture the discontinuity in the continuum covariance dynamics as correlations lengths tend to zero. Our analytical and numerical results show that this leads to significant, spurious variance loss in certain regions, and gain in others. The results suggest that developing local covariance propagation methods designed specifically to capture covariance evolution near the diagonal may prove a useful alternative to current methods of covariance propagation.

covariance propagation↗
DOE OSTI2026

Three-dimensional high-content imaging of unstained soft tissue with subcellular resolution using a laboratory-based X-ray microscope

With increasing interest in studying biological systems across spatial scales—from centimeters down to nanometers—histology continues to be the gold standard for tissue imaging at cellular resolution, providing an essential bridge between macroscopic and nanoscopic analysis. However, its inherently destructive and two-dimensional nature limits its ability to capture the full three-dimensional complexity of tissue architecture. Here, we show that phase-contrast X-ray microscopy can enable three-dimensional virtual histology with subcellular resolution. This technique provides direct quantification of electron density without restrictive assumptions, allowing for direct characterization of cellular nuclei in a standard laboratory setting. By combining high spatial resolution and soft tissue contrast, with automated segmentation of cell nuclei, we demonstrated virtual Hematoxylin and Eosin (H&E) staining using machine learning-based style transfer, yielding volumetric datasets compatible with existing histopathological analysis tools. Furthermore, by integrating electron density and the sensitivity to nanometric features of the dark field contrast channel, we achieve stain-free, high-content imaging capable of distinguishing nuclei and extracellular matrix.

3D virtual histology↗
NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD API, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return (MSR) mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency (ESA) with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Hume L. Peabody↗
NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD API, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return (MSR) mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency (ESA) with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Thermal desktop↗
NASA NTRSDate not supplied

Implementation of a Hybrid Edge Node-Centroid Node Approach for the Generation of Reduced Thermal Models

Reduced thermal models are often required for delivery to organizations that manage observatory or launch models at the highest levels of assembly. However, the effort to generate reduced models, and verify against their detailed counterparts, is a challenge that has not yet been conclusively solved. Higher level organizations often place a limit on the number of nodes for delivered models with the assumption that smaller models generally result in less computation time. However, the burden of producing and verifying the accuracy of the reduced models is placed primarily on the lower-level organizations, which in turn consumes resources needed to produce these models. Limiting the total number of allowable nodes may also prevent users from taking full advantage of software capabilities that allow for faster generation of models, such as finite elements, which generally require more nodes than older centroid based models. A methodology using Thermal Desktop was described in 2010 which used: (1) finite elements and edge nodes for a conduction matrix, (2) centroid nodes for capacitance and radiative computations, and (3) the super network feature to produce a conduction matrix based only on the centroid nodes. At that time, the methodology was clear, but the implementation would have had to be done manually; however, with the inclusion of the OpenTD Application Programming Interface, this methodology can now be implemented programmatically and for the first time, be a viable approach for the generation of reduced models. The approach was implemented and developed at the NASA Goddard Space Flight Center (GSFC) for the Capture, Containment, and Return System (CCRS) payload on the Earth Return Orbiter (ERO) as part of the Mars Sample Return mission, resulting in the TCYEE tool. ERO features a spacecraft bus provided by Airbus through the European Space Agency with node limitations on the delivered CCRS model provided by GSFC. TCYEE was used to generate the reduced model for this delivery and the predictions compared favorably to the detailed model currently in use for the thermal performance evaluation. Furthermore, TCYEE is being explored for potential use on the Roman Space Telescope (RST) for the generation of reduced models for delivery to the launch provider, which also has node limit requirements on the RST observatory model for use in launch simulation analyses. This paper describes the methodology, its implementation, and compares the performance of reduced models generated by TCYEE to their detailed counterparts.

Thermal Desktop↗
NASA NTRS2019

High Velocity Impact Performance of a Dual Layer Thermal Protection System for the Mars Sample Return Earth Entry Vehicle

The Mars Sample Return (MSR) Earth Entry Vehicle (EEV) is currently planned on being released from its Micro-Meteorite/Orbital Debri (MM/OD) shielding housing about two days before the Earth entry phase. This leaves the EEV exposed to incoming MM/OD impacts, potentially damaging the heat shield and compromising its Entry, Decent, and Landing (EDL) integrity. Currently, two materials are proposed to comprise the MSR-EEV heat shield, a dual layer material Heat-shield for Extreme Entry Environment Technology (HEEET) and Phenolic Impregnated Carbon Ablator (PICA). PICA has been well characterized for OD class impacts, ~7km/s high mass impacts, from previous testing done in the Orion program, but hasn't under-gone extensive MM impact testing. HEEET is a rela-tively new material with minimal prior testing in re-gards to High Velocity Impacts (HVI). In order to in-form selectability of a material, it is crucial to under-stand the material performance when faced with an HVI, directly affecting mission success probability.The current measure of a Thermal Protection Sys-tem's (TPS) performance against an HVI is evaluating a thermally sized material against its derived Balistic Limit Equation (BLE). A BLE is generated empirically from multiple shots of HVI testing, and is used as a first order method in evaluating TPS's performance against the expected MM/OD environment. This method has proved useful for previous uniform densi-ty TPS materials, but has never been validated against a dual layer recession material such as HEEET. Testing in the MSR program for FY19 has a re-quirement to assess the effects of a dual layer TPS by testing various thicknesses of HEEET's Recession Layer, seen in Figure 1. From this data, BLE's will be derived for the individual thickness ratio samples, as well as the material as a whole to evaluate if a heritage form of the BLE can capture the complex physics associated with a dual layer system. Crater morpholo-gy will also be assessed with post-test Non-Destructive Evaluation (NDE) methods such as CT scanning to visualize if a BLE can well predict the associated pene-tration depths, since a BLE assumes full disinigration of the impacting particle and is generally only used to size a spherical crater ? disregarding any shrapenel effects from a high density impactor.To inform this analysis, expected MM environ-ments from the Meteoroid Engineering Model (MEM) and analytical equations for mass flux of incoming MM were evaluated against the notional MSR-EEV trajectory [1]. Using those dispersions, a monte-carlo was run to determine the most probable particle pa-rameters, as well as the riskiest in terms of full bondline penetration. From these probabilities, a test matrix was designed to test against bounding cases for the various parameters of the BLE: projectile density, projectile mass, projectile velocity, and the impact angle.This presentation will discuss the performance of the dual layer TPS material HEEET against a wide range of impact kinetic energies and densitites, as well as the comparison of HEEET to PICA in terms of MM/OD performance and selectability criteria.

Libben, B. J.↗
DOE OSTI2022

Comparison of Fission Product Release Predictions using PARFUME and BISON with Results from the AGR-3/4 Irradiation Experiment

The PARFUME (PARticle Fuel ModEl) fuel performance modeling code and the BISON nuclear fuel performance application built on the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element library were used to predict the fission product release from tristructural isotropic (TRISO) coated fuel particles and compacts during the third and fourth irradiation experiment of the Advanced Gas Reactor (AGR-3/4) Fuel Development and Qualification Program. The fuel performance modeling codes PARFUME and BISON modeled the AGR-3/4 irradiation experiment using the fuel compact time-averaged volume averaged (TAVA) daily temperatures for a total irradiation duration of 369.1 effective full power days (EFPD) to predict the release fraction of the fission product silver (Ag-110m) from a representative TRISO-coated fuel particle from AGR-3/4 compacts. Post-irradiation examination (PIE) measurements provided data on the release of these fission products in the compacts outside of the silicon carbide (SIC) layer. The PARFUME and BISON results were then compared to the silver release measured from compact gamma scanning. The results showed good agreement between PARFUME and BISON but both codes under-predicted the silver release fraction for all the compacts. In addition, BISON was used to model and predict the fission product concentration radial profile outside of the compacts in capsules’ inner and outer rings. These rings were either comprised of matrix and/or structural graphite. To obtain the concentration profiles of silver, cesium, and strontium, a sorption isotherm model was developed in BISON to capture the effects of fission product transport across the gaps between the concentric rings. The general shape of the concentration radial profiles as calculated by BISON were similar in the inner ring (IR) but varied in the outer ring (OR) depending on the fission product of interest or capsule temperature. Using this methodology and model, BISON now has the capability to aid in developing new fission product diffusion coefficients for matrix or structural graphite materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗
NASA NTRS2014

The Tethered Moon

A reasonable initial condition on Earth after the Moonforming impact is that it begins as a hot global magma ocean1,2. We therefore begin our study with the mantle as a liquid ocean with a surface temperature on the order of 3000- 4000 K at a time some 100-1000 years after the impact, by which point we can hope that early transients have settled down. A 2nd initial condition is a substantial atmosphere, 100-1000 bars of H2O and CO2, supplemented by smaller amounts of CO, H2, N2, various sulfur-containing gases, and a suite of geochemical volatiles evaporated from the magma. Third, we start the Moon with its current mass at the relevant Roche limit. The 4th initial condition is the angular momentum of the Earth-Moon system. Canonical models hold this constant, whilst some recent models begin with considerably more angular momentum than is present today. Here we present a ruthlessly simplified model of Earth's cooling magmasphere based on a full-featured atmosphere and including tidal heating by the newborn Moon. Thermal blanketing by H2O-CO2 atmospheres slows cooling of a magma ocean. Geochemical volatiles - chiefly S, Na, and Cl - raise the opacity of the magma ocean's atmosphere and slow cooling still more. We assume a uniform mantle with a single internal (potential) temperature and a global viscosity. The important "freezing point" is the sharp rheological transition between a fluid carrying suspended crystals and a solid matrix through which fluids percolate. Most tidal heating takes place at this "freezing point" in a gel that is both pliable and viscous. Parameterized convection links the cooling rate to the temperature and heat generation inside the Earth. Tidal heating is a major effect. Tidal dissipation in the magma ocean is described by viscosity. The Moon is entwined with Earth by the negative feedback between thermal blanketing and tidal heating that comes from the temperature-dependent viscosity of the magma ocean. Because of this feedback, the rate that the Moon's orbit evolves is limited by the modest radiative cooling rate of Earth's atmosphere, which in effect tethers the Moon to the Earth. Consequently the Moon's orbit evolves orders of magnitude more slowly than in conventional models. Slow orbital evolution promotes capture by orbital resonances that may have been important in the Earth-Moon system

Earth-Moon system↗
DOE OSTI2025

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗
DOE OSTI2025

Accelerated Fuel Qualification of Fast Modular Reactor Fuel in a Thermal Reactor: Modeling and Simulation Paired with Irradiation Testing

The accelerated fuel qualification (AFQ) methodology is applied by simulating accelerated fuel tests of the General Atomics Electromagnetic Systems’ fuel system for its 44-MW(electric) gas-cooled, fast-spectrum fast modular reactor (FMR). This fuel is comprised of UO 2 pellets in SiGA® cladding, a silicon carbide ceramic matrix composite. Fast reactors, like the FMR, offer many benefits, including high fuel utilization and flexibility, but may require a lengthy material design process if tests are performed using fast neutron irradiation alone. A thermal neutron irradiation can instead be used to rapidly test how well key components of the current material models extend to high burnup. Thermal neutrons produce a different radial power distribution within the pin than fast neutrons. However, the temperature and burnup values for the two neutron types are comparable, and the differences between the simulated fuel responses are relatively small, demonstrating the weak sensitivity of the physics-based fuel model calculations on the neutron type and the irradiation rate. Furthermore, the deformation of the SiGA cladding saturates after about 1 displacement per atom for both neutron spectra. In an accelerated fuel test, the irradiation time required to reach the target fuel burnup can be reduced by a factor of 3 by using a small rodlet with a 45% smaller pellet diameter while maintaining the same linear power. Therefore, the time for data collection up to high burnup can be significantly reduced while maintaining the same temperature profile, which largely determines the material response. Tests of fuel rodlets of standard and compact size will be carried out in the Idaho National Laboratory’s Advanced Test Reactor (ATR), including full size and compact rodlets with varying gap sizes. By applying physics-based mechanistic modeling and simulation in accordance with the AFQ methodology, this type of compact rodlet testing in a thermal test reactor captures the necessary phenomena to test fuel material models up to high burnup and to simulate the expected impact of fast neutron radiation on the fuel in FMR operations. Finally, this approach to testing fast reactor fuels in existing thermal test reactors, paired with advanced physics-based mechanistic modeling and simulation, is expected to be applicable to a range of advanced fuels and will decrease the overall fuel qualification timeframe from decades to years.

Advanced test reactor (ATR)↗