Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion constraints”

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

Physics-informed Karhunen-Loeve and Neural Network Approximations for Solving Inverse Differential Equation Problems

Here we present the PI-CKL-NN method for parameter estimation in differential equation (DE) models given sparse measurements of the parameters and states. In the proposed approach, the space- or time-dependent parameters are approximated by Karhunen-Loeve (KL) expansions that are conditioned on the parameters’ measurements, and the states are approximated by deep neural networks (DNNs). The unknown weights in the KL expansions and DNNs are found my minimizing the cost function that enforces the measurements of the states the DE constraint. Regularization is achieved by adding the l2 norm of the conditional KL coefficients into the loss function. Our approach assumes that the parameter fields are correlated in space or time and enforces the statistical knowledge (the mean and the covariance function) in addition to the DE constraints and measurements as opposed to the physics-informed neural network (PINN) and other similar physics-informed machine learning methods where only DE constraints and data are used for parameter estimation. We use the PI-CKL-NN method for parameter estimation in an ordinary differential equation with an unknown time-dependent parameter and the one- and two-dimensional partial differential diffusion equations with unknown space-dependent diffusion coefficients. We also demonstrate that PI-CKL-NN is more accurate than the PINN method, especially when the observations of the parameters are very sparse

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Understanding protein-complex assembly through grand canonical maximum entropy modeling

Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded systems where the abundance of these proteins vary with cell types. While some protein complexes form putative structures that can be visualized with imaging, there are far more protein complexes that are yet to be solved because of their dynamic associations with one another. Nevertheless, it is possible to infer these protein complexes through a physical model. However, it is often not clear to physicists what kind of data from biology is necessary for such a modeling endeavor. Here, we aim to model these clusters of coarse-grained protein assemblies from multiple subunits through the constraints of interactions among the subunits and the chemical potential of each subunit. We obtained the constraints on the interactions among subunits from the known protein structures. We inferred the chemical potential that dictates the particle number distribution of each protein subunit from the knowledge of protein abundance from experimental data. Guided by the maximum entropy principle, we formulated an inverse statistical mechanical method to infer the distribution of particle numbers from the data of protein abundance as chemical potentials for a grand canonical multicomponent mixture. Using grand canonical Monte Carlo simulations, we captured a distribution of high-order clusters in a protein complex of succinate dehydrogenase with four known subunits. The complexity of hierarchical clusters varies with the relative protein abundance of each subunit in distinctive cell types such as lung, heart, and brain. When the crowding content increases, we observed that crowding stabilizes emergent clusters that do not exist in dilute conditions. We, therefore, proposed a testable hypothesis that the hierarchical complexity of protein clusters on a molecular scale is a plausible biomarker of predicting the phenotypes of a cell.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

INTEGRATE - Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements

