Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mathematics and Computing, Geosciences”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Inference of Rock Flow and Mechanical Properties from Injection-Induced Microseismic Events During Geologic CO 2 Storage

Monitoring microseismic activities during CO 2 injection into geologic formations is important for ensuring the safety of the storage operations. The resulting data provide insight into the response of the storage formation to CO 2 injection and can be used to infer the underlying rock flow and mechanical properties. In this paper, assimilation of microseismic data is performed for dynamic characterization of the storage formation by using a stochastic simulation model to forecast the microseismic response of a geologic formation during CO 2 injection. Two modeling approaches are adopted to predict the space-time distribution of the injection-induced microseismicity. The first model is based on pore pressure relaxation assumption, while the second model uses coupled flow and geomechanics simulation to establish the complex physical relation between the storage formation properties and the corresponding microseismic responses during CO 2 injection. The stochastic predictive models in each case are used in ensemble data assimilation frameworks to estimate rock properties from the observed microseismic data. Two data assimilation methods are considered: (i) a new ensemble-based stochastic point process filter (EnPPF) that can directly integrate discrete microseismic events, and (ii) a variant of ensemble smoother, known as the ensemble smoother with multiple data assimilation (ES-MDA), which requires continuous representation of microseismic events for assimilation. The two methods are successfully applied to a geologically realistic model of the Farnsworth Field in Texas, with complex geologic flow units and interacting fault systems.

42 ENGINEERING↗
DOE OSTI2024

Advancing spatiotemporal forecasts of CO 2 plume migration using deep learning networks with transfer learning and interpretation analysis

Accurate and timely forecasts of CO 2 plume distribution throughout the injection and post-injection phases are crucial for detecting plume migration, assessing leakage risks, and supporting operational decisions in geologic carbon storage (GCS). Current convolutional neural network-based approaches primarily focus on spatial information and overlook temporal dependencies in plume distributions, thus limiting their ability to capture dynamic movement effects and provide accurate predictions of plume migration. In this work, we propose two deep learning models, Auto-Encoder (AE)-LSTM and Encoder-Decoder (ED)-ConvLSTM, each uniquely designed to capture both spatial and temporal features. We apply the proposed methods to forecast the dynamic distribution of CO 2 plumes based on 108 reservoir simulations over a 30-year injection and a 30-year post-injection period. The results indicate that the ED-ConvLSTM model outperforms the AE-LSTM model in accurately predicting the spatiotemporal dynamics of CO 2 plume migration, achieving R 2 values above 0.99. To provide a deeper understanding of these model predictions, we employ a gradient-based explanation method on the trained models. This approach provides insights into the influence of input variables on plume migration forecasts and uncovers the underlying prediction mechanisms of the proposed models. Furthermore, we introduce a transfer learning technique, enabling fast and accurate plume migration forecasting in the post-injection phase by leveraging the trained model during the injection phase. This reduces the necessity for extensive data collection or re-training. In conclusion, the methods proposed in our work enhances the performance and interpretability of CO 2 plume migration forecasts, thereby facilitating informed decision-making throughout the entire lifecycle of GCS applications.

58 GEOSCIENCES↗
DOE OSTI2023

Comparative verification of hydro-mechanical fracture behavior: Task G of international research project DECOVALEX–2023

Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, in this study, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.

58 GEOSCIENCES↗
DOE OSTI2020

Enhanced multiscale restriction-smoothed basis (MsRSB) preconditioning with applications to porous media flow and geomechanics

A novel method to enable application of the Multiscale Restricted Smoothed Basis (MsRSB) method to non M-matrices is presented. The original MsRSB method is enhanced with a filtering strategy enforcing M-matrix properties to enable the robust application of MsRSB as a preconditioner. Through applications to porous media flow and linear elastic geomechanics, the method is proven to be effective for scalar and vector problems with multipoint finite volume (FV) and finite element (FE) discretization schemes, respectively. Realistic complex (un)structured two- and three-dimensional test cases are considered to illustrate the method's performance.

58 GEOSCIENCES↗
DOE OSTI2022

Scalable preconditioning for the stabilized contact mechanics problem

We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.

58 GEOSCIENCES↗
DOE OSTI2024

Multilevel well modeling in aggregation-based nonlinear multigrid for multiphase flow in porous media

A full approximation scheme (FAS) nonlinear multigrid solver for two-phase flow and transport problems driven by wells with multiple perforations is developed here. It is an extension to our previous work on FAS solvers for diffusion and transport problems. The solver is applicable to discrete problems defined on unstructured grids as the coarsening algorithm is aggregation-based and algebraic. To construct coarse basis that can better capture the radial flow near wells, coarse grids in which perforated well cells are not near the coarse-element interface are desired. This is achieved by an aggregation algorithm proposed in this paper that makes use of the location of well cells in the cell-connectivity graph. Numerical examples in which the FAS solver is compared against Newton's method on benchmark problems are given. In particular, for a refined version of the SAIGUP model, the FAS solver is at least 35% faster than Newton's method for time steps with a CFL number greater than 10.

