Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space physics”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

DuctGPT: A Generative Transformer for Forward Screening of Ductile Refractory Multi-Principal Element Alloys

Designing ductile materials for extreme environments such as fusion reactors requires a deep understanding of the complex interplay between electronic structure, mechanical stability, and wide compositional space. Here, in this work, we introduce DuctGPT, a physics-informed, GPT-powered machine learning platform that enables rapid and accurate prediction of ductility across a wide range of refractory multi-principal element alloys (MPEAs). Trained on both experimental and high-fidelity computational data, DuctGPT integrates descriptors such as density of states at the Fermi level, elastic constants, and valence electron concentration to capture the fundamental mechanisms governing ductile versus brittle behavior. Using this framework, we screen over 1000 compositions in of body-centered cubic (BCC) MPEAs, including two new alloy classes, i.e., NbTa-rich (NbTa $>$ 50 at.%) NbTa-Ti-V and W-rich ($>$ 50 at.%) W-Ti-V MPEAs, to rapidly identify promising alloy compositions with enhanced ductility. Validation against experimental data confirms the model's ability to predict ductility with high fidelity and low uncertainty. By leveraging conversational AI and robust physical modeling, DuctGPT provides a blueprint for the next generation of alloy design assistants, enabling human-AI collaboration in the accelerated discovery of ductile, high-performance materials for fusion, aerospace, and advanced manufacturing.

AI/ML↗
DOE OSTI2020

A multiscale cohesive law for carbon fiber networks

Better predictive models of mechanical failure in low-weight heat shield composites would aid material certification for missions with aggressive atmospheric entry conditions. In this study, we develop such a model for the rapid engineering analysis of the failure limits of phenolic impregnated carbon ablator (PICA) - a leading heat shield material whose structural component is a carbon fiber network. We hypothesize inelastic deformation failure mechanisms and model their behavior using molecular dynamics simulations to calculate the binding energy. We then upscale this binding energy to the macroscale using a renormalization argument. The approach delivers insightful and reasonably accurate macroscale predictions that compare favorably to experiments. In application, the model is validated for a particular variety of PICA by comparison to experiment and would then be used to study design scenarios in different entry conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Crystal Growth of Quaternary Rare Earth Selenosilicates by Using the Flux-Assisted Boron Chalcogen Mixture Method: Investigation of Their Magnetic and Optical Properties

A series of rare earth magnesium selenosilicates, RE3Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd, Tb, Dy) were obtained as single crystal using the flux assisted boron chalcogen mixture (BCM) method. The structures of the crystals were determined by single-crystal X-ray diffraction. The RE3Mg 0.5 SiSe 7 series crystallizes in the hexagonal crystal system in the space group P6 3 . Polycrystalline powders were synthesized to perform physical property measurements. Magnetic measurements over the 2–300 K temperature range reveal that Ce 3 Mg 0.5 SiSe 7 and Gd 3 Mg 0.5 SiSe 7 exhibit paramagnetic behavior with negative Weiss constants (θ W = −14.50, θ W = −6.13 K, respectively). The optical properties of RE 3 Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd) were measured by ultraviolet–visible (UV–vis) diffuse reflectance. Density functional theory (DFT) electronic structure calculations were performed. Furthermore, a second harmonic generation measurement was performed on a polycrystalline powder of Ce 3 Mg 0.5 SiSe 7 and was found to be SHG active with an efficiency of 0.11 times the standard potassium dihydrogen phosphate (KDP).

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

Crystal Growth of Quaternary RE 2 EuSi 2 S 8 ( RE = Ce–Nd, Sm, Gd, Tb) Using Flux-Assisted Boron Chalcogen Mixture (BCM) Method: Investigation of Magnetic and Luminescence Properties

A series of quaternary rare-earth containing thiosilicates with the general formula RE 2 EuSi 2 S 8 (RE = Ce–Nd, Sm, Gd, Tb) has been synthesized via the flux-assisted boron chalcogen mixture (BCM) crystal growth method. High-quality single crystals were obtained, and their crystal structures were determined by single-crystal X-ray diffraction. The RE 2 EuSi 2 S 8 series crystallizes in the trigonal system, adopting the space group R-3c. Polycrystalline samples were employed for physical property measurements, including magnetic susceptibility measurements, UV–visible diffuse reflectance, and photoluminescent response. Magnetic data of RE 2 EuSi 2 S 8 ( RE = Ce, Nd, and Gd) were collected over the 2–300 K temperature range. The samples were paramagnetic behavior with negative Weiss constants (θ W = −11.05, −10.55, and −1.35 K respectively). Their thermal stability was investigated using thermogravimetric analysis (TGA). Optical band gaps, estimated from diffuse reflectance spectra, were determined to be 2.2(1) eV for Ce 2 EuSi 2 S 8 , 1.8(1) eV for Nd 2 EuSi 2 S 8 , and 1.7(1) eV for Gd 2 EuSi 2 S 8 respectively. Finally, photoluminescence measurements were collected on Ce 2 EuSi 2 S 8 and Tb 2 EuSi 2 S 8 single crystals.

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