The INTEGRATE (Inverse Network Transformations for Efficient Generation of Robust Airfoil and Turbine Enhancements) project is developing a new inverse-design capability for the aerodynamic design of wind turbine rotors using invertible neural networks. This AI-based design technology can capture complex non-linear aerodynamic effects while being 100 times faster than design approaches based on computational fluid dynamics. This project enables innovation in wind turbine design by accelerating time to market through higher-accuracy early design iterations to reduce the levelized cost of energy. INVERTIBLE NEURAL NETWORKS Researchers are leveraging a specialized invertible neural network (INN) architecture along with the novel dimension-reduction methods and airfoil/blade shape representations developed by collaborators at the National Institute of Standards and Technology (NIST) learns complex relationships between airfoil or blade shapes and their associated aerodynamic and structural properties. This INN architecture will accelerate designs by providing a cost-effective alternative to current industrial aerodynamic design processes, including: - Blade element momentum (BEM) theory models: limited effectiveness for design of offshore rotors with large, flexible blades where nonlinear aerodynamic effects dominate - Direct design using computational fluid dynamics (CFD): cost-prohibitive - Inverse-design models based on deep neural networks (DNNs): attractive alternative to CFD for 2D design problems, but quickly overwhelmed by the increased number of design variables in 3D problems AUTOMATED COMPUTATIONAL FLUID DYNAMICS FOR TRAINING DATA GENERATION - MERCURY FRAMEWORK The INN is trained on data obtained using the University of Marylands (UMD) Mercury Framework, which has with robust automated mesh generation capabilities and advanced turbulence and transition models validated for wind energy applications. Mercury is a multi-mesh paradigm, heterogeneous CPU-GPU framework. The framework incorporates three flow solvers at UMD, 1) OverTURNS, a structured solver on CPUs, 2) HAMSTR, a line based unstructured solver on CPUs, and 3) GARFIELD, a structured solver on GPUs. The framework is based on Python, that is often used to wrap C or Fortran codes for interoperability with other solvers. Communication between multiple solvers is accomplished with a Topology Independent Overset Grid Assembler (TIOGA). NOVEL AIRFOIL SHAPE REPRESENTATIONS USING GRASSMAN SPACES We developed a novel representation of shapes which decouples affine-style deformations from a rich set of data-driven deformations over a submanifold of the Grassmannian. The Grassmannian representation as an analytic generative model, informed by a database of physically relevant airfoils, offers (i) a rich set of novel 2D airfoil deformations not previously captured in the data , (ii) improved low-dimensional parameter domain for inferential statistics informing design/manufacturing, and (iii) consistent 3D blade representation and perturbation over a sequence of nominal shapes. TECHNOLOGY TRANSFER DEMONSTRATION - COUPLING WITH NREL WISDEM Researchers have integrated the inverse-design tool for 2D airfoils (INN-Airfoil) into WISDEM (Wind Plant Integrated Systems Design and Engineering Model), a multidisciplinary design and optimization framework for assessing the cost of energy, as part of tech-transfer demonstration. The integration of INN-Airfoil into WISDEM allows for the design of airfoils along with the blades that meet the dynamic design constraints on cost of energy, annual energy production, and the capital costs. Through preliminary studies, researchers have shown that the coupled INN-Airfoil + WISDEM approach reduces the cost of energy by around 1% compared to the conventional design approach. This page will serve as a place to easily access all the publications from this work and the repositories for the software developed and released through this pr...

aerodynamics↗
DOE OSTI2022

Stochastic Learning Approach for Binary Optimization: Application to Bayesian Optimal Design of Experiments

Here, we present a novel stochastic approach to binary optimization suited for optimal experimental design (OED) for Bayesian inverse problems governed by mathematical models such as partial differential equations. The OED utility function, namely, the regularized optimality criterion, is cast into a stochastic objective function in the form of an expectation over a multivariate Bernoulli distribution. The probabilistic objective is then solved by using a stochastic optimization routine to find an optimal observational policy. This formulation (a) is generally applicable to binary optimization problems with soft constraints and is ideal for OED and sensor placement problems; (b) does not require differentiability of the original objective function (e.g., a utility function in OED applications) with respect to the design variable, and thus it enables direct employment of sparsity-enforcing penalty functions such as $\ell_0$, without needing to utilize a continuation procedure or apply a rounding technique; (c) exhibits much lower computational cost than traditional gradient-based relaxation approaches; and (d) can be applied to both linear and nonlinear OED problems with proper choice of the utility function. The proposed approach is analyzed from an optimization perspective with detailed convergence analysis of the optimization approach and is also analyzed from a machine learning perspective with correspondence to policy gradient reinforcement learning. The approach is demonstrated numerically by using an idealized two-dimensional Bayesian linear inverse problem and validated by extensive numerical experiments carried out for sensor placement in a parameter identification setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nuclear level density and γ -ray strength function of Ni 67 and the impact on the i process

Proton- γ coincidences from ( d , p ) reactions between a Ni 66 beam and a deuterated polyethylene target have been analyzed with the inverse-Oslo method to find the nuclear level density (NLD) and γ -ray strength function ( γ SF ) of Ni 67 . The Ni 66 ( n , γ ) capture cross section has been calculated using the Hauser-Feshbach model in using the measured NLD and γ SF as constraints. The results confirm that the Ni 66 ( n , γ ) reaction acts as a bottleneck when relying on one-zone nucleosynthesis calculations. However, the impact of this reaction is strongly dampened in multizone models of low-metallicity AGB stars experiencing i -process nucleosynthesis. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The significant impact of aerosol vertical structure on lower atmosphere stability and its critical role in aerosol–planetary boundary layer (PBL) interactions

