Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse problems”

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

Isospin 0 and 2 two-pion scattering at physical pion mass using all-to-all propagators with periodic boundary conditions in lattice QCD

A study of two-pion scattering for the isospin channels, 𝐼 = 0 and 𝐼 = 2, using lattice QCD is presented. Möbius domain-wall fermions, on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions, are used for the lattice computations, which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD Collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of 𝐿 = 4.63 and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states—two pions with a nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the 𝐼 = 0 and 𝐼 = 2 channels are consistent with the Roy equation and chiral perturbation theory, though at this preliminary stage our errors for 𝐼 = 0 are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying 𝐾 → 𝜋⁢𝜋 decay, on physical-mass ensembles using the GEVP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Near-optimal ground state preparation

Preparing the ground state of a given Hamiltonian and estimating its ground energy are important but computationally hard tasks. However, given some additional information, these problems can be solved efficiently on a quantum computer. We assume that an initial state with non-trivial overlap with the ground state can be efficiently prepared, and the spectral gap between the ground energy and the first excited energy is bounded from below. With these assumptions we design an algorithm that prepares the ground state when an upper bound of the ground energy is known, whose runtime has a logarithmic dependence on the inverse error. When such an upper bound is not known, we propose a hybrid quantum-classical algorithm to estimate the ground energy, where the dependence of the number of queries to the initial state on the desired precision is exponentially improved compared to the current state-of-the-art algorithm proposed in [Ge et al. 2019]. These two algorithms can then be combined to prepare a ground state without knowing an upper bound of the ground energy. We also prove that our algorithms reach the complexity lower bounds by applying it to the unstructured search problem and the quantum approximate counting problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A three-dimensional laser ray-tracing methodology for radiation-hydrodynamics simulations

We report on a methodology for performing laser ray-tracing in three spatial dimensions for radiation-hydrodynamics simulation codes. Our method, which is an extension of that developed in Haines et al., Comput. Fluids 201, 104478 (2020), utilizes an automatically generated separate mesh for the laser ray-tracing from the radiation-hydrodynamics mesh. This enables the laser mesh to be tailored to minimize ray noise with significantly fewer rays than would be required when the ray-tracing is performed on the radiation-hydrodynamics mesh, primarily by allowing the use of high-aspect-ratio cells that are not suitable for hydrodynamics solvers. For a planar target, we show that our method provides a ≈ 100× reduction in computational expense to achieve a fixed level of ray noise relative to ray-tracing directly on the radiation-hydrodynamics mesh. The relatively low ray requirement also enables efficient computation of cross-beam energy transfer. Each cell in the logically cubic laser mesh is a non-convex dodecahedron with triangular sides, and numerical integration of the ray trajectories and inverse bremsstrahlung is performed by mapping each cell to the unit cube. We will describe our methodology in detail as well as its implementation in the xRAGE radiation-hydrodynamics code, discuss performance, and present the results from applying the methodology to test problems with analytic solutions for laser ray-tracing through a quadratic density gradient with an analytic solution as well as for a laser-driven heat front. In 3D radiation-hydrodynamics simulations of laser-driven experiments performed on the National Ignition Facility, laser ray-tracing with our methodology uses less than 1% of total computational time while introducing acceptably low levels of ray noise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Temperature of a tungsten spallation target

We discuss the heat load inside a neutron spallation target made of tungsten hit by an 800 MeV proton beam. This problem is important for the Neutron Target project, where a spallation target will be surrounded by an extended graphite or beryllium cube to create a neutron target. The idea of a free-neutron target enables the measurement of neutron-induced reactions in inverse kinematics. This idea is part of the LANSCE strategy to stay a worldwide leader for neutron-induced research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

AutoPhaseNN: unsupervised physics-aware deep learning of 3D nanoscale Bragg coherent diffraction imaging

