Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion methods”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Calorimeter with Bayesian unfolding of spectra of high-flux broadband x rays

We report the development of a multipurpose differential x-ray calorimeter with a broad energy bandwidth. The absorber architecture is combined with a Bayesian unfolding algorithm to unfold high energy x-ray spectra generated in high-intensity laser–matter interactions. Particularly, we show how to extract absolute energy spectra and how our unfolding algorithm can reconstruct features not included in the initial guess. The performance of the calorimeter is evaluated via Monte Carlo generated data. The method accuracy to reconstruct electron temperatures from bremsstrahlung is shown to be 5% for electron temperatures from 1 to 50 MeV. We study bremsstrahlung generated in solid target interaction showing an electron temperature of 0.56 ± 0.04 MeV for a 700 μm Ti titanium target and 0.53 ± 0.03 MeV for a 50 μm target. We investigate bremsstrahlung from a target irradiated by laser-wakefield accelerated electrons showing an endpoint energy of 551 ± 5 MeV, inverse Compton generated x rays with a peak energy of 1.1 MeV, and calibrated radioactive sources. The total energy range covered by all these sources ranges from 10 keV to 551 MeV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Estimation of constituent properties of concrete materials with an artificial neural network based method

Multi-scale models are developed for heterogeneous concrete materials to estimate their macroscopic mechanical properties in terms of micro-structural data. One crucial challenge of those models is the identification of local properties of constituent phases. In this paper, we present an efficient method based on Artificial Neural Networks (ANN). Typical concrete materials are taken as example. A macroscopic analytical strength criterion is established from three steps of nonlinear homogenization procedure. The macroscopic strength of materials is determined as a function of the frictional coefficient and cohesion of solid cement particles at nanometer scale, intra-particle pores, inter-particle pores and aggregates (inclusions). The objective is to identify the nanoscopic frictional coefficient and cohesion of cement particle from measured macroscopic values of uniaxial compression and tensile strengths. For this purpose, a numerical method based on the ANN is developed. With the analytical macroscopic strength criterion, sensitivity studies are first realized to identify the most important micro-structural parameters influencing the macroscopic strength of concrete. A simplified analytical macroscopic strength criterion is then proposed. A large dataset is further constructed through the inversion of the analytical strength criterion by using the aggregates volume fraction, porosity, macroscopic uniaxial tensile and compressive strengths as input variables and the frictional coefficient and cohesion of cement particles as output unknowns. An ANN model containing four hidden layers and 100 neurons in each layer is constructed and trained by using this dataset. Various types of validation of the ANN model are performed. It is found that the proposed ANN based model can effectively predict the frictional coefficient and cohesion of porous cement paste at the microscopic scale with a very good accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Source size of x rays from self-modulated laser wakefield accelerators

A comparative study of x-ray sources generated with different mechanisms from self-modulated laser wakefield acceleration (SM-LWFA) electrons was performed to compare the source size or spatial resolution for use in high energy density science applications. We examine the source size of betatron, inverse Compton scattering, and bremsstrahlung radiation with a Fresnel diffraction based formalism and a modified x-ray ray tracing model. We observe the dependence of source size on the radiation generation process, laser parameters, and compare to what is possible in other regimes of LWFA, as well as current methods. This information is significant as we begin to explore the use of light sources driven by SM-LWFA for use as a diagnostic at large-scale laser facilities where blowout regime LWFA is not possible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Structure studies of 13 Be from the 12 Be ⁢(𝑑,𝑝) reaction in inverse kinematics on a solid deuteron target

The low-lying structure of 13 Be has remained an enigma for decades. Despite numerous experimental and theoretical studies, large inconsistencies remain. Being both unbound and one neutron away from 14 Be, the heaviest bound beryllium nucleus, 13 Be is difficult to study through simple reactions with weak radioactive-ion beams or more complex reactions with stable-ion beams. Here, we present the results of a study using the 12 Be (𝑑,𝑝)⁢ 13 Be reaction in inverse kinematics using a 9.5 MeV per nucleon 12 Be beam from the ISAC-II facility. The solid deuteron target of IRIS was used to achieve an increased areal thickness compared to conventional deuterated polyethylene targets. The 𝑄-value spectrum below −4.4 MeV was analyzed using a Bayesian method with geant4 simulations. A three-point angular distribution with the same 𝑄-value gate was fit with a mixture of 𝑠- and 𝑝-wave, 𝑠- and 𝑑-wave, or pure 𝑝-wave transfer. The 𝑄-value spectrum was also compared with geant simulations obtained using the energies and widths of states reported in four previous works. As a result, it was found that our results are incompatible with works that revealed a wide 5/2 + resonance but shows better agreement with ones that reported a narrower width.

6 ≤ A ≤ 19↗
DOE OSTI2025

