Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “3D Modeling”

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

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

At least 1,297 records · Page 72

Feature‐shared adaptive‐boost deep learning for invasiveness classification of pulmonary subsolid nodules in CT images

Purpose In clinical practice, invasiveness is an important reference indicator for differentiating the malignant degree of subsolid pulmonary nodules. These nodules can be classified as atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA), or invasive adenocarcinoma (IAC). The automatic determination of a nodule's invasiveness based on chest CT scans can guide treatment planning. However, it is challenging, owing to the insufficiency of training data and their interclass similarity and intraclass variation. To address these challenges, we propose a two‐stage deep learning strategy for this task: prior‐feature learning followed by adaptive‐boost deep learning. Methods The adaptive‐boost deep learning is proposed to train a strong classifier for invasiveness classification of subsolid nodules in chest CT images, using multiple 3D convolutional neural network (CNN)‐based weak classifiers. Because ensembles of multiple deep 3D CNN models have a huge number of parameters and require large computing resources along with more training and testing time, the prior‐feature learning is proposed to reduce the computations by sharing the CNN layers between all weak classifiers. Using this strategy, all weak classifiers can be integrated into a single network. Results Tenfold cross validation of binary classification was conducted on a total of 1357 nodules, including 765 noninvasive (AAH and AIS) and 592 invasive nodules (MIA and IAC). Ablation experimental results indicated that the proposed binary classifier achieved an accuracy of with an AUC of 81.3 . These results are superior compared to those achieved by three experienced chest imaging specialists who achieved an accuracy of , , and , respectively. About 200 additional nodules were also collected. These nodules covered 50 cases for each category (AAH, AIS, MIA, and IAC, respectively). Both binary and multiple classifications were performed on these data and the results demonstrated that the proposed method definitely achieves better performance than the performance achieved by nonensemble deep learning methods. Conclusions It can be concluded that the proposed adaptive‐boost deep learning can significantly improve the performance of invasiveness classification of pulmonary subsolid nodules in CT images, while the prior‐feature learning can significantly reduce the total size of deep models. The promising results on clinical data show that the trained models can be used as an effective lung cancer screening tool in hospitals. Moreover, the proposed strategy can be easily extended to other similar classification tasks in 3D medical images.

Wang, Jun↗

A naturally grown three-dimensional nonlinear photonic crystal

Nonlinear frequency conversion via three-dimensional (3D) quasi-phase matching (QPM) process is experimentally realized based on a Ba{sub 0.77}Ca{sub 0.23}TiO{sub 3} (BCT) crystal. The ferroelectric domains in BCT crystal are observed, and the results reveal that the antiparallel domains distribute in three dimensions and can provide 3D reciprocal lattice vectors for QPM processes. Broadband petal-like second-harmonic patterns are achieved, which are well consistent with the theoretical quasi-cubic model of 3D nonlinear photonic crystals. Our work not only promotes the development of QPM technique but also builds a platform for 3D nonlinear optics and quantum optics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

QUIC-Fire: 3D Fire-Atmosphere Feedback Model for Wildland Fire Management

Los Alamos National Laboratory, U.S. Forest Service, and Tall Timbers Research Station developed QUIC-Fire, the first laptop-capable, real-time wildland fire prediction software that captures the critical influences of 3D vegetation structure, interactions between multiple fires, variable winds, and complex topography at meter-scale resolutions. Prescribed burns are needed to prevent catastrophic wildfires in the U.S. QUIC-Fire transforms a prescribed fire manager’s ability to assess risk, optimize fuel treatments, and plan prescribed burns. Los Alamos, U.S. Forest Service, Tall Timbers Research Station, U.S. Fish and Wildlife Service, Spatial Information Group, and Applied Interagency Prescribed Training Center have used QUIC-Fire.

54 ENVIRONMENTAL SCIENCES↗

Examining Runner’s Outdoor Heat Exposure Using Urban Microclimate Modeling and GPS Trajectory Mining

