Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable predictive control”

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 91 records · Page 5

Measurement of 𝑑 2⁢ 𝜎/𝑑⁢|$\vec{q}$|⁢𝑑⁢𝐸 avail in charged current 𝜈 𝜇 -nucleus interactions at ⟨𝐸 𝜈 ⟩=1.86 GeV using the NOvA Near Detector

Double- and single-differential cross sections for inclusive charged-current 𝜈 𝜇 -nucleus scattering are reported for the kinematic domain 0 to 2 GeV/𝑐 in three-momentum transfer and 0 to 2 GeV in available energy, at a mean 𝜈 𝜇 energy of 1.86 GeV. The measurements are based on an estimated 995,760 𝜈 𝜇 charged-current (CC) interactions in the scintillator medium of the NOvA Near Detector. The subdomain populated by 2-particle-2-hole (2p2h) reactions is identified by the cross section excess relative to predictions for 𝜈 𝜇 -nucleus scattering that are constrained by a data control sample. Models for 2-particle-2-hole processes are rated by 𝜒 2 comparisons of the predicted-versus-measured 𝜈 𝜇 CC inclusive cross section over the full phase space and in the restricted subdomain. Shortfalls are observed in neutrino generator predictions obtained using the theory-based València and SuSAv2 2p2h models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Scientific Discovery with Physics-Informed System Identification (Abbreviated Report)

My fellowship research focused on making physics-based simulations faster and more useful through machine learning. Many problems in science and engineering are governed by partial differential equations, but high-fidelity simulations are often too expensive to run repeatedly. I worked on improving Latent Space Dynamics Identification (LaSDI), a reduced-order modeling framework that compresses large simulation data sets into a smaller representation and then learns how that representation evolves over time. The motivation was to develop reduced models that remain accurate for more challenging systems, especially when predictions must remain reliable over long time intervals or when the underlying dynamics are more complicated than standard methods can easily handle. I also contributed to related work on Quandary, a high-performance software effort for simulation and control of open quantum systems, before focusing primarily on Latent Space Dynamics Identification methods. The main outcomes of the fellowship were two new algorithms (both of which were published), Rollout-LaSDI and Higher-Order LaSDI, together with supporting work on multi-stage Latent Space Dynamics Identification. Rollout-LaSDI improved long-term prediction by training the model to stay accurate over extended time horizons, and Higher-Order LaSDI broadened the method so it could model systems with higher-order time dynamics. My contributions to multistage Latent Space Dynamics Identification also helped show that its later training stages could be simplified without losing effectiveness, and that this behavior held across different model architectures and training strategies. Taken together, these advances improved the accuracy, flexibility, and practical value of reduced-order modeling tools for computational science.

97 MATHEMATICS AND COMPUTING↗

Four-dimensional phase-space reconstruction of flat and magnetized beams using neural networks and differentiable simulations

Beams with cross-plane coupling or extreme asymmetries between the two transverse phase spaces are often encountered in particle accelerators. Flat beams with large transverse-emittance ratios are critical for future linear colliders. Similarly, magnetized beams with significant cross-plane coupling are expected to enhance the performance of electron cooling in hadron beams. Preparing these beams requires precise control and characterization of the four-dimensional transverse phase space. In this study, we employ generative phase-space reconstruction techniques to rapidly characterize magnetized and flat-beam phase-space distributions using a conventional quadrupole-scan method. The reconstruction technique is experimentally demonstrated on an electron beam produced at the Argonne Wakefield Accelerator and successfully benchmarked against conventional diagnostics techniques. Specifically, we show that predicted beam parameters from the reconstructed phase-space distributions (e.g., as magnetization and flat-beam emittances) are in excellent agreement with those measured from the conventional diagnostic methods. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Commentary: Duckweeds as model organisms for metabolic studies

