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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Elastic nucleon-pion scattering at m π = 200 MeV from lattice QCD

Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with N f = 2 + 1 dynamical quark flavors and m π = 200 MeV. The s-wave scattering lengths with both total isospins I = 1/2 and I = 3/2 are inferred from the finite-volume spectrum below the inelastic threshold together with the I = 3/2 p-wave containing the Δ(1232) resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels, enabling a determination of the I = 3/2 scattering length with statistical errors below 5%, while the I = 1/2 scattering length is somewhat less precisely evaluated. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a step toward future computations down to physical light quark masses with multiple lattice spacings and volumes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Generalized Wigner theorem for noninvertible symmetries

In this article, we establish the conditions under which a conservation law associated with a non-invertible operator may be realized as a symmetry in quantum physics. As established by Wigner, all quantum symmetries must be represented by either unitary or antiunitary transformations. Relinquishing an implicit assumption of invertibility, we demonstrate that the fundamental invariance of quantum transition probabilities under the application of symmetries mandates that all non-invertible symmetries may only correspond to projective unitary or antiunitary transformations, i.e., partial isometries. This extends the notion of physical states beyond conventional rays in Hilbert space to equivalence classes in an extended, gauged Hilbert space, thereby broadening the traditional understanding of symmetry transformations in quantum theory. Our generalized theorem applies irrespective of the origin of the (non)invertible symmetry, holds in arbitrary spatial dimensions, and is independent of the Hamiltonian or action. We explore its physical consequences and, using simple model systems, illustrate how the distinction between invertible and non-invertible symmetries can sometimes be tied to the choice of boundary conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Ionically Active MXene Nanopore Actuators

Reversible electrochemical intercalation of cations into the interlayer space of 2D materials induces tunable physical and chemical properties in them. In MXenes, a large class of recently developed 2D carbides and nitrides, low intercalation energy, high storage capacitance, and reversible intercalation of various cations have led to their improved performance in sensing and energy storage applications. In this work, a coupled nanopore-actuator system where an ultrathin free-standing MXene film serves as a nanopore support membrane and ionically active actuator is reported. In this system, the contactless MXene membrane in the electric field affects the cation movement in the field through their (de)intercalation between individual MXene flakes. This results in reversible swelling and contraction of the membrane monitored by ionic conductance through the nanopore. This unique nanopore coupled to a mechanical actuation system could provide new insights into designing single-molecule biosensing platforms at the nanoscale.

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

Spatial‐Uniformity–Driven Bayesian Optimization for Rapid Development of Printed Perovskite Solar Cells

Printed metal halide perovskites can enable rapid, roll-to-roll manufacturing of a broad class of optoelectronics—flexible solar cells and imagers among them—while promising cost and speed advantages over incumbent silicon. However, though current methods offer high throughput and patterning capabilities, perovskite films’ spatial heterogeneity remains a challenge for large-area devices. Here, a spatial-uniformity-driven Bayesian optimization (BO) approach is leveraged to accelerate the development of printed perovskite solar cells and improve large-area device performance. Using a BO surrogate model, a 6D design space of ink chemistry and printing physics is explored via extensive iterative experimentation (≈100) informed by an objective function capturing spatial photoluminescence (PL) variance. It is discovered that optimizing for uniformity drives rapid advances in photovoltaic performance, yielding ≈20% power conversion efficiency (PCE) for small area (0.134 cm 2 ) devices and > 16% for large area (1 cm 2 ) devices. This machine-learning approach simultaneously enables rheological comparison of ink formulations that accelerate the leveling of Saffman-Taylor artifacts and improve film uniformity. Here, this showcases uniformity-driven BO as an efficient approach for uncovering the key printing physics and mitigating spatial heterogeneity to enable device scaling beyond small cell areas.

14 SOLAR ENERGY↗
DOE OSTI2025

The no boundary density matrix

