Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constitutive models”

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

Planar, longitudinal, compressive waves in solids: Thermodynamics and uniaxial strain restrictions

A new tri-axial pressure-based constitutive expression has been found using Cauchy's stress tensor. This stress state emphasizes pressure and shear stress. The description is a pressure plus an effective shear stress allowing for a constitutive law based on atomic solid-state phase changes in crystalline cells due to pressure plus shear-based dislocation motion commonly associated with plasticity. Pressure has a new role in the material's constitutive response as it is separated from plasticity. The thermo-mechanical system describes third-order Gibbs’ expressions without specific volume restrictions placed upon the material. Isothermally, the ratio of heat to shear work in elastic copper is shown to approach zero at a very low temperature and become larger than one as temperature approaches melting. Wave compression models investigated are elastic and plastic: in fully elastic materials, the planar wave is restricted by Poisson's effect although plastic shear changes this constraint. Plastic deformation, dominated by dissipative shear stresses in uniaxial strain, heats the material while excluding phase changes from hydrostatic pressures. The material properties per se across Hugoniot shocks are described with entropy concepts. Shock waves are exceedingly complex since the constitutive laws are linked at extreme temperatures, pressures, and shear stresses. Isothermal, isentropic, isochoric, and iso-shear conditions are used throughout with Jacobian algebra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The Impact of Wettability on the Co-moving Velocity of Two-Fluid Flow in Porous Media

Abstract The impact of wettability on the co-moving velocity of two-fluid flow in porous media is analyzed herein. The co-moving velocity, developed by Roy et al. (Front Phys 8:4, 2022), is a novel representation of the flow behavior of two fluids through porous media. Our study aims to better understand the behavior of the co-moving velocity by analyzing simulation data under various wetting conditions. We analyzed 46 relative permeability curves based on the Lattice–Boltzmann color fluid model and two experimentally determined relative permeability curves. The analysis of the relative permeability data followed the methodology proposed by Roy et al. (Front Phys 8:4, 2022) to reconstruct a constitutive equation for the co-moving velocity. Surprisingly, the coefficients of the constitutive equation were found to be nearly the same for all wetting conditions. On the basis of these results, a simple approach was proposed to reconstruct the relative permeability of the oil phase using only the co-moving velocity relationship and the relative permeability of the water phase. This proposed method provides new information on the interdependence of the relative permeability curves, which has implications for the history matching of production data and the solution of the associated inverse problem. The research findings contribute to a better understanding of the impact of wettability on fluid flow in porous media and provide a practical approach for estimating relative permeability based on the co-moving velocity relationship, which has never been shown before.

Engineering↗
DOE OSTI2021

The exo-β- N -acetylmuramidase NamZ from Bacillus subtilis is the founding member of a family of exo-lytic peptidoglycan hexosaminidases

Endo-β-N-acetylmuramidases, commonly known as lysozymes, are well-characterized antimicrobial enzymes that catalyze an endo-lytic cleavage of peptidoglycan; i.e., they hydrolyze the β-1,4-glycosidic bonds connecting N-acetylmuramic acid (MurNAc) and N-acetylglucosamine (GlcNAc). In contrast, little is known about exo-β-N-acetylmuramidases, which catalyze an exo-lytic cleavage of β-1,4-MurNAc entities from the non-reducing ends of peptidoglycan chains. Such an enzyme was identified earlier in the bacterium Bacillus subtilis, but the corresponding gene has remained unknown so far. We now report that ybbC of B. subtilis, renamed namZ, encodes the reported exo-β-N-acetylmuramidase. A ΔnamZ mutant accumulated specific cell wall fragments and showed growth defects under starvation conditions, indicating a role of NamZ in cell wall turnover and recycling. Recombinant NamZ protein specifically hydrolyzed the artificial substrate para-nitrophenyl β-MurNAc and the peptidoglycan-derived disaccharide MurNAc-β-1,4-GlcNAc. Together with the exo-β-N-acetylglucosaminidase NagZ and the exo-muramoyl-L-alanine amidase AmiE, NamZ degraded intact peptidoglycan by sequential hydrolysis from the non-reducing ends. A structure model of NamZ, built on the basis of two crystal structures of putative orthologs from Bacteroides fragilis, revealed a two-domain structure including a Rossmann-fold-like domain that constitutes a unique glycosidase fold. Thus, NamZ, a member of the DUF1343 protein family of unknown function, is now classified as the founding member of a new family of glycosidases (CAZy GH171). NamZ-like peptidoglycan hexosaminidases are mainly present in the phylum Bacteroidetes and less frequently found in individual genomes within Firmicutes (Bacilli, Clostridia), Actinobacteria, and γ-proteobacteria.

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

