SEARCH · Engineering Papers

Results for “inverse modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗

DOE OSTI2021

Coexisting normal and intruder configurations in 32 Mg

Situated in the so-called “island of inversion,” the nucleus 32 Mg is considered as an archetypal example of the disappearance of magicity at N = 20. We report on high statistics in-beam spectroscopy of 32 Mg with a unique approach, in that two direct reaction probes with different sensitivities to the underlying nuclear structure are employed at the same time. More specifically, states in 32 Mg were populated by knockout reactions starting from 33 Mg and 34 Si, lying inside and outside the island of inversion, respectively. The momentum distributions of the reaction residues and the cross sections leading to the individual final states were confronted with eikonal-based reaction calculations, yielding a significantly updated level scheme for 32 Mg and spin-parity assignments. By fully exploiting observables obtained in this measurement, a variety of structures coexisting in 32 Mg was unraveled. Comparisons with theoretical predictions based on shell-model overlaps allowed for clear discrimination between different structural models, revealing that the complete theoretical description of this key nucleus is yet to be achieved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗

DOE OSTI2026

Synthetic modeling of soft x-ray emissivity for magnetic island analysis in LTX-β

We present a synthetic soft x-ray (SXR) forward-modeling framework that characterizes the emissivity structure of rotating magnetic islands in the Lithium Tokamak eXperiment-β across space and time. Magnetic islands associated with tearing modes produce modulations in line-integrated SXR brightness. Traditional tomographic methods struggle to resolve these structures in devices with limited sightlines, resulting in an under-determined inversion problem, or otherwise necessitate equilibrium reconstruction data, preventing their use in active control. Our approach removes this challenge by combining a fully three-dimensional ray-tracing model with a time-dependent emissivity prescription of a $m/n$ = 2/1 magnetic island. The model is validated through helical island geometry, rotation measured by magnetic diagnostics, and equilibrium constraints from PSI-Tri reconstructions. Synthetic brightness signals are generated for each photodiode sightline from a modeled emissivity profile and directly compared with experimental data from a tangential midplane SXR array. By fitting the synthetic diagnostic output to the observed brightness across time, we simultaneously infer key geometric parameters-such as island radial location, island width, and rotation frequency-without requiring full tomographic reconstruction, all with an average absolute deviation of less than 3%. This work demonstrates that forward modeling with a single tangential array can extract key island parameters in a spherical tokamak, providing a computationally efficient alternative to conventional SXR tomography and providing a pathway toward real-time magnetic-island characterization in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Precision bootstrap for the $\mathcal{N}$ = 1 super-Ising model

In this note we report an improved determination of the scaling dimensions and OPE coefficients of the minimal supersymmetric extension of the 3d Ising model using the conformal bootstrap. We also show how this data can be used as input to the Lorentzian inversion formula, finding good agreement between analytic calculations and numerical extremal spectra once mixing effects are resolved.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Crustal Structure of the Mesopotamian Plain, East of Iraq

The crustal structure of Iraq was investigated through analyzing teleseismic data from 12 new seismic stations. Three seismic stations are located within the Zagros Fold-Thrust Belt, and nine are within the Mesopotamian Plain. Joint inversion of P wave receiver function and Rayleigh wave dispersion data were employed to resolve S wave velocity structure models beneath each station. Combining these models with available Moho depths within and around the study area reveals that the Moho depth is smoothly increasing from the Arabian platform toward the Zagros Mountains. Further, an exceptional deeper root is observed near the eastern edge of the Mesopotamian Plain where sedimentary pile is the thickest. This root was interpreted as a structure inherited from the successive Mesozoic rifting of the NE Arabian platform enhanced by sedimentary loading and progressive subsidence. A low-velocity uppermost mantle was observed beneath the Arabian Foredeep attesting the existence of a low-strength lithospheric mantle beneath the region southwest from the Zagros deformation front. The weak rheology of the uppermost mantle may have allowed local sinking of the crust and deepening of the Moho boundary due to vertical loading and the Late Miocene lateral contraction.

58 GEOSCIENCES↗

DOE OSTI2022

Excited quantum anomalous and spin Hall effect: dissociation of flat-bands-enabled excitonic insulator state

Quantum anomalous Hall effect (QAHE) and quantum spin Hall effect (QSHE) are two interesting physical manifestations of 2D materials that have an intrinsic nontrivial band topology. In principle, they are ground-state equilibrium properties characterized by Fermi level lying in a topological gap, below which all the occupied bands are summed to a non-zero topological invariant. In this work, we propose theoretical concepts and models of ‘excited’ QAHE (EQAHE) and EQSHE generated by dissociation of an excitonic insulator (EI) state with complete population inversion (CPI), a unique many-body ground state enabled by two yin-yang flat bands (FBs) of opposite chirality hosted in a diatomic Kagome lattice. The two FBs have a trivial gap in between, i.e. the system is a trivial insulator in the single-particle ground-state, but nontrivial gaps above and below, so that upon photoexcitation the quasi-Fermi levels of both electrons and holes will lie in a nontrivial gap achieved by the CPI-EI state, as demonstrated by exact diagonalization calculations. Then dissociation of singlet and triplet EI state will lead to EQAHE and EQSHE, respectively. Realizations of yin-yang FBs in real materials are also discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2022