It is important to quantify human heat exposure in order to evaluate and mitigate the negative impacts of heat on human well-being in the context of global warming. This study proposed a human-centric framework to examine human personal heat exposure based on anonymous GPS trajectories data mining and urban microclimate modeling. The mean radiant temperature (Tmrt) that represents the human body’s energy balance was used to indicate human heat exposure. The meteorological data and high-resolution 3D urban model generated from multispectral remotely sensed images and LiDAR data were used as inputs in urban microclimate modeling to map the spatio-temporal distribution of the Tmrt in the Boston metropolitan area. The anonymous human GPS trajectory data collected from fitness Apps was used to map the spatiotemporal distribution of human outdoor activities. By overlaying the anonymous GPS trajectories on the generated spatio-temporal maps of Tmrt, this study further examined the heat exposure of runners in different age-gender groups in the Boston area. Results show that there is no significant difference in terms of heat exposure for female and male runners. The female runners in the age of 45-54 are exposed to more heat than female runners of 18-24 and 25-34, while there is no significant difference among male runners. This study proposed a novel method to estimate human heat exposure, which would shed new light on mitigating the negative impacts of heat on human health.

Personal heat exposure↗

Modeling of Hydraulic Fracture Propagation at the kISMET Site Using a Fully Coupled 3D Network-Flow and Quasi- Static Discrete Element Model

Aimed at supporting the design of hydraulic fracturing experiments at the kISMET site, ~1500 m below ground in a deep mine, we performed pre-experimental hydraulic fracturing simulations in order to estimate the breakdown pressure, propagation pressure, fracture geometry, and the magnitude of induced seismicity using a newly developed fully coupled three-dimensional (3D) network flow and quasi-static discrete element model (DEM). The quasi-static DEM model, which is constructed by Delaunay tessellation of the rock volume, considers rock fabric heterogeneities by using the “disordered” DEM mesh and adding random perturbations to the stiffness and tensile/shear strengths of individual DEM elements and the elastic beams between them. A conjugate 3D flow network based on the DEM lattice is constructed to calculate the fluid flow in both the fracture and porous matrix. One distinctive advantage of the model is that fracturing is naturally described by the breakage of elastic beams between DEM elements. It is also extremely convenient to introduce mechanical anisotropy into the model by simply assigning orientation-dependent tensile/shear strengths to the elastic beams. In this paper, the 3D hydraulic fracturing model was verified against the analytic solution for a penny-shaped crack model. We applied the model to simulate fracture propagation from a vertical open borehole based on initial estimates of rock mechanical properties and in-situ stress conditions. The breakdown pressure and propagation pressure are directly obtained from the simulation. In addition, the released elastic strain energies of individual fracturing events were calculated and used as a conservative estimate for the magnitudes of the potential induced seismic activities associated with fracturing. The comparisons between model predictions and experimental results are still ongoing.

15 GEOTHERMAL ENERGY↗

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗

Antineutrino analysis for continuous monitoring of nuclear reactors: Sensitivity study

This paper explores the various contributors to uncertainty on predictions of the antineutrino source term which is used for reactor antineutrino experiments and is proposed as a safeguard mechanism for future reactor installations. The errors introduced during simulation of the reactor burnup cycle from variation in nuclear reaction cross sections, operating power, and other factors are combined with those from experimental and predicted antineutrino yields, resulting from fissions, evaluated, and compared. The most significant contributor to uncertainty on the reactor antineutrino source term when the reactor was modeled in 3D fidelity with assembly-level heterogeneity was found to be the uncertainty on the antineutrino yields. Using the reactor simulation uncertainty data, the dedicated observation of a rigorously modeled small, fast reactor by a few-ton near-field detector was estimated to offer reduction of uncertainty on antineutrino yields in the 3.0–6.5 MeV range to a few percent for the primary power-producing fuel isotopes, even with zero prior knowledge of the yields.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

79 ASTRONOMY AND ASTROPHYSICS↗

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗

Integrating Analytical Solutions and U-Net Model for Predicting Groundwater Contaminant Plumes in Pump-and-Treat Systems

Pump-and-treat (P&T) is a common technique for groundwater remediation involving the extraction and treatment of contaminated water above ground. Optimizing the design and operation of the P&T well network is essential for maximizing the system’s effectiveness and efficiency. However, this optimization often necessitates many model evaluations, leading to computationally demanding tasks. This study introduces a novel approach that integrates analytical solutions for groundwater dynamics with the U-Net (Ronneberger et al., 2015) deep learning framework to predict groundwater contaminant plume migration under dynamic pumping conditions. By incorporating the Thiem equation (Thiem, 1906) into the input preprocessing, the U-Net model transforms sparse well data into a continuous spatial field that captures the hydraulic impacts of pumping activities. This integration enables the model to leverage both deep learning capabilities and classical physics-based groundwater theories, enhancing prediction accuracy and computational efficiency. These advancements can facilitate rapid, large-scale evaluations of P&T optimization simulations, allowing for timely and effective decision-making in well placement and system management. We demonstrate the model's robust performance across both simplified transient 2D models and a more complex 3D heterogeneous site model at the 200 West P&T facility at the Hanford Site. The U-Net-based model offers substantial computational advantages, reducing simulation times significantly compared to full physics-based models and providing a powerful tool for rapid site evaluation and P&T system optimization, such as evaluating alternative P&T well network designs. Our findings highlight the potential of advanced machine learning models to significantly enhance the efficiency and sustainability of groundwater remediation efforts, offering a novel application of U-Net architecture in environmental science.

Pump-and-treat↗

Structural Modeling Using "Scanning and Mapping" Technique

Supported by NASA Glenn Center, we are in the process developing a structural damage diagnostic and monitoring system for rocket engines, which consists of five modules: Structural Modeling, Measurement Data Pre-Processor, Structural System Identification, Damage Detection Criterion, and Computer Visualization. The function of the system is to detect damage as it is incurred by the engine structures. The scientific principle to identify damage is to utilize the changes in the vibrational properties between the pre-damaged and post-damaged structures. The vibrational properties of the pre-damaged structure can be obtained based on an analytic computer model of the structure. Thus, as the first stage of the whole research plan, we currently focus on the first module - Structural Modeling. Three computer software packages are selected, and will be integrated for this purpose. They are PhotoModeler-Pro, AutoCAD-R14, and MSC/NASTRAN. AutoCAD is the most popular PC-CAD system currently available in the market. For our purpose, it plays like an interface to generate structural models of any particular engine parts or assembly, which is then passed to MSC/NASTRAN for extracting structural dynamic properties. Although AutoCAD is a powerful structural modeling tool, the complexity of engine components requires a further improvement in structural modeling techniques. We are working on a so-called "scanning and mapping" technique, which is a relatively new technique. The basic idea is to producing a full and accurate 3D structural model by tracing on multiple overlapping photographs taken from different angles. There is no need to input point positions, angles, distances or axes. Photographs can be taken by any types of cameras with different lenses. With the integration of such a modeling technique, the capability of structural modeling will be enhanced. The prototypes of any complex structural components will be produced by PhotoModeler first based on existing similar components, then passed to AutoCAD for modification and correction of any discrepancies seen in the Photomodeler version of the 3Dmodel. These three software packages are fully compatible. The DXF file can be used to transfer drawings among those packages. To begin this entire process, we are using a small replica of an actual engine blade as a test object. This paper introduces the accomplishment of our recent work.

Amos, Courtney L.↗

Circulating miRNA Signature Predicts Health Risks Associated with Cancer and Spaceflight

