Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large displacement”

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 325 records · Page 18

Atomic-scale insights into topotactic transformations in an extra-large-pore zeolite using time-resolved 3D electron diffraction

Understanding the atomic-scale structural dynamics of phase transformations is crucial for developing materials and tailoring their properties. However, many materials are obtained as polycrystalline powders with large unit cells and/or complex structures, making it challenging to investigate detailed structural changes using conventional X-ray diffraction techniques. Here we employ time-resolved three-dimensional electron diffraction to reveal the topotactic reactions and transformations that convert the extra-large-pore silicate zeolite ECNU-45 into ECNU-46. ECNU-45 features three-dimensional interconnecting 24 × 10 × 10-ring channels, while ECNU-46 consists of one-dimensional 24-ring channels connected to 10-ring pockets. ECNU-45 and ECNU-46 are both examples of pure silicate zeolites with pore openings larger than 22-ring. Our findings indicate changes at six distinct tetrahedral silicon sites, involving atom displacement, addition and removal of framework atoms through bond breakage and formation. This work presents the synthesis of zeolites and also provides atomic-level insights into the dynamic processes of topotactic reactions. Our results have implications for advancing materials engineering and understanding complex solid-state reactions at an atomic scale.

Inorganic chemistry↗

Accelerated Fuel Qualification of Fast Modular Reactor Fuel in a Thermal Reactor: Modeling and Simulation Paired with Irradiation Testing

The accelerated fuel qualification (AFQ) methodology is applied by simulating accelerated fuel tests of the General Atomics Electromagnetic Systems’ fuel system for its 44-MW(electric) gas-cooled, fast-spectrum fast modular reactor (FMR). This fuel is comprised of UO 2 pellets in SiGA® cladding, a silicon carbide ceramic matrix composite. Fast reactors, like the FMR, offer many benefits, including high fuel utilization and flexibility, but may require a lengthy material design process if tests are performed using fast neutron irradiation alone. A thermal neutron irradiation can instead be used to rapidly test how well key components of the current material models extend to high burnup. Thermal neutrons produce a different radial power distribution within the pin than fast neutrons. However, the temperature and burnup values for the two neutron types are comparable, and the differences between the simulated fuel responses are relatively small, demonstrating the weak sensitivity of the physics-based fuel model calculations on the neutron type and the irradiation rate. Furthermore, the deformation of the SiGA cladding saturates after about 1 displacement per atom for both neutron spectra. In an accelerated fuel test, the irradiation time required to reach the target fuel burnup can be reduced by a factor of 3 by using a small rodlet with a 45% smaller pellet diameter while maintaining the same linear power. Therefore, the time for data collection up to high burnup can be significantly reduced while maintaining the same temperature profile, which largely determines the material response. Tests of fuel rodlets of standard and compact size will be carried out in the Idaho National Laboratory’s Advanced Test Reactor (ATR), including full size and compact rodlets with varying gap sizes. By applying physics-based mechanistic modeling and simulation in accordance with the AFQ methodology, this type of compact rodlet testing in a thermal test reactor captures the necessary phenomena to test fuel material models up to high burnup and to simulate the expected impact of fast neutron radiation on the fuel in FMR operations. Finally, this approach to testing fast reactor fuels in existing thermal test reactors, paired with advanced physics-based mechanistic modeling and simulation, is expected to be applicable to a range of advanced fuels and will decrease the overall fuel qualification timeframe from decades to years.

Advanced test reactor (ATR)↗

Dynamic Rupture Modeling of the M7.1, 2019 Ridgecrest, California, Earthquake

An important part of the simulation technique implemented in the strong motion prediction recipe developed in Japan, originally proposed by Irikura and Miyake (2011) is the earthquake rupture characterization that is based on established empirical scaling relationships between source parameters and seismic moment. These scaling relationships and rupture characterization are derived from recorded data and are tested by simulations. Recent analysis of results obtained with the Recipe for the 2016 M7.1 Kumamoto earthquake performed by several researchers, indicate that in order to reproduce near-fault long-period ground motions with permanent displacements, typical for earthquakes with surface rupturing, the recipe needs modifications. In particular, improvements are required in characterizing the rupture kinematics at shallow depths, above the seismogenic zone (Irikura et al., 2019). Modeling of rupture dynamics for recorded earthquakes, and numerical simulations for different rupture scenarios, provides physical constrains on slip evolution in time and space, which can be used in adjustments of empirical source scaling relationships for large crustal earthquakes. In this study we mainly focused on rupture dynamics modeling of the 2019 M7 Ridgecrest, California earthquake. The characterization of the kinematic rupture model derived from the dynamic rupture models of this earthquake can inform further refinements of rupture modeling proposed by the Recipe. Through different rupture scenarios we investigated potential relationships between characteristics of Long Period Ground Motion Generation Area (LMGA), used in the Recipe, and asperity depth and rupture initiation.

58 GEOSCIENCES↗

Geophysical measurements of a demolition at NRF

On May 31 st , 2025, at approximately 14:40:00 UTC, the former Submarine First Westinghouse (SW1) structure housing a large crane was demolished at the Naval Reactors Facility (NRF) located within the perimeter of the Idaho National Laboratory (INL). Several targeted explosives were used to facilitate this demolition. The INL Seismic Monitoring Program (SMP) installed four temporary seismic instrument suites around the site of the demolition in order to monitor ground motion including measuring peak ground acceleration (PGA), peak ground velocity (PGV), and to calculate peak ground displacement (PGD). The demolition produced three types of waves: body waves, surface waves, and acoustic waves. The acoustic wave represents the largest signal received. The PGA, measured on the eastern component of the instrument located approximately 95m away from the demolition site, was 0.187g. The PGV measured was 0.0004 m/s, and the PGD was approximately 0.000006m. Rayleigh waves incident on an instrument located at the Advanced Test Reactor, about 7.6km away, were used to estimate a magnitude of 0.3 at the ATR, representing a yield of energy equivalent to approximately 3 grams of TNT.

58 - GEOSCIENCES↗

Development and experimental qualification of novel disruption prevention techniques on DIII-D

Novel disruption prevention solutions spanning a range of control regimes are being developed and tested on DIII-D to enable ITER success. First, a new real-time control algorithm has been developed and tested for regulating nearness to stability limits and maintaining safety-margins. Its first application has been for reliable prevention of vertical displacement events (VDEs) by adjusting plasma elongation (κ) and the inner-gap between the plasma and inner-wall in response to real- time open-loop VDE growth rate (γ) estimators. VDEs were robustly prevented up to average open-loop growth rates of 800 rad/s with initial tunings, with only applying shape modification when near safety limits. Second, the disruption risk during fast, emergency shutdown after large tearing and locked modes can be significantly improved by transitioning to a limited topology during shutdown. More than 50% of emergency limited shutdowns after locked modes reach a final normalized current I_N < 0.3 before terminating, scaling to the 3 MA ITER requirement. Furthermore, this is in contrast to diverted shutdowns, the majority of which disrupt at I_N > 0.8. Despite improvements, these results highlight the critical importance of early prevention. Third, a novel emergency shut down method has been developed which excites instabilities to form a warm, helical core post-thermal quench. The current quench extends to ~100ms and avoids VDEs and runaway electron generation. Novel real-time machine learning disruption prediction has been integrated with the DIII-D proximity controller, and a real- time compatible multi-mode MHD spectroscopy technique has been developed. Results presented here were enabled by a focused effort, the Disruption Free Protocol, in DIII-D’s 2019-20 campaign to complement disruption prevention experiments with a large piggy-back program. In addition to testing novel techniques, it is estimated to have helped avoid 32 potential disruptions in piggyback operations with rapid, early shutdowns after large rotating n=1 or locked modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development and experimental qualification of novel disruption prevention techniques on DIII-D

Novel disruption prevention solutions spanning a range of control regimes are being developed and tested on DIII-D to enable ITER success. First, a new real-time control algorithm has been developed and tested for regulating nearness to stability limits and maintaining safety-margins. Its first application has been for reliable prevention of vertical displacement events (VDEs) by adjusting plasma elongation (κ) and the inner-gap between the plasma and inner-wall in response to real- time open-loop VDE growth rate (γ) estimators. VDEs were robustly prevented up to average open-loop growth rates of 800 rad/s with initial tunings, with only applying shape modification when near safety limits. Second, the disruption risk during fast, emergency shutdown after large tearing and locked modes can be significantly improved by transitioning to a limited topology during shutdown. More than 50% of emergency limited shutdowns after locked modes reach a final normalized current I N < 0.3 before terminating, scaling to the 3 MA ITER requirement. This is in contrast to diverted shutdowns, the majority of which disrupt at I N > 0.8. Despite improvements, these results highlight the critical importance of early prevention. Third, a novel emergency shut down method has been developed which excites instabilities to form a warm, helical core post-thermal quench. The current quench extends to ~100ms and avoids VDEs and runaway electron generation. Novel real-time machine learning disruption prediction has been integrated with the DIII-D proximity controller, and a real- time compatible multi-mode MHD spectroscopy technique has been developed. Results presented here were enabled by a focused effort, the Disruption Free Protocol, in DIII-D’s 2019-20 campaign to complement disruption prevention experiments with a large piggy-back program. In addition to testing novel techniques, it is estimated to have helped avoid 32 potential disruptions in piggyback operations with rapid, early shutdowns after large rotating n=1 or locked modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bioproducts from high-protein algal biomass: An economic and environmental sustainability review and risk analysis

This report has High-protein algal biomass as an important bio-commodity that has the potential to provide a new source of sustainable protein products. Herein is a critical review that identifies 1) the most relevant sustainability findings related to the processing of proteinaceous algal biomass to higher value protein products and 2) the potential pathways to improve life cycle assessment (LCA) and techno-economic analysis (TEA) metrics, including life-cycle carbon dioxide equivalent (CO 2 eq), life cycle energy, and minimum selling price (MSP) of these products. The critical review of the literature revealed a large variation in model input parameters relating to these metrics. Therefore, a Monte Carlo analysis was conducted to assess the risk associated with these input variations. To understand the uncertainties that propagate into high-protein algae to products’ systems, we reviewed more than 20 state-of-the-art unit operations for algal biomass processing., including cell disruption, protein solubilization, protein precipitation and purification, and protein concentration. We evaluated displacement of proteinaceous products by algal-bioproducts, including ruminant feed, aquaculture feed, protein tablets, and biopolymers and biopolyesters, with prices in the market ranging from 1.9 to 120 $\$$ kg -1 protein. This review realized that the MSP of ruminant and non-ruminant feed ranges from 0.65±0.56 to 2.9±1.1 $\$$ kg -1 protein, and bioplastics’ MSP ranges from 0.97 to 7.0 $\$$ kg -1 protein. Regarding LCA metrics, there is limited research on life cycle energy in proteinaceous biomass concentration and bioproduct systems, reported at 32.7 MJ kg protein -1 , for animal feed displacement. Animal feed emissions in the literature report negative fluxes, representing environmental benefits, as low as -3.7 kgCO 2 eq kg -1 protein and positive fluxes, i.e., global warming potential, as high as 12.8 kgCO 2 eq kg -1 protein. There is limited research on bioplastics life cycle emissions reported at 0.6 kgCO 2 eq kg -1 protein. In general, the studies to date of algae-derived protein bioproducts showed similar life cycle emissions to soybean meals, nylon, polymers, and polystyrenes. Our risk analysis realized that more than 50% of scenarios can result in negative-net life cycle CO 2 eq emissions. This review and risk analysis assess and demonstrate the scenarios that improve economic and environmental sustainability metrics in high-protein algal bioproducts systems.

09 BIOMASS FUELS↗

Search for heavy neutral leptons in decays of W bosons using leptonic and semi-leptonic displaced vertices in $\sqrt{s}$ = 13 TeV pp collisions with the ATLAS detector

A search is performed for long-lived heavy neutral leptons (HNLs), produced through the decay of a W boson along with a muon or electron. Two channels are explored: a leptonic channel, in which the HNL decays into two leptons and a neutrino, and a semi-leptonic channel, in which the HNL decays into a lepton and a charged pion. The search is performed with 140 fb −1 of $\sqrt{s}$ = 13 TeV proton-proton collision data collected by ATLAS during Run 2 of the Large Hadron Collider. No excess of events is observed; Dirac-like and Majorana-like HNLs with masses below 14.5 GeV and mixing coefficients as small as 10 −7 are excluded at the 95% confidence level. The results are interpreted under different assumptions on the flavour of the leptons from the HNL decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The Role of Land‐Atmosphere Feedbacks in Midlatitude Wintertime Surface Temperature Variability

Accurately representing synoptic near-surface temperature variability is crucial to predict weather extremes, yet models remain biased. Existing studies primarily attribute wintertime midlatitude near-surface temperature variability to tropospheric large-scale advection, assuming minimal land influence. However, nudging the model's circulation toward observations yields little improvement in wintertime temperature variance over Northern Hemisphere land, suggesting that land-atmosphere interactions also warrant attention. We introduce a new scaling framework for temperature variance that incorporates local land-atmosphere feedbacks. Comparing our framework to the mixing length approach—which links temperature variance to the meridional temperature gradient and air parcel displacement (mixing length)—shows that land-atmosphere feedbacks are inherently embedded in the mixing length, a connection previously overlooked. Roles of land–atmosphere feedbacks are evaluated via model experiments with perturbed circulation, land, or both. We find that longwave radiative damping dominates temperature variance responses over meridional temperature gradient when both land and circulation are perturbed.

atmosheric science↗

Superionic-like diffusion in yttrium dihydride

For the next-generation high temperature microreactors, yttrium dihydride (YH 2 ) is an attractive solid state neutron moderator. Despite a number of recent investigations, the mechanism of hydrogen transport remains poorly understood. Experimental evaluations of diffusivity are inconclusive with large variations in diffusivities and activation energies. In this work, we perform ab initio molecular dynamics (AIMD) simulations on YH 2 for temperatures spanning 300 K to 1200 K. Our main finding is that YH 2 shows a superionic-like behavior with hydrogen atoms hopping from one native site to another above a characteristic temperature of 800 K. This correlated motion results in quasi-one-dimensional string-like displacements that enable the hydrogen atoms to diffuse rapidly. We confirm that the octahedral sites are mostly unoccupied, although channeling through them is the most favored pathway between lattice hops above 800 K. At the highest temperature of 1200 K, the string relaxation time is merely of the order of a few picoseconds, which indicates a liquid-like diffusive behavior. Based on the formation of spontaneous thermal vacancies, an order-disorder crossover temperature T α ~ 800 K is established for YH 2 with an activation energy of 0.83 eV for hydrogen diffusion in the superionic-like state.

Superionic-like Diffusion↗

Investigating the linkage between simulated precipitation climatology and ENSO-related precipitation anomaly based on multi-model and perturbed parameter ensembles

It remains a challenge for climate models to correctly capture the relationship between precipitation and ENSO. This studyexamines the linkage between the simulated precipitation climatology and ENSO-related precipitation anomaly during borealwinter based on the multi-model ensemble from the Atmospheric Model Intercomparison Project Phase 5 (AMIP5) and perturbedparameter ensemble (PPE) with the Beijing Climate Center (BCC) atmospheric model. The AMIP5 models have large biases insimulating the tropical precipitation anomaly during El Niño, such as the shifts of the Inter-Tropical Convergence Zone (ITCZ)and South Pacific Convergence Zone (SPCZ). The inter-model differences show that the precipitation change in response to seasurface temperature (SST) change increases with enhanced precipitation climatology. The ENSO-related precipitation anomalycan also be related to the spatial distribution of the mean-state precipitation. The simulated ITCZ/SPCZ displacements aresignificantly correlated with the spatial precipitation-SST relationship in the mean state. Models with stronger mean-stateprecipitation-SST relationship also produce stronger SPCZ/ITCZ displacements. The mean-state tropical precipitation also hasstrong impacts on the ENSO-related precipitation anomalies over East Asia. In the BCC PPE, the connection between the mean-state precipitation and ENSO-related precipitation anomaly is overall consistent with that in AMIP5. Parameters associated withthe low-cloud and deep convection processes are the most influential ones for the precipitation simulations in BCC. Comparedwith the version in AMIP5, the new BCC model can better simulate the precipitation climatology and the relationship betweenENSO and precipitation over southern China. These results have important implications for model development and model-basedclimate predictions.

Yang, Ben↗

Using high pressure solutions of polyfluoroacrylate and CO 2 to Seal cement cracks for improved wellbore integrity

Polyfluoroacrylate (PFA) is a hydrophobic and oleophobic polymer that is soluble in high pressure carbon dioxide (CO 2 ). In this study, the ability of PFA-CO 2 solutions to greatly reduce the apparent permeability of split or cracked Portland cement cylindrical samples is assessed. The apparent permeability values of confined samples were determined before and after treatment with PFA-CO 2 solutions. In four tests, PFA-CO 2 solutions were continuously displacing pure CO 2 from the cracked cement and the decrease in apparent permeability due to PFA adsorption and wettability alteration was monitored. The lowest apparent permeability cracked cement sample (81 nD) was completely sealed with a very small amount of solution. Samples with initial apparent permeabilities of 89 μD and 29.4 mD exhibited 92% and 99% reductions in permeability, respectively, before the experiments had to be stopped because of the excessively large increase in pressure drop. A 50% reduction in apparent permeability was observed with a 3.80 mD sample. Four other split cement samples (bound together with tape) with an initial apparent permeability in the 9.0–70 mD range were removed from the core holder and immersed in a PFA-CO 2 solution for 24 h to allow for PFA adsorption. Then the PFA-CO 2 solution was depressurized, allowing for the deposition of additional PFA from the solution within the crack as the pressure fell below the cloud point pressure of the PFA-CO 2 solution. These four samples were then confined again in a core holder and apparent permeability reductions of 29–93% were observed. Finally, results from these eight experiments indicates that the more substantial reductions in the nD – mD apparent permeability of the cracked cement correlated to lower initial crack permeability, higher PFA concentration, and slower injection rate of the PFA-CO 2 solution into the crack.

58 GEOSCIENCES↗

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nearly hyperuniform, nonhyperuniform, and antihyperuniform density fluctuations in two-dimensional transition metal dichalcogenides with defects

Hyperuniform many-body systems in d-dimensional Euclidean space R d are characterized by completely suppressed (normalized) infinite-wavelength density fluctuations, and appear to be endowed with novel exotic physical properties. Recently, hyperuniform systems of disordered varieties have been observed in the context of various atomic-scale two-dimensional (2D) materials. Here, we analyze the effects of localized defects on the density fluctuations across length scales and on the hyperuniformity property of experimental samples of two-dimensional transition metal dichalcogenides. In particular, we extract atomic coordinates from time series annular dark field-scanning transmission electron microscopy (ADF-STEM) imaging data of 2D tungsten chalcogenides with the 2H structure (Te-doped 2H-WSe 2 ) showing continuous development and evolution of electron-beam induced defects, and construct the corresponding chemical-bonding informed coordination networks between the atoms. We then compute a variety of pair statistics and bond-orientational statistics to characterize the samples. At low defect concentrations, the corresponding materials are nearly hyperuniform, characterized by significantly suppressed scattering at the zero wave-number limit (omitting forward/ballistic scattering). As more defects are introduced, the (approximate) hyperuniformity of the materials is gradually destroyed, and the system becomes non-hyperuniform even when the material still contains a significant amount of crystalline regions. At high defect concentrations, the structures become antihyperuniform with diverging (normalized) large-scale density fluctuations, mimicking those typically observed at the thermal critical points associated with phase transitions. Overall, the defected materials possess varying degrees of orientation order, and there is apparently no intermediate hexatic phase emerging. To understand the observed nearly hyperuniform density fluctuations in the slightly defected materials, we construct a minimalist structural model and demonstrate that the experimental samples can be essentially viewed as perturbed honeycomb crystals with small correlated displacements and double chalcogen vacancies. Moreover, the small correlated displacements alone can significantly degrade hyperuniformity of the perfect honeycomb structure. Therefore, even a small amount of vacancies, when coupled with correlated displacements, can completely destroy hyperuniformity of the system.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Multiaxial Plastic Deformation of Zircaloy-4 Nuclear Fuel Cladding Tubes

Here, this work is motivated by the desire to devise an internal pressure test that can mimic a displacement-controlled loading scenario and demonstrate how to apply the multiaxial stress and strain data from the test to develop an elastic/plastic constitutive model for a thin-walled tubular component. This is achieved by conducting simultaneous measurements of tangential and axial strain during the pressure test and integrating these strain measures into a feedback loop with the pressure controller. It is shown how data from such a test can be used to develop a large mechanical property data set relevant to biaxial loading conditions. The data obtained have high confidence evidenced by their low variability and alignment with other literature studies. Additionally, data from these internal pressure tests combined with full-tube axial tensile tests allow for the derivation of the Hill anisotropic yield function. The developed Hill yield function is validated by comparing the plastic strain ratios from the full-tube tension tests and by comparing the predicted yield stress in the tangential direction with measured values from ring tension tests in a previous study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Design of a Level-1 Long-Lived Particle Trigger in the CMS Hadron Calorimeter

The Hadron Calorimeter (HCAL) in the Compact Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC) was recently upgraded for Run 3 to introduce depth segmentation and online timing measurements, expanding the physics capabilities. In particular, the augmentation of the calorimeter information at the hardware trigger level enables quick identification and recording of long-lived particle decays using lower thresholds on energy-based event quantities. The depth segmentation and online timing capabilities are utilized in novel HCAL-based hardware-level triggers to identify displaced and delayed long-lived particles (LLPs), either decaying inside the calorimeter volume or arriving at a delayed time. This increases the sensitivity to LLP decays occurring up to almost 6 m from the collision point. This two-pronged calorimeter trigger approach leverages the new capabilities of the CMS HCAL to expand the phase space accessible in ongoing LLP searches. These triggers were deployed for Run 3 of the LHC, beginning data-taking in 2022, and we review the trigger implementation and calibration.

Kopp, Gillian Baron [Princeton University, NJ (Uni↗

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗

Growth-induced Donnan exclusion influences swelling kinetics in highly charged dynamic polymerization hydrogels

Polymeric gels crosslinked by DNA sequences can exploit DNA strand-displacement reactions to promote swelling through dynamic polymerization. The degree of swelling and the rate of swelling must be directly tunable to achieve the promise of programmable soft matter. Though the kinetics of the strand-displacement reaction provide insertion rates up to 10 4 /Molar/second as measured in bulk solution, DNA hydrogel swelling can take upwards of 30 h to complete. Computational modeling of the reaction-induced swelling of these gels with our recently-developed reactive electrochemomechanical theory (Zimmerman et al., 2024) suggests that their extraordinarily slow swelling is partly due to a scaling mismatch between the addition of charge and the addition of fluid volume, leading to a large transient increase in the fixed charge density. The significant increase in the gel’s fixed charge density, due to the binding of negatively charged DNA, sharply restricts the concentration of mobile hairpins through the phenomenon of Donnan charge exclusion, an effect commonly exploited in nanofiltration applications using polymeric membranes. The scaling problem is overcome when the mean additional swelling provided to the hydrogel by addition of a crosslink is above a critical value, thus the swelling outpaces the charge accumulation, leading the fixed charge density to drop and significantly accelerating the swelling process. This study shows that Donnan exclusion can explain the kinetics of DNA hydrogel swelling, and studies ways to modulate the reaction speed by either modifying the salt concentration or increasing or decreasing the number of base pairs in each DNA sequence.

42 ENGINEERING↗