Enabling GPU accelerated computing in the SUNDIALS time integration library
Not Available
97 MATHEMATICS AND COMPUTING↗
SEARCH · Engineering Papers
Results for “numerical computation”
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.
Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals
In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Shock-driven kinetic and diffusive mix in high-Z pusher ICF designs
Revolver and Double Shell Inertial Confinement Fusion capsule designs hope to achieve a robust volumetric thermonuclear burn via the use of a high-Z pusher shell filled with a cryogenic D–T fuel. Unfortunately, mix of the pusher material into the fuel (gas) may adversely impact the burn performance. Hydrodynamic instability of the metal/gas interface as the mix source is an obvious concern, but 1D effects may also be detrimental. Such effects include plasma diffusion at material interfaces, which has been the subject of numerous theoretical, computational, and experimental investigations. However, other 1D mix mechanisms may exist, which have yet to be thoroughly explored. In particular, plasma kinetic effects may drive the mix when a shock breaks out of the metal/gas interface. Using the state-of-the-art, hybrid (kinetic-ion/fluid electron), multi-ion Vlasov–Fokker–Planck code, iFP, we show herein that shock-driven kinetic effects can reconfigure the interface and the interfacial width subsequently grows diffusively. Finally, we consider any implications for high-Z pusher designs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
MASTODON: An Open-Source Software for Seismic Analysis and Risk Assessment of Critical Infrastructure
Seismic analysis and risk assessment of safety-critical infrastructure like hospitals, nuclear power plants, dams, and facilities handling radioactive materials involve computationally intensive numerical models and coupled multiphysics scenarios. They are also performed in a strict regulatory environment that requires high software quality assurance standards, and in the case of safety-related nuclear facilities, a conformance to the American Society of Mechanical Engineers Nuclear Quality Assurance (NQA-1) standard. This paper introduces the open-source finite-element software, MASTODON (Multi-hazard Analysis of Stochastic Time-Domain Phenomena), which implements state-of-the-art seismic analysis and risk assessment tools in a quality-controlled environment. MASTODON is built on MOOSE (Multi-physics Object-Oriented Simulation Environment), which is a highly parallelizable, NQA-1 conforming, coupled multiphysics, finite-element framework developed at Idaho National Laboratory. MASTODON is capable of fault rupture and source-to-site wave propagation using the domain reduction method, nonlinear site response, and soil-structure interaction analysis, implicit and explicit time integration, automated stochastic simulations, and seismic probabilistic risk assessment. When coupled with other MOOSE applications, MASTODON can also solve strongly and weakly coupled multiphysics problems. This paper presents a summary of the capabilities of MASTODON and some demonstrative examples.
42 ENGINEERING↗
Parametric Finite Element Analysis of Naturally Corroded Steel Specimens Using 3D Surface Laser Scans
Corrosion is considered a uniform thickness reduction design guideline of the maritime industry. However, additionally, the corroded and irregular morphology of the surface affects the steel's load-bearing capacity and its impact on the strength and elongation behaviour of the steel is not yet fully understood. These effects on the local behaviour of steel structures under tensile loading were investigated with tensile tests on naturally corroded steel specimens and nonlinear finite element simulations including the corroded surface morphology with a uniform surface idealation. The models also include the deformed specimen shape. The developed approach led to highly accurate parametric finite element models predicting the ultimate tensile strength and longitudinal position of fracture. The results show that all included aspects are essential for accurate simulations, while solely the maximum available surface resolution was not as decisive.
corrosion↗
Window convolution of the galaxy clustering bispectrum
In galaxy survey analysis, the observed clustering statistics do not directly match theoretical predictions but rather have been processed by a window function that arises from the survey geometry including the sky footprint, redshift-dependent background number density and systematic weights. While window convolution of the power spectrum is well studied, for the bispectrum with a larger number of degrees of freedom, it poses a significant numerical and computational challenge. In this work, we consider the effect of the survey window in the tripolar spherical harmonic decomposition of the bispectrum and lay down a formal procedure for their convolution via a series expansion of configuration-space three-point correlation functions, which was first proposed by Sugiyama et al. (2019). We then provide a linear algebra formulation of the full window convolution, where an unwindowed bispectrum model vector can be directly premultiplied by a window matrix specific to each survey geometry. To validate the pipeline, we focus on the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) luminous red galaxy (LRG) sample in the South Galactic Cap (SGC) in the redshift bin 0.4 ≤ z ≤ 0.6. We first perform convergence checks on the measurement of the window function from discrete random catalogues, and then investigate the convergence of the window convolution series expansion truncated at a finite of number of terms as well as the performance of the window matrix. This work highlights the differences in window convolution between the power spectrum and bispectrum, and provides a streamlined pipeline for the latter for current surveys such as DESI and the Euclid mission.
79 ASTRONOMY AND ASTROPHYSICS↗
Asynchronous quadratic control for constrained hidden markov jump linear systems with incomplete MTPM and MOCPM
Abstract This paper investigates the quadratic optimal control problem for constrained Markov jump linear systems with incomplete mode transition probability matrix (MTPM). Considering original system mode is not accessible, observed mode is utilized for asynchronous controller design where mode observation conditional probability matrix (MOCPM), which characterizes the emission between original modes and observed modes is assumed to be partially known. An LMI optimization problem is formulated for such constrained hidden Markov jump linear systems with incomplete MTPM and MOCPM. Based on this, a feasible state-feedback controller can be designed with the application of free-connection weighting matrix method. The desired controller, dependent on observed mode, is an asynchronous one which can minimize the upper bound of quadratic cost and satisfy restrictions on system states and control variables. Furthermore, clustering observation where observed modes recast into several clusters, is explored for simplifying the computational complexity. Numerical examples are provided to illustrate the validity.
Zhu, Jin↗
Fast calculation of diffraction patterns from an ensemble of aligned molecules
We report an algorithm to calculate electron diffraction patterns for molecules with anisotropic angular distribution, which is significantly faster than existing methods. The algorithm uses a transform to convert the molecular orientation distribution, which is a function of three Euler angles, to the atom-pair distribution functions which depend on the polar and azimuthal angles. The diffraction signal can then be calculated from the atom-pair distributions. We demonstrate the computation method numerically by calculating electron diffraction patterns for a symmetric top molecule (trifluoroiodomethane) and an asymmetric top molecule (formaldehyde) and show that it reduces the calculation time by approximately two orders of magnitude compared to the standard brute-force method. Here, the method can also be applied to the calculation of x-ray diffraction patterns.
74 ATOMIC AND MOLECULAR PHYSICS↗
High-temperature superconductivity in LaH 10
The recent discovery of a high critical temperature $T_c$ in compressed H 3 S has been followed by the prediction of Liu et al. of $T_c$ ≈ 250 K in the clathrate LaH 10 structure. This report has been confirmed experimentally by Somayazulu et al. and Drozdov et al. Additional theoretical work by Wang et al. and Quan et al. further established the mechanism of electron-phonon interaction and the dominant role of hydrogen. In the present Rapid Communication we follow the classic McMillan paper, which separates the electron and phonon contributions to the electron-phonon coupling λ. We first compute the numerator of McMillan's expression, the Hopfield parameter η, using the theory of Gaspari and Gyorffy (GG), and obtain the force constants in the denominator from Wang et al. and Quan et al. The resulting λ is used in the Allen-Dynes equation to calculate $T_c$. The value of $T_c$ reaches a maximum in the range of 236–263 K at pressures of 255 GPa and decreases for smaller or larger pressures. Further, we provide a thorough analysis of the different terms of the GG equation and draw the conclusion that the $\textit{sp}$ channel of the hydrogen is the most important contribution to obtain high values of $T_c$ in this material. Consistent with Wang et al., we find large values of λ that decrease with increasing pressure. In addition, we find that the hydrogen sites are the largest contributors to the total value of the coupling constant λ. That is, the acoustic mode associated with La contributes only 2% to the total λ, while the optic modes associated with H contribute 18% for the H1 site and 80% for the H2 site. These relative contributions to λ are consistent with those given by Wang et al. and by Quan et al. Thus, our results strongly support the view that LaH 10 is a metallic hydrogen superconductor.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Mean-field entanglement transitions in random tree tensor networks
Entanglement phase transitions in quantum chaotic systems subject to projective measurements and in random tensor networks have emerged as a new class of critical points separating phases with different entanglement scaling. We propose a mean-field theory of such transitions by studying the entanglement properties of random tree tensor networks. As a function of bond dimension, we find a phase transition separating area-law from logarithmic scaling of the entanglement entropy. Using a mapping onto a replica statistical mechanics model defined on a Cayley tree and the cavity method, we analyze the scaling properties of such transitions. Our approach provides a tractable, mean-field-like example of an entanglement transition. Furthermore, we verify our predictions numerically by computing directly the entanglement of random tree tensor network states.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Basis light-front quantization for a chiral nucleon-pion Lagrangian
We present the first application of the basis light-front quantization method to a simple chiral model of the nucleon-pion system as a relativistic bound state for the physical proton. In this work, the light-front mass-squared matrix of the nucleon-pion system is obtained within a truncated basis. The mass and the corresponding light-front wave function (LFWF) of the proton are computed by numerical diagonalization of the resulting mass-squared matrix. With the boost-invariant LFWF, we calculate the probability density distribution of the pion's longitudinal momentum fraction and the Dirac form factor of the proton.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Elastic dijet production in electron scattering on a longitudinally polarized proton at small x : A portal to orbital angular momentum distributions
We calculate the elastic production of dijets from electron collisions with a longitudinally polarized proton target at small values of the Bjorken x variable. Building on the pioneering proposals of Hatta [Gluon orbital angular momentum at small- x , .] and Bhattacharya [Signature of the gluon orbital angular momentum, ., Probing the quark orbital angular momentum at electron-ion colliders using exclusive π 0 production, ., and Exploring orbital angular momentum and spin-orbit correlation for gluons at the Electron-Ion Collider, .] for measuring the quark and gluon orbital angular momentum (OAM) distributions, our focus is on both the longitudinal double spin asymmetry (DSA) and longitudinal single spin asymmetry (SSA). We compute the numerators of these asymmetries in the small- x formalism of the light-cone operator treatment. Utilizing the small- x expressions for the OAM distributions derived earlier in Kovchegov and Manley [Orbital angular momentum at small x revisited, .], we demonstrate that the DSA provides a robust probe for both the quark and gluon OAM distributions within the proton. In contrast, we find that while the SSA is also sensitive to the OAM distributions, extraction of the latter from the SSA would require new developments in small- x theory and phenomenology, and is probably not feasible at this point in time. These findings highlight the potential of DSA measurements in elastic dijet production at the future Electron-Ion Collider to provide the first-ever direct access to the quark and gluon OAM distributions at small x , paving the way for new insights into the proton spin puzzle. Published by the American Physical Society 2025
Kovchegov, Yuri V. (ORCID:0000000169909173)↗
Theory of Trotter Error with Commutator Scaling
The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.
quantum computing, numerical analysis↗
First-principles thermodynamics of Al 10 V: An analytical treatment of localized anharmonic modes
Many complex intermetallic structures possess cagelike environments that can host additional guest atoms. In Al 10 V, these atoms give rise to low-frequency, localized vibrations (Einstein modes) that dominate the thermodynamic response at low temperature. They become imaginary under volume expansion as temperature rises, invalidating the harmonic approximation. Here, we develop a framework to incorporate these strongly anharmonic vibrational modes into first-principles thermodynamic calculations. By explicitly modeling the cage potential and solving the associated Schrödinger equation numerically, we compute the full anharmonic free energy contribution and demonstrate its impact on the thermodynamic behavior of Al 10 V. This allows us to examine structures with different cage fillings and construct the Al-V phase diagram in the relevant composition range. Our results reproduce key experimental signatures, including the anomalous rise in the thermal expansion coefficient and heat capacity at low temperatures, and reveal that the presence and the extent of cage filling by guest atoms is essential to stabilizing the Al 10 V phase at elevated temperatures.
anharmonic lattice dynamics↗
Enhanced Planar Antenna Efficiency Through Magnetic Thin-Films
This work proposes to use magnetic material as the substrate of planar antennas to overcome the platform effect caused by the conducting ground plane. The upper bound of the radiation efficiency of an electric-current-driven low-profile antenna is theoretically derived, which is inversely proportional to the Gilbert damping factor of the magnetic material. Meanwhile, the improvement of radiation due to the use of magnetic material is demonstrated by a three-dimensional (3D) multiphysics and multiscale time-domain model. The simulation results match the theoretical derivation, showing 25% radiation efficiency from a planar antenna backed by a FeGaB thin film with 2.56 m thickness. Furthermore, for conductive ferromagnetic materials, it is shown that the eddy current loss can be well suppressed by laminating the thin film into multiple layers. Here, the radiation efficiency of the modeled antenna with a conductive ferromagnetic substrate is improved from 2.2% to 11.8% by dividing the substrate into 10 layers, with a ferromagnetic material fill factor of 93%.
42 ENGINEERING↗
Multicode Signaling in a Filter Bank Multicarrier Spread Spectrum System and Its Application to HF Communications
The use of filter banks for implementing multicarrier spread spectrum systems leads to a class of effective waveforms that are highly resilient to partial-band interferers. Such waveforms can be also designed to keep the peak-to-average power ratio (PAPR) of the resulting signal at a minimum level. The use of multiple spreading gain vectors (known as multicodes), on the other hand, is an effective method for increasing the data rate in spread spectrum systems, in general. This paper presents a detailed analysis of a class of filter bank multicarrier spread spectrum (FBMC-SS) waveforms and demonstrates an effective receiver implementation of them when multicodes are applied. Application of the developed multicode waveform for communications over high-frequency (HF) skywave channels is also explored, and the benefits that it provides are studied both numerically, through computer simulations, and experimentally, by examining the receivers performance over a variety of skywave links.
42 ENGINEERING↗
METRICS FOR INTERCOMPARISON OF REMAPPING ALGORITHMS (MIRA)
The MIRA project comprises of a set of Python drivers that enable rigorous intercomparison of remapping algorithms by providing the infrastructure for computing the numerical metrics of interest for ESM. This work is funded by the SciDAC Coupling Approaches for Next-Generation Architectures (CANGA) project.
SciDAC↗
A Brief Overview of Reactive Transport Codes Used in CO 2 Applications
Reactive transport modelling consists of computational and numerical models that describe the coupled physical, chemical, mechanical, and biological processes interacting with each other over a broad range of spatial and temporal scales. This report reviews some of the reactive transport codes available in the literature for carbon dioxide (CO 2 ) applications and can assist the scientific community with the applicability of each code to solve a specific problem. Approaches implemented in these codes can be categorized in three groups: (i) continuum-scale, (ii) pore-scale, and (iii) hybrid-scale approaches. Theoretical foundations, numerical implementations, and application examples using the models are described. In the continuum approach, flow and transport are formulated in terms of a representative control volume of the medium and coupled with geochemical reactions. The medium is characterized by bulk parameters such as porosity, permeability, or reactive surface area. In the pore-scale approach, each point of space in the pore network is occupied by either a fluid or solid phase. The pore-scale approach requires an exact knowledge of the spatial and temporal phase distribution. In the hybrid-scale approach, media are described by multiple characteristic length-scales, with some regions using a pore-scale approach while others used a continuum approach.
01 COAL, LIGNITE, AND PEAT↗