Biological risks associated with space radiation and microgravity are major concerns for long-term space travel. Through a Systems Biology approach, our previous NASA work has shown both TGF signaling pathways and miRNAs have a critical impact on defining health risks with and without space irradiation. We hypothesize that circulating microRNA (miRNA) signatures are driving microvascular disease and muscle degeneration associated with accelerating aging and will be enhanced by exposure to the space environment (radiation and microgravity). We are investigating this hypothesis with both in vivo and in vitro models to test novel antagonist therapies to these miRNA signatures as countermeasures to reduce space radiation-induced health risks. A comprehensive Systems Biology approach is utilized to examine the influence by high atomic number by high (H) atomic number (Z) and energy (E) (HZE) irradiation. To simulate low-dose exposure due to galactic cosmic rays (GCR), we used ions, energy, and doses determined by a NASA consensus formula of 7 different ions to represent GCR (referred to as GCR sim model). To similate high-dose radiation exposure due to solar particle events (SPE), we used a solar particle event (SPE) sim model which gave a total dose of 1Gy protons with energy ranges from 50MeV to 150MeV. C57BL/6 wild-type female mice were utilized for the irradiations with our established simulated microgravity model (hindlimb suspension model) and an in vitro 3D microvasculature tissue model under simulated microgravity (clinostat) conditions was also irradiated. To expand on the circulating miRNA signature determined from our preliminary data, we determined a group of conserved miRNAs which are commonly being regulated in the majority of the organs and tissues throughout the host using our established techniques. MiRNA-sequencing was done on serum (at time of sacrifice), liver, heart, and muscle (soleus muscle) tissue for all radiation groups. Additional validation of the key miRNAs was performed by droplet digital PCR (ddPCR). This revealed a key circulating miRNA signature (consisting of multiple miRNAs) impacting cardiovascular and muscular disease risk. Further in vitro experiments with CRISPR/Cas9 system to knockout the key miRNA signatures, novel self-delivering antagomirs, overexpression of the miRNAs test the functional impact of the miRNA signatures on both microvascular disease and muscle degeneration due to space irradiation. The current work has started to allow the possible development of a novel minimally invasive miRNA based radioprotector to be used as a countermeasure for space radiation. Collectively, understanding of how whole body space radiation impacts microvascular and tissue degeneration through circulating miRNAs will greatly enhance health risk prognostication and provide possible new mechanisms for protection against space radiation. This work is supported by the Translational Research Institute through NASA Cooperative Agreement NNX16AO69A (T-0404) awarded to AB.

Beheshti, Afshin↗

Explicit Cloud Representation in the Atmos 1D Climate Model for Earth and Rocky Planet Applications

1D climate models are less sophisticated than 3D global circulation models (GCMs), however their computational time is much less expensive, allowing a large number of runs in a short period of time to explore a wide parameter space. Exploring parameter space is particularly important for predicting the observable properties of exoplanets, for which few parameters are known with certainty. Therefore, 1D climate models are still very useful tools for planetary studies. In most of these 1D models, clouds are not physically represented in the atmosphere, despite having a well-known, significant impact on a planetary radiative budget. This impact is simulated by artificially raising surface albedo, in order to reproduce the observed-averaged surface temperature (i.e. 288 K for modern Earth) and a radiative balance at the top of the atmosphere. This non-physical representation of clouds, causes atmospheric long wave and shortwaves fluxes to not match observational data. Additionally, this technique represents a parameter that is highly-tuned to modern Earth’s climate, and may not be appropriate for planets that deviate from modern Earth’s climate conditions. In this paper, we present an update to the climate model within the Atmos 1­D atmospheric modeling package with a physical representation of clouds. We show that this physical representation of clouds in the atmosphere allows both long wave and shortwave fluxes to match observational data. This improvement will allow us to study the energy fluxes for a variety of cloudy rocky planets, and increase our confidence in future simulations of temperature profile and net energy balance.

Thomas Fauchez↗

Coupled Neutronics - Thermal-Hydraulics - Thermomechanics modeling of Molten Chloride Reactor

