Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix approximation”

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 487 records · Page 27

Diffusion-Model-Assisted Supervised Learning of Generative Models for Density Estimation

Here, we present a supervised learning framework of training generative models for density estimation. Generative models, including generative adversarial networks (GANs), normalizing flows, and variational auto-encoders (VAEs), are usually considered as unsupervised learning models, because labeled data are usually unavailable for training. Despite the success of the generative models, there are several issues with the unsupervised training, e.g., requirement of reversible architectures, vanishing gradients, and training instability. To enable supervised learning in generative models, we utilize the score-based diffusion model to generate labeled data. Unlike existing diffusion models that train neural networks to learn the score function, we develop a training-free score estimation method. This approach uses mini-batch-based Monte Carlo estimators to directly approximate the score function at any spatial-temporal location in solving an ordinary differential equation (ODE), corresponding to the reverse-time stochastic differential equation (SDE). This approach can offer both high accuracy and substantial time savings in neural network training. Once the labeled data are generated, we can train a simple, fully connected neural network to learn the generative model in the supervised manner. Compared with existing normalizing flow models, our method does not require the use of reversible neural networks and avoids the computation of the Jacobian matrix. Compared with existing diffusion models, our method does not need to solve the reverse-time SDE to generate new samples. As a result, the sampling efficiency is significantly improved. We demonstrate the performance of our method by applying it to a set of 2D datasets as well as real data from the University of California Irvine (UCI) repository.

97 MATHEMATICS AND COMPUTING↗

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗

QCD factorization with multihadron fragmentation functions

Important aspects of quantum chromodynamics (QCD) factorization theorems are the properties of the objects involved that can be identified as universal. One example is that the definitions of parton densities and fragmentation functions for different types of hadrons differ only in the identity of the nonperturbative states that form the matrix elements, but are otherwise the same. This leads to independence of perturbative calculations on nonperturbative details of external states. It also lends support to interpretations of correlation functions as encapsulations of intrinsic nonperturbative properties. These characteristics have usually been presumed to still hold true in fragmentation functions even when the observed nonperturbative state is a small-mass cluster of n hadrons rather than simply a single isolated hadron. However, the multidifferential aspect of cross sections that rely on these latter types of fragmentation functions complicates the treatment of kinematical approximations in factorization derivations. That has led to recent claims that the operator definitions for fragmentation functions need to be modified from the single hadron case with nonuniversal prefactors. With such concerns as our motivation, we retrace the steps for factorizing the unpolarized semi-inclusive e + e − annihilation cross section and confirm that they do apply without modification to the case of a small-mass multihadron observed in the final state. In particular, we verify that the standard operator definition from single hadron fragmentation, with its usual prefactor, remains equally valid for the small-mass n -hadron case with the same hard parts and evolution kernels, whereas the more recently proposed definitions with nonuniversal prefactors do not. Our results reaffirm the reliability of most past phenomenological applications of dihadron fragmentation functions. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

When Do Extended Physics-Informed Neural Networks (XPINNs) Improve Generalization?

Physics-informed neural networks (PINNs) have become a popular choice for solving high-dimensional partial differential equations (PDEs) due to their excellent approximation power and generalization ability. Recently, extended PINNs (XPINNs) based on domain decomposition methods have attracted considerable attention due to their effectiveness in modeling multiscale and multiphysics problems and their parallelization. However, theoretical understanding of their convergence and generalization properties remains unexplored. In this study, we take an initial step towards understanding how and when XPINNs outperform PINNs. Specifically, for general multilayer PINNs and XPINNs, we first provide a prior generalization bound via the complexity of the target functions in the PDE problem and a posterior generalization bound via the posterior matrix norms of the networks after optimization. Moreover, based on our bounds, we analyze the conditions under which XPINNs improve generalization. Concretely, our theory shows that the key building block of XPINN, namely, the domain decomposition, introduces a tradeoff for generalization. On the one hand, XPINNs decompose the complex PDE solution into several simple parts, which decreases the complexity needed to learn each part and boosts generalization. On the other hand, decomposition leads to less training data being available in each subdomain, and hence such a model is typically prone to overfitting and may become less generalizable. Empirically, we choose five PDEs to show when XPINNs perform better than, similar to, or worse than PINNs, hence demonstrating and justifying our new theory.