Quantitative Nonlinear Optical Polarimetry with High Spatial Resolution

Nonlinear optical microscopy such as in the optical second-harmonic generation (SHG) modality has become a popular tool today for probing materials in the physical and biological sciences. While imaging and spectroscopy are widely used in the microscopy mode, nonlinear polarimetry, which can shed light on materials’ symmetry and microstructure, is relatively underdeveloped. This is partly because quantitative analytical modeling of the optical SHG response for anisotropic crystals and films largely assumes low-numerical aperture (NA) focusing of light, where the plane-wave approximation is sufficient. Tight focusing provides unique benefits in revealing out-of-plane polarization responses, which cannot be detected by near-plane-wave illumination at normal incidence. Here, we outline a method for quantitatively analyzing SHG polarimetry measurements obtained under high-NA focusing within a microscope geometry. Experiments and simulations of a variety of standard samples, from single crystals to thin films, are in good agreement, including measured and simulated spatial SHG maps of ferroelectric domains. A solution to the inverse problem is demonstrated, where the spatial distribution of an SHG tensor with unknown tensor coefficient magnitudes is determined by experimentally measured polarimetry. The ability to extract the out-of-plane component of the nonlinear polarization in normal incidence is demonstrated, which can be valuable for high-resolution polarimetry of 2D materials, thin films, heterostructures, and uniaxial crystals with a strong out-of-plane response.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects

Researchers at New Mexico Tech, Los Alamos National Laboratory (LANL), Sandia National Laboratory (SNL), and other collaborators are developing a methodology to improve characterization of stress in the subsurface by means of a model-based inversion of reservoir engineering data, time lapse seismic measurements, and microseismicity. Historically, various geophysical techniques have been used in efforts to understand the state of stress in the subsurface through direct and indirect imaging of stress sensitive features (faults and fractures) and observations of transient stress related observations (time variant elastic moduli and microseismicity). Standard direct and indirect techniques for seismic fault and fracture imaging suffer from detectability limits due to reliance on high data quality and multiplicity. Variations in effective elastic properties from inversion of high quality seismic (Vertical Seismic Profile (VSP)) data may be used to infer stress near the wellbore through integration with independent experimental characterization of the stress-velocity relationship. Given a sufficiently robust observation network, microseismic emissions may be inverted to characterize focal mechanisms which, together with supporting assumptions and constraints, inform estimates of in-situ stress. While each of these techniques contributes in part to the characterization of stress within limited spatial and temporal domains, no one method provides an unambiguous stress measurement or a predictive capability over a site scale spatial extent. Challenges associated with solution non-uniqueness, measurement ambiguity, and irregular sampling may be greatly minimized through combination of one or more of independent measurements within a common framework in which realistic geological, hydrodynamic, geomechanical, and seismological constituent process models may act as constraints.

02 PETROLEUM↗
DOE OSTI2021

Evolutionary inverse design of defects at graphene 2D lateral interfaces

Grain boundaries (GBs) in two-dimensional (2D) materials often have a profound impact on various material properties from mechanical to optical to electronic, yet predicting all possible GB formations is a challenge. Here, we introduce a workflow based on an evolutionary algorithm for exploring possible GBs formed at a lateral 2D interface. In a departure from conventional genetic algorithm based structure optimization methods, we perform genetic operations in the near interface region that allow us to be computationally efficient. We benchmark our method using graphene, which is a well-studied 2D material with a wide range of point defects. An empirical potential was used as the surrogate of the evolutionary search. More than 11.5 × 106 structures in total were evaluated for 128 GB orientations, and for each orientation, the ten best structures are recorded. A subset of low energy GBs predicted by empirical potential based search was relaxed by first-principles calculations and used to validate the energetic rank order. With the validated formation energy, we rank-ordered the best 128 GB structures and performed a detailed statistical analysis of primitive rings to find the correlation between the ring distribution and the formation energy. We found that for low energy GBs (below 0.5eV/Å), Stone–Wales defects will dominate, while structures with a higher energy (0.5–1.1eV/Å) show an increasing population of heptagons and nine-membered rings to form seven-nine defect pairs. For structures with energy higher than 1.1eV/Å, the percentage of octagons and nine-membered rings increases, which indicates that these two types of rings are not energetically favorable. Our proposed methodology is broadly applicable to explore defective low dimensional materials and represents a powerful tool that enables a systematic search of GBs of lateral interfaces for 2D materials.

Zhang, Jianan↗
DOE OSTI2021

Semiclassical shell-structure micro-macroscopic approach for the level density