We discuss a no-boundary proposal for a subregion of the universe. In the classical approximation, this density matrix involves finding a specific classical solution of the equations of motion with no boundary. Beyond the usual no boundary condition at early times, we also have another no boundary condition in the region we trace out. We can find the prescription by starting from the usual Hartle-Hawking proposal for the wavefunction on a full slice and tracing out the unobserved region in the classical approximation. We discuss some specific subregions and compute the corresponding solutions. These geometries lead to phenomenologically unacceptable probabilities, as expected.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Estimating building occupancy: a machine learning system for day, night, and episodic events

Building occupancy research increasingly emphasizes understanding the social and physical dynamics of how people occupy space. Opportunities in the open source domain including social media, Volunteered Geographic Information, crowdsourcing, and sensor data have proliferated, resulting in the exploration of building occupancy dynamics at varying spatiotemporal scales. At Oak Ridge National Laboratory, research into building occupancies through the development of a global learning framework that accommodates exploitation of open source authoritative sources, including governmental census and surveys, journal articles, real estate databases, and more, to report national and subnational building occupancies across the world continues through the Population Density Tables (PDT) project. This probabilistic learning system accommodates expert knowledge, experience, and open-source data to capture local, socioeconomic, and cultural information about human activity. It does so through a systematic process of data harmonization techniques in the development of observation models for over 50 building types to dynamically update baseline estimates and report probabilistic diurnal and episodic building occupancy estimates. This discussion will explore how PDT is implemented at scale and expanded based on the development of observation model classes and will explain how to interpret and spatially apply the reported probability occupancy estimates and uncertainty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

TTDFT: A GPU accelerated Tucker tensor DFT code for large-scale Kohn-Sham DFT calculations

We present the Tucker tensor DFT (TTDFT) code which uses a tensor-structured algorithm with graphic processing unit (GPU) acceleration for conducting ground-state DFT calculations on large-scale systems. The Tucker tensor DFT algorithm uses a localized Tucker tensor basis computed from an additive separable approximation to the Kohn-Sham Hamiltonian. The discrete Kohn-Sham problem is solved using Chebyshev filtered subspace iteration method that relies on matrix-matrix multiplications of a sparse symmetric Hamiltonian matrix and a dense wavefunction matrix, expressed in the localized Tucker tensor basis. These matrix-matrix multiplication operations, which constitute the most computationally intensive step of the solution procedure, are GPU accelerated providing ~8-fold GPU-CPU speedup for these operations on the largest systems studied. In conclusion, the computational performance of the TTDFT code is presented using benchmark studies on aluminum nano-particles and silicon quantum dots with system sizes ranging up to ~7,000 atoms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Flux and estimated spectra from a low-intensity laser-driven X-ray source

Laser-driven X-rays as probes for high-energy-density physics spans an extremely large parameter space with laser intensities varying by 8 orders of magnitude. We have built and characterized a soft X-ray source driven by a modest intensity laser of 4 × 10 13 W/cm 2 . Emitted X-rays were measured by diamond radiation detectors and a filtered soft X-ray camera. A material-dependence study on Al, Ti, stainless steel alloy 304, Fe, Cu and Sn targets indicated 5-μm-thick Cu foils produced the highest X-ray yield. X-ray emission in the laser direction and emission in the reverse direction depend strongly on the foil material and the thickness due to the opacity and hydrodynamic disassembly time. The time-varying X-ray signals are used to measure the material thinning rate and is found to be ~1.5 μm/ns for the materials tested implying thermal temperature around 0.6 eV. The X-ray spectra from Cu targets peaks at ~2 keV with no emission >4 keV and was estimated using images with eight different foil filters. One-dimensional hydrodynamic and spectral calculations using HELIOS-CR provide qualitative agreement with experimental results. Modest intensity lasers can be an excellent source for nanosecond bursts of soft X-rays.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

From calibration to parameter learning: Harnessing the scaling effects of big data in geoscientific modeling