Abstract The problem of phase retrieval underlies various imaging methods from astronomy to nanoscale imaging. Traditional phase retrieval methods are iterative and are therefore computationally expensive. Deep learning (DL) models have been developed to either provide learned priors or completely replace phase retrieval. However, such models require vast amounts of labeled data, which can only be obtained through simulation or performing computationally prohibitive phase retrieval on experimental datasets. Using 3D X-ray Bragg coherent diffraction imaging (BCDI) as a representative technique, we demonstrate AutoPhaseNN, a DL-based approach which learns to solve the phase problem without labeled data. By incorporating the imaging physics into the DL model during training, AutoPhaseNN learns to invert 3D BCDI data in a single shot without ever being shown real space images. Once trained, AutoPhaseNN can be effectively used in the 3D BCDI data inversion about 100× faster than iterative phase retrieval methods while providing comparable image quality.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An FFT-based approach for Bloch wave analysis: application to polycrystals

A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Improved surrogates in inertial confinement fusion with manifold and cycle consistencies

Neural networks have become the method of choice in surrogate modeling because of their ability to characterize arbitrary, high-dimensional functions in a data-driven fashion. This paper advocates for the training of surrogates that are 1) consistent with the physical manifold, resulting in physically meaningful predictions, and 2) cyclically consistent with a jointly trained inverse model; i.e., backmapping predictions through the inverse results in the original input parameters. We find that these two consistencies lead to surrogates that are superior in terms of predictive performance, are more resilient to sampling artifacts, and tend to be more data efficient. Using inertial confinement fusion (ICF) as a test-bed problem, we model a one-dimensional semianalytic numerical simulator and demonstrate the effectiveness of our approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A new process-based and scale-aware desert dust emission scheme for global climate models – Part I: Description and evaluation against inverse modeling emissions

Desert dust accounts for most of the atmosphere's aerosol burden by mass and produces numerous important impacts on the Earth system. However, current global climate models (GCMs) and land-surface models (LSMs) struggle to accurately represent key dust emission processes, in part because of inadequate representations of soil particle sizes that affect the dust emission threshold, surface roughness elements that absorb wind momentum, and boundary-layer characteristics that control wind fluctuations. Furthermore, because dust emission is driven by small-scale (~ 1 km or smaller) processes, simulating the global cycle of desert dust in GCMs with coarse horizontal resolutions (~ 100 km) presents a fundamental challenge. This representation problem is exacerbated by dust emission fluxes scaling nonlinearly with wind speed above a threshold wind speed that is sensitive to land-surface characteristics. Here, we address these fundamental problems underlying the simulation of dust emissions in GCMs and LSMs by developing improved descriptions of (1) the effect of soil texture on the dust emission threshold, (2) the effects of nonerodible roughness elements (both rocks and green vegetation) on the surface wind stress, and (3) the effects of boundary-layer turbulence on driving intermittent dust emissions. We then use the resulting revised dust emission parameterization to simulate global dust emissions in a standalone model forced by reanalysis meteorology and land-surface fields. We further propose (4) a simple methodology to rescale lower-resolution dust emission simulations to match the spatial variability of higher-resolution emission simulations in GCMs. The resulting dust emission simulation shows substantially improved agreement against regional dust emissions observationally constrained by inverse modeling. We thus find that our revised dust emission parameterization can substantially improve dust emission simulations in GCMs and LSMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Enhanced relaxed physical factorization preconditioner for coupled poromechanics

The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled poromechanics. For its application, however, it is necessary to invert blocks with the algebraic form C^ = (C + βFF T ), where C is a symmetric positive definite matrix, FF T a rank-deficient term, and β a real non-negative coefficient. The inversion of C^, performed in an inexact way, can become unstable for large values of β, as it usually occurs at some stages of a full poromechanical simulation. In this work, we propose a family of algebraic techniques to stabilize the inexact solve with C^. This strategy can prove useful in other problems as well where such an issue might arise, such as augmented Lagrangian preconditioning techniques for Navier-Stokes or incompressible elasticity. First, we introduce an iterative scheme obtained by a natural splitting of matrix C^. Second, we develop a technique based on the use of a proper projection operator annihilating the near-kernel modes of C^. Both approaches give rise to a novel class of preconditioners denoted as Enhanced RPF (ERPF). Furthermore, effectiveness and robustness of the proposed algorithms are demonstrated in both theoretical benchmarks and real-world large-size applications, outperforming the native RPF preconditioner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tomographic optical emission spectroscopy of an atmospheric pressure plasma jet and surface ionization waves on planar and structured surfaces

Here in this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N$_2^+$, and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N$_2^+$ and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N$_2^+$ that spatially converged to form the 'plasma bullet' or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.

