SEARCH · Engineering Papers

Results for “Physics Regularization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Multi-head physics-informed neural networks for learning functional priors and uncertainty quantification

In numerous applications, the integration of prior knowledge and historical information is essential, particularly for tasks requiring the solution of ordinary or partial differential equations (ODEs/PDEs) in data-sparse or noisy environments. For instance, achieving accurate solutions to time-dependent PDEs with limited initial condition measurements necessitates an effective strategy for embedding prior knowledge. Hard-parameter sharing architectures in neural networks (NNs) have demonstrated success in both traditional and scientific machine learning domains, facilitating the learning of informative representations. Here, in this study, we introduce a novel, yet efficient, method to enhance physics-informed neural networks (PINNs) by incorporating a multi-head structure that enables the learning of functional priors from both empirical data and governing physical laws. This prior information can then be used to address data sparsity and high-level noise in solving ODE/PDE problems with uncertainty quantification (UQ). The approach, termed Multi-Head PINN (MH-PINN), consists of a shared body NN and multiple head NNs, each corresponding to an individual PINN instance. Our framework for functional prior learning is carried out in two stages: (1) training the MH-PINNs to develop a shared body NN alongside multiple head NNs, and (2) employing these trained head NNs to estimate a prior distribution through a normalizing flow-based density estimator. The learned functional prior can then be applied as a regularization mechanism in deterministic contexts or as an informative prior within a Bayesian inference framework, aiding in the resolution of subsequent ODE/PDE tasks. We evaluate the efficacy of MH-PINNs across five benchmark problems, including a high-dimensional parametric PDE, all characterized by data sparsity or substantial noise levels. Our findings reveal that MH-PINNs deliver accurate solutions and robust UQ, demonstrating adaptability across a range of complex and challenging scenarios.

Bayesian inference↗

DOE OSTI2024

On the connection between least squares, regularization, and classical shadows

Classical shadows (CS) offer a resource-efficient means to estimate quantum observables, circumventing the need for exhaustive state tomography. Here, we clarify and explore the connection between CS techniques and least squares (LS) and regularized least squares (RLS) methods commonly used in machine learning and data analysis. By formal identification of LS and RLS ``shadows'' completely analogous to those in CS---namely, point estimators calculated from the empirical frequencies of single measurements---we show that both RLS and CS can be viewed as regularizers for the underdetermined regime, replacing the pseudoinverse with invertible alternatives. Through numerical simulations, we evaluate RLS and CS from three distinct angles: the tradeoff in bias and variance, mismatch between the expected and actual measurement distributions, and the interplay between the number of measurements and number of shots per measurement. Compared to CS, RLS attains lower variance at the expense of bias, is robust to distribution mismatch, and is more sensitive to the number of shots for a fixed number of state copies---differences that can be understood from the distinct approaches taken to regularization. Conceptually, our integration of LS, RLS, and CS under a unifying ``shadow'' umbrella aids in advancing the overall picture of CS techniques, while practically our results highlight the tradeoffs intrinsic to these measurement approaches, illuminating the circumstances under which either RLS or CS would be preferred, such as unverified randomness for the former or unbiased estimation for the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Resolving extreme jet substructure

We study the effectiveness of theoretically-motivated high-level jet observables in the extreme context of jets with a large number of hard sub-jets (up to N = 8). Previous studies indicate that high-level observables are powerful, interpretable tools to probe jet substructure for N ≤ 3 hard sub-jets, but that deep neural networks trained on low-level jet constituents match or slightly exceed their performance. We extend this work for up to N = 8 hard sub-jets, using deep particle-flow networks (PFNs) and Transformer based networks to estimate a loose upper bound on the classification performance. A fully-connected neural network operating on a standard set of high-level jet observables, 135 N-subjetiness observables and jet mass, reach classification accuracy of 86.90%, but fall short of the PFN and Transformer models, which reach classification accuracies of 89.19% and 91.27% respectively, suggesting that the constituent networks utilize information not captured by the set of high-level observables. We then identify additional high-level observables which are able to narrow this gap, and utilize LASSO regularization for feature selection to identify and rank the most relevant observables and provide further insights into the learning strategies used by the constituent-based neural networks. The final model contains only 31 high-level observables and is able to match the performance of the PFN and approximate the performance of the Transformer model to within 2%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Nonergodic dynamics in disorder-free potentials