Illuminating Ligand Field Contributions to Molecular Qubit Spin Relaxation via T 1 Anisotropy

Electron spin relaxation in paramagnetic transition metal complexes constitutes a key limitation on the growth of molecular quantum information science. However, there exist very few experimental observables for probing spin relaxation mechanisms, leading to a proliferation of inconsistent theoretical models. Here we demonstrate that spin relaxation anisotropy in pulsed electron paramagnetic resonance is a powerful spectroscopic probe for molecular spin dynamics across a library of highly coherent Cu(II) and V(IV) complexes. Here, neither the static spin Hamiltonian anisotropy nor contemporary computational models of spin relaxation are able to account for the experimental T 1 anisotropy. Through analysis of the spin-orbit coupled wavefunctions, we derive an analytical theory for the T 1 anisotropy that accurately reproduces the average experimental anisotropy of 2.5. Furthermore, compound-by-compound deviations from the average anisotropy provide a promising approach for observing specific ligand field and vibronic excited state coupling effects on spin relaxation. Finally, we present a simple density functional theory workflow for computationally predicting T 1 anisotropy. Analysis of spin relaxation anisotropy leads to deeper fundamental understanding of spin-phonon coupling and relaxation mechanisms, promising to complement temperature-dependent relaxation rates as a key metric for understanding molecular spin qubits.

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

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Chaotic Dynamics Driven by Particle-Core Interactions

High-intensity beams in modern linacs are frequently encircled by diffuse halos, which drive sustained particle losses and result in gradual degradation of accelerating structures. In large part, the growth of halos is facilitated by internal space-charge forces within the beams, and detailed characterization of this process constitutes an active area of ongoing research. A partial understanding of dynamics that ensue within space-charge dominated beams is presented by the particle-core interaction paradigm – a mathematical model wherein single particle dynamics, subject to the collective potential of the core, are treated as a proxy for the broader behavior of the beam. In this work, we investigate the conditions for the onset of large-scale chaos within the framework of this model, and demonstrate that the propensity towards stochastic evolution is strongly dependent upon the charge distribution of the beam. In particular, we show that while particle motion within a uniformly charged beam is dominantly regular, rapid deterministic chaos readily arises within space-charge dominated Gaussian beams. Importantly, we find that for sufficiently high values of the beam’s space charge and beam pulsation amplitude, enhanced chaotic mixing between the core and the halo can lead to an enhanced radial diffusion of charged particles. We explain our results from analytic grounds by demonstrating that chaotic motion is driven by the intersection of two principal resonances of the system, and derive the relevant overlap conditions. Additionally, our analysis illuminates a close connection between the mathematical formulation of the particle-core interaction model and the Andoyer family of integrable Hamiltonians

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Catastrogenesis: DM, GWs, and PBHs from ALP string-wall networks

Axion-like particles (ALPs), a compelling candidate for dark matter (DM), are the pseudo Nambu-Goldstone bosons of a spontaneously and explicitly broken global U(1) symmetry. When the symmetry breaking happens after inflation, the ALP cosmology predicts the formation of a string-wall network which must annihilate early enough, producing gravitational waves (GWs) and primordial black holes (PBHs), as well as non-relativistic ALPs. We call this processcatastrogenesis. We show that, under the generic assumption that the potential has several degenerate minima, GWs from string-wall annihilation at temperatures below 100 eV could be detected by future CMB and astrometry probes, for ALPs with mass from 10 -16 to 10 6 eV. In this case, structure formation could limit ALPs to constitute a fraction of the DM and the annihilation would produce mostly "stupendously large" PBHs. For larger annihilation temperatures, ALPs can constitute 100% of DM, and the annihilation could produce supermassive black holes with a mass of up to 10 9 M ⊙ as found at the center of large galaxies. Furthermore our model can solve two mysteries, the nature of the DM and the origin of these black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Rapidly evolving aerosol emissions are a dangerous omission from near-term climate risk assessments