The aerosol–planetary boundary layer (PBL) interaction was proposed as an important mechanism to stabilize the atmosphere and exacerbate surface air pollution. Despite the tremendous progress made in understanding this process, its magnitude and significance still have large uncertainties and vary largely with aerosol distribution and meteorological conditions. In this study, we focus on the role of aerosol vertical distribution in thermodynamic stability and PBL development by jointly using micropulse lidar, sun photometer, and radiosonde measurements taken in Beijing. Despite the complexity of aerosol vertical distributions, cloud-free aerosol structures can be largely classified into three types: well-mixed, decreasing with height, and inverse structures. The aerosol–PBL relationship and diurnal cycles of the PBL height and PM 2.5 associated with these different aerosol vertical structures show distinct characteristics. The vertical distribution of aerosol radiative forcing differs drastically among the three types, with strong heating in the lower, middle, and upper PBL, respectively. Such a discrepancy in the heating rate affects the atmospheric buoyancy and stability differently in the three distinct aerosol structures. Absorbing aerosols have a weaker effect of stabilizing the lower atmosphere under the decreasing structure than under the inverse structure. As a result, the aerosol–PBL interaction can be strengthened by the inverse aerosol structure and can be potentially neutralized by the decreasing structure. Moreover, aerosols can both enhance and suppress PBL stability, leading to both positive and negative feedback loops. This study attempts to improve our understanding of the aerosol–PBL interaction, showing the importance of the observational constraint of aerosol vertical distribution for simulating this interaction and consequent feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Gaussian processes for inferring parton distributions

The extraction of parton distribution functions (PDFs) from experimental or lattice QCD data is an ill-posed inverse problem, where regularization strongly impacts both systematic uncertainties and the reliability of the results. We study a framework based on Gaussian Process Regression (GPR) to reconstruct PDFs from lattice QCD matrix elements. Within a Bayesian framework, Gaussian processes serve as flexible priors that encode uncertainties, correlations, and constraints without imposing rigid functional forms. We investigate a wide range of kernel choices, mean functions, and hyperparameter treatments. We quantify information gained from the data using the Kullback-Leibler divergence. Synthetic data tests demonstrate the consistency and robustness of the method. Our study establishes GPR as a systematic and non-parametric approach to PDF reconstruction, offering controlled uncertainty estimates and reduced model bias in lattice QCD analyses.

hadronic spectroscopy↗
DOE OSTI2021

Jet fragmentation transverse momentum distributions in pp and p-Pb collisions at $ \sqrt{s}$, $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV

Jet fragmentation transverse momentum ( j T ) distributions are measured in proton-proton (pp) and proton-lead (p-Pb) collisions at s NN = 5 . 02 TeV with the ALICE experiment at the LHC. Jets are reconstructed with the ALICE tracking detectors and electromagnetic calorimeter using the anti- k T algorithm with resolution parameter R = 0 . 4 in the pseudorapidity range |η| < 0 . 25. The j T values are calculated for charged particles inside a fixed cone with a radius R = 0 . 4 around the reconstructed jet axis. The measured j T distributions are compared with a variety of parton-shower models. Herwig and P ythia 8 based models describe the data well for the higher j T region, while they underestimate the lower j T region. The j T distributions are further characterised by fitting them with a function composed of an inverse gamma function for higher j T values (called the “wide component”), related to the perturbative component of the fragmentation process, and with a Gaussian for lower j T values (called the “narrow component”), predominantly connected to the hadronisation process. The width of the Gaussian has only a weak dependence on jet transverse momentum, while that of the inverse gamma function increases with increasing jet transverse momentum. For the narrow component, the measured trends are successfully described by all models except for Herwig. For the wide component, Herwig and PYTHIA 8 based models slightly underestimate the data for the higher jet transverse momentum region. These measurements set constraints on models of jet fragmentation and hadronisation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Development of a coarse-grained molecular dynamics model for poly(dimethyl- co -diphenyl)siloxane

Polydimethylsiloxane is an important polymeric material with a wide range of applications. However, environmental effects like low temperature can induce crystallization in this material with resulting changes in its structural and dynamic properties. The incorporation of phenyl-siloxane components, e.g., as in a poly(dimethyl-co-diphenyl)siloxane random copolymer, is known to suppress such crystallization. Molecular dynamics (MD) simulations can be a powerful tool to understand such effects in atomistic detail. Unfortunately, all-atomistic molecular dynamics (AAMD) is limited in both spatial dimensions and simulation times it can probe. Here, to overcome such constraints and to extend to more useful length- and time-scales, we systematically develop a coarse-grained molecular dynamics (CGMD) model for the poly(dimethyl-co-diphenyl)siloxane system with bonded and non-bonded interactions determined from all-atomistic simulations by the iterative Boltzmann inversion (IBI) method. Additionally, we propose a lever rule that can be useful to generate non-bonded potentials for such systems without reference to the all-atomistic ground truth. Our model captures the structural and dynamic properties of the copolymer material with quantitative accuracy and is useful to study long-time dynamics of highly-entangled systems, sequence-dependent properties, phase behaviour, etc.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Search for the 1/2 + intruder state in 35 P