Highlights: • Dynamics in tilted lattices with harmonic confinement is reviewed. • Both spinless and spinful fermions are analyzed. • Local effective tilt is the relevant parameter for nonergodicity. • Neel or domain wall initial states are studied. • Domains with similar spin orientation act as transport and communication barriers. We review the dynamics of interacting particles in disorder-free potentials concentrating on a combination of a harmonic binding with a constant tilt. We show that a simple picture of an effective local tilt describes a variety of cases. Our examples include spinless fermions (as modeled by Heisenberg spin chain in a magnetic field), spinful fermions as well as bosons that enjoy a larger local on-site Hilbert space. We also discuss the domain-wall dynamics that reveals nonergodic features even for a relatively weak tilt as suggested by Doggen et al. (2020). By adding a harmonic potential on top of the static field we confirm that the surprising regular dynamics is not entirely due to Hilbert space shuttering. It seems better explained by the inhibited transport within the domains of identically oriented spins. Once the spin-1/2 restrictions are lifted as, e.g., for bosons, the dynamics involve stronger entanglement generation. Again for domain wall melting, the effect of the harmonic potential is shown to lead mainly to an effective local tilt.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Calculation of the Thermal Neutron Scattering Cross-Section of Solids Using OCLIMAX

The thermal neutron scattering cross-section of a solid depends on the energy (or wavelength) of the incident neutrons. Devising a method to calculate the energy dependence from first principles, without the approximations built in the scattering theory, has been a major undertaking in nuclear engineering. Here, we demonstrate such a calculation method using the program OCLIMAX. In this work, our approach eliminates various approximations and limitations involved in a regular calculation with the LEAPR module of NJOY code, and the results are compared with available experimental and theoretical data. It is also demonstrated how additional insight can be obtained from the calculated full dynamical structure factor. The results reported here show the great potential and excellent platform provided by OCLIMAX for future development in the study of neutron thermalization in solid materials for different applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Electronic, magnetic, and structural properties of CoVMnSb: Ab initio study

We present computational results on electronic, magnetic, and structural properties of CoVMnSb, a quaternary Heusler alloy. Our calculations indicate that this material may crystallize in two energetically close structural phases: inverted and regular cubic. The inverted cubic phase is the ground state, with ferrimagnetic alignment, and around 80% spin polarization. Despite having a relatively large bandgap in the minority-spin channel close to the Fermi level, this phase does not undergo a half-metallic transition under pressure. This is explained by the “pinning” of the Fermi level at the minority-spin states at the Γ point. At the same time, the regular cubic phase is half-metallic and retains its perfect spin polarization under a wide range of mechanical strain. Transition to a regular cubic phase may be attained by applying uniform pressure (but not biaxial strain). In practice, this pressure may be realized by an atomic substitution of non-magnetic atoms (Sb) with another non-magnetic atom (Si) of a smaller radius. Our calculations indicate that 25% substitution of Sb with Si results in a half-metallic regular cubic phase being the ground state. In addition, CoVMnSb0.5Si0.5 retains its half-metallic properties under a considerable range of mechanical pressure, as well as exhibits thermodynamic stability, thus making this alloy attractive for potential spintronic applications. We hope that the presented results will stimulate experimental efforts to synthesize this compound.

Physics↗

DOE OSTI2024

The reverse quantum limit and its implications for unconventional quantum oscillations in YbB 12

The quantum limit in a Fermi liquid, realized when a single Landau level is occupied in strong magnetic fields, gives rise to unconventional states, including the fractional quantum Hall effect and excitonic insulators. Stronger interactions in metals with nearly localized $f$-electron degrees of freedom increase the likelihood of these unconventional states. However, access to the quantum limit is typically impeded by the tendency of $f$-electrons to polarize in a strong magnetic field, consequently weakening the interactions. In this study, we propose that the quantum limit in such systems must be approached in reverse, starting from an insulating state at zero magnetic field. In this scenario, Landau levels fill in the reverse order compared to regular metals and are closely linked to a field-induced insulator-to-metal transition. We identify YbB 12 as a prime candidate for observing this effect and propose the presence of an excitonic insulator state near this transition.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Modifiers versus Channels: Creating Shape-Selective Catalysis of Metal Nanoparticles/Porous Nanomaterials