Duckweeds have many practical applications, for example in human nutrition, as animal feed, in the production of bioplastics or vaccines, and phytoremediation (Acosta et al., 2021). Under most conditions, they reproduce asexually which provides genetically uniform material with predictable patterns of growth that make them ideal as sentinel organisms for phytotoxicity testing (Park et al., 2021). Asexual growth also results in high biomass production which makes duckweeds promising candidates as biofuel feedstocks (Acosta et al., 2021; Liang et al., 2023). In addition, duckweed species like Lemna minor and Spirodela polyrhiza are also reemerging as model organisms in plant biology as high-quality full genome assemblies and other genomic resources become available (Chang et al., 2016; Acosta et al., 2021). We argue that duckweed species are particularly of interest for the study of primary plant metabolism. Primary metabolism concerns the part of metabolism that is directly involved in the growth and development of plants, and which tends to be highly conserved among plant species. What makes duckweeds particularly attractive is that when grown on liquid media more precise control of physiological conditions can be attained relative to growth of plants in soil. Also, due to their relatively simple anatomical structure and asexual reproduction of fronds by budding, precise characterization of the physiological state under study is possible through one simple metric, i.e., the specific growth rate (rate of dry weight increase per existing dry weight), which can be incorporated relatively easily into metabolic models. This is not possible for land plants, such as Arabidopsis, because over the course of their life cycle, they go through multiple growth stages and phases of anatomical differentiation, which are much more complex to quantify. Furthermore, duckweeds can grow on organic substrates under heterotrophic or photomixotrophic conditions that facilitate isotope tracer studies. For example, in a previous study on duckweed by one of the authors, Lemna gibba (L). was grown on glucose with a position-specific 13 C-label that can be detected and resolved by Mass Spectrometry or Nuclear Magnetic Resonance spectrometry. Using this approach, the 13 C-label was traced into biomass compounds formed from glucose, particularly isoprenoid compounds. Some of the resulting labeling patterns were in apparent disagreement with predictions based on known metabolic pathways for the biosynthesis of isopentenyl pyrophosphate, the universal building block for isoprenoids, (Lichtenthaler et al., 1997). From this data it was deduced that isoprenoid compounds such as carotenoids and isoprenoid chains of phytol and plastoquinone, synthesized in the chloroplast, are produced via a previously unreported plant metabolic pathway, now known as the methylerythitol/deoxyxylulose-5-phosphate pathway (Lichtenthaler et al., 1997).

59 BASIC BIOLOGICAL SCIENCES↗

Probabilistic error estimation for non-intrusive reduced models learned from data of systems governed by linear parabolic partial differential equations

This work derives a residual-based a posteriori error estimator for reduced models learned with non-intrusive model reduction from data of high-dimensional systems governed by linear parabolic partial differential equations with control inputs. It is shown that quantities that are necessary for the error estimator can be either obtained exactly as the solutions of least-squares problems in a non-intrusive way from data such as initial conditions, control inputs, and high-dimensional solution trajectories or bounded in a probabilistic sense. Here, the computational procedure follows an offline/online decomposition. In the offline (training) phase, the high-dimensional system is judiciously solved in a black-box fashion to generate data and to set up the error estimator. In the online phase, the estimator is used to bound the error of the reduced-model predictions for new initial conditions and new control inputs without recourse to the high-dimensional system. Numerical results demonstrate the workflow of the proposed approach from data to reduced models to certified predictions.

97 MATHEMATICS AND COMPUTING↗

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗

Interpreting T-History Data on Supercooling and Phase Change Enthalpy of a Salt-Hydrogel Complex for Building Thermal Management