The excitation energy of deformed intruder states (specifically the 2p2h bandhead) as a function of proton number Z along N = 20 is of interest both in terms of better understanding the evolution of nuclear structure between spherical Ca 40 and the Island of Inversion nuclei, and for benchmarking theoretical descriptions in this region. At the center of the N = 20 Island of Inversion, the n p n h (where n = 2 , 4 , 6 ) neutron excitations across a diminished N = 20 gap result in deformed and collective ground states, as observed in Mg 32 . In heavier isotones, n p n h excitations do not dominate in the ground states but are present in the relatively low-lying level schemes. With the aim of identifying the expected 2 p 2 h x s 1 / 2 + state in P 35 , the only N = 20 isotone for which the neutron 2p2h excitation bandhead has not yet been identified, the S 36 ( d , He 3 ) P 35 reaction has been revisited in inverse kinematics with the HELical Orbit Spectrometer (HELIOS) at the Argonne Tandem Linac Accelerator System (ATLAS). While a candidate state has not been located, an upper limit for the transfer reaction cross section to populate such a configuration within a 2.5 to 3.6 MeV energy range provides a stringent constraint on the wave function compositions in both S 36 and P 35 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

CLASSY III. The Properties of Starburst-driven Warm Ionized Outflows

We report the results of analyses of galactic outflows in a sample of 45 low-redshift starburst galaxies in the COS Legacy Archive Spectroscopic SurveY (CLASSY), augmented by five additional similar starbursts with Cosmic Origins Spectrograph (COS) data. The outflows are traced by blueshifted absorption lines of metals spanning a wide range of ionization potential. The high quality and broad spectral coverage of CLASSY data enable us to disentangle the absorption due to the static interstellar medium (ISM) from that due to outflows. We further use different line multiplets and doublets to determine the covering fraction, column density, and ionization state as a function of velocity for each outflow. We measure the outflow's mean velocity and velocity width, and find that both correlate in a highly significant way with the star formation rate, galaxy mass, and circular velocity over ranges of four orders of magnitude for the first two properties. We also estimate outflow rates of metals, mass, momentum, and kinetic energy. We find that, at most, only about 20% of silicon created and ejected by supernovae in the starburst is carried out in the warm phase we observe. The outflows' mass-loading factor increases steeply and inversely with both circular and outflow velocity (log–log slope ~–1.6), and reaches ~10 for dwarf galaxies. We find that the outflows typically carry about 10%–100% of the momentum injected by massive stars and about 1%–20% of the kinetic energy. We show that these results place interesting constraints on, and new insights into, models and simulations of galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Application of physics-informed neural networks (PINNs) solution to coupled thermal and hydraulic processes in silty sands

Abstract The accurate modeling of water and heat transport in soils is crucial for both geo-environmental and geothermal engineering. Traditional modeling methods are problematic because they require well-defined boundaries and initial conditions. Recently, physics-informed neural networks (PINNs), which incorporate partial differential equations (PDEs) to solve forward and inverse problems, have attracted increasing attention in machine learning research. In this study, we applied PINNs to tackle hydraulic and thermal transport coupling forward problems in silty sands. A fully connected deep neural network was utilized for training. This neural network model leverages automatic differentiation to apply the governing equations as constraints, based on the mathematical approximations established by the neural network itself. We conducted forward problems and compared the solutions derived from PINNs with those from Finite Element Method (FEM) simulations. The forward problem results demonstrate the PINNs model’s capability in predicting hydraulic transport, heat transport, and thermal–hydraulic coupling in silty sands under various boundary conditions. The PINNs exhibited great performance in simulating the thermal–hydraulic coupling problem. The accuracy of the PINNs solutions shows its potential for simulation in geotechnical engineering.

Feng, Yuan↗
DOE OSTI2026

Input specific neural networks

