Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “metric geometry”

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 163 records · Page 9

Network Uncertainty Quantification for Analysis of Multi-Component Systems

To impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.

42 ENGINEERING↗

Preliminary neutronic and thermal analysis of Inertial Fusion Energy systems with a focus on thick liquid wall concepts

This paper provides an assessment of the neutronic behavior of Inertial Fusion Energy systems using the Serpent 2 Monte Carlo code. The study explores different types of reactor concepts, including dry-wall and thick liquid wall systems. In particular, it characterizes the tritium breeding and energy multiplication performance, as well as the effect of long-term irradiation on structural materials. For this purpose, the displacement per atom (DPA) and the gas production rate in structural materials are calculated. A simplified geometry with varying geometrical parameters is used as a representation for the different kinds of devices. Additionally, the case of the HYLIFE-II concept is used to assess modeling inaccuracies arising from certain approximations. A comparison between the simplified geometry and a realistic one is performed, thanks to Serpent’s capability to run transport calculations on unstructured meshes. The study then addresses the impact of these discrepancies on the temperature distribution in the first wall using the GeN-Foam code for conjugate heat transfer analysis. The results highlight the necessity to consider the heterogeneous nature of the liquid wall, as simplified models tend to underestimate radiation damage metrics and heat loads.

Conjugate heat transfer↗

Development and validation of the cavitation-induced erosion risk assessment tool

This work presents the development of a cavitation-induced erosion risk assessment (CIERA) tool that links multiphase flow simulation predictions with the progress towards material erosion. To develop a robust erosion modeling tool, the cavitation and erosion predictions for pressurized diesel fuel flow within a channel geometry were validated over a range of Reynolds and cavitation number conditions in two different aluminum channel geometries, one featuring a rounded inlet corner and the other with a sharp inlet corner. The multiphase flow development within the channel was modeled using a compressible mixture model, where phase change was represented with the homogeneous relaxation model and the turbulent flow evolution was modeled using a dynamic structure approach for Large Eddy Simulations. To improve representation of the incubation period before material rupture over existing approaches, a physics-based metric was derived based on the cumulative energy absorbed by the solid material from repeated hydrodynamic impacts. When the average peak pressure was related to the incubation period, the incubation period and its sensitivity to changes in flow conditions were found to be overpredicted. In contrast, predictions from CIERA provided a more accurate means to qualitatively and quantitatively predict the influence of flow conditions on the incubation period before material erosion. When the predicted stored energy was related to the solid material properties to estimate the incubation period, multiphase flow simulations demonstrated accurate representation of the sensitivity of erosion severity to changes in flow conditions. The use of CIERA led to quantitative agreement of the predicted incubation period within 5% of the experimentally measured incubation period.

42 ENGINEERING↗

Modeling of Supercritical CO2 Shell-and-Tube Heat Exchangers Under Extreme Conditions: Part II: Heat Exchanger Model

Abstract Heat exchangers play a critical role in supercritical CO2 Brayton cycles by providing necessary waste heat recovery. Supercritical CO2 thermal cycles potentially achieve higher energy density and thermal efficiency operating at elevated temperatures and pressures. Accurate and computationally efficient estimation of heat exchanger performance metrics at these conditions is important for the design and optimization of sCO2 systems and thermal cycles. In this paper (Part II), a computationally efficient and accurate numerical model is developed to predict the performance of shell-and-tube heat exchangers (STHXs). Highly accurate correlations reported in Part I of this study are utilized to improve the accuracy of performance predictions, and the concept of volume averaging is used to abstract the geometry and reduce computation time. The numerical model is validated by comparison with computational fluid dynamics (CFD) simulations and provides high accuracy and significantly lower computation time compared to existing numerical models. A preliminary optimization study is conducted and the advantage of using supercritical CO2 as a working fluid for energy systems is demonstrated.

Engineering↗

Probing electromagnetic nonreciprocity with quantum geometry of photonic states