Phase change materials (PCMs) are promising materials for storing and discharging latent heat in buildings to regulate the thermal environment. Measurement of salt hydrate PCM thermal performance is often challenging because of supercooling and phase segregation. Though differential scanning calorimetry (DSC) has been widely used to measure thermal properties of various phase change materials, DSC uses milligrams of sample and is unrealistic for predicting large-scale materials performance. The T-history method [1], in contrast, handles higher volumes but the analysis typically requires a lumped capacitance assumption. Here, we report and analyze T-History data on a novel salt-hydrogel complex with high viscosity, that renders low Biot number sample preparation extremely challenging. We re-examine the standard T-history method and compared against DSC data in this complex material. We report details of an experimental setup, built with necessary additions to a controlled environmental chamber. We carefully examine the validity of the lumped capacitance assumption and suggest a modified analysis using a computational approach for reliable data. Experimental parameters such as ramp rate and temperature range of the environmental chamber are considered via computational modeling to understand the effect of ramp rate on supercooling [2]. Finally, we report thermal cycling experiments on the salt-hydrogel complex. Our modified approach to measuring supercooling and enthalpy of fusion at large scales is important in better understanding the performance of phase-change materials at scale for building thermal storage. References: [1] Marin, Jose M, Belen Zalba, Luisa F Cabeza, and Harald Mehling. “Determination of Enthalpy Temperature Curves of Phase Change Materials with the Temperature-History Method: Improvement to Temperature Dependent Properties.” Measurement Science and Technology 14, no. 2 (February 1, 2003): 184–89. https://doi.org/10.1088/0957-0233/14/2/305. [2] Safari, A., R. Saidur, F.A. Sulaiman, Yan Xu, and Joe Dong. “A Review on Supercooling of Phase Change Materials in Thermal Energy Storage Systems.” Renewable and Sustainable Energy Reviews 70 (April 2017): 905–19. https://doi.org/10.1016/j.rser.2016.11.272

thermodynamics↗

Differential effects of drought on nonstructural carbohydrate storage in seedlings and mature trees of four species in a subtropical forest

Nonstructural carbohydrates (NSC) play important roles in forest vulnerability to climate change, especially under increasing drought intensity and frequency. Understanding NSC dynamics is essential for accurately predicting the resistance and resilience of forests in response to drought. However, our knowledge of NSC re-sponses to drought is still limited due to the lack of research in trees of di?erent life stages. In this study, we conducted a throughfall exclusion experiment (TFE) with four subtropical evergreen tree species to examine drought e?ects on NSC in mature trees and understory seedlings. Our results showed the differential effects of drought on NSC dynamics of understory seedlings and mature trees. In the TFE experiment, mature trees of all four species were relatively homoeostatic with the insigni?cantly changed NSC pools, photosynthesis, and growth under the drought treatment compared to the control. In contrast, understory seedlings displayed sig-ni?cant decreases in total NSC and soluble sugars (-14.70 ± 3.66% and -16.93 ± 3.85%, respectively) with the exception of Lithocarpus glaber (L. glaber). The seedlings of L. glaber with the highest hydraulic resistance maintained or slightly increased NSC and its components in response to drought. Our study highlights the im-portance of life stage in assessing drought e?ects of trees on NSC storage and then forest C cycling, which could be incorporated into the dynamic global vegetation models (DGVMs) to better understand drought effects on forest C balance in the future.

Zhang, Peipei↗

Epitaxial stabilization and oxygen vacancy control of EuNiO 3 thin films

Rare-earth nickelates exhibit valuable behavior for neuromorphic computing at low temperature: Building blocks for biologically inspired microelectronic neurons like electrically driven insulator–metal transitions (IMTs), negative differential resistance, and self-oscillations have been shown up to 230 K for SmNiO 3 and NdNiO 3 . EuNiO 3 raises the IMT far above room temperature (460 K) but high-quality thin films are challenging to synthesize. Here, we explore the epitaxial stabilization of EuNiO 3 using pulsed laser deposition. X-ray diffraction reciprocal space maps, x-ray absorption spectroscopy, and transmission electron microscopy show that higher growth temperature (800 °C) reduces oxygen vacancy concentrations in EuNiO 3 . Pseudomorphic EuNiO 3 is demonstrated on both SrLaAlO 4 and NdGaO 3 substrates, and LaNiO 3 buffer layers are incorporated to facilitate future vertical device fabrication. In contrast to bulk thermodynamic predictions, the greater oxidation and crystallinity at higher temperature we observe indicates that epitaxial substrates can stabilize EuNiO 3 at O 2 pressures less than 1 atm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Molecular-scale Insights into Cooperativity Switching of x TAB Adsorption on Gold Nanoparticles