Shape-selective catalysis is of great research interest for its key role in chemical synthesis. Porous nanomaterials with uniform pore structures, have been proven as ideal supports for metal nanoparticles (MNPs) to generate efficient shape-selective catalysis. However, many widely used commercial irregular porous nanomaterials usually face the challenge to realize satisfactory shape selectivity due to the lack of molecular sieving structures. In this work, we report a concept of creating shape selectivity in MNPs/porous-MNPs using appropriate modifiers. The remained MNPs within the porous nanomaterials can cooperate with the channels to generate selectivity. Such strategy not only applies to regular porous nanomaterials (metal-organic frameworks, zeolites, etc.), but also can be extended to irregular porous nanomaterials (active carbon, P25, etc.). Potentially, the matching among different MNPs, corresponding modifiers and porous nanomaterials makes our strategy promising for wide applications in selective catalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Efficient numerical integration of thermal interaction rates

In many problems in particle cosmology, interaction rates are dominated by 2 ↔ 2 scatterings, or get a substantial contribution from them, given that 1 ↔ 2 and 1 ↔ 3 reactions are phase-space suppressed. We describe an algorithm to represent, regularize, and evaluate a class of thermal 2 ↔ 2 and 1 ↔ 3 interaction rates for general momenta, masses, chemical potentials, and helicity projections. A key ingredient is an automated inclusion of virtual corrections to 1 ↔ 2 scatterings, which eliminate logarithmic and double-logarithmic IR divergences from the real 2 ↔ 2 and 1 ↔ 3 processes. We also review thermal and chemical potential induced contributions that require resummation if plasma particles are ultrarelativistic.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

The reverse quantum limit and its implications for unconventional quantum oscillations in YbB12

Abstract The quantum limit in a Fermi liquid, realized when a single Landau level is occupied in strong magnetic fields, gives rise to unconventional states, including the fractional quantum Hall effect and excitonic insulators. Stronger interactions in metals with nearly localized f -electron degrees of freedom increase the likelihood of these unconventional states. However, access to the quantum limit is typically impeded by the tendency of f -electrons to polarize in a strong magnetic field, consequently weakening the interactions. In this study, we propose that the quantum limit in such systems must be approached in reverse, starting from an insulating state at zero magnetic field. In this scenario, Landau levels fill in the reverse order compared to regular metals and are closely linked to a field-induced insulator-to-metal transition. We identify YbB 12 as a prime candidate for observing this effect and propose the presence of an excitonic insulator state near this transition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Lie-algebraic Kähler sigma models with U(1) isotropy

We discuss various questions that emerge in connection with the Lie-algebraic deformation of the sigma model in two dimensions. First, we supersymmetrize the original model endowing it with the minimal and extended supersymmetries. Then we derive the general hypercurrent anomaly in both cases. In the latter case this anomaly is one-loop but is somewhat different from the standard expressions one can find in the literature because the target manifold is nonsymmetric. We also show how to introduce the twisted masses and the term, and study the Bogomol’nyi–Prasad–Sommerfield equation for instantons, in particular the value of the topological charge. Then we demonstrate that the second loop in the function of the nonsupersymmetric Lie-algebraic sigma model is due to an infrared effect. To this end we use a supersymmetric regularization. We also conjecture that the above statement is valid for higher loops too, similar to the parallel phenomenon in four-dimensional super-Yang-Mills. In the second part of the paper we develop a special dimensional reduction—namely, starting from the two-dimensional Lie-algebraic model we arrive at a quasi-exactly solvable quantum-mechanical problem of the Lamé type.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

A tailored deep eutectic solvent for high-yield conversion of poplar residues to bio-based building blocks at mild conditions

The generation of condensed compounds in a deep eutectic solvent (DES) pretreatment has become a paramount factor that inhibits the further conversion of pretreated cellulose-rich substrate. Here, this study proposed a novel three-constituent DES system for the value-added utilization of poplar residues by efficiently suppressing hemicellulose derived condensation reactions. The results showed that adding ethanol into the DES system could effectively deconstruct the recalcitrant structure of poplar residues, simultaneously hindering the condensed composition formation. Under optimum ethanol addition (30%), the maximum xylan and lignin removal of 87.94% and 90.99% could be achieved, with more than 90.47% of glucan being preserved. With this, the glucan enzymatic hydrolysis efficiency was remarkably improved to 100% under the pretreatment conditions of 100 °C for 60 min, which was 84.69 times higher than the regular DES-pretreated poplar residues (1.17%). The working mechanism of relieving condensation reactions was characterized by Py-GC–MS, NMR, SEM, FTIR, and others. The results demonstrated that the addition of ethanol could regulate the degradation/conversion pathway of hemicellulose and lignin precipitation, thus inhibit the formation of condensed compounds during the pretreatment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Full treatment of the thrust distribution in single inclusive e + e – → h X processes