plasma interactions with complex surfaces↗
DOE OSTI2023

Geometrically frustrated interactions drive structural complexity in amorphous calcium carbonate

Amorphous calcium carbonate is an important precursor for biomineralization in marine organisms. Key outstanding problems include understanding the structure of amorphous calcium carbonate and rationalizing its metastability as an amorphous phase. Here we report high-quality atomistic models of amorphous calcium carbonate generated using state-of-the-art interatomic potentials to help guide fits to X-ray total scattering data. Exploiting a recently developed inversion approach, we extract from these models the effective Ca∙∙∙Ca interaction potential governing the structure. This potential contains minima at two competing distances, corresponding to the two different ways that carbonate ions bridge Ca 2+ -ion pairs. We reveal an unexpected mapping to the Lennard-Jones–Gauss model normally studied in the context of computational soft matter. The empirical model parameters for amorphous calcium carbonate take values known to promote structural complexity. We thus show that both the complex structure and its resilience to crystallization are actually encoded in the geometrically frustrated effective interactions between Ca 2+ ions.

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

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. For this work, we have implemented CSP in a C++ open-source library CSPlib using the Kokkos parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Beam focusing and consequences for Doppler backscattering measurements

The phenomenon of focusing of microwave beams in a plasma near a turning-point caustic is discussed by exploiting the analytical solution to the Gaussian beam-tracing equations in the two-dimensional (2-D) linear-layer problem. The location of maximum beam focusing and the beam width at that location are studied in terms of the beam initial conditions. This focusing must be taken into account to interpret Doppler backscattering (DBS) measurements. We find that the filter function that characterises the scattering intensity contribution along the beam path through the plasma is inversely proportional to the beam width, predicting enhanced scattering from the beam focusing region. We show that the DBS signal enhancement for decreasing incident angles between the beam path and the density gradient is due to beam focusing and not due to forward scattering, as was originally proposed by (Gusakov et al., (Plasma Phys. Contr. Fusion, vol. 56, 2014, p. 0250092014, 2017); Plasma Phys. Rep. vol. 43(6), 2017, pp. 605–613). The analytic beam model is used to predict the measurement of the k y density-fluctuation wavenumber power spectrum via DBS, showing that, in an NSTX-inspired example, the spectral exponent of the turbulent, intermediate-to-high k y density-fluctuation spectrum might be quantitatively measurable via DBS, but not the spectral peak corresponding to the driving scale of the turbulent cascade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Polynomial-time preparation of low-temperature Gibbs states for two-dimensional toric code

In this work, we propose a polynomial-time algorithm for preparing the Gibbs state of the two-dimensional toric code Hamiltonian at any temperature, starting from any initial state, significantly improving upon prior estimates that suggested exponential scaling with inverse temperature. We prove that fast mixing at low temperature for the two-dimensional toric code can be achieved by augmenting local jump operators with simple global jump operators, which enable efficient transitions between logical sectors. To establish tight lower bounds on the spectral gap, we introduce a new reduction method that eventually maps the problem to estimating the spectral gap of a perturbed graph Laplacian on a stair graph. Our proof also shows that the Lindblad dynamics with a digitally implemented low-temperature local Davies generator is able to efficiently drive the quantum state toward the ground state manifold.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Future Atmospheric Rivers and Impacts on Precipitation: Overview of the ARTMIP Tier 2 High–Resolution Global Warming Experiment

Atmospheric rivers (ARs) are long, narrow synoptic scale weather features important for Earth’s hydrological cycle typically transporting water vapor poleward, delivering precipitation important for local climates. Understanding ARs in a warming climate is problematic because the AR response to climate change is tied to how the feature is defined. The Atmospheric River Tracking Method Intercomparison Project (ARTMIP) provides insights into this problem by comparing 16 atmospheric river detection tools (ARDTs) to a common data set consisting of high resolution climate change simulations from a global atmospheric general circulation model. ARDTs mostly show increases in frequency and intensity, but the scale of the response is largely dependent on algorithmic criteria. Across ARDTs, bulk characteristics suggest intensity and spatial footprint are inversely correlated, and most focus regions experience increases in precipitation volume coming from extreme ARs. The spread of the AR precipitation response under climate change is large and dependent on ARDT selection.

54 ENVIRONMENTAL SCIENCES↗