Quantifying adsorption behaviors is crucial for various applications such as catalysis, separation, and sensing, yet it is generally challenging to access in solution. Here, we report a combined experimental and computational study of the adsorption behaviors of alkyl-trimethylammonium bromides (xTAB), a class of ligands important for colloidal nanoparticle stabilization and shape control, with various alkyl chain lengths x on Au nanoparticles. We use density functional theory (DFT) to calculate xTAB binding energies on Au{111} and Au{110} surfaces with standing-up and lying-down configurations, which provides insights into the adsorption affinity and cooperativity differences of xTAB on these two facets. We demonstrate the key role of van der Waals interactions in determining the xTAB adsorption behavior. These computational results predict and explain the experimental discovery of xTAB’s adsorption behavior switch from stronger affinity, negative cooperativity to weaker affinity, positive cooperativity when the concentration of xTAB increases in solution. We also show that in the standing-up configuration, bilayer adsorption may occur on both facets, which can lead to different differential binding energies and consequently adsorption crossover between the two facets when the ligand concentration increases. Our combined experimental and computational approaches demonstrate a paradigm for gaining molecular-scale insights into adsorbate–surface interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluation of Digital Twin Modeling and Simulation

A digital twin has intelligent modules that continuously monitor the condition of the individual components and the whole of a system. Digital twins can provide nuclear power plants (NPP) operators an unprecedented level of monitoring, control, supervision, and security by contributing a greater volume of data for more comprehensive data analysis and increased accuracy of insights and predictions for decision making throughout the entire NPP lifecycle. NPP operators and managers have historically relied on limited, second hand or incomplete data. With proper implementation, digital twins can provide a central hub of all intel that allows for a multidisciplinary view of an NPP. This equips operators and managers with the ability to have more information, context, and intel that can be used for greater granularity during planning and decision making. Digital twins can be used in many activities as the technology has many different concepts surrounding it. From the various definitions of a digital twin within the industry, digital twins can be differentiated by levels of integration/automation. The three main models include digital model, digital shadow, and digital twin. Digital twins offer many potential advancements to the nuclear industry that could reduce costs, improve designs, provide safer operation, and improve their overall security.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Convergence of Eigenvector Continuation

Eigenvector continuation is a computational method that finds the extremal eigenvalues and eigenvectors of a Hamiltonian matrix with one or more control parameters. It does this by projection onto a subspace of eigenvectors corresponding to selected training values of the control parameters. The method has proven to be very efficient and accurate for interpolating and extrapolating eigenvectors. However, almost nothing is known about how the method converges, and its rapid convergence properties have remained mysterious. In this Letter, we present the first study of the convergence of eigenvector continuation. In order to perform the mathematical analysis, we introduce a new variant of eigenvector continuation that we call vector continuation. We first prove that eigenvector continuation and vector continuation have identical convergence properties and then analyze the convergence of vector continuation. Our analysis shows that, in general, eigenvector continuation converges more rapidly than perturbation theory. The faster convergence is achieved by eliminating a phenomenon that we call differential folding, the interference between nonorthogonal vectors appearing at different orders in perturbation theory. From our analysis we can predict how eigenvector continuation converges both inside and outside the radius of convergence of perturbation theory. Further, while eigenvector continuation is a nonperturbative method, we show that its rate of convergence can be deduced from power series expansions of the eigenvectors. Our results also yield new insights into the nature of divergences in perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CyDER: A Cyber Physical Co-simulation Platform for Distributed Energy Resources in Smartgrids