Level density ρ(E,A) is derived for a one-component nucleon system with a given energy E and particle number A within the mean-field semiclassical periodic-orbit theory beyond the saddle-point method of the Fermi gas model. We obtain ρ∝I ν (S)/S ν , with I ν (S) being the modified Bessel function of the entropy S. Within the micro-macro-canonical approximation (MMA), for a small thermal excitation energy U, with respect to rotational excitations E rot , one obtains ν = 3/2 for ρ(E,A). In the case of excitation energy U larger than E rot but smaller than the neutron separation energy, one finds a larger value of ν = 5/2. A role of the fixed spin variables for rotating nuclei is discussed. The MMA level density ρ reaches the well-known grand-canonical ensemble limit (Fermi gas asymptote) for large S related to large excitation energies, and also reaches the finite micro-canonical limit for small combinatorial entropy S at low excitation energies (the constant “temperature” model). Fitting the ρ(E,A) of the MMA to the experimental data for low excitation energies, taking into account shell and, qualitatively, pairing effects, one obtains for the inverse level density parameter K a value which differs essentially from that parameter derived from data on neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Deep Learning the Electromagnetic Properties of Metamaterials—A Comprehensive Review

Deep neural networks (DNNs) are empirically derived systems that have transformed traditional research methods, and are driving scientific discovery. Here, artificial electromagnetic materials (AEMs)—including electromagnetic metamaterials, photonic crystals, and plasmonics—are research fields where DNN results valorize the data driven approach; especially in cases where conventional methods have failed. In view of the great potential of deep learning for the future of artificial electromagnetic materials research, the status of the field with a focus on recent advances, key limitations, and future directions is reviewed. Strategies, guidance, evaluation, and limits of using deep networks for both forward and inverse AEM problems are presented.

14 SOLAR ENERGY↗
DOE OSTI2021

Bedrock architecture, soil texture, and hyporheic zone characterization combining electrical resistivity and induced polarization imaging

The shallow (<15 m depth) subsurface environment of a short reach along a bedrock stream was investigated with electrical resistivity and induced polarization (IP) to map details of bedrock and soil sediments. The bedrock is mostly comprised of limestone, and is generally resistive compared to the overlying soil. The soil-bedrock interface was determined through a trial and error approach using a sharp boundary feature in the inversion model. The inferred bedrock surface determined from the inversions exhibited undulating patterns with troughs and ridges. A near continuous trough ran alongside the stream within the floodplain, and is suggestive of a paleochannel among other interpretations of this feature. The structure of the electrical resistivity above the bedrock showed small-scale elongated features. The chargeability from the IP method showed larger scale features. High values of chargeability were associated with the sediments in the floodplain, and low values were associated with bedrock, stream, and soil on the elevated banks above the stream. If the chargeability is associated with membrane polarization characteristic of clayey soils, then IP seems to highlight the mere existence of clay, while the resistivity may be more discernable of the relative proportion of clay. Vegetation differences may also explain the chargeability distribution, where parts of the survey with high chargeability had dense pine with no understory making the soils more organically rich.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Quantifying uncertainty in analysis of shockless dynamic compression experiments on platinum. I. Inverse Lagrangian analysis

Absolute measurements of solid-material compressibility by magnetically driven shockless dynamic compression experiments to multi-megabar pressures have the potential to greatly improve the accuracy and precision of pressure calibration standards for use in diamond anvil cell experiments. Here, to this end, we apply characteristics-based inverse Lagrangian analysis (ILA) to 11 sets of ramp-compression data on pure platinum (Pt) metal and then reduce the resulting weighted-mean stress–strain curve to the principal isentrope and room-temperature isotherm using simple models for yield stress and Grüneisen parameter. We introduce several improvements to methods for ILA and quasi-isentrope reduction, the latter including calculation of corrections in wave speed instead of stress and pressure to render results largely independent of initial yield stress while enforcing thermodynamic consistency near zero pressure. More importantly, we quantify in detail the propagation of experimental uncertainty through ILA and model uncertainty through quasi-isentrope reduction, considering all potential sources of error except the electrode and window material models used in ILA. Compared to previous approaches, we find larger uncertainty in longitudinal stress. Monte Carlo analysis demonstrates that uncertainty in the yield-stress model constitutes by far the largest contribution to uncertainty in quasi-isentrope reduction corrections. We present a new room-temperature isotherm for Pt up to 444 GPa, with 1-sigma uncertainty at that pressure of just under ±1.2%; the latter is about a factor of three smaller than uncertainty previously reported for multi-megabar ramp-compression experiments on Pt. The result is well represented by a Vinet-form compression curve with (isothermal) bulk modulus K 0 = 270.3 ± 3.8 GPa, pressure derivative K$^{'}_{0}$= 5.66 ± 0.10, and correlation coefficient $R_{K_{0},K^{'}_{0}}$= –0.843.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2020

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging. Part II