Reciprocal and nonreciprocal effects in dielectric and magnetic materials provide crucial information about the microscopic properties of electrons. However, experimentally distinguishing the two has proven to be challenging, especially when the associated effects are extremely small. To this end, we propose a contactless detection using a cross-cavity device where a material of interest is placed at its center. We show that the optical properties of the material, such as Kerr and Faraday rotation, or birefringence, manifest in the coupling between the cavity's electromagnetic modes and in the shift of their resonant frequencies. By calculating the dynamics of a geometrical photonic state, we formulate a measurement protocol based on the quantum metric and quantum process tomography that isolates the individual components of the material's complex refractive index and minimizes the quantum mechanical Cramér-Rao bound on the variance of the associated parameter estimation. Our approach is expected to be applicable across a broad spectrum of experimental platforms including Fock states in optical cavities, or coherent states in microwave and THz resonators. Published by the American Physical Society 2025

Petrides, Ioannis (ORCID:0000000174842444)↗

Welding of Haynes 282 to Steels to Enable Modular Rotors for Advanced Ultra Super-Critical Steam Turbines

Steam turbines for an Advanced Ultra Super Critical (AUSC) fossil fired power plant will operate at temperatures well above those of current commercial steam cycles with inlet temperatures more than 760 C , which is beyond the capabilities of alloy steels presently used for steam turbine applications and requires advanced materials. Large components such as steam turbine rotors may be made of nickel based super alloys such as Haynes 282 (H282). However, monolithic forgings of superalloys in the sizes required for large steam turbines rotors can be prohibitively expensive besides many technical challenges. To minimize cost and alleviate the related technical challenges, superalloy use needs to be limited to locations on the steam turbine rotor where strength and temperature requirements cannot be met by conventional steels. This is possible if nickel-based superalloys can be successfully welded to steels and the related technical challenges - machining parts made of dissimilar welded materials, non-destructive examination of such welds for flaw detection; and, material properties of such hybrid components – are sufficiently addressed. In this technology development project, we successfully welded H282 to plates up to ~ 75 mm (~ 3 inches) to a 3.5NiCrMoV steel of similar thickness. Advanced ultrasonic inspection technique called Phased Array Ultrasonic Testing (PAUT) was employed to examine the dissimilar H282-Steel welds, into which flat bottomed side drilled holes (SDH) of various diameters were introduced, to determine the minimum detectable feature sizes; it was shown that with PAUT SDH of dia. down to 0.5 mm could be detected in the base alloys and SDH with dia. down to 2.4 mm could be detected in the weld metal under multiple orientations successfully. An autonomous machining process monitoring system was developed and demonstrated whereby the forces acting on the cutting tool could be actively monitored as the cutting tool transitioned from H282 to Steel across the weld using which the machining parameters can be potentially altered without interruption to extend tool life. This project successfully achieved its objectives of - i. Developing a welding methodology and viable welding geometries, to successfully join H282 to steel 3.5CrMoNiV steel to enable manufacture of modular steam turbine rotors for AUSC applications (conditions of at least 760 °Celsius and 3,100 psia (pounds per square inch absolute pressure) and evaluate the material properties of the welded specimen. ii. Employ the advanced ultrasonic inspection technique to the dissimilar weld metal joint and determine the minimum detectable feature sizes. iii. Develop effective machining techniques to machine such hybrid structures with online tool force monitoring and effect machine state metrics for optimal results.

20 FOSSIL-FUELED POWER PLANTS↗

New superstrata from three-dimensional supergravity

We find a two-parameter family of generalized superstrata that emerge as smooth, supersymmetric solutions in three-dimensional gauged supergravity coupled to additional scalar fields. This new family of generalized superstrata are smooth microstate geometries and may be thought of as supersymmetric Coulomb-branch extensions of the original superstrata in which the underlying supertube undergoes an elliptical deformation. These solutions had already been obtained numerically, and as series solutions, to the equations of motion, and some of them were conjectured to be supersymmetric. Here we prove the supersymmetry of an entire two-parameter family and we obtain a highly non-trivial analytic and smooth solution for a one-parameter limit in which the global symmetry of the metric is enhanced to SO(3). We also confirm that the other known families of microstrata are not supersymmetric. We conclude with a cursory analysis of some of the singular brane distributions that can be accessed from three-dimensional gauged supergravity while preserving the same supersymmetries as the superstratum, and therefore of the three-charge black hole.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Convolutional double copy in (anti) de Sitter space

The double copy is a remarkable relationship between gauge theory and gravity that has been explored in a number of contexts, most notably scattering amplitudes and classical solutions. The convolutional double copy provides a straightforward method to bridge the two theories via a precise map for the fields and symmetries at the linearised level. This method has been thoroughly investigated in flat space, offering a comprehensive dictionary both with and without fixing the gauge degrees of freedom. In this paper, we extend this to curved space with an (anti) de Sitter background metric. We work in the temporal gauge, and employ a modified convolution that involves the Mellin transformation in the time direction. As an example, we show that the point-like charge in gauge theory double copies to the (dS-) Schwarzschild black hole solution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electroslag additive manufacturing: A pathway for high throughput near net shape production

Electroslag Additive Manufacturing (ESAM), a new high-throughput additive manufacturing (AM) method that combines Electroslag Strip Cladding (ESC) and wire arc AM (WAAM) is introduced. This combination enables the high deposition rate of ESC (more than 20 kg/h with a 60 mm strip electrode) to benefit from the precise geometric control of WAAM. As a precursor to ESAM, the ESC process is investigated in an AM context independently by evaluating both direct and staggered bead-stacking strategies and analyzing the microstructural and mechanical properties of each. This is followed by an ESAM demonstration producing an annular geometry by pairing ESC with gas tungsten arc welding (GTAW), wherein GTAW is utilized to construct annular walls that are subsequently infilled via ESC. The microstructure and mechanical properties of ESC-only AM are compared with that of the ESAM method and it is shown that printed integral retaining walls do not impact the resulting mechanical properties of ESAM. Furthermore, results indicate that ESAM-produced Alloy 625 parts exhibit tensile properties on par with cast counterparts, supporting the method’s scalability to components exceeding one metric ton, and possibly making ESAM a viable future manufacturing approach for competitive production of large-scale components currently manufactured by casting and forging.

Additive manufacturing↗

Data, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗

Optimization of the Second Target Station cold source moderators using an automated workflow

The Second Target Station (STS) at the US Department of Energy’s Oak Ridge National Laboratory is designed to become the world’s highest peak-brightness spallation source of cold neutrons. Successful completion of the STS, which is currently in the preliminary design phase, will provide transformative new capabilities to examine novel materials for future technologies. At STS, neutrons will be generated by spallation reactions in a solid tungsten target. They will be moderated and thermalized in two cold (20 K) para-hydrogen moderators. Careful optimization of these moderators is essential to the project’s success. To find optimal moderator designs, an advanced optimization workflow integrates high-fidelity neutronics calculations using the Monte Carlo N-Particle (MCNP) transport code MCNP6.2 with state-of-the-art optimization algorithms in the Dakota optimization toolkit. For each design iteration, a parametrized solid CAD geometry is generated in Creo and automatically converted into an unstructured mesh geometry by Attila 4MC for the neutronics calculation with MCNP. Iterations repeat until optimal designs are found. Herein this paper presents the results of a sensitivity and optimization study for the cylindrical and tube moderators. Both moderators can be optimized for maximum peak brightness, maximum time-integrated brightness, or any combination between these extremes. Maximum peak brightness is achieved by using smaller optimal dimensions of the moderators, whereas maximum time-integrated brightness is achieved by using larger dimensions. A Pareto front details the designs that optimally balance both brightness metrics. The Pareto front can be found in only 40–110 iterations with 4–5 design parameters when using the efficient global and Pareto-set optimization algorithms in Dakota. Additionally, important engineering constraints can be taken into account, such as the coupling between the cylindrical moderator radius and aluminum vessel wall thicknesses required to ensure structural integrity of the vessels. This interaction has a significant impact on the resulting optimal designs. Our new, highly efficient, fully automated optimization workflow will be used to optimize additional STS components in the future and can be adopted for design and optimization studies at other experimental neutron and accelerator facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Mesh-based multiphysics coupling acceleration for fusion neutronics through clustering for fusion blanket applications

Accurate modeling of particle transport within fusion blankets is essential for predicting performance metrics such as heat deposition and the tritium breeding ratio (TBR). However, high-fidelity coupling of thermal fluids from computational fluid dynamics (CFD) to neutronics simulations often incurs significant computational costs due to the complexity of surface intersection calculations in Monte Carlo codes. This paper presents an accelerated multiphysics coupling method for neutronics that utilizes hierarchical agglomerative clustering to map complex material property distributions to a neutronics model. Implemented within the fusion reactor design and assessment (FREDA) framework, the method leverages existing Python packages to automate the creation of clustered geometries for OpenMC. The approach is demonstrated on a sector model of an ARC-class tokamak with an immersion molten salt blanket, and an simple geometry with varying isotopic concentrations. Results show that the clustering method significantly reduces computational burden without compromising fidelity, providing a foundation for agile iteration of neutronics simulations involving multiple coupled material properties.

Bae, Jin Whan [ORNL] (ORCID:0000000326548907)↗

Concavity-based local erosion and sphere-size-based local dilation applied to lithium-ion battery electrode microstructures for particle identification

Performance metrics of lithium-ion batteries can be extracted from the analysis of electrode microstructures nanoscale imaging. The characterization workflow can involve a challenging particle identification, or instance segmentation, step. In this work, we propose a new identification method based on an original transformation: a sphere-size-based local dilation followed by a concavity-based local erosion, that is local morphology closing. The new transformation is much more efficient than the global morphology closing, with correct identification achieved with only 1.7 % dilation volume and 2.6 % erosion volume on a test geometry, compared to 39.2 % and more than 50 %, respectively, with its global counterpart. The new method has been then benchmarked versus other identification algorithms (watershed and pseudo coulomb repulsive field) on a real electrode microstructure with equal or better segmentation achieved.

25 ENERGY STORAGE↗

Evaluating the Importance of Conformers for Understanding the Vacuum-Ultraviolet Spectra of Oxiranes: Experiment and Theory

Vacuum-ultraviolet (VUV) absorption spectroscopy enables electronic transitions that offer the unambiguous identification of molecules. As target molecules become more complex, multifunctional species present a great challenge to both experimental and computational spectroscopy. This research reports both experimental and theoretical studies of oxiranes. Computationally, the nuclear ensemble approach has been used to accurately predict experimental spectra for a variety of molecules. However, this approach incurs great computational cost, as ensembles generally consist of thousands of geometries. The present study aims to drastically reduce the ensemble by evaluating the significance of the conformers to the predicted spectra. This approach was applied to 11 substituted oxiranes using the Conformer Rotamer Ensemble Sampling Tool (CREST) of Grimme to generate an ensemble of unique conformers determined by their Boltzmann populations. Five TD-DFT functionals (BMK, CAM-B3LYP, M06-2X, MN15, ωB97X-D) and EOM-CCSD were used to simulate the spectrum of each substituted oxirane ensemble. Computed spectra were then compared to the experiment using both qualitative and quantitative metrics. Based on these metrics, it was observed that certain conformers may not be necessary to characterize this set of oxiranes despite the temperature (323 K) of the experiment. A single conformer can then be used with TD-DFT and EOM-CCSD to replicate the experimental spectra of these medium-sized combustion species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Minimizing Ground Risk in Cellular-Connected Drone Corridors With mmWave Links

Unmanned Aircraft Systems (UASs) have been receiving significant interest and support from academia, industry, and regulatory bodies over the past decade due to their various use cases. To safely integrate UAS operations into the national airspace, particularly overpopulated regions, the risk posed to ground users, buildings, and vehicles due to unmanned aerial vehicle (UAV) flight should be minimized. This risk can be represented by a numerical metric, which we refer to in this article as the “ground risk.” Many UAS applications also depend on the presence of a reliable wireless communication link between the UAV and a control station for the transmission of UAV position, surveillance video, UAV payload commands, and other mission-related data. Such wireless communication requirements also need to be considered in the design of UAS operations. In this article, we consider both these aspects and study the design of nonintersecting trajectories for UAS operations to minimize ground risk, subject to constraints on the wireless signal strength and geometry of the trajectory, specified in terms of: 1) an enclosing cylinder within which the trajectory must lie and 2) an integrated angular change along the UAV's trajectory. The performance of a computationally expensive optimal algorithm is compared with that of a computationally faster heuristic approach within the dense urban environment of Manhattan, NY, USA. Performance evaluation using ray-tracing simulations shows that the heuristic approach performs close to the optimal algorithm at a reduced computation cost. In conclusion, this research can be utilized to make UAS operations safe and reliable and accelerate their adoption.

99 GENERAL AND MISCELLANEOUS↗

Soft particles and infinite-dimensional geometry

Abstract In the sigma model, soft insertions of moduli scalars enact parallel transport ofS-matrix elements about the finite-dimensional moduli space of vacua, and the antisymmetric double-soft theorem calculates the curvature of the vacuum manifold. We explore the analogs of these statements in gauge theory and gravity in asymptotically flat spacetimes, where the relevant moduli spaces are infinite-dimensional. These models have spaces of vacua parameterized by (trivial) flat connections on the celestial sphere, and soft insertions of photons, gluons, and gravitons parallel transportS-matrix elements about these infinite-dimensional manifolds. We argue that the antisymmetric double-soft gluon theorem ind + 2 bulk dimensions computes the curvature of a connection on the infinite-dimensional space Map ( S d , G ) / G , whereGis the global part of the gauge group. The analogous metrics in abelian gauge theory and gravity are flat, as indicated by the vanishing of the antisymmetric double-soft theorems in those models. In other words, Feynman diagram calculations not only know about the vacuum manifold of Yang–Mills theory, they can also be used to compute its curvature. The results have interesting implications for flat space holography.

Astronomy & Astrophysics↗

Ab-Initio Simulation Studies of Cr Solvation in Fluoride Molten Salts

Understanding molten salt chemistry is essential in ongoing research of the molten salt nuclear reactor (MSR). In this context, detailed understanding of the mechanisms underlying selective oxidation of metal species, such as Cr, is required to guide the design of effective corrosion mitigation strategies in molten salts. An important starting point for such mechanistic understanding is knowledge of the solvation structure and its role in controlling metal speciation. In this work, we use ab initio molecular dynamics simulations to study the short-range (on the scale of the nearest-neighbor bond lengths) and medium-range (over length scales of several neighbor spacings) structure in three different fluoride melts with and without Cr addition; namely, 2KF-NaF, 2LiF-BeF 2 , and 3LiF-AlF3. We find that Cr 0 ,Cr 2+ ,Cr 3+ can each be coordinated by different numbers of F-, with the variance in coordination number decreasing as oxidation state increases, and that these coordination geometries are largely independent of solvent. The manner by which Cr changes the medium-range structure, however, is found to be solvent-dependent. While 2KF-NaF melts show short and medium range order that is highly dynamic, 2LiF-BeF 2 and 3LiF-AlF 3 are characterized by molecular associates that are relatively long-lived that organize into oligomer structures on larger length scales. Rather than being solvated by F- ions alone, we find that Cr can incorporate into and be solvated within this oligomer structure. Fluoroacidity, alone, may therefore prove too simple a metric for assessing the corrosivity of molten fluorides. As our work suggests, the ability of Cr to solvate must be understood in the context of the short- and medium-range structure of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