The CyDER project aimed at developing an open-source, modular and scalable co-simulation platform for power grids with large shares of Distributed Energy Resources (DERs). The project partners are the Lawrence Berkeley National Lab (LBNL), Lawrence Livermore National Lab (LLNL), PG&E, SolarCity, and ChargePoint. The prime recipient is LBNL; SolarCity and ChargePoint were partners for the project’s first two years. Increased DER integration introduces a number of challenges in power grid operation including a more dynamic interaction between the transmission grid and distribution grids, and increased modeling complexity. Although specialized software exists to precisely model different components of the power system, it is far from trivial to integrate all various models and perform a holistic simulation. Instead of replicating all models in a common simulation program, a commonly accepted approach to tackle this model diversity is to couple third-party simulators and models through a co-simulation platform that coordinates information exchange among the various components. Following this line of research, this project’s objective was to develop a co-simulation platform based on a widely accepted industrial standard called Functional Mockup Interface (FMI). Within this process, the project developed models compliant with the FMI standard, called Functional Mockup Units (FMUs), and used them to perform various operational and planning power system analyses. Relying and building upon an industrial standard is the main differentiation of this project compared with previous or parallel efforts in the co-simulation area. Particular emphasis was put on delivering software utilities to facilitate setting up and running co-simulations by end-users. Furthermore, a strong aspect of this project is demonstrating that co-simulation techniques can be used to perform Hardware-in-the-Loop (HIL) simulations that couple software components (e.g., simulated models) with hardware components (e.g., real devices such PV systems and batteries). The long-term goal of CyDER project is to help establish FMI as a powerful standard for co-simulation and promote adoption by electric utilities and other interested stakeholders. The main accomplishments of the project include the development of several FMUs including distribution and transmission grid models, PV inverters with Volt/Var/Watt controllers, batteries, and predictive optimal controllers. Additionally, a unique software package was developed, called SimulatorToFMU, which is capable of exporting any Python-driven simulator or Python script as an FMU. This is an important contribution towards establishing FMI as one of the main co-simulation standards, because more and more third-party programs for sub-system modeling and simulation are delivered with Python APIs. The CyDER platform was used to perform PV hosting capacity analyses in real utility feeders with and without smart inverter controls, battery storage, and EV charging. Smart inverter controls include conventional Volt/Var/Watt controls for reactive power support and active power curtailment, but also predictive controls that optimize the charging and discharging profile of the battery connected on the DC side in order to minimize the customer’s economic benefit. Finally, an important result of this project is delivering an experimental setup that consists of residential-scale PV inverters with battery storage, a real-time grid simulator with an ideal voltage source as grid emulator, and micro Phasor Measurement Units (PMUs). All these components and additional software modules are coupled to one another using the FMI standard and can be co-simulated with the CyDER platform.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Formation kinetics and thermodynamic stability of the Aurivillius compounds in Bi 4 Ti 3 O 12 –BiFeO 3 system

Abstract The Aurivillius compounds in the Bi 2 O 3 –Fe 2 O 3 –TiO 2 system, combining ferroelectric, semiconducting, and ferromagnetic properties, have attracted particular interest. Formation kinetics and thermodynamic stability are the fundamental knowledge needed for modeling and predicting the temporal microstructure and property evolution during materials processing but have not yet been addressed by quantitative experimental measurement. This article focuses on the Bi n +1 Fe n –3 Ti 3 O 3 n +3 Aurivillius compounds on the Bi 4 Ti 3 O 12 –BiFeO 3 tie‐line to elucidate the mechanisms and thermodynamic controls responsible for phase formation of compounds with various perovskite‐like layers. Five high‐purity Aurivillius compounds Bi 4 Ti 3 O 12 , Bi 5 FeTi 3 O 15 , Bi 6 Fe 2 Ti 3 O 18 , Bi 7 Fe 3 Ti 3 O 21 , and Bi 8 Fe 4 Ti 3 O 24 with integer n = 3–7 values were synthesized and their phase transformation properties and enthalpies of formation were studied by X‐ray diffraction in situ, high temperature differential scanning calorimetry, and high temperature oxide melt solution calorimetry. Thermodynamic stability of the compounds decreases with increasing n , and formation kinetics gradually slow down, demonstrating the inherent difficulty to synthesize pure Aurivillius compounds with n larger than 8. This difficulty was confirmed by an impurity phase coexisting with Bi 9 Fe 5 Ti 3 O 27 .

Gong, Weiping↗

Taylor–Couette flow for astrophysical purposes