The Illumination of Thunderclouds by Lightning: 1. The Extent and Altitude of Optical Lightning Sources

Abstract Optical space‐based lightning sensors including NOAA's Geostationary Lightning Mapper (GLM) detect lightning though its transient illumination of the surrounding clouds. What space‐based optical lightning sensors measure is influenced by the physical attributes of the light source, the location of the source within the cloud scene, and the spatial variations in cloud composition. We focus on the lightning channels that serve as optical sources for GLM groups and flashes in this first part of our thundercloud illumination study. We match Lightning Mapping Array (LMA) sources with GLM groups and flashes during two thunderstorms to examine channel segments that are active during optical emission. We find that in each storm, the LMA sources matched with GLM groups are small (median: 3 km) compared to GLM pixels (nominal: 8 km), and preferentially come from high altitudes in the cloud (>8–10 km). The detection advantage for high‐altitude sources permits GLM to resolve faint optical pulses near the cloud top that might be missed from lower altitudes. However, the most energetic groups can be detected from all altitudes, and the largest groups largely originate at low altitudes. The relationship between group energy and illuminated area depends on flash development within the cloud medium, and flash development into different cloud regions can be identified by tracking GLM metrics of cloud illumination over time.

58 GEOSCIENCES↗
DOE OSTI2020

Inverse methods for design of soft materials

Functional soft materials, comprising colloidal and molecular building blocks that self-organize into complex structures as a result of their tunable interactions, enable a wide array of technological applications. Inverse methods provide a systematic means for navigating their inherently high-dimensional design spaces to create materials with targeted properties. Furthermore, while multiple physically motivated inverse strategies have been successfully implemented in silico, their translation to guiding experimental materials discovery has thus far been limited to a handful of proof-of-concept studies. In this perspective, we discuss recent advances in inverse methods for design of soft materials that address two challenges: (1) methodological limitations that prevent such approaches from satisfying design constraints and (2) computational challenges that limit the size and complexity of systems that can be addressed. Strategies that leverage machine learning have proven particularly effective, including methods to discover order parameters that characterize complex structural motifs and schemes to efficiently compute macroscopic properties from the underlying structure. We also highlight promising opportunities to improve the experimental realizability of materials designed computationally, including discovery of materials with functionality at multiple thermodynamic states, design of externally directed assembly protocols that are simple to implement in experiments, and strategies to improve the accuracy and computational efficiency of experimentally relevant models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3–4 orders of magnitude fewer Hamiltonian evaluations and 2–3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2–3/1–2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ~1.3×10 9 parameters to obtain accurate dynamics.

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

PolyODENet: Deriving mass-action rate equations from incomplete transient kinetics data

Kinetics of a reaction network that follows mass-action rate laws can be described with a system of ordinary differential equations (ODEs) with polynomial right-hand side. However, it is challenging to derive such kinetic differential equations from transient kinetic data without knowing the reaction network, especially when the data are incomplete due to experimental limitations. We introduce a program, PolyODENet, toward this goal. Based on the machine-learning method Neural ODE, PolyODENet defines a generative model and predicts concentrations at arbitrary time. As such, it is possible to include unmeasurable intermediate species in the kinetic equations. Importantly, we have implemented various measures to apply physical constraints and chemical knowledge in the training to regularize the solution space. Here, using simple catalytic reaction models, we demonstrate that PolyODENet can predict reaction profiles of unknown species and doing so even reveal hidden parts of reaction mechanisms.

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

Efficient many-jet event generation with flow matching

We apply for the first time, to the best of our knowledge, the flow matching method to the problem of phase-space sampling for event generation in high-energy collider physics. By training the model to remap the random numbers used to generate the momenta and helicities of the scattering matrix elements as implemented in the portable partonic event generator pepper, we find substantial efficiency improvements in the studied processes. We focus our study on the highest final-state multiplicities in Drell-Yan and top-antitop pair production used in simulated samples for the Large Hadron Collider, which computationally are the most relevant ones. We find that the unweighting efficiencies improve by factors of 184 and 25, respectively, when compared to the standard approach of using a vegas-based optimization. We also compare continuous normalizing flows trained with flow matching against the previously studied normalizing flows based on coupling layers and find that the former leads to better results, faster training and a better scaling behavior across the studied multiplicity range, while the latter evaluate faster. When combining the advantages of both methods using the regflow approach, we find parton-level unweighted event generation walltime gains of about a factor of 10 at the highest final-state multiplicities.