X-ray Bragg coherent diffraction imaging (BCDI) has been demonstrated as a powerful 3D microscopy approach for the investigation of sub-micrometre-scale crystalline particles. The approach is based on the measurement of a series of coherent Bragg diffraction intensity patterns that are numerically inverted to retrieve an image of the spatial distribution of the relative phase and amplitude of the Bragg structure factor of the diffracting sample. This 3D information, which is collected through an angular rotation of the sample, is necessarily obtained in a non-orthogonal frame in Fourier space that must be eventually reconciled. To deal with this, the approach currently favored by practitioners (detailed in Part I) is to perform the entire inversion in conjugate non-orthogonal real- and Fourier-space frames, and to transform the 3D sample image into an orthogonal frame as a post-processing step for result analysis. In this article, which is a direct follow-up of Part I, two different transformation strategies are demonstrated, which enable the entire inversion procedure of the measured data set to be performed in an orthogonal frame. The new approaches described here build mathematical and numerical frameworks that apply to the cases of evenly and non-evenly sampled data along the direction of sample rotation ( i.e. the rocking curve). The value of these methods is that they rely on the experimental geometry, and they incorporate significantly more information about that geometry into the design of the phase-retrieval Fourier transformation than the strategy presented in Part I. Finally, two important outcomes are (1) that the resulting sample image is correctly interpreted in a shear-free frame and (2) physically realistic constraints of BCDI phase retrieval that are difficult to implement with current methods are easily incorporated. Computing scripts are also given to aid readers in the implementation of the proposed formalisms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Design, Analysis, and Experimental Testing of Hydrogen Lean Direct Injection Nozzles at Elevated Pressure

Abstract There are many challenges of commissioning a hydrogen combustor into future gas turbine engines; especially regarding achieving emissions goals. Previously, Escudero et al. and Tran et al. conducted a study to adapt the liquid fuel Lean Direct Injection (LDI) concept from Jet-A to gaseous natural gas-hydrogen blends and pure hydrogen [1], [2]. Experimental data was collected at atmospheric conditions using a Box Behnken design of experiments. The design of experiments suggested that biasing the air split in favor of the inner air circuit and increasing the swirl strength of this inner air passage resulted in improved NOx emissions, while the inverse was true for stability, which was quantified by studying the lean blowoff point (LBO) [1], [2]. The trends revealed by the original experiment [1], [2] provided a design direction for further iterations of the experimental hardware. The study presented herein describes the further investigation of such LDI injectors through experimental methods and computational fluid dynamic (CFD) simulations at atmospheric conditions, which were used to identify potential flow behaviors driving enhanced emissions performance. Further evaluation of select injectors from both studies was then conducted at elevated pressures up to 6 atmospheres. The results from both experiments are presented in this study, which include flame observations, emissions measurements, and operational challenges. NOx emissions results are reported on a volume basis in ppmvd corrected to 15% O2 and corrected for fuel. A predictive model for relating NOx emissions to test conditions at atmospheric conditions show high significance to adiabatic flame temperature while little to no significance to fuel composition for the best performing configurations. The results illustrate the connection between atmospheric testing and testing elevated pressures. The design direction indicated by the initial tests and CFD results in promising configurations for implementation into a Multi-point LDI array.

08 HYDROGEN↗
DOE OSTI2025

Self-supervised physics-informed generative networks for phase retrieval from a single X-ray hologram

X-ray phase contrast imaging significantly improves the visualization of structures with weak or uniform absorption, broadening its applications across a wide range of scientific disciplines. Propagation-based phase contrast is particularly suitable for time- or dose-critical in vivo/in situ/operando (tomography) experiments because it requires only a single intensity measurement. However, the phase information of the wave field is lost during the measurement and must be recovered. Conventional algebraic and iterative methods often rely on specific approximations or boundary conditions that may not be met by many samples or experimental setups. In addition, they require manual tuning of reconstruction parameters by experts, making them less adaptable for complex or variable conditions. Here we present a self-learning approach for solving the inverse problem of phase retrieval in the near-field regime of Fresnel theory using a single intensity measurement (hologram). A physics-informed generative adversarial network is employed to reconstruct both the phase and absorbance of the unpropagated wave field in the sample plane from a single hologram. Unlike most state-of-the-art deep learning approaches for phase retrieval, our approach does not require paired, unpaired, or simulated training data. This significantly broadens the applicability of our approach, as acquiring or generating suitable training data remains a major challenge due to the wide variability in sample types and experimental configurations. The algorithm demonstrates robust and consistent performance across diverse imaging conditions and sample types, delivering quantitative, high-quality reconstructions for both simulated data and experimental datasets acquired at beamline P05 at PETRA III (DESY, Hamburg), operated by Helmholtz-Zentrum Hereon. Furthermore, it enables the simultaneous retrieval of both phase and absorption information.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