IGR J18249−3243: a new GeV-emitting FR II and the emerging population of high-energy radio galaxies

ABSTRACT The advent of new all-sky radio surveys such as the VLA Sky Survey and the Rapid ASKAP Continuum Survey, performed with the latest generation radio telescopes, is opening new possibilities on the classification and study of extragalactic γ-ray sources, specially the underrepresented ones like radio galaxies. In particular, the enhanced sensitivity (sub-mJy level) and resolution (a few arcsec) provides a better morphological and spectral classification. In this work, we present the reclassification of a Fermi/Large Area Telescope (LAT) source as a new Fanaroff–Riley II radio galaxy from the International Gamma-Ray Astrophysics Laboratory sample found to emit at GeV energies. Through a broad-band spectral fitting from radio to γ-ray, we find that the commonly invoked jet contribution is not sufficient to account for the observed γ-ray flux. Our modelling suggests that the observed emission could mainly originate in the lobes (rather than in the radio core) by inverse Compton scattering of radio-emitting electrons off the ambient photon fields. In addition, we cross-correlated the latest generation radio surveys with a list of Fermi/LAT candidate misaligned AGN from the literature, finding four new radio galaxies with a double-lobed morphology. Additional four objects could be classified as such thanks to previous studies in the literature, for a total of nine new radio galaxies with GeV emission presented in this work. We foresee that further objects of this class might be found in the near future with the advent of the Square Kilometer Array, populating the GeV sky.

Bruni, G. (ORCID:0000000251826289)↗

DOE OSTI2022

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Coherent High-Frequency Axial Oscillations in a Partially Magnetized Direct Current Magnetron Discharge

High-frequency oscillations are observed in a neon plasma of a direct current magnetron discharge. At low discharge currents, we see highly coherent 60 MHz fluctuations. Above a distinct current threshold, secondary 5–10 MHz fluctuations emerge in addition to turbulent fluctuations in the 60–100 MHz range. The oscillations in the total discharge current suggest axial wave propagation. A lower-hybrid wave theory is invoked to model the high-frequency oscillations. Here, we attribute the low-frequency modes to a turbulence-driven inverse cascade process, as suggested by recent simulations.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

The Artificial Scientist: in-Transit Machine Learning of Plasma Simulations

Large-scale simulations or scientific experiments produce petabytes of data per run. This poses massive challenges for I/O and storage when scientific analysis workflows are run manually offline. Unsupervised deep learning-based techniques to extract patterns and non-linear relations from these large amounts of data provide a way to build scientific understanding from raw data, reducing the need for manual pre-selection of analysis steps, but require exascale compute and memory to process the full dataset available. In this paper, we demonstrate a heterogeneous streaming workflow in which plasma simulation data is streamed directly to a Machine Learning (ML) application training a model on the simulation data in-transit, completely circumventing the capacity-constrained filesystem bottleneck. This workflow employs openPMD to provide a high level interface to describe scientific data and also uses ADIOS2, to transfer volumes of data that exceed the capabilities of the filesystem. We employ experience replay to avoid catastrophic forgetting in learning from this non-steady state process in a continual manner and adapt it to improve model convergence while learning in-transit. As a proof-of-concept, we approach the ill-posed inverse problem of predicting particle dynamics from radiation in a particle-incell (PIConGPU) simulation of the Kelvin-Helmholtz instability (KHI). We detail hardware-software co-design challenges as we scale PIConGPU to full Frontier, the Top-1 system as of June 2024 Top500 list.

Kelling, Jeffrey [Helmholtz-Zentrum Dresden Rossen↗

DOE OSTI2021

Universal dwell time optimization for deterministic optics fabrication

Computer-Controlled Optical Surfacing (CCOS) has been greatly developed and widely used for precision optical fabrication in the past three decades. It relies on robust dwell time solutions to determine how long the polishing tools must dwell at certain points over the surfaces to achieve the expected forms. However, as dwell time calculations are modeled as ill-posed deconvolution, it is always non-trivial to reach a reliable solution that 1) is non-negative, since CCOS systems are not capable of adding materials, 2) minimizes the residual in the clear aperture 3) minimizes the total dwell time to guarantee the stability and efficiency of CCOS processes, 4) can be flexibly adapted to different tool paths, 5) the parameter tuning of the algorithm is simple, and 6) the computational cost is reasonable. In this study, we propose a novel Universal Dwell time Optimization (UDO) model that universally satisfies these criteria. First, the matrix-based discretization of the convolutional polishing model is employed so that dwell time can be flexibly calculated for arbitrary dwell points. Second, UDO simplifies the inverse deconvolution as a forward scalar optimization for the first time, which drastically increases the solution stability and the computational efficiency. Finally, the dwell time solution is improved by a robust iterative refinement and a total dwell time reduction scheme. The superiority and general applicability of the proposed algorithm are verified on the simulations of different CCOS processes. A real application of UDO in improving a synchrotron X-ray mirror using Ion Beam Figuring (IBF) is then demonstrated. The simulation indicates that the estimated residual in the 92.3 mm × 15.7 mm CA can be reduced from 6.32 nm Root Mean Square (RMS) to 0.20 nm RMS in 3.37 min. After one IBF process, the measured residual in the CA converges to 0.19 nm RMS, which coincides with the simulation.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Turbulent Model of Crab Nebula Radiation