Abstract Anthropogenic aerosol emissions are expected to change rapidly over the coming decades, driving strong, spatially complex trends in temperature, hydroclimate, and extreme events both near and far from emission sources. Under-resourced, highly populated regions often bear the brunt of aerosols’ climate and air quality effects, amplifying risk through heightened exposure and vulnerability. However, many policy-facing evaluations of near-term climate risk, including those in the latest Intergovernmental Panel on Climate Change assessment report, underrepresent aerosols’ complex and regionally diverse climate effects, reducing them to a globally averaged offset to greenhouse gas warming. We argue that this constitutes a major missing element in society’s ability to prepare for future climate change. We outline a pathway towards progress and call for greater interaction between the aerosol research, impact modeling, scenario development, and risk assessment communities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Rain Drop Size Distributions Estimated from NOAA Snow-Level Radar Data

Using NOAA’s S-band High-Power Snow-Level Radar (HPSLR), a technique for estimating the rain drop size distribution (DSD) above the radar is presented. This technique assumes the DSD can be described by a four parameter, generalized gamma distribution (GGD). Using the radar’s measured average Doppler velocity spectrum and a value (assumed, measured, or estimated) of the vertical air motion w, an estimate of the GGD is obtained. Four different methods can be used to obtain w. One method that estimates a mean mass-weighted raindrop diameter D m from the measured reflectivity Z produces realistic DSDs compared to prior literature examples. These estimated DSDs provide evidence that the radar can retrieve the smaller drop sizes constituting the “drizzle” mode part of the DSD. Here, this estimation technique was applied to 19 h of observations from Hankins, North Carolina. Results support the concept that DSDs can be modeled using GGDs with a limited range of parameters. Further work is needed to validate the described technique for estimating DSDs in more varied precipitation types and to verify the vertical air motion estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Community Based Data of Uranium Adsorption onto Quartz