58 GEOSCIENCES↗
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Intercomparison of Thermal Regime Algorithms in 1-D Lake Models

Lakes are an important component of the global weather and climate system, but the modeling of their thermal regimes has shown large uncertainties due to the highly diverse lake properties and model configurations. Here we evaluate the algorithms of four key lake thermal processes including turbulent heat fluxes, wind-driven mixing, light extinction, and snow density, using a highly diverse lake dataset provided by the Inter-sectoral Impact Model Intercomparison Project (ISIMIP) 2a lake sector. Algorithm codes are configured and run separately within the same parent model to rule out any interference from factors apart from the algorithms examined. Evaluations are based on both simulation accuracy and recalibration complexity for application to global lakes. For turbulent heat fluxes, the non-Monin-Obukhov similarity (MOS) based, more simplified algorithms perform better in predicting lake epilimnion temperatures and achieve high convergence in the values of the calibrated parameters. For wind-driven mixing, a two-algorithm strategy considering lake shape and season is suggested with the regular mixing algorithm used for spring and earlier summer and the mixing-enhanced algorithm for summer steady stratification and fall overturn periods. There are no evident differences in the simulated thermocline depths using different light extinction algorithms or the observation. Finally, for lake ice phenology, a constant snow density at around 110 kg m-3 is found to be sufficient for most northern lakes while the Arctic lakes require a higher value. Our study provides highly practical guides for improving 1-D lake models and feasible parameterization strategies to better simulate global lake thermal regimes.

58 GEOSCIENCES↗
DOE OSTI2021

Enabling the Application of Large Footprint Open-Bottom Permeameters through New Shape Factors

Open-bottom permeameters (OBP) are commonly used for estimating the hydraulic properties of sediment beds under water bodies and fluxes across the sediment-water interface. The specific shape factor for an OBP flow situation is necessary for data interpretation. It is determined by the ratio of the permeameter tube radius to the tube's penetration depth within the sediment. Previous permeameter designs favored small radii which allowed for deeper penetration and stability of the tube. The shape factor is provided here for large footprint OBPs that allows for shallow penetration and greater radius. Through a numerical and analytical modeling study of the OBP hydrodynamics, previous knowledge on the OBP shape factor is generalized to the entire range of radius-to-penetration ratios that might be practically applied. The results reveal the advantages of the large footprint design, including better sensitivity due to the ability to impose larger discharge over a shorter test duration.

42 ENGINEERING↗
DOE OSTI2022

Improving Seasonal Forecast Using Probabilistic Deep Learning

The path toward realizing the potential of seasonal forecasting and its socioeconomic benefits relies on improving general circulation model (GCM) based dynamical forecast systems. To improve dynamical seasonal forecasts, it is crucial to set up forecast benchmarks, and clarify forecast limitations posed by model initialization errors, formulation deficiencies, and internal climate variability. With huge costs in generating large forecast ensembles, and limited observations for forecast verification, the seasonal forecast benchmarking and diagnosing task proves challenging. Here, we develop a probabilistic deep learning-based statistical forecast methodology, drawing on a wealth of climate simulations to enhance seasonal forecast capability and forecast diagnosis. By explicitly modeling the internal climate variability and GCM formulation differences, the proposed Conditional Generative Forecasting (CGF) methodology enables bypassing crucial barriers in dynamical forecast, and offers a top-down viewpoint to examine how complicated GCMs encode the seasonal predictability information. We apply the CGF methodology for global seasonal forecast of precipitation and 2 m air temperature, based on a unique data set consisting 52,201 years of climate simulation. Results show that the CGF methodology can faithfully represent the seasonal predictability information encoded in GCMs. We successfully apply this learned relationship in real-world seasonal forecast, achieving competitive performance compared to dynamical forecasts. Using this CGF as benchmark, we reveal the impact of insufficient forecast spread sampling that limits the skill of the considered dynamical forecast system. Finally, we introduce different strategies for composing ensembles using the CGF methodology, highlighting the potential for leveraging the strengths of multiple GCMs to achieve advantgeous seasonal forecast.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Combining Deep Learning With Physics Based Features in Explosion-Earthquake Discrimination

This paper combines the power of deep-learning with the generalizability of physics-based features, to present an advanced method for seismic discrimination between earthquakes and explosions. The proposed method contains two branches: a deep learning branch operating directly on seismic waveforms or spectrograms, and a second branch operating on physics-based parametric features. These features are high-frequency P/S amplitude ratios and the difference between local magnitude (ML) and coda duration magnitude (MC). The combination achieves better generalization performance when applied to new regions than models that are developed solely with deep learning. Further, we also examined which parts of the waveform data dominate deep learning decisions (i.e., via Grad-CAM). Such visualization provides a window into the black-box nature of the machine-learning models and offers new insight into how the deep learning derived models use data to make decisions.

58 GEOSCIENCES↗
DOE OSTI2023