Abstract The behaviors and skills of models in many geosciences (e.g., hydrology and ecosystem sciences) strongly depend on spatially-varying parameters that need calibration. A well-calibrated model can reasonably propagate information from observations to unobserved variables via model physics, but traditional calibration is highly inefficient and results in non-unique solutions. Here we propose a novel differentiable parameter learning (dPL) framework that efficiently learns a global mapping between inputs (and optionally responses) and parameters. Crucially, dPL exhibits beneficial scaling curves not previously demonstrated to geoscientists: as training data increases, dPL achieves better performance, more physical coherence, and better generalizability (across space and uncalibrated variables), all with orders-of-magnitude lower computational cost. We demonstrate examples that learned from soil moisture and streamflow, where dPL drastically outperformed existing evolutionary and regionalization methods, or required only ~12.5% of the training data to achieve similar performance. The generic scheme promotes the integration of deep learning and process-based models, without mandating reimplementation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Highly sensitive 2D X-ray absorption spectroscopy via physics informed machine learning

Abstract Improving the spatial and spectral resolution of 2D X-ray near-edge absorption structure (XANES) has been a decade-long pursuit to probe local chemical reactions at the nanoscale. However, the poor signal-to-noise ratio in the measured images poses significant challenges in quantitative analysis, especially when the element of interest is at a low concentration. In this work, we developed a post-imaging processing method using deep neural network to reliably improve the signal-to-noise ratio in the XANES images. The proposed neural network model could be trained to adapt to new datasets by incorporating the physical features inherent in the latent space of the XANES images and self-supervised to detect new features in the images and achieve self-consistency. Two examples are presented in this work to illustrate the model’s robustness in determining the valence states of Ni and Co in the LiNi x Mn y Co 1-x-y O 2 systems with high confidence.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A modal wave-packet model for the multi-mode Richtmyer–Meshkov instability

In this work, a model for multimode perturbations subject to the Richtmyer–Meshkov (RM) instability is presented and compared with simulations and experiments for conditions relevant to inertial confinement fusion. The model utilizes the single mode response to the RM impulse whereby its amplitude h(k, t) first grows with an initial velocity V 0 ∝ kh(k, 0) that eventually decays in time as 1/kV 0 t. Both the growth and saturation stages are subject to nonlinearities since they depend explicitly on the initial amplitude. However, rather than using the individual mode amplitude h(k, t), nonlinearity is taken to occur when the root-mean-square amplitude h rms (k, t) of a wave-packet within wavenumbers k ± δk becomes comparable to 1/k. This is done because nearby sidebands can act in unison for an auto-correlation distance 1/δk beyond nonlinearity as observed in the beam-plasma instability. Thus, the nonlinear saturation amplitude for each mode is reduced from the usual 1/k by a phase space factor that depends on the physical dimensionality, as in the Haan model for the Rayleigh–Taylor instability. In addition, for RM, the average value of kh rms for the initial spectrum is used to calculate a nonlinear factor F NL that reduces V 0 , as observed for single modes. For broadband perturbations, the model describes self-similar growth ∝t θ as successively longer wavelength modes reach saturation. The growing and saturated modes must be discerned because only the former promote θ and are enhanced by reshock and spherical convergence. All of these flows are described here by the model in good agreement with simulations and experiments.

42 ENGINEERING↗
DOE OSTI2024

Perspectives and progress on wurtzite ferroelectrics: Synthesis, characterization, theory, and device applications