97 MATHEMATICS AND COMPUTING↗

Validation of neutron-induced reactions on natural carbon using an active target at neutron energies up to 22 MeV at LANSCE

A single crystal chemical vapor deposited (sCVD) diamond detector is used as an active target to measure neutron-induced reactions on natural carbon using the neutrons produced by spallation, with a broad energy spectrum at LANSCE. Additionally, the neutron-induced reactions are detected in the diamond as low as E n = 400 keV and up to approximately 100 MeV. Relative cross sections for C 12 ( n , α 0 ) , C 12 ( n , p 0 ) , C 12 ( n , d 0 + p 1 ) , and C 13 ( n , α 0 ) are reported up to E n = 22 MeV and comparisons on detected pulse-height spectra and detector response of scattering reactions are made with GEANT4 simulations using the ENDF/B-VIII.0 evaluated nuclear data library up to 20 MeV. The results are compared with past experimental data, including other works that incorporate diamond detectors as an active carbon target. In addition, R-matrix calculations for the C 13 + n system are presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Lattice quantum chromodynamics at large isospin density

We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of n-pion correlation functions for n ϵ {1, … , N} using a single N × N matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 charged pions using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass m π = 170 MeV and spacetime volumes of (4.4 3 × 8.8) fm 4 and (5.8 3 × 11.6) fm 4 . We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well-approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential μ I ~ 1.5m π , signaling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound, c s , in the medium is seen to exceed the ideal-gas (conformal) limit ($c^{2}_{s} ≤ 1/3)$ over a wide range of chemical potential before falling towards the asymptotic expectation at μ I ~ 15m π . These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The response of accident tolerant fuel cladding to LOCA burst testing: A comparative study of leading concepts

Accident tolerant fuel claddings seek to improve safety margins during loss-of-coolant accident (LOCA) scenarios by reducing cladding steam oxidation rates, hydrogen production, delaying or mitigating cladding burst, and reducing cladding deformation. After rapid and extensive progress in recent years, Cr coated Zr-based alloys, ferritic FeCrAl, and SiC/SiC ceramic matrix composites (CMCs) have emerged as the most promising cladding concepts, but there are few studies comparing the accident tolerance of all three concepts simultaneously. Here, in this study, LOCA burst testing has been conducted on Cr high impulse power magnetron sputtered Zry-4 (Cr/Zry-4), C26M (FeCrAl), and SiC/SiC SiGA (CMC). Experimental observations indicate Cr/Zry-4 burst at higher temperatures than Zry-4 however the effectiveness of the coating on mitigating ZrO 2 formation via Cr 2 O 3 formation depended on physical proximity to the burst opening. No dramatic improvement in diametrical strain nor opening geometry was observed. C26M burst at higher temperatures than Zry-4 at lower pressures, exhibiting small amounts of strain and smaller openings, and formed slow-growing Al 2 O 3 regardless of location. The CMCs did not burst or show any macroscopic signs of deformation, all the way up to an internal cladding pressure of 15.2 MPa. High temperature expanding Nb plug testing indicated that CMC room temperature strength may persist up to 1650°C, with an approximately 40% loss of strength at 1900°C. CMC steam oxidation was also compared to reference chemically vapor deposited SiC up to 1700°C, and CMCs were found to oxidize at slightly accelerated rates that were still orders of magnitude lower than bare Zry-4. Projected H 2 generation during LOCA burst testing showed that when accounting for cladding burst exposing inner diameters to steam, Cr coated Zry-4 reduced H 2 generation by a factor of two compared with bare Zry-4, while both C26M and SiC/SiC reduced H 2 generation by orders of magnitude relative to bare Zry-4.

36 MATERIALS SCIENCE↗

Binary operations on neuromorphic hardware with application to linear algebraic operations and stochastic equations