This data upload includes a compilation of experiments for quantifying uranium adsorption onto quartz. The provided .csv file has been compiled from the literature in a findable, accessible, interoperable, reusable (FAIR) data format. This was accomplished using the Lawrence Livermore National Laboratory Surface Complexation Database Converter (SCDC) code written in the R programming language (free licensing available at https://ipo.llnl.gov/technologies/software/llnl-surface-complexation-database-converter-scdc). This constitutes all current data compiled on uranium-quartz interactions in the L-SCIE (LLNL Surface Complexation/Ion Exchange) database (as of 02/28/2022). This data was used to develop surface complexation models that fit the global community dataset (https://doi.org/10.1021/acs.est.1c07109). The FAIR-formatted dataset also enables the implementation of alternative machine-learning approaches that can be explored in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Identifying Characteristics of Wildfire Towers and Troughs

Wildfire behavior is dictated by the complex interaction of numerous physical phenomena including dynamic ambient and fire-induced winds, heat transfer, aerodynamic drag on the wind by the fuel and combustion. These phenomena create complex feedback effects between the fire and its surroundings. In this study, we aim to study the mechanisms by which buoyant flame dynamics along with vortical motions and instabilities control wildfire propagation. Specifically, this study employs a suite of simulations conducted with the physics-based coupled fire-atmosphere behavior model (FIRETEC). The simulations are initialized with a fire line and the fires are allowed to propagate on a grass bed, where the fuel heights and wind conditions are varied systematically. Flow variables are extracted to identify the characteristics of the alternating counter-rotational vortices, called towers and troughs, that drive convective heat transfer and fire spread. These vortices have previously been observed in wildfires and laboratory fires, and have also been observed to arise spontaneously in FIRETEC due to the fundamental physics incorporated in the model. However, these past observations have been qualitative in nature and no quantitative studies can be found in the literature which connected these coherent structures fundamental to fire behavior with the constitutive flow variables. To that end, a variety of state variables are examined in the context of these coherent structures under various wind profile and grass height conditions. Identification of various correlated signatures and fire-atmosphere feedbacks in simulations provides a hypothesis that can be tested in future observational or experimental efforts, potentially assisting experimental design, and can aid in the interpretation of data from in situ detectors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Polycrystal thermo-elasticity revisited: theory and applications

The self-consistent (SC) theory is the most commonly used mean-field homogenization method to estimate the mechanical response behavior of polycrystals based on the knowledge of the properties and orientation distribution of constituent single-crystal grains. The original elastic SC method can be extended to thermo-elasticity by adding a stress-free strain to an elastic constitutive relation that expresses stress as a linear function of strain. With the addition of this independent term, the problem remains linear. Although the thermo-elastic self-consistent (TESC) model has important theoretical implications for the development of self-consistent homogenization of non-linear polycrystals, in this paper, we focus on TESC applications to actual thermo-elastic problems involving non-cubic (i.e. thermally anisotropic) materials. To achieve this aim, here we provide a thorough description of the TESC theory, which is followed by illustrative examples involving cooling of polycrystalline non-cubic metals. The TESC model allows studying the effect of crystallographic texture and single-crystal elastic and thermal anisotropy on the effective thermo-elastic response of the aggregate and on the internal stresses that develop at the local level.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

New insights into fuel blending effects: Intermolecular chemical kinetic interactions affecting autoignition times and intermediate-temperature heat release

Fuel blending effects on chemically-dominated fuel properties, such as gasoline anti-knock quality, are influenced by fundamental chemical kinetic interactions between the blending agent and the base fuel. Historically, quantification of such interactions has focused on changes to the radical pool, including $\dot{O}$H and HO 2 , while intermolecular interactions pertaining to carbonated, non-fuel-specific intermediates are typically overlooked. In this regard, this work aims to derive new insight into intermolecular chemical kinetic interactions that are intrinsic to fuel blending effects via a case study on blends of 0–30% ethanol (by volume) into FGF-LLNL (a multi-component gasoline surrogate for FACE-F research gasoline) using a rapid compression machine at a diluted/stoichiometric fuel loading, compressed pressure of 40 bar and low- to intermediate-temperature regimes that are representative of boosted SI engine operation. Ethanol blending effects on the intermediate temperature heat release (ITHR) of FGF-LLNL are characterized using experimental measurements, where ethanol is found to promote the extent of ITHR and suppress the transition from ITHR to main ignition. Chemical kinetic modeling is undertaken using a recently updated gasoline surrogate model. Sensitivity analyses on ITHR characteristics further corroborate the ethanol blending effects, and highlight the significant dependence of ITHR on both fuel-specific and non-fuel-specific reactions. An approach allowing comprehensive characterization of the complex intermolecular chemical kinetic interactions between constitutes in a fuel blend is then proposed. Application of the approach to FGF-LLNL/E0–E30 reveals that ethanol perturbs the heat release and autoignition characteristics of FGF-LLNL not only by directly changing the $\dot{O}$H and HO 2 radical pools via fuel-specific reactions, but also through intermolecular interactions where participating intermediates can be produced and consumed by various sub-chemistries. Disabling the intermolecular interactions in carbonated species between ethanol and FGF-LLNL sub-chemistries leads to somewhat slower ITHR evolution and lower ignition reactivity. The role of the individual intermolecular interaction is also characterized using the proposed approach. Lastly, implications of intermolecular interactions for future studies that aim to improve model performance are highlighted, where it is found that further investigations on the core C0-C4 chemistry are needed for developing highly accurate chemistry models for complex fuel blends.

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

Observable signatures of initial state momentum anisotropies in nuclear collisions

In this work, we show that the correlation between the elliptic momentum anisotropy v 2 and the average transverse momentum [ p T ] at fixed multiplicity in small system nuclear collisions carries information on the origin of the observed momentum anisotropy. A calculation using a hybrid IP-Glasma+ music +UrQMD model that includes contributions from final state response to the initial geometry as well as initial state momentum anisotropies of the color glass condensate predicts a characteristic sign change of the correlator ρ ^ ( v 2 2 , [ p T ] ) as a function of charged particle multiplicity in p + Au and d + Au collisions at s = 200 GeV , and p + Pb collisions at s = 5.02 TeV . This sign change is absent in calculations without initial state momentum anisotropies. The model further predicts a qualitative difference between the centrality dependence of ρ ^ ( v 2 2 , [ p T ] ) in Au + Au collisions at s = 200 GeV and Pb + Pb collisions at s = 5.02 TeV , with only the latter showing a sign change in peripheral events. Predictions for O + O collisions at different collision energy show a similar behavior. Experimental observation of these distinct qualitative features of ρ ^ ( v 2 2 , [ p T ] ) in small and large systems would constitute strong evidence for the presence and importance of initial state momentum anisotropies predicted by the color glass condensate effective theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Tropical High Cloud Feedback Relationships to Climate Sensitivity

Clouds constitute a large portion of uncertainty in predictions of equilibrium climate sensitivity (ECS). While low cloud feedbacks have been the focus of intermodel studies due to their high variability among global climate models, tropical high cloud feedbacks also exhibit considerable uncertainty. Here, we apply the cloud radiative kernel technique of Zelinka et al. to 22 models across the CMIP5 and CMIP6 ensembles to survey tropical high cloud feedbacks and analyze their relationship to ECS. We find that the net high cloud feedback and its altitude and optical depth feedback components are significantly positively correlated with ECS in the tropical mean. On the other hand, the tropical mean high cloud amount feedback is not correlated with ECS. These relationships are most pronounced outside of areas of strong climatological ascent, suggesting the importance of thin cirrus feedbacks. Finally, we explore connections between high cloud feedbacks, climate sensitivity, and mean state high cloud properties. In general, high ECS models are cloudier in the upper troposphere but have a thinner high cloud population. Furthermore, we find that having more thin cirrus in the mean state relates to more positive high cloud altitude and optical depth feedbacks, and it either amplifies or dampens the high cloud amount feedback depending on the large-scale dynamical regime (amplifying in descent and dampening in ascent). In summary, our analysis highlights the importance of tropical high cloud feedbacks for driving intermodel spread in ECS and suggests that mean state high cloud characteristics might provide a unique opportunity for observationally constraining high cloud feedbacks.

atmosphere↗
DOE OSTI2022

Not all types of secondary organic aerosol mix: two phases observed when mixing different secondary organic aerosol types

Abstract. Secondary organic aerosol (SOA) constitutes a large fraction of atmospheric aerosol. To assess its impacts on climate and air pollution, knowledge of the number of phases in internal mixtures of different SOA types is required. Atmospheric models often assume that different SOA types form a single phase when mixed. Here, we present visual observations of the number of phases formed after mixing different anthropogenic and biogenic SOA types. Mixing SOA types generated in environmental chambers with oxygen-to-carbon (O/C) ratios between 0.34 and 1.05, we found 6 out of 15 mixtures of two SOA types to result in two phase particles. We demonstrate that the number of phases depends on the difference in the average O/C ratio between the two SOA types (Δ(O/C)). Using a threshold Δ(O/C) of 0.47, we can predict the phase behavior of over 90 % of our mixtures, with one- and two-phase particles predicted for Δ(O/C)<0.47 and Δ(O/C)≥0.47, respectively. This threshold ΔO/C value provides a simple parameter to predict whether mixtures of fresh and aged SOA form one- or two-phase particles in the atmosphere. In addition, we show that phase-separated SOA particles form when mixtures of volatile organic compounds emitted from real trees are oxidized.

Mahrt, Fabian (ORCID:0000000270596765)↗
DOE OSTI2020

Coupled Time-Lapse Full-Waveform Inversion for Subsurface Flow Problems Using Intrusive Automatic Differentiation

We describe a novel framework for estimating subsurface properties, such as rock permeability and porosity, from time-lapse observed seismic data by coupling full-waveform inversion (FWI), subsurface flow processes, and rock physics models. For the inverse modeling, we handle the back propagation of gradients by an intrusive automatic differentiation strategy that offers three levels of user control: (1) At the wave physics level, we adopted the discrete adjoint method in order to use our existing high-performance FWI code; (2) at the rock physics level, we used built-in automatic differentiation operators from the TensorFlow backend; (3) at the flow physics level, we implemented customized partial differential equation (PDE) operators for the multiphase flow equations. The three-level coupled inversion strategy strikes a good balance between computational efficiency and programming efforts, and when the gradients are chained together, it constitutes a coupled inverse system. Our numerical experiments demonstrate that the three-level coupled inverse problem is superior in terms of accuracy to a traditional decoupled inversion strategy. Additionally, our method is able to simultaneously invert for parameters in empirical relationships such as the rock physics models. Our proposed inverted model can be used for reservoir performance prediction and reservoir management/optimization purposes.

54 ENVIRONMENTAL SCIENCES↗