Bothmann, E. [CERN; Gottingen U.] (ORCID:000000016↗
DOE OSTI2021

Superconductinglike response in a driven gapped bosonic system

In this study, we propose a mechanism for the photoinduced superconductinglike response recently reported in cuprates and other strongly correlated materials. This mechanism relies on quantum-fluctuating bosons consisting of electron pairs. With periodic drive, the electron pairs and vacancies of pairs form a coherent nonequilibrium condensate, different from conventional superconductors, yet showing a superconductinglike response in some regime even with dissipation. Unlike the case of driven fermionic bands, which results in the familiar Floquet bands with hybridization gaps, for driven bosons the “gap” opens up in the momentum direction, resulting in a resonant region in momentum space where the eigenvalues are complex. We give a simple physical argument for why this picture leads to a “perfect conductor” which exhibits superconductinglike frequency-dependent conductivity but no Meissner response. While our model is quite general, in the case of cuprates, a quantum-fluctuating pair density wave in the pseudogap region may serve as the origin of the quantum-fluctuating electron pairs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Light-quark connected intermediate-window contributions to the muon 𝑔 − 2 hadronic vacuum polarization from lattice QCD

We present a lattice-QCD calculation of the light-quark connected contribution to window observables associated with the leading-order hadronic vacuum polarization contribution to the anomalous magnetic moment of the muon, $a^{HVP,LO}_μ$. We employ the MILC Collaboration’s isospin-symmetric QCD gauge-field ensembles, which contain four flavors of dynamical highly improved staggered quarks with four lattice spacings between $a ≈ 0.06-0.15$ fm and close-to-physical quark masses. We consider several effective-field-theory-based schemes for finite volume and other lattice corrections and combine the results via Bayesian model averaging to obtain robust estimates of the associated systematic uncertainties. After unblinding, our final results for the intermediate and "W2" windows are $a^{U,W}_μ$(conn) = 206.6(1.0) x 10 -10 and $a^{U,W2}_μ$ (conn) = 100.7(3.2) x 10 -10 , respectively

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Decaying sterile neutrinos at short baselines

Long-standing anomalous experimental results from short-baseline neutrino experiments have persisted for decades. These results, when interpreted with one or more light sterile neutrinos, are inconsistent with numerous null results experimentally. However, if the sterile neutrino decays en route to the detector, this can mimic ν μ → ν e oscillation signals while avoiding many of these external constraints. We revisit this solution to the MiniBooNE and LSND puzzles in view of new data from the MicroBooNE experiment at Fermilab. Using MicroBooNE’s liquid-argon time-projection chamber search for an excess of ν e in the Booster beam, we derive new limits in two models’ parameter spaces of interest; where the sterile neutrino decays (I) via mixing with the active neutrinos, or (II) via higher-dimensional operators. We also provide an updated, comprehensive fit to the MiniBooNE neutrino- and antineutrino-beam data, including appearance ( ν e ) and disappearance ( ν μ ) channels. Despite alleviating the tension with muon neutrino disappearance experiments, we find that the latest MicroBooNE analysis rules out the decaying sterile neutrino solution in a large portion of the parameter space at more than 99% CL. Published by the American Physical Society 2024

Hostert, Matheus (ORCID:0000000295848877)↗
DOE OSTI2022

ExaAM: Metal additive manufacturing simulation at the fidelity of the microstructure

Additive manufacturing (AM), or 3D printing, of metals is transforming the fabrication of components, in part by dramatically expanding the design space, allowing optimization of shape and topology. However, although the physical processes involved in AM are similar to those of welding, a field with decades of experimental, modeling, simulation, and characterization experience, qualification of AM parts remains a challenge. The availability of exascale computational systems, particularly when combined with data-driven approaches such as machine learning, enables topology and shape optimization as well as accelerated qualification by providing process-aware, locally accurate microstructure and mechanical property models. We describe the physics components comprising the Exascale Additive Manufacturing simulation environment and report progress using highly resolved melt pool simulations to inform part-scale finite element thermomechanics simulations, drive microstructure evolution, and determine constitutive mechanical property relationships based on those microstructures using polycrystal plasticity. We report on implementation of these components for exascale computing architectures, as well as the multi-stage simulation workflow that provides a unique high-fidelity model of process–structure–property relationships for AM parts. In addition, we discuss verification and validation through collaboration with efforts such as AM-Bench, a set of benchmark test problems under development by a team led by the National Institute of Standards and Technology.

3D printing↗
DOE OSTI2023

The battle of clean and dirty qubits in the era of partial error correction

When error correction becomes possible it will be necessary to dedicate a large number of physical qubits to each logical qubit. Error correction allows for deeper circuits to be run, but each additional physical qubit can potentially contribute an exponential increase in computational space, so there is a trade-off between using qubits for error correction or using them as noisy qubits. In this work, we look at the effects of using noisy qubits in conjunction with noiseless qubits (an idealized model for error-corrected qubits), which we call the “clean and dirty” setup. We employ analytical models and numerical simulations to characterize this setup. Numerically we show the appearance of Noise-Induced Barren Plateaus (NIBPs), i.e., an exponential concentration of observables caused by noise, in an Ising model Hamiltonian variational ansatz circuit. We observe this even if only a single qubit is noisy and given a deep enough circuit, suggesting that NIBPs cannot be fully overcome simply by error-correcting a subset of the qubits. On the positive side, we find that for every noiseless qubit in the circuit, there is an exponential suppression in the concentration of gradient observables, showing the benefit of partial error correction. Finally, our analytical models corroborate these findings by showing that observables concentrate with a scaling in the exponent related to the ratio of dirty-to-total qubits.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Optical Test Bench (SYSM 5620 Final Report)

The National Ignition Facility (NIF) is a research facility that focuses on High Energy Density Science and Inertial Confinement Fusion. In the NIF, 192 laser beams are directed into a miniscule Target Capsule to generate an implosion – this is called a NIF Shot. Many NIF Users conduct experiments using different Targets and laser inputs to generate reactions that are of interest to the scientific community. Those reactions are monitored using various Target Diagnostic instruments, including Shot Time Photography – which provides crucial images of the Target as the implosion occurs. A subset of these NIF Users is less interested in the reaction itself and more interested in subjecting materials or devices to the intense environment that these reactions can create. Using the Target as a source of x-rays, neutrons, or blast loading for particular materials has been significant to the shock physics community and to those who work in extreme design spaces. However, the Material and Radiation Effects team – that performs these experiments on behalf of the NIF Users – has not been able to make good use of the existing Shot Time Photography Diagnostics. The STP cameras are focused on the Targets rather than the MRE Customers’ materials, the current system is not oriented to have a good view of the MRE Customers’ materials, and the cameras are limited in their capabilities (black and white, poor resolution, etc.). As MRE experiments increase in quantity and importance to the scientific community, the desire for visual monitoring of the experiments to substantiate models and inform other data is also rising.

42 ENGINEERING↗
DOE OSTI2021

Muon to electron conversion search in the presence of Al nuclei at the Fermilab Mu2e experiment: Motivation, Design and Progress

The Mu2e experiment aims to find charged lepton flavor violation(CLFV) by measuring the monochromatic electrons from $\mu^- N\rightarrow e^- N$ conversion with an unprecedented single event sensitivity of $3 \times 10^{-17}$. When completed the experiment will improve the current limit by four orders of magnitude and make a previously unexplored phase space available for the search for beyond the standard model physics. Physics motivation for the experiment is discussed and design choices for the detectors and various sub systems are explained in this document. The construction of the Mu2e experiment is well underway and commissioning will begin in late 2022 with first data taking in late 2025 and 4-5 years of data taking thereafter to reach our goal.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

The battle of clean and dirty qubits in the era of partial error correction

When error correction becomes possible it will be necessary to dedicate a large number of physical qubits to each logical qubit. Error correction allows for deeper circuits to be run, but each additional physical qubit can potentially contribute an exponential increase in computational space, so there is a trade-off between using qubits for error correction or using them as noisy qubits. In this work we look at the effects of using noisy qubits in conjunction with noiseless qubits (an idealized model for error-corrected qubits), which we call the "clean and dirty" setup. We employ analytical models and numerical simulations to characterize this setup. Numerically we show the appearance of Noise-Induced Barren Plateaus (NIBPs), i.e., an exponential concentration of observables caused by noise, in an Ising model Hamiltonian variational ansatz circuit. We observe this even if only a single qubit is noisy and given a deep enough circuit, suggesting that NIBPs cannot be fully overcome simply by error-correcting a subset of the qubits. On the positive side, we find that for every noiseless qubit in the circuit, there is an exponential suppression in concentration of gradient observables, showing the benefit of partial error correction. Finally, our analytical models corroborate these findings by showing that observables concentrate with a scaling in the exponent related to the ratio of dirty-to-total qubits.

97 MATHEMATICS AND COMPUTING↗