SEARCH · Engineering Papers
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.
Demonstration of Teko preconditioning capability in Aria. Proposed future development work to aid preconditioner selection via simple heuristic algorithms is presented. Some highlight applications of Teko are included. These contain brief physics descriptions, solver performance, and solver convergence information details. Notably, no geometric information or otherwise sensitive information is provided.
Atmospheric convection spans a wide range of spatial and temporal scales and involves complex interactions with the surrounding dynamic and thermodynamic environment, particularly over tropical continental regions. These processes remain a major source of uncertainty in weather and climate models, including persistent biases in the diurnal cycle of convective precipitation that directly affect estimates of climate sensitivity. Addressing these challenges requires the combined use of high-resolution observations and cloud-resolving modeling frameworks. In this context, the DOE Atmospheric Radiation Measurement (ARM) program’s Large-Eddy Simulation ARM Symbiotic Simulation and Observation (LASSO) activity provides a powerful platform that pairs comprehensive observations with numerical simulations to enable process-level understanding of atmospheric convection. Within this context, this Research and Development Partnership Pilot (RDPP) project was designed to initiate and expand DOE ARM/ASR research capacity at minority-serving institutions, while advancing scientific understanding of convective processes over the Amazon rainforest. Consistent with the RDPP mission, the project emphasized partnership development, training, and workforce capacity building alongside exploratory research activities. On the scientific side, the project produced two peer-reviewed journal articles, and one manuscript currently under review (see list in section 3.1). Together, these studies combine long-term ARM observations and cloud-resolving and convection-permitting modeling to investigate the environmental controls on the shallow-to-deep convective transition during the Amazon wet season. The results demonstrate the central role of early-day moisture preconditioning and large-scale dynamical forcing in regulating isolated deep convection, provide mechanistic insight into convective evolution, and establish physically informed modeling frameworks for future sensitivity experiments. These scientific outcomes are described in sections 2.1 to 2.3 and were disseminated in 8 conference presentations (see section 3.2) and 5 invited talks (see section 3.3), reflecting broad engagement with our community. Equally important, the project achieved its RDPP capacity-building objectives (see section 2.4). A sustained research partnership was established among the University of Maryland, Baltimore County (UMBC), Morgan State University (MSU), and Howard University (HU), and extended to include collaboration with Pacific Northwest National Laboratory (PNNL). The project organized multiple multi-day training events focused on ARM data, LASSO simulations, and quantitative analysis methods, directly engaging students, postdoctoral researchers, and faculty across institutions. These activities broadened participation in ASR research and led to independent adoption of LASSO workflows by students beyond the immediate project team. Finally, the project successfully positioned the participating institutions to pursue future DOE research. Preliminary scientific results, coupled with strengthened partnerships and technical capacity, enabled the submission of follow-on proposals to DOE ASR funding opportunities. In this way, the project fulfilled the RDPP goal of seeding durable research capacity and laying the foundation for larger-scale, sustained engagement with DOE ARM and ASR programs.
Lattice gauge theory research, with its computationally intensive simulations and complex multi‑stage workflows, is well positioned to benefit from agentic AI systems. We demonstrate how such tools can support key components of lattice gauge theory research, including novel simulation code development using standard LQCD frameworks, HPC job orchestration, simulation data analysis, and expert‑guided tuning of algorithmic parameters such as Hasenbusch mass preconditioning and multigrid solvers. Our results show that agentic AI can reduce manual effort, improve productivity, and accelerate the research cycle while maintaining essential human oversight.
This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.
This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.
Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.
Abstract We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass (10 3 M ⊙ ) black holes. The length scales required to resolve the initiation mechanism are not easily reached in 3D, so instead we have devised 2D proxy models, which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 M ⊙ white dwarf stars parameterized by tidal strengths in the range β = 4–23. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find that many diagnostics, such as temperature, total released energy, and iron-group products, are fairly well converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate-mass transients like calcium, which remain uncertain up to factors of 2, even at 1 km resolution.
Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.
Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.
As the perovskite technology is ramping up into commercialization, reliable and accurate power rating of large-size perovskite modules becomes a prominent aspect for its future deployment in the PV market. It is known that the performance calibration of perovskite PV devices is very challenging due to its complex dynamic response during a conventional current-voltage (IV) measurement. PV researchers have previously proposed several steady-state performance calibration methods to reliably extract PV efficiencies, but mostly focus on small area research-type cells. In this paper, we emphasize the importance of reliable performance calibration on large-size perovskite modules. Extending the NREL Cell and Module Performance (CMP) group’s steady-state performance calibration protocol (i.e., Asymptotic PMAX Scan) for perovskite cells to modules, we justify the necessity of reporting steady-state efficiencies for perovskite cells and discuss the challenges of applying this protocol to modules. We also present our protocol for Maximum Power Point Tracking (MPPT), which is a technique often used for performance calibration of perovskite cells and modules, and show a comparison between MPPT and Asymptotic PMAX. Using MPPT we demonstrate the interplay between metastability and degradation in perovskite modules, and emphasize the necessity to develop preconditioning protocols for stabilizing these devices. Our aim is to promote development of consensus protocols for performance calibration of perovskite modules, and to advance their credible power ratings, which will be beneficial to the growth of perovskite technology in the PV market.
Gauss-Seidel (GS) relaxation is often employed as a preconditioner for a Krylov solver or as a smoother for Algebraic Multigrid (AMG). However, the requisite sparse triangular solve is difficult to parallelize on many-core architectures such as graphics processing units (GPUs). In the present study, the performance of the sequential GS relaxation based on a triangular solve is compared with two-stage variants, replacing the direct triangular solve with a fixed number of inner Jacobi-Richardson (JR) iterations. When a small number of inner iterations is sufficient to maintain the Krylov convergence rate, the two-stage GS (GS2) often outperforms the sequential algorithm on many-core architectures. The GS2 algorithm is also compared with JR. When they perform the same number of ops for SpMV (e.g. three JR sweeps compared to two GS sweeps with one inner JR sweep), the GS2 iterations, and the Krylov solver preconditioned with GS2, may converge faster than the JR iterations. Moreover, for some problems (e.g. elasticity), it was found that JR may diverge with a damping factor of one, whereas two-stage GS may improve the convergence with more inner iterations. Finally, to study the performance of the two-stage smoother and preconditioner for a practical problem, these were applied to incompressible uid ow simulations on GPUs.
Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.
As part of the U.S. DOE EVs@Scale consortium Next-Gen Profiles (NGP) project, the profile capture of production electric vehicles undergoing high power charging (HPC) is conducted over a wide range of conditions to explore variance and performance. Charge session parameters are collected from both the electric vehicle (EV) and electric vehicle supply equipment (EVSE) at a rate of 10Hz and entered into a time-series database for analysis. These charge profiles are captured under nominal and off-nominal conditions, exploring the impact of battery state of charge (SOC), battery temperature, vehicle condition, smart charge management (SCM), and EVSE limitations. Nominal conditions are defined to be ideal conditions that should transfer the maximum allowable energy in the minimum possible amount of time. Nominal condition profiles are compared across EVs to characterize state-of-the-art EV charging performance against one another. Off-nominal condition profiles are compared against its nominal condition profile counterpart to highlight the variance across less desirable starting conditions within a single EV. Under nominal conditions, most EVs can achieve the original equipment manufacturer (OEM) rated peak-power and charge times. Peak power across EVs has variance due to the vehicle design, battery topology, charging strategy, etc. These 10-100% nominal preconditioned charge profiles across 13 EVs (11 light-duty (LD), 2 heavy-duty (HD)) were used for analysis in power curve analysis, power distribution, SOC and range comparison, thermal impacts of current draw, battery pack size and energy charged, and ramp rates. Power curve analysis show the uniqueness in power vs time curves across all EVs, breaking down the features of a typical profile and where variation is typically seen. Power distribution results showed that over 50% of charge time is spent below 50kW and only 12.1% is spent above 200kW. SOC and range performance yielded different top performing EVs when exploring goals of performance from SOC and range gained after 10-min (EV2) and 20-min (EV8), and time-to-achieve 80% SOC (EV8) and 200 miles of range (EV1). Current draw from 400-volt EVs had a higher thermal impact to cable/connector temperatures when compared to 800-volt EVs, due to a higher current requirement to achieve similar power levels. There was high variance in C rating, a useful metric when comparing the relationship between peak/average charge session power and relative battery pack size, across EVs under test. Ramp rates during initial power transfer were examined, fastest and slowest speeds ranging from 192.5kW/second to 2.6kW/second respectively. A similar analysis of ramp rates was also conducted for OCPP curtailment testing, where EVs underwent a 2-minute 65A curtailment request before returning to full-power charge. Under off-nominal conditions, most EVs experienced variation from the nominal condition profiles. EVs that underwent the full set of NGP defined testing conditions were compared, analysis of which was categorized by 800-volt and 400-volt EVs. It should be noted that the 800-volt EVs had significantly higher peak power ratings than the 400-volt EVs, and thus were more prone to variance. Initial state of charge, battery temperature, and vehicle condition had a considerable impact on peak power levels and charge time for the 800-volt vehicles under test. EVSE limited tests lowered achievable power down to 200kW, greatly impacted 800-volt charging power, and had little to no effect on 400-volt charging power. Adapter testing was performed on a single 400-volt EV, where power and current limitations were found but overall charge time was not significantly impacted. Adapter and boost converter testing was performed on a single 800-volt EV, where both charge power and charge time were greatly impacted. Charge profiles are unique, and comparing such requires analysis that examines starting conditions, vehicle and battery topologies,
The treatment of radioactive wastewater from the operating nuclear industry (energy production, reprocessing of spent fuel, decommissioning of shut-down facilities, etc.) or in case of incidental or accidental situation (such as Fukushima Daichii NPP accident) is an important issue requiring constant improvement of the technologies and techniques used to prevent radioactive releases into the environment and/or to reduce the production of secondary waste. CEA, Orano and CTI are collaborating since 2014, with the support of the French government's program 'investments for the future', to develop Sorbmatech{sup R}Cs, an innovative and patented selective sorbent (Grandjean and Barre 2014) to deal with the removal of radioactive {sup 137}Cs, one of the most hazardous isotope to be managed in nuclear waste, highly efficient either in agitated reactor or in column process. Sorbmatech{sup R}Cs sorbent consists of mesoporous silica particles loaded with potassium-copper hexacyanoferrate nanoparticles (KCu-HCF). This peculiar structure enables much faster sorption kinetics compared to denser crystalline structures, such as Silicotitanate or Zeolite sorbents, whose drawback consists in hindering the diffusion process through the material, and whose larger sorbent particles slow down the ion exchange process. The equilibrium of Sorbmatech{sup R}Cs is reached in few minutes, even at high flow rates in a fixed-bed column process (up to 10 m.h{sup -1} Darcy velocity) with low drop pressure. Moreover, Sorbmatech{sup R}Cs can be used over a wide pH range from 2 to 10 and has a very high selectivity towards {sup 137}Cs even in brine solution (saline, mineral or sea water). This new sorbent also has a K{sub d} higher than 10{sup 5} mL.g{sup -1} depending on the salinity of the wastewater, enabling a highly efficient removal of Cs trace in radioactive wastewater compared to Cs selective sorbent such as Zeolites and Silicotitanates. Its radiological resistance enables the treatment of high, medium and low activity wastewater. The CEA-Orano-CTI/ALSYS partnership leads to the development of a simple and low-cost synthetic route, which demonstrated large-scale manufacturing of the product on a batch of 100 kg. Two products have been developed and are ready for short run market introduction so far with grain size well adapted either for column processes or agitated reactors.. Finally, the industrial implementation has been studied and largely benefits from the properties of the sorbent. Indeed, its high sorption kinetics enable high flow rates resulting in a decrease of operation times, unit sizes, and volumes of secondary waste produced. Moreover, Sorbmatech{sup R}Cs is designed to be combined into a full and integrated solution, ranging from preliminary optimization of operating conditions through numerical dimensioning model, industrial implementation and treatment of the wastewater (through cartridge or preconditioned columns systems developed thanks to Orano's extensive expertise), to long-term stabilization of generated waste (under progress). (authors)
A decoupled finite element algorithm is developed for simulating the vortex dynamics on an elastic superconductor which couples the time-dependent Ginzburg- Landau equation with the complex-valued superconducting order parameter and the vector-valued magnetic potential, and the elasticity equation. We present an iterative algorithm for the decoupled system arising from the time and spatial discretization using a combination of preconditioner, algebraic multigrid method (AMG) and preconditioned conjugate gradient method (PCG). The iterative algorithm allows us to perform large-scale three-dimensional simulations of mesoscale pattern formation during superconducting phase transitions with arbitrary elastic boundary conditions. Here, the performance and efficiency of the algorithm are numerically verified by several benchmark problems, exhibiting up to two orders of magnitude improvement depending on the scale of discrete system compared to the exact solver.
Upon application of a sufficiently strong electric field, electrons break away from thermal equilibrium and approach relativistic speeds. These highly energetic ‘runaway’ electrons (~ MeV) play a significant role in tokamak disruption physics, and therefore their accurate understanding is essential to develop reliable mitigation strategies. As such, we have developed a fully implicit solver for the 0D-2P (i.e., including two momenta coordinates) relativistic nonlinear Fokker–Planck equation (rFP). As in earlier implicit rFP studies (NORSE, CQL3D), electron–ion interactions are modeled using the Lorentz operator, and synchrotron damping using the Abraham–Lorentz–Dirac reaction term. However, our implementation improves on these earlier studies by (1) ensuring exact conservation properties for electron collisions, (2) strictly preserving positivity, and (3) being scalable algorithmically and in parallel. Key to our proposed approach is an efficient multigrid preconditioner for the linearized rFP equation, a multigrid elliptic solver for the Braams–Karney potentials, and a novel adaptive technique to determine the associated boundary values. We verify the accuracy and efficiency of the proposed scheme with numerical results ranging from small electric-field electrical conductivity measurements to the accurate reproduction of runaway tail dynamics when strong electric fields are applied.
We focus on the fully implicit solution of the linear systems arising from a three-field mixed finite element approximation of Biot’s poroleasticity equations. The objective is to develop algebraic block preconditioners for the efficient solution of such systems by Krylov subspace methods. In this work, we investigate the use of approximate inverse-based techniques to decouple the native system of equations and obtain explicit sparse approximations of the Schur complements related to the physics-based partitioning of the unknowns by field type. Here, the proposed methods are tested in various numerical experiments including real-world applications dealing with petroleum and geotechnical engineering.