Abstract Non-von Neumann computational hardware, based on neuron-inspired, non-linear elements connected via linear, weighted synapses—so-called neuromorphic systems—is a viable computational substrate. Since neuromorphic systems have been shown to use less power than CPUs for many applications, they are of potential use in autonomous systems such as robots, drones, and satellites, for which power resources are at a premium. The power used by neuromorphic systems is approximately proportional to the number of spiking events produced by neurons on-chip. However, typical information encoding on these chips is in the form of firing rates that unarily encode information. That is, the number of spikes generated by a neuron is meant to be proportional to an encoded value used in a computation or algorithm. Unary encoding is less efficient (produces more spikes) than binary encoding. For this reason, here we present neuromorphic computational mechanisms for implementing binary two’s complement operations. We use the mechanisms to construct a neuromorphic, binary matrix multiplication algorithm that may be used as a primitive for linear differential equation integration, deep networks, and other standard calculations. We also construct a random walk circuit and apply it in Brownian motion simulations. We study how both algorithms scale in circuit size and iteration time.

97 MATHEMATICS AND COMPUTING↗

Phenomenological R-Matrix parameterization of direct, doorway, and compound nuclear reactions [Slides]

The presented R-matrix formalism could improve evaluations of resolved resonance nuclear cross section data by parameterizing direct reaction processes via an orthonormal channel mixing matrix and doorway state processes. Together, they can parameterize direct-semidirect capture. A formal similarity among the R-, K-, T- matrix formalisms was leveraged to parameterize direct and doorway processes. The advantage of this R-matrix parameterization is that it is easily extensible due to its application of projection operators directly upon the (H-E) on p. 6.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spin-spin entanglement in diffractive heavy-quark production

We calculate the spin density matrix of a heavy-quark-antiquark pair ($b\bar{⁢b}$, $c\bar{⁢c}$ or $s\bar{⁢s}$) diffractively produced in deep inelastic scattering and ultraperipheral collisions. We show that the Pomeron exchange leaves characteristic imprints on the entanglement pattern between the quark and the antiquark. For the longitudinally polarized virtual photon, the pair always exhibits maximal entanglement and maximal violation of the Bell-Clauser-Horne-Shimony-Holt inequality. For the transversely polarized photon, the pair is always entangled and Bell violating, reaching maximal entanglement and maximal violation simultaneously when the transverse momentum approximately equals the quark mass.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Optimization-Based Azeotropic Distillation System Synthesis Using Geometric Insights

The synthesis of azeotropic distillation systems is challenging due to the existence of compartments in the residue curve map coupled with the combinatorial aspects from numerous possible system configurations. In this work, an optimization-based approach is introduced to synthesize homogeneous azeotropic distillation systems. The approach employs a network-based representation generated via a matrix method. To design the distillation columns, the modified Underwood equations are adopted, in which pseudocomponent-based compositions are used. For cases where separatrices are significantly curved, corrections via piece-wise linear functions and collinearity properties are implemented. Two examples are presented to illustrate the proposed approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Trojan horse method as an indirect approach to study resonant reactions in nuclear astrophysics

The primary goal of the Trojan horse method (THM) is to analyze resonant rearrangement reactions when the density of the resonance levels is low and statistical models cannot be applied. The main difficulty of the analysis is related with the facts that in the final state the THM reaction involves three particles and that the intermediate particle, which is transferred from the Trojan horse particle to the target nucleus to form a resonance state, is virtual. Another difficulty is associated with the Coulomb interaction between the particles, especially, taking into account that the goal of the THM is to study resonant rearrangement reactions at very low energies important for nuclear astrophysics. The exact theory of such reactions with three charged particles is very complicated and is not available. This is why different approximations are used to analyze THM reactions. In this review paper we describe a new approach based on a few-body formalism that provides a solid basis for deriving the THM reaction amplitude taking into account rescattering of the particles in the initial, intermediate and final states of the THM reaction. Since the THM uses a two-step reaction in which the first step is the transfer reaction populating a resonance state, we address the theory of the transfer reactions. The theory is based on the surface-integral approach and R-matrix formalism. We also discuss application of the THM to resonant reactions populating both resonances located on the second energy sheet and subthreshold resonances, which are subthreshold bound states located at negative energies close to thresholds. We consider the application of the THM to determine the astrophysical factors of resonant radiative-capture reactions at energies so low that direct measurements can hardly be performed due to the negligibly small penetrability factor in the entry channel of the reaction. Here, we elucidated the main ideas of the THM and outline necessary conditions to perform the THM experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A high‐order discontinuous Galerkin approach for physics‐based thermospheric modeling