Wurtzite ferroelectrics are an emerging material class that expands the functionality and application space of wide bandgap semiconductors. Promising physical properties of binary wurtzite semiconductors include a large, reorientable spontaneous polarization, direct band gaps that span from the infrared to ultraviolet, large thermal conductivities and acoustic wave velocities, high mobility electron and hole channels, and low optical losses. The ability to reverse the polarization in ternary wurtzite semiconductors at room temperature enables memory and analog type functionality and quasi-phase matching in optical devices and boosts the ecosystem of wurtzite semiconductors, provided the appropriate combination of properties can be achieved for any given application. In this article, advances in the design, synthesis, and characterization of wurtzite ferroelectric materials and devices are discussed. Highlights include: the direct and quantitative observation of polarization reversal of ~135 μC/cm 2 charge in Al 1-x B x N via electron microscopy, Al 1-x B x N ferroelectric domain patterns poled down to 400 nm in width via scanning probe microscopy, and full polarization retention after over 1000 h of 200 °C baking and a 2× enhancement relative to ZnO in the nonlinear optical response of Zn 1-x Mg x O. In conclusion, the main tradeoffs, challenges, and opportunities in thin film deposition, heterostructure design and characterization, and device fabrication are overviewed.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2025

A novel closed-form inversion of the convection–diffusion equation for rapid convection, diffusion, and source profile estimation

To simplify and routinize particle transport analysis in fusion devices, a novel closed form linear inversion of the 1-D convection diffusion equation to estimate diffusion and convection profiles D(r ⃗ ), v(r ⃗ ) and source distribution s(r ⃗ ), of a single species from measured data is derived and demonstrated on synthetic data. Profile estimates of D(r ⃗ ), v(r ⃗ ), s(r ⃗ ) and their uncertainties are given as a matrix expression constructed directly from the incoming density data of the transported species in space and time, as well as physics assumptions such as particle conservation and experimental geometry. The derived matrix expression can be applied to a pumped or non-pumped recycling species, or a non-recycling species that is effectively “pumped” by plasma-facing surfaces.

Hinson, Edward [ORNL] (ORCID:000000019713140X)↗
DOE OSTI2022

Efficient computation of N -point correlation functions in D dimensions

We present efficient algorithms for computing the N-point correlation functions (NPCFs) of random fields in arbitrary D-dimensional homogeneous and isotropic spaces. Such statistics appear throughout the physical sciences and provide a natural tool to describe stochastic processes. Typically, algorithms for computing the NPCF components have $\mathscr O$(n N ) complexity (for a dataset containing n particles); their application is thus computationally infeasible unless N is small. By projecting the statistic onto a suitably defined angular basis, we show that the estimators can be written in a separable form, with complexity $\mathscr O$(n 2 ) or $\mathscr O$(n g log n g ) if evaluated using a Fast Fourier Transform on a grid of size n g . Our decomposition is built upon the D-dimensional hyperspherical harmonics; these form a complete basis on the (D – 1) sphere and are intrinsically related to angular momentum operators. Concatenation of (N – 1) such harmonics gives states of definite combined angular momentum, forming a natural separable basis for the NPCF. As N and D grow, the number of basis components quickly becomes large, providing a practical limitation to this (and all other) approaches: However, the dimensionality is greatly reduced in the presence of symmetries; for example, isotropic correlation functions require only states of zero combined angular momentum. We provide a Julia package implementing our estimators and show how they can be applied to a variety of scenarios within cosmology and fluid dynamics. The efficiency of such estimators will allow higher-order correlators to become a standard tool in the analysis of random fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Fast emulation of quantum three-body scattering

Here, we develop a class of emulators for solving quantum three-body scattering problems. They are based on combining the variational method for scattering observables and the recently proposed eigenvector continuation concept. The emulators are first trained by the exact scattering solutions of the governing Hamiltonian at a small number of points in its parameter space, and then employed to make interpolations and extrapolations in that space. Through a schematic nuclear-physics model with finite-range two and three-body interactions, we demonstrate the emulators to be extremely accurate and efficient. The computing time for emulation is on the scale of milliseconds (on a laptop), with relative errors ranging from 10 –13 to 10 –4 depending on the case. The emulators also require little memory. We argue that these emulators can be generalized to even more challenging scattering problems. Furthermore, this general strategy may be applicable for building the same type of emulators in other fields, wherever variational methods can be developed for evaluating physical models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