A concise review is given of astrophysically motivated experimental and theoretical research on Taylor–Couette flow. The flows of interest rotate differentially with the inner cylinder faster than the outer, but are linearly stable against Rayleigh’s inviscid centrifugal instability. At shear Reynolds numbers as large as 10 6 , hydrodynamic flows of this type (quasi-Keplerian) appear to be nonlinearly stable: no turbulence is seen that cannot be attributed to interaction with the axial boundaries, rather than the radial shear itself. Direct numerical simulations agree, although they cannot yet reach such high Reynolds numbers. This result indicates that accretion-disc turbulence is not purely hydrodynamic in origin, at least insofar as it is driven by radial shear. Theory, however, predicts linear magnetohydrodynamic (MHD) instabilities in astrophysical discs: in particular, the standard magnetorotational instability (SMRI). MHD Taylor–Couette experiments aimed at SMRI are challenged by the low magnetic Prandtl numbers of liquid metals. High fluid Reynolds numbers and careful control of the axial boundaries are required. The quest for laboratory SMRI has been rewarded with the discovery of some interesting inductionless cousins of SMRI, and with the recently reported success in demonstrating SMRI itself using conducting axial boundaries. Some outstanding questions and near-future prospects are discussed, especially in connection with astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Surface Carboxylate Sensitivity to Electron and Hole Relaxation in Photoexcited Cadmium Sulfide Nanocrystals

Understanding how passivating surface ligands couple to excitonic states in nanocrystal photocatalysts is crucial for controlling nonradiative relaxation pathways which compete with interfacial charge transfer. Here, we report femtosecond transient infrared (IR) spectroscopy to resolve ∼100 fs ligand-exciton coupling between 1S exciton states in oleate-capped cadmium sulfide (CdS) nanocrystals and vibrational modes of surface carboxylates. Differential mid-IR spectra show distinct negative amplitude and positive photoinduced absorption signals at ∼1540 cm –1 (carboxylate asymmetric stretch) and ∼1440 cm –1 (carboxylate symmetric stretch), respectively. Fluence-dependent transient IR measurements reveal that the symmetric stretch is uniquely sensitive to picosecond Auger recombination, while the asymmetric stretch shows no analogous decay. Our results provide direct measurement of femtosecond ligand-exciton coupling in CdS nanocrystals and demonstrate how surface-bound carboxylate ligands serve as carrier-specific reporters of nanocrystal photophysics. Furthermore, these findings offer critical insights for designing and developing predictive models for ligand-mediated strategies in next-generation nanocrystal photocatalysts.

Cadmium sulfide↗

Validation of prediction capability of operating space for plasma initiation in MAST-U

DYON is a plasma initiation modelling code that solves the differential equation system of the full circuit equations (plasma current, active coil currents and eddy currents in full passive structures) and 0D global energy and particle balance equations (Kim 2022 Nucl. Fusion 62 126012). In order to test the capability of the full electromagnetic plasma initiation model to predict individual discharges in experiments and thus the operating space in the device, a dedicated experimental database was built in MAST-U by scanning the prefilled gas pressure p 0 and the induced loop voltage V loop . In the experimental operating space of p 0 and V loop the lower and the upper limits of p 0 are determined by the plasma breakdown failure and the plasma burn-through failure, respectively. The lower limit of V loop is determined by the plasma burn-through failure. By directly reading the control room data used in each discharge (i.e. currents in the solenoid, poloidal field coils, and toroidal field coils, p 0 , and gas puffing rate), the full electromagnetic DYON consistently predicted the failed breakdown, failed burn-through, and successful plasma initiation discharges in the experimental database, demonstrating its capability to predict the operating space for inductive plasma initiation. The Paschen curve calculated with the effective connection length in MAST-U indicates a much higher p 0 required for plasma breakdown than the experimental data, indicating that individual field line evaluation is necessary to calculate the quantitative requirements for Townsend breakdown. The demonstration in this paper shows that the full electromagnetic DYON could be a useful simulation tool to assess the feasibility of inductive plasma initiation and to optimise operating scenarios in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