We construct a turbulent model of the Crab Nebula's nonthermal emission. The present model resolves a number of long-standing problems of the Kennel–Coroniti model: (i) the sigma problem, (ii) the hard spectrum of radio electrons, (iii) the high peak energy of gamma-ray flares, (iv) and the spatial evolution of the infrared (IR) emission. The Nebula contains two populations of injected particles: Component-I, accelerated at the wind termination shock via the Fermi-I mechanism; and Component-II, accelerated in reconnecting turbulence in highly magnetized (σ >> 1) plasma in the central part of the Crab Nebula. The reconnecting turbulence in Component-II extends from radio to gamma-rays: it accelerates radio electrons with a hard spectrum, destroys the large-scale magnetic flux (and thus resolves the sigma problem), and occasionally produces gamma-ray flares (from the largest-scale reconnection events). The model reproduces the broadband spectrum of the Crab Nebula, from low-frequency synchrotron emission in radio to inverse-Compton emission at TeV energies, as well as the spatially resolved evolution of the spectral indices in the IR and optical bands.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Highly-scalable, physics-informed GANs for learning solutions of stochastic PDEs

Uncertainty quantification for forward and inverse problems is a central challenge across physical and biomedical disciplines. We address this challenge for the problem of modeling subsurface flow at the Hanford Site by combining stochastic computational models with observational data using physics-informed GAN models. The geographic extent, spatial heterogeneity, and multiple correlation length scales of the Hanford Site require training a computationally intensive GAN model to thousands of dimensions. We develop a highly optimized implementation that scales to {\herogpucount} NVIDIA Volta GPUs. We develop a hierarchical scheme based on a multi-player game-theoretic approach for exploiting domain parallelism, map discriminators and generators to multiple GPUs, and employ efficient communication schemes to ensure training stability and convergence. Our implementation scales to {\heronodecount} nodes on the Summit supercomputer with a \weakscalingeff\% scaling efficiency, achieving peak and sustained half-precision rates of {\heropeakperf} PF/s and {\herosustainedperf} PF/s.

physics informed ML, machine learning, generative ↗

DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗

DOE OSTI2026

Verification of an energy-conserving semi-implicit electrostatic particle-in-cell scheme for modeling high-density plasma at scale

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Cyclotron resonance↗

DOE OSTI2023

Field evaluation of semi‐automated moisture estimation from geophysics using machine learning

Geophysical methods can provide three-dimensional (3D), spatially continuous estimates of soil moisture. However, point-to-point comparisons of geophysical properties to measure soil moisture data are frequently unsatisfactory, resulting in geophysics being used for qualitative purposes only. This is because (1) geophysics requires models that relate geophysical signals to soil moisture, (2) geophysical methods have potential uncertainties resulting from smoothing and artifacts introduced from processing and inversion, and (3) results from multiple geophysical methods are not easily combined within a single soil moisture estimation framework. To investigate these potential limitations, an irrigation experiment was performed wherein soil moisture was monitored through time, and several surface geophysical datasets indirectly sensitive to soil moisture were collected before and after irrigation: ground penetrating radar, electrical resistivity tomography (ERT), and frequency domain electromagnetics (FDEM). Data were exported in both raw and processed form, and then snapped to a common 3D grid to facilitate moisture prediction by standard calibration techniques, multivariate regression, and machine learning. A combination of inverted ERT data, raw FDEM, and inverted FDEM data was most informative for predicting soil moisture using a random regression forest model (one-thousand 60/40 training/test cross-validation folds produced root mean squared errors ranging from 0.025–0.046 cm 3 /cm 3 ). This cross-validated model was further supported by a separate evaluation using a test set from a physically separate portion of the study area. Machine learning was conducive to a semi-automated model-selection process that could be used for other sites and datasets to locally improve accuracy.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A Regional Phase Amplitude Model of 2-D Attenuation for North America

We analyzed seismic attenuation patterns across the North America using over 70,000,000 Lg wave amplitudes recorded by the various networks at frequencies of 0.1-32 Hz. Our inversion solved for laterally varying attenuation, site terms, moments, and apparent stress following Phillips et al., (2016). The inversion was anchored by independently constrained: corner frequencies (via coda spectral ratios) to control the attenuation-stress tradeoff and moment measurements of teleseismic (GCMT, USGS) and regional (St. Louis University and UC Berkeley) earthquakes to provided absolute scaling. The quality factor (Q) shows clear regional patterns: low values in coastal, volcanic, and tectonically active regions, and high values in stable areas like the Great Plains and major plateaus throughout North America. These 2-D Q models enable improved regional source characterization, magnitude estimation, and yield determination.

58 GEOSCIENCES↗