Neural networks have emerged as powerful tools for mapping between inputs and outputs. However, their black-box nature limits the ability to encode or impose specific structural relationships between inputs and outputs. Many scientific and engineering problems, such as constitutive modeling in solid mechanics, require networks that can enforce convexity, monotonicity, or other structural constraints to ensure physical consistency. Here, we introduce the Input Specific Neural Network (ISNN), a new architecture that enables multiple, distinct constraints to be imposed on different input subsets for scalar-valued outputs. This framework unifies convex, monotone–convex, monotone, and arbitrary mappings within a single network for the first time. Two ISNN architectures with analytical first- and second-order derivatives are developed. We demonstrate the performance on synthetic toy problems, inverse problems in isotropic hyperelasticity, and finite element simulations. ISNNs achieve improved extrapolation behavior, require fewer invariant inputs than standard input convex networks for polyconvex potentials, and enable significant computational savings via manual differentiation. We also show how ISNNs can be used to learn structural relationships between inputs and outputs via a binary gating mechanism. Particularly, ISNNs are employed to model a homogenized anisotropic free energy potential in a decoupled multiscale setting. The network learns whether or not the potential should be modeled as polyconvex and retains only the relevant layers while using the minimum number of inputs. ISNNs provide a flexible foundation for embedding structural priors into neural networks, enhancing both interpretability and stability. They are broadly applicable across computational mechanics and other scientific domains requiring constrained functional relationships.

Jadoon, Asghar A. [Univ. of Texas, Austin, TX (Uni↗
DOE OSTI2024

Quantum limits of superconducting-photonic links and their extension to millimeter waves

Photonic addressing of superconducting circuits has been proposed to overcome wiring complexity and heat-load challenges. However, such superconducting-photonic links suffer from an efficiency-noise trade-off that limits scalability. This trade-off arises because increasing power conversion efficiency entails reducing optical power, which makes the converted signal susceptible to shot noise. We analyze this trade-off and find that the laser-driven qubit gate infidelity scales inversely with the number of photons used. While methods such as nonlinear detection or squeezed light could mitigate this effect, we consider generating higher-frequency electrical signals, such as millimeter waves (100 GHz), using laser light. At these higher frequencies, circuits have higher operating temperatures and cooling power budgets and alleviate constraints posed by the trade-off. Therefore, we demonstrate an optically driven cryogenic millimeter-wave source with a maximum power efficiency of 8 ×10 −5 that can generate a maximum of 0.7 μ⁢W of 80 GHz power, with a 1200-thermal-photon equivalent of added noise at 4K. Using this source, we perform frequency-domain spectroscopy of superconducting NbTiN resonators at 80–90 GHz. Our results show a promising approach to lessen the efficiency-noise constraints on superconducting-photonic links, while leveraging the benefits of photonic signal delivery. Further optimization of power efficiency and noise at high frequencies could make scalable photonic control of superconducting qubits viable at temperatures exceeding 1K.

Optoelectronics↗
DOE OSTI2024

β -delayed neutron emission of Mn 64 , Cr 62 , and Fe 65

The β-decay properties of nuclei near the second nuclear “island of inversion” around neutron rich nuclei with neutron number 40 are important tests of nuclear structure models and interactions. In particular, the β-delayed neutron emission branch (Ρ n ), is useful for investigating β-strength and neutron-γ competition above the neutron separation energies of the daughter nuclei. We report new constraints for Ρ n values for three nuclei in the region: 62 Cr (Ρ n <1%), 64 Mn (Ρ n = 1.5⁢(6)%), and 65 Fe (Ρ n < 1%), measured with the Neutron Emission Ratio Observer (NERO) neutron long counter system and the Beta Counting Station (BCS) at the National Superconducting Cyclotron Laboratory (NSCL). Our results resolve the large discrepancy between previous direct and indirect measurements for 64 Mn and confirm the predictions of global theoretical models when a statistical treatment of the γ and neutron decays of the daughter states is included. Here we also obtain improved half-lives for 62 Cr [206(5) ms] and the short-lived isomer in the 62 Fe daughter [112(7) ms] from β-delayed γ emission data obtained in the same experiment with the Summing NaI (SuN) total absorption spectrometer. Finally, we use γ emission data to obtain a new upper limit for the 62 Cr β-decay population of the long-lived isomeric state in 62 Mn.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 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↗