Simulating the Transport and Rupture of Pollen in the Atmosphere

Pollen, one type of primary biological aerosol particle (PBAP), is emitted from the terrestrial biosphere and can undergo physical changes in the atmosphere via particle rupture. To examine the fate of pollen and its atmospheric processing, a pollen emission and transport scheme is coupled to the Weather Research and Forecasting Model with Chemistry (WRF-Chem). We simulate the emission of pollen and its impacts on the cloud properties and precipitation in the Southern Great Plains from 12 to 19 April 2013, a period with both high pollen emissions and convective activity. We conduct a suite of ensemble runs that simulate primary pollen and three different pollen rupture mechanisms that generate subpollen particles, including (a) high humidity-induced surface rupture, (b) high humidity-induced in-atmosphere plus surface rupture, and (c) lightning-induced rupture, where in-cloud and cloud-to-ground lightning strikes trigger pollen rupture events. When relative humidity is high (>80%), coarse primary pollen (~1 μg m –3 ) is converted into fine subpollen particles (~1.2e –4 μg m –3 ), which produces 80% more subpollen particles than lightning-induced rupture. The in-atmosphere humidity-driven rupture predominantly produces subpollen particles, which is further enhanced during a frontal thunderstorm. During strong convection, vertical updrafts lift primary pollen and subpollen particles (~0.5e –4 μg m –3 ) to the upper troposphere (~12 km) and laterally transports the ruptured pollen in the anvil top outflow. In regions of high pollen and strong convection, ruptured pollen can influence warm cloud formation by decreasing low cloud (<4 km) cloud water mixing ratios and increasing ice phase hydrometeors aloft (>10 km)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Influence of Climate Variability and Future Climate Change on Atlantic Hurricane Season Length

Abstract Atlantic hurricane season length is important for emergency management preparation, motivating the need to understand its variability and change. We investigated the influence of ocean variability on Atlantic hurricane season length in observations and a future climate simulated by the Energy Exascale Earth System Model (E3SM). We found that multiple factors influence hurricane season length, through their influence on season start and end. Warm western subtropical Atlantic sea‐surface temperature anomalies (SSTAs) during boreal spring (before the official hurricane season start) drive early starts to the hurricane season, and vice versa for cool SSTAs. Meanwhile, La Niña in autumn (before the official hurricane season end) drives late ends to the hurricane season, and vice versa for El Niño. E3SM projects a 27‐day increase in future Atlantic hurricane season length given La Niña and warm northern tropical Atlantic SSTAs. This research documents sources of predictability for Atlantic hurricane season length.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Effects of Pollen on Hydrometeors and Precipitation in a Convective System

Anemophilous (wind-driven) pollen is one type of primary biological aerosol particle, which can rupture under high humidity conditions and form smaller sub-pollen particles (SPPs). Both pollen and SPPs can reach the upper troposphere under convective conditions, acting as cloud condensation nuclei (CCN) and ice nucleating particles (INPs), thus influencing cloud formation and precipitation. However, the impacts of these biological aerosols on cold cloud formation and local climate remain unclear as there are large uncertainties on their emission flux and ice nucleating abilities. Here, we incorporate pollen emission and rupture processes in the Weather Research and Forecasting Model with Chemistry (WRF-Chem) simulations and update the Morrison microphysics scheme within WRF-Chem using aerosol-aware INP parameterizations to account for pollen in addition to other anthropogenic and biogenic aerosol. INP parameterizations for pollen and SPP are derived from laboratory experiments. When including pollen rupture rates as observed in a series of chamber studies, SPP concentrations increase, leading to an increase of cloud ice and water by up to 50% and potentially extending the duration of the convective system. Among all simulated hydrometeors, graupel and raindrops exhibit the largest enhancements from the inclusion of SPPs, with intensifying precipitation at the backside of the convective system and a greater spatial extent. Sensitivity simulations indicate that SPPs have a greater effect on cloud microphysical processes than whole pollen grains, and further observational evidence is needed to constrain these processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗
DOE OSTI2026

Antarctic Meltwater Accelerates Southern Ocean Evolution Under Projected Atmospheric Warming

Increasing basal meltwater from Antarctic ice shelves may impact the Southern Ocean properties that feed back on the rate of melting. We investigate this feedback in a high‐emissions scenario using an Earth‐system model with interactive ice‐shelf basal melting, an improvement on previous studies that did not have the capability to evolve melt rates and the ocean state self‐consistently. We find that when interactive melt increases, it primarily accelerates the evolution of a spatial pattern of continental shelf warming and cooling that is initiated by freshening and sea‐ice formation decline due to projected atmospheric warming. The competition between enhanced warming at depth from reduced ventilation and enhanced continental shelf cooling from reduced dense water export leads to net ~35% reduction in ice‐shelf meltwater input into the Southern Ocean over the 21st century. Omitting this feedback introduces a bias in the timing of projected ocean‐melt‐driven ice loss from Antarctica.

58 GEOSCIENCES↗