Abstract The accurate prediction of aerodynamic drag on satellites orbiting in the upper atmosphere is critical to the operational success of modern space technologies, such as satellite‐based communication or navigation systems, which have become increasingly popular in the last few years due to the deployment of constellations of satellites in low‐Earth orbit. As a result, physics‐based models of the ionosphere and thermosphere have emerged as a necessary tool for the prediction of atmospheric outputs under highly variable space weather conditions. This paper proposes a high‐fidelity approach for physics‐based space weather modeling based on the solution of the Navier–Stokes equations using a high‐order discontinuous Galerkin method, combined with a matrix‐free strategy suitable for high‐performance computing on GPU architectures. The approach consists of a thermospheric model that describes a chemically frozen neutral atmosphere in nonhydrostatic equilibrium driven by the external excitation of the Sun. A novel set of variables is considered to treat the low densities present in the upper atmosphere and to accommodate the wide range of scales present in the problem. At the same time, and unlike most existing approaches, radial and angular directions are treated in a nonsegregated approach. The study presents a set of numerical examples that demonstrate the accuracy of the approximation and validate the current approach against observational data along a satellite orbit, including estimates of established empirical and physics‐based models of the ionosphere‐thermosphere system. Finally, a one‐dimensional radial derivation of the physics‐based model is presented and utilized for conducting a parametric study of the main thermal quantities under various solar conditions.

Engineering↗

Topological structure enhanced nanostructure of high temperature polymer exhibiting more than ten times enhancement of dipolar response

Nanofillers in polymer composites will induce interface regions in the polymer with non-uniform nanostructures which can be approximated as multilayer shells surrounding nanofillers. This paper studies topological nanostructures in polymers that interface with one-dimensional (1-D) and zero-dimensional (0-D) nanofillers. This new class of dielectric polymer nanocomposites involves nanofillers at ultra-low volume loading (< 1 vol percent) that generate large dielectric enhancement. The results show that the dipolar response of polyetherimide, a high temperature dielectric polymer, is increased by more than ten times with ultra-low (0.75 vol %) loading of 1-D nanofillers. The results reveal that 1-D nanofillers induce cylindrical shell-topology in the polymer matrix which is more effective in enhancing the high dielectric constant interfacial region than the spherical shell-topology generated by 0-D fillers. Furthermore, the cylindrical shells generated by 1-D fillers provide much higher dielectric enhancement over a broader interfacial region in the polymer matrix, as compared to the spherical shells induced by 0-D nanofillers.

42 ENGINEERING↗

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗

Parameters of fast and high-yield InAs/GaAs quantum dot semiconductor scintillator

InAs quantum dots (QDs) embedded into a waveguiding GaAs semiconductor matrix may produce scintillation detectors with exceptional speed and yield, making them valuable for nuclear security, medical imaging, and high energy physics applications. In this work, we developed thick (~25 μm) epitaxial heterostructres with high luminescence efficiency composed of self-assembled nano-engineered InAs QDs grown by molecular beam epitaxy. Here, the bulk GaAs acts as a stopping material for incident particles and as a waveguide when layer-transferred onto a low-index substrate. Waveguiding and self-absorption (< 1 cm –1 ) were studied using photoluminescence with scanning laser excitation and modeled with ray optics approximation and geometrical coupling of high-index waveguide to a collection fiber. Scintillating signals from α-particles were analyzed with an external photodiode (PD) and an integrated PD which provided an improved optical coupling. The mean charge collected by the integrated PD corresponded to 3 × 10 4 photoelectrons per 1 MeV of deposited energy, or ~13% of the theoretically achievable light yield. Combined with the previously measured QD scintillation time of 0.3-0.6 ns, this makes the InAs/GaAs QD heterostructures the fastest high yield scintillation material reported.

36 MATERIALS SCIENCE↗