Molten Salt Reactors (MSRs) are innovative Generation-IV nuclear reactors known for their enhanced safety features, fuel economy, and fuel cycle advancements. However, designing and studying these reactors present unique safety and modeling challenges. The Department of Energy (DOE)'s Nuclear Energy Advanced Modeling and Simulation program addresses these challenges through the development of flexible multi-fidelity, multiphysics simulation tools. This study exemplifies NEAMS tools' capabilities by modeling a chloride-based micro-MSR: Griffin for neutronics simulations, Pronghorn for thermal-hydraulics simulations, and BISON for thermomechanics simulations. These tools are integrated into a comprehensive 3D multiphysics model, successfully capturing steady-state and transient reactor operations. A loss of flow accident is examined as an illustrative example of transient behavior. Additionally, this paper contributes to demonstrating the theoretical feasibility of micro-MSR concepts, which have received less attention compared to other designs. The subsequent sections provide a brief description of the MSR concept studied, an overview of the computational tools used, and present results for both steady-state operation and a loss of flow accident.

42 - ENGINEERING↗

A thermodynamically compatible splitting procedure in hyperelasticity

A material is hyperelastic if the stress tensor is obtained by variation of the stored energy function. The corresponding 3D mathematical model of hyperelasticity written in the Eulerian coordinates represents a system of 14 conservative partial differential equations submitted to stationary differential constraints. A classical approach for numerical solving of such a 3D system is a geometrical splitting: the 3D system is split into three 1D systems along each spatial direction and solved then by using a Godunov type scheme. Each 1D system has 7 independent eigenfields corresponding to contact discontinuity, longitudinal waves and shear waves. The construction of the corresponding Riemann solvers is not an easy task even in the case of isotropic solids. Indeed, for a given specific energy it is extremely difficult, if not impossible, to check its rank-one convexity which is a necessary and sufficient condition for hyperbolicity of the governing equations. In this paper, we consider a particular case where the specific energy is a sum of two terms. The first term is the hydrodynamic energy depending only on the density and the entropy, and the second term is the shear energy which is unaffected by the volume change. In this case a very simple criterion of hyperbolicity can be formulated. We propose then a new splitting procedure which allows us to find a numerical solution of each 1D system by solving successively three 1D sub-systems. Each sub-system is hyperbolic, if the full system is hyperbolic. Moreover, each sub-system has only three waves instead of seven, and the velocities of these waves are given in explicit form. The last property allows us to construct reliable Riemann solvers. Numerical 1D tests confirm the advantage of the new approach. A multi-dimensional extension of the splitting procedure is also proposed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

Stellar Activity Effects on Moist Habitable Terrestrial Atmospheres Around M Dwarfs

Transit spectroscopy of terrestrial planets around nearby M dwarfs will be a primary goal of space missions in coming decades. Three-dimensional climate modeling has shown that slow-synchronous rotating terrestrial planets may develop thick clouds at the substellar point, increasing the albedo. For M dwarfs with T(eff)>3000 K, such planets at the inner habitable zone (IHZ) have been shown to retain moist greenhouse conditions, with enhanced stratospheric water vapor ( fH2O>10(exp −3)) and low Earth-like surface temperatures. However, M dwarfs also possess strong UV activity, which may effectively photolyze stratospheric H2O. Prior modeling efforts have not included the impact of high stellar UV activity on the H2O. Here, we employ a 1D photochemical model with varied stellar UV, to assess whether H2O destruction driven by high stellar UV would affect its detectability in transmission spectroscopy. Temperature and water vapor profiles are taken from published 3D climate model simulations for an IHZ Earth-sized planet around a 3300 K M dwarf with an N2–H2O atmosphere; they serve as self-consistent input profiles for the 1D model. We explore additional chemical complexity within the 1D model by introducing other species into the atmosphere. We find that as long as the atmosphere is well-mixed up to 1 mbar, UV activity appears to not impact detectability of H2O in the transmission spectrum. The strongest H2O features occur in the James Webb Space Telescope MIRI instrument wavelength range and are comparable to the estimated systematic noise floor of ∼50 ppm.

Badhan, Mahmuda Afrin↗