Extending the transverse momentum dependent factorization to thrust dependent observables entails a series of difficulties, ultimately associated to the behavior of soft radiation. As a consequence, the definition of the transverse momentum dependent functions has to be revised, while preserving (and possibly extending) their universality properties. Moreover, the regularization of the rapidity divergences generates non trivial correlations between rapidity and thrust. In this paper, we show how to deal with these correlations in a consistent treatment of the thrust dependence of e + e – → h X cross section, where the hadron transverse momentum is measured with respect to the thrust axis. In this framework all results obtained in the past few years properly fit together, leading to a remarkable phenomenological description of the experimental measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Three-dimensional reconstruction of x-ray emission volumes in magnetized liner inertial fusion from sparse projection data using a learned basis

The ability to visualize x-ray and neutron emission from fusion plasmas in 3D is critical to understand the origin of the complex shapes of the plasmas in experiments. Unfortunately, this remains challenging in experiments that study a fusion concept known as Magnetized Liner Inertial Fusion (MagLIF) due to a small number of available diagnostic views. Here, we present a basis function-expansion approach to reconstruct MagLIF stagnation plasmas from a sparse set of x-ray emission images. A set of natural basis functions is “learned” from training volumes containing quasi-helical structures whose projections are qualitatively similar to those observed in experimental images. Tests on several known volumes demonstrate that the learned basis outperforms both a cylindrical harmonic basis and a simple voxel basis with additional regularization, according to several metrics. Two-view reconstructions with the learned basis can estimate emission volumes to within 11% and those with three views recover morphology to a high degree of accuracy. The technique is applied to experimental data, producing the first 3D reconstruction of a MagLIF stagnation column from multiple views, providing additional indications of liner instabilities imprinting onto the emitting plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

MHD-induced SOL filaments and divertor heat flux striations in NSTX

Magnetohydrodynamic (MHD) activity in NSTX has been found to induce striated heat flux patterns at the divertor. These striated heat fluxes were previously interpreted as strike-point splitting caused by stochastic magnetic lobes near the X-point. Recent analysis indicates that MHD can directly generate filaments in the scrape-off layer (SOL), and the MHD-related filaments are consistent with the rotation and the frequency of the MHD. Magnetic field line tracing from the mid-plane of the MHD-related filaments to the divertor is consistent with the location of striated heat fluxes on the divertor. Additionally, gas puff imaging (GPI) diagnostic data shows the generation of regular blobs following each edge harmonic oscillation (EHO) toroidal mode number n = 1 perturbation, which indicates that the MHD perturbation directly drives the blobs. These blobs are responsible for the divertor heat flux striations. Here, this finding provides a new framework for understanding MHD-induced edge transport and has important implications for heat-load control in NSTX-U and future spherical tokamaks.

MHD↗

DOE OSTI2022

Improvements to the 115g/115m Cd Nuclear Data

Repeated measurements of the isotope 115m Cd by Los Alamos National Laboratory (LANL), Pacific Northwest National Laboratory (PNNL), and Atomic Weapons Establishment (AWE) have revealed a consistent discrepancy between beta and gamma counting. This suggests that a significant disparity exists between the true value of the isotope’s gamma-ray branching ratios and the values reported by the National Nuclear Data Center (NNDC). Enriched 114 Cd was irradiated using thermal neutrons to produce a high-purity 115m/g Cd source for counting and analysis using traditional singles gamma-ray spectroscopy, liquid scintillation counting, gas-proportional counting, and analysis using the new Gamma-Alpha-Beta-Gamma (GABγ) and Gamma-Alpha-Beta Radio-Isotope EvaLuator (GABRIEL) coincidence detection systems. Using the D-T fusion generated neutrons, reduces the wait time required to let the short-lived 115Cd ground state decay away by leveraging the 1:1 production ratio of the ground and metastable state of 115 Cd from the 115 In (n,p) reaction relative to the 10:1 production ratio of the ground and metastable isotopes of 115 Cd by thermal neutron capture on 114 Cd. However, the much higher flux in from the Washington State University TRIGA reactor is a more rapid method of production of 115m/gCd but requires access to enriched 114 Cd. High purity 115m/g Cd samples were analyzed on regular intervals over a period of 140-days. The activity of 115m Cd measured by beta and gamma show a clear discrepancy using the current best known gamma-ray decay branching ratios. Based on these observations and measurements, a set of new branching ratio recommendations have been produced using conventional counting techniques and the advanced GAB? and GABRIEL instruments.

07 ISOTOPE AND RADIATION SOURCES↗