Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 541 records · Page 30

Mechanistic verification of empirical UO 2 fuel fracture models

Standard UO 2 fuel pellets used in light-water reactors fracture during irradiation due to the large thermal gradient in the radial direction. Over the decades, numerous researchers have explored fuel cracking from experimental and modeling points of view. To date, there have been both empirical and mechanistic approaches to predict the number of fragments that form in UO 2 . The empirical models only consider maximum power and burnup as inputs. Existing mechanistic approaches for normal operation have not accounted for irradiation effects. Here, this work employs a mechanistic fuel cracking model using the extended finite element method to explore radial crack formation while including a sensitivity analysis that accounts for the randomization of tensile strength within the fuel, the strength randomization criteria (uniform or volume-weighted Weibull), power ramping rates, computational mesh density, maximum power level, and irradiation (burnup) effects. The results indicate that the uncertainty in this mechanistic modeling approach envelopes the predicted values from three different empirical correlations in almost all cases. This means that, for computationally intensive analyses involving UO 2 fragmentation, the empirical correlations can be used. However, since the mechanistic calculations bound those of the empirical correlations, there is confidence in the applicability of the mechanistic approach developed in this work to generate a correlation for fuel types where limited data exists (e.g., doped-UO 2 , U 3 Si 2 ).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Alloying effects on deformation induced microstructure evolution in copper

In this work, we investigated the effects of alloying elements on plastic deformation and microstructure evolution in polycrystalline copper (Cu) and Cu alloyed with 1 wt. % lead (Cu-1 %Pb). These materials were selected due to the size mismatch between Cu and Pb, with the latter forming precipitates at grain boundaries. Multi-modal characterization techniques, including neutron diffraction, electron backscatter diffraction (EBSD), and transmission electron microscopy (TEM), along with finite element simulations were employed to study the deformation behavior across multiple length scales. While both Cu and Cu-1%Pb exhibited similar macroscale response and final deformation textures, both dislocation line profile analysis and TEM revealed increased dislocation density in deformed Cu-1%Pb specimens. The presence of lead precipitates also significantly affected local plastic deformation during compression, with their influence diminishing with increasing strain. These results demonstrate the complex relationships between alloying elements, plastic deformation, microstructural evolution, and material behavior under load. The insights gained from this multi-scale and multi-technique approach contribute to the fundamental understanding of microstructural evolution in immiscible alloys and are valuable for tailoring the properties of structural materials for specific engineering applications.

36 MATERIALS SCIENCE↗

Can local white dwarfs constrain dark matter interactions?

We investigate whether nearby white dwarfs (WDs) can constrain dark matter (DM) interactions with ordinary matter. As experimental sensitivity improves, driven by the Gaia mission, the sample volume of nearby WDs has been increasing over recent years. We carefully select a sample of ten cold, isolated, non-magnetic WDs within 13 pc of the Sun. We model their carbon-oxygen interior using a finite temperature relativistic equation of state and model atmospheres to infer their core temperatures. This enables us to perform a thorough estimation of the DM capture rate and evaporation mass using actual astrophysical observations. Given the low local DM density, we focus on DM that annihilates into long-lived mediators, which escape the WD and later decay into photons. While Fermi-LAT data shows no significant gamma-ray excess, future telescopes, CTA North & South, LHAASO, SWGO, could probe DM-nucleon cross sections down to ∼ 10 -41 cm 2 for DM masses above the TeV scale. Our results are competitive with current direct detection bounds (e.g., LZ) in the multi-TeV regime. This work underscores the importance of systematic WD studies in the broader landscape of DM detection and demonstrates the synergy between astrophysical and terrestrial searches in exploring DM interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling hydraulic fracture opening and closure with proppant transport and settlement

Hydraulic fracturing is a widely used reservoir stimulation technique for improving fluid circulation in rock formations with extremely low permeability, particularly in enhanced geothermal systems (EGS). To better understand the complex processes involved and improve hydraulic stimulation performance, we have developed ELK (ELectrical fracKing), a MOOSE-based 3D finite element application designed to model the behavior of proppant-fluid mixtures in propagating fractures. ELK integrates both the fluid and proppant components, incorporating particle-driven processes such as gravity settling, particle-particle interactions, and strong density and viscosity contracts, in addition to conventional fluid-driven fracture propagation. In this contribution, we extend ELK to model propped fracture closure, which occurs after the injection phase due to a dramatic drop in the effective stress on the fracture plane. During the shut-in, flowback, and production periods, the fracture width decreases, with the closure behavior depending on proppant concentration.. At low concentrations, closure follows a nonlinear joint law linked to the stiffness of asperities in the fracture walls. While at high concentrations, it is controlled by the properties of packed proppant bed. The extended ELK application is validated against several benchmark examples, including the propagation of an inclined frictional crack, fracture opening and sliding in response to fluid injection, and flowback analysis. We believe that ELK’s enhanced capabilities can serve as a valuable tool for the design and optimization of EGS deployment.

15 - GEOTHERMAL ENERGY↗

Adaptive Uncertainty Quantification for Stochastic Hyperbolic Conservation Laws

Here, we propose a predictor-corrector adaptive method for the study of hyperbolic partial differential equations (PDEs) under uncertainty. Constructed around the framework of stochastic finite volume (SFV) methods, our approach circumvents sampling schemes or simulation ensembles while also preserving fundamental properties, in particular hyperbolicity of the resulting systems and conservation of the discrete solutions. Furthermore, we augment the existing SFV theory with a priori convergence results for statistical quantities, in particular push-forward densities, which we demonstrate through numerical experiments. By linking refinement indicators to regions of the physical and stochastic spaces, we drive anisotropic refinements of the discretizations, introducing new degrees of freedom where deemed profitable. To illustrate our proposed method, we consider a series of numerical examples for nonlinear hyperbolic PDEs based on Burgers’ and Euler’s equations.

97 MATHEMATICS AND COMPUTING↗

Chirality reversal at finite magnetic impurity strength and local signatures of a topological phase transition

Here, we study the honeycomb lattice with a single magnetic impurity modeled by adding imaginary next-nearest-neighbor hopping 𝑖⁢ℎ on a single hexagon. This Haldane defect gives a topological mass term to the gapless Dirac cones and generates chirality. For a small density of defects, Neehus et al. [Phys. Rev. Lett. 135, 126604 (2025)] found that the system's chirality reverses at a critical ℎ 𝑐 ≈ 0.95 associated with an unexpected tricritical point of Dirac fermions at zero defect density. We investigate this zero-density limit by analyzing a single defect and computing two experimentally relevant measures of chirality: (1) orbital magnetization via local Chern marker, a bulk probe of all occupied states; and (2) electronic currents of low-energy states. Both probes show a chirality reversal at a critical ℎ 𝑐 ≈ 0.9–1.0. Motivated by this consistency, we propose a defect-scale toy model whose low-energy states reverse their chirality at ℎ$^{'}_{c}$ ≈ 0.87. Remarkably, the same pair of zero-energy bound states also generates the critical point ℎ 𝑐 in the full impurity projected T-matrix. Our results show how the chirality reversal produced by an impurity can be observed either in local probes or in the global topology, and suggest a possible role of the microscopic defect structure at the critical point.

Chern insulators↗

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.

Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗

Study of a High-speed Outer Rotor Burst Containing Enclosure

The increasing demand for high-speed electric machines in many applications pushes the development of high-power density electric motors. Outer rotor motors, which can be designed with a larger airgap diameter than inner rotor motors for the same overall diameter and, therefore, provide higher torque, distinguish themselves in space-constrained but high torque applications. However, increasing the motor speed and diameter results in significant centrifugal loads for the outer rotor motor due to its higher airgap diameter. Additional attention must be paid at the design stage to consider the extreme situation of possible mechanical failure, such as the sudden burst of the high-speed rotating parts in the motor, and measures must be taken to prevent their damage to other system components in advance. Through the finite element dynamic impact analysis, we studied the damage caused by broken parts on the inner wall of the electric motor housing when the motor rotates at 20,000 RPM. The proposed method is expected to provide reasonable recommendations of the enclosure material and its thickness required to protect other components in the same powertrain system when the high-power density motor fails unexpectedly.

Lin, Lianshan↗

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coherent spin states and emergent de Sitter quasinormal modes

As a toy model for the microscopic description of matter in de Sitter space, we consider a Hamiltonian acting on the spin-j representation of SU(2). This is a model with a finite-dimensional Hilbert space, from which quasinormal modes emerge in the large-spin limit. The path integral over coherent spin states can be evaluated at the semiclassical level and from it we find the single-particle de Sitter density of states, including 1/j corrections. Along the way, we discuss the use of quasinormal modes in quantum mechanics, starting from the paradigmatic upside-down harmonic oscillator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Flat band separation and robust spin Berry curvature in bilayer kagome metals

Kagome materials have emerged as a setting for emergent electronic phenomena that encompass different aspects of symmetry and topology. It is debated whether the XV 6 Sn 6 kagome family (where X is a rare-earth element), a recently discovered family of bilayer kagome metals, hosts a topologically non-trivial ground state resulting from the opening of spin–orbit coupling gaps. These states would carry a finite spin Berry curvature, and topological surface states. Here we investigate the spin and electronic structure of the XV 6 Sn 6 kagome family. We obtain evidence for a finite spin Berry curvature contribution at the centre of the Brillouin zone, where the nearly flat band detaches from the dispersing Dirac band because of spin–orbit coupling. In addition, the spin Berry curvature is further investigated in the charge density wave regime of ScV 6 Sn 6 and it is found to be robust against the onset of the temperature-driven ordered phase. Utilizing the sensitivity of angle-resolved photoemission spectroscopy to the spin and orbital angular momentum, our work unveils the spin Berry curvature of topological kagome metals and helps to define its spectroscopic fingerprint.

36 MATERIALS SCIENCE↗

Building materials genome from ground‐state configuration to engineering advance

Individual phases are commonly considered as the building blocks of materials. However, the accurate theoretical prediction of properties of individual phases remains elusive. The top-down approach by decoding genomic building blocks of individual phases from experimental observations is nonunique. The density functional theory (DFT), as a state-of-the-art solution of quantum mechanics, prescribes the existence of a ground-state configuration at 0 K for a given system. It is self-evident that the ground-state configuration alone is insufficient to describe a phase at finite temperatures as symmetry-breaking non-ground-state configurations are excited statistically at temperatures above 0 K. Our multiscale entropy approach (recently terms as Zentropy theory) postulates that the entropy of a phase is composed of the sum of the entropy of each configuration weighted by its probability plus the configurational entropy among all configurations. Consequently, the partition function of each configuration in statistical mechanics needs to be evaluated by its free energy rather than total energy. The combination of the ground-state and symmetry-breaking non-ground-state configurations represents the building blocks of materials and can be used to quantitatively predict free energy of individual phases with the free energy of each configuration predicted from DFT as well as all properties derived from free energy of individual phases.

CALPHAD↗

Soft fragmentation on the celestial sphere

We develop two approaches to the problem of soft fragmentation of hadrons in a gauge theory for high energy processes. The first approach directly adapts the standard resummation of the parton distribution function’s anomalous dimension (that of twist-two local operators) in the forward scattering regime, using k T -factorization and BFKL theory, to the case of the fragmentation function by exploiting the mapping between the dynamics of eikonal lines on transverse-plane to the celestial-sphere. Critically, to correctly resum the anomalous dimension of the fragmentation function under this mapping, one must pay careful attention to the role of regularization, despite the manifest collinear or infra- red finiteness of the BFKL equation. The anomalous dependence on energy in the celestial case, arising due to the mismatch of dimensionality between positions and angles, drives the differences between the space-like and time-like anomalous dimension of parton densities, even in a conformal theory. The second approach adapts an angular-ordered evolution equation, but working in 4 – 2ϵ dimensions at all angles. The two approaches are united by demanding that the anomalous dimension in 4 – 2ϵ dimensions for the parton distribution function determines the kernel for the angular-ordered evolution to all orders.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Efficient and high-capacity storage materials are indispensable for a hydrogen-based economy. In silico tools can accelerate the process of discovery of new adsorbent materials with optimal hydrogen adsorption enthalpies. Density functional theory is well-poised to become a very useful tool for enabling high-throughput screening of potential materials. In this work, we have identified density functional approximations that provide good performance for hydrogen binding applications following a two-pronged approach. First, we have compiled a data set (H2Bind275) that comprehensively represents the hydrogen binding problem capturing the chemical and mechanistic diversity in the binding sites encountered in hydrogen storage materials. We have also computed reference interaction energies for this data set using coupled-cluster theory. Second, we have assessed the performance of 55 density functional approximations for predicting H 2 interaction energies and have identified two hybrid density functionals (ωB97X-V and ωB97M-V), two double hybrid density functionals (DSD-PBEPBE-D3(BJ) and PBE0-DH), and one semilocal density functional (B97M-V) as the best performing ones. We have recommended the addition of empirical dispersion corrections to systematically underbinding density functionals such as revPBE, BLYP, and B3LYP for improvements in performance at negligible additional cost. We have also recommended the usage of the def2-TZVPP basis set as it represents a good compromise between accuracy and cost, limiting the finite basis set errors to less than 1 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fast and flexible neutrino decoupling. Part I. The Standard Model

Cosmological determinations of the number of relativistic neutrino species, N eff , are becoming increasingly accurate, and further improvements are expected both from CMB and BBN data. Given this context, we update the evaluation of N eff and the current entropy density via the momentum-averaged approach. This allows for a numerically fast description of neutrino decoupling, easily portable to an array of new physics scenarios. We revisit all aspects of this approach, including collision terms with full electron mass dependence, finite temperature QED corrections to the equation of state, neutrino oscillations, and the modelling of neutrino ensembles with effective chemical potentials. For integrated observables, our results differ by less than 0.04% from the solution of the momentum-dependent evolution equation. We outline how to extend the approach to BSM settings, and will highlight its power in Part II. To facilitate the practical implementation, we release a Mathematica and python code within nudec_BSM_v2, easily linkable to BBN codes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