Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BUBBLE”

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

Favorable morphology and electronic conductivity of functional sublayers for highly efficient water splitting electrodes

Low electronic conductivities and improper morphologies of anode electrodes greatly limit the reaction area, catalyst utilization and efficiency in proton exchange membrane water electrolyzers. In this study, conductive sublayers with different conductivities and morphologies were introduced into anode electrodes in membrane-based water electrolyzers. In-situ and ex-situ investigation results showed that conductive sublayers (Au mesh and carbon nanotube (CNT) film) augmented the sheet conductivity of anode electrodes by up to 4000 times (from 2000 to 0.5 ohm square -1 ), and the ohmic resistance of water electrolyzers was reduced to 1/3 when inserting conductive sublayers. In addition, CNT film provided a higher electrochemical active area than Au mesh, because of favorable morphologies (large porosity and surface area) of CNT fibers on CNT films. Therefore, the current density of water splitting was increased by 3 times (from 4.55 mA cm -2 to 14.83 mA cm -2 ) at 2.5 V compared to a conventional anode electrode. Visualizations on bubble dynamics showed improved performances with conductive sublayers; this was mainly due to greatly increased number of reaction sites, highly spread reaction area (from 50 to 1000 um), and reduced activation overpotential. Therefore, a balance between high electronic conductivity and nanoporous morphology is essential to the anode electrode for larger reaction sites and areas in highly efficient water electrolyzers.

25 ENERGY STORAGE↗
DOE OSTI2022

Model development for flash boiling spray and validations with isooctane, hexane, ethanol and their binary mixtures

The objective of this study is to develop a comprehensive model for spray atomization under the consideration of flash boiling, which includes nucleation, bubble growth, and droplet breakup. Although flash boiling has been previously studied, there is still no well-established analytical and numerical model of their mechanisms, or general guidelines on the implementation of this atomization strategy. In this study, detailed formulations are elaborated for the unified spray atomization model under both non-flash and flash boiling conditions. Next, the model was implemented into CFD software and 3-D simulations were performed for model validations. Two fuel blends are tested: ethanol-isooctane blends and hexane-isooctane blends. Simulation results were compared with the experimental measurements conducted using two methodologies, DBI (Diffuse Backlight Imaging) system for spray images captured by a high-speed camera, and a PDA (Phase Doppler Anemometry) system for droplet size measurement during the spray injection process. It was found that the current flash boiling model can accurately capture the spray characteristics under different environment conditions, such as the spray plume interactions and liquid penetration. The averaged Sauter mean diameter (SMD) calculations match the measurement sufficiently well. The spray structures for both single- and multi-component fuels during the injection process are well represented by current model.

42 ENGINEERING↗
DOE OSTI2023

Methane foam performance in oil-wet unconsolidated porous media: A systematic experimental investigation at reservoir conditions

Methane foam performance was systematically evaluated using high-pressure and high-temperature experiments performed on oil-wet proppant packs. The sensitivities of foam performance to various key factors, including surfactant concentration, foam quality, total injection rate, salinity, operating pressure, and initial oil saturation were examined. The results showed that even though the foamability and foam strength in oil-wet proppant packs were adversely affected by wettability and presence of oil, the steady-state pressure drops across the proppant packs were insensitive to the quantity of oil initially present. The presence of oil above the critical oil saturation impeded the onset of foam generation. This suggested that the foamability of surfactants depends on the synergetic effects of surfactant’s ionic nature and foam parameters, such as foam quality and total injection rate, on their ability to reduce the oil saturation to a threshold value below which favorable bubble generation sites could be created. It was observed that, for low injection rate, anionic surfactant was ineffective at low salinity conditions, while amphoteric agent performed well at all salinities. Additionally, an increase in salinity resulted in enhanced pressure drop across the proppant packs for both anionic and amphoteric surfactants. Further, the results also established that too high operating pressures are somewhat detrimental to the foam strength and stability, irrespective of the ionic nature of foaming agents and the wettability of the porous medium. For the anionic surfactant, the increase in the total flow rate was found to mitigate the adverse effects of low salt concentration and operating pressure. Results from quality-sensitivity foam tests identified two regions for the foam behavior, named low-quality and high-quality regimes, separated by the transition foam quality. Finally, we discuss the implications of the results for the design and effective implementation of enhanced oil recovery (EOR) schemes that deploy hydrocarbon-based foams.

02 PETROLEUM↗
DOE OSTI2023

Effect of He seeding on hydrogen isotope permeation in tungsten by H-D mixed plasma exposure

Hydrogen (H) and deuterium (D) permeation fluxes for tungsten (W) by H-D mixed plasma were studied to clarify the hydrogen isotope effect on permeation by a newly designed plasma driven permeation (PDP) device at Shizuoka University. A mixed H-D plasma was produced by introducing H and D gasses with various ratios. The atomic ratio and permeation flux in the plasma were measured during the experiments. It was found that the H:D atomic ratio in the plasma was almost proportional to the introduced H:D gas ratio. When single species of H or D was irradiated, the H permeation rate showed approximately 1.6 times as high as that of D. Under H-D mixed plasma irradiation, the highest HD permeation rate was observed when the atomic ratio of H:D in the plasma was 50:50. In addition, the effect of He seeding on H-D mixed plasma driven permeation behavior was also evaluated. The recycling of hydrogen isotopes on the upstream side was enhanced due to the formation of He bubbles near the surface, which would prevent the dissolution of hydrogen isotopes into bulk W and enhance the reflection process during plasma exposure. H and D permeation fluxes were clearly reduced by seeding He in H-D mixed plasma. At lower temperature (~723 K), the ratio of steady permeation fluxes of H and D was about 75 : 25, confirming the high permeation flux of H . These results have demonstrated a significant hydrogen isotope effect, even if He was seeded or not. As the temperature increased, H and D permeation fluxes approached unity at 923 K. Therefore,this indicates that the hydrogen isotope effect was determined by the combination of recycling / reflection at or beneath the surface, which will be controlled by the temperature during plasma exposure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗
DOE OSTI2023

Outgassing behavior and heat treatment optimization of JSC-1A lunar regolith simulant

As NASA strives towards a long duration presence on the Moon, it has become increasingly important to learn how to better utilize resources from the lunar surface for everything from habitats, vehicle infrastructure, and chemical extraction. To that end, a variety of lunar simulants have been sourced from terrestrially available volcanic minerals and glass as Apollo regolith is unavailable for experimentation needing large masses. However, while mineralogy and chemical composition can approach that of lunar material in these simulants, there are still distinct non-lunar phases such as hydrates, carbonates, sulfates, and clays that can cause simulants to behave distinctly non-lunar in a variety of processing conditions that maybe applied in-situ to lunar material. Notably, severe glassy bubbling has been documented in a variety of vacuum sintering experiments on JSC-1A lunar mare simulant heated via microwaves. The origins of this outgassing have not been well understood but are normally attributed to the decomposition of non-lunar contaminates intrinsic to virtually all terrestrially sourced simulants. As such, a series of controlled environmental tests were performed to ascertain the origins of the high temperature outgassing and to develop heat treatments that can drive JSC-1A closer to lunar composition and behavior. It was found that in JSC-1A at elevated temperatures distinct gas evolutions of water, carbon dioxide, and sulfur dioxide occur in both inert gas and vacuum. Additionally, the presence of hydrogen during heat treatments was shown to dramatically change gas evolutions, leading to distinctly more lunar-like composition and behavior from JSC-1A simulant.

58 GEOSCIENCES↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2026

Predicting critical heat flux using localized sensing at invisible vapor-liquid interfaces

Predicting critical heat flux (CHF) in two-phase electronics cooling systems remains a significant challenge due to the sudden onset of boiling crisis and the difficulty in directly visualizing vapor-liquid interfaces. Existing sensing methods rely on lagging temperature measurements, optically accessible systems, or spatially averaged signals that cannot pinpoint CHF initiation at localized high-heat-flux regions. Here, we report a planar capacitive sensing approach that enables real-time, localized detection of vapor-liquid interface dynamics for CHF prediction in boiling heat transfer. The capacitive sensor exploits the dielectric constant difference between liquid and vapor phases to capture bubble nucleation, growth, and departure dynamics with a temporal resolution down to 2 ms. The capacitive sensing reveals distinct signals across boiling regimes: from high-frequency fluctuations during strong nucleate boiling to low-frequency fluctuations with increased amplitudes when approaching CHF. The multi-sensor array experiments demonstrate real-time localized sensing, where each sensor responds exclusively to boiling in its immediate vicinity without crosstalk from neighboring regions. This non-intrusive sensing approach provides predictive rather than lagging sensing signals of CHF occurrence, offering predictive diagnosis of two-phase liquid cooling for the thermal management of high-power-density electronics.

CHF↗
DOE OSTI2025

Droplet breakup by multimodal nonlinear Rayleigh Taylor instability

A droplet impacted by a shock wave will undergo a process of fragmentation due to the development of interfacial hydrodynamic instabilities. The interface experiences variable acceleration and shear that result in the development of both inertial (Rayleigh–Taylor) and shear (Kelvin–Helmholtz) instabilities. These perturbations grow in time and drive the fragmentation and breakup of the deformed droplet. Experiments are performed on nominally 0.86 mm water droplet subjected to a Mach 7.6 detonation wave, resulting in a high Weber number (~36,000) breakup event. Perturbation growth is measured from a series of high-speed (> MHz) shadowgraph images. It is proposed that, given the size of the large-scale perturbations observed in experiments, these instabilities are growing in the non-linear regime and can be described by bubble-merger models for nonlinear mixing. Calculations are performed for the growth rates and size of these instabilities using deformation and external flow models to establish the time-dependent boundary conditions. The concurrence of the measured perturbation widths and the predictions of the simple model lend credence to the theory. This novel approach serves to open a new avenue in the characterization of droplet breakup via hydrodynamic instabilities.

Droplet breakup↗
DOE OSTI2025

Analysis of anion exchange membrane water electrolyzer performance and its evolution over time

Understanding water, evolved gas, and ionic transport in membrane-electrode-assemblies (MEAs) is essential for the development of high performance and durable anion exchange membrane water electrolyzers (AEMWEs). This study evaluates the MEA conditioning process, operating conditions, and short-term stability in a 1 M potassium hydroxide (KOH) electrolyte, focusing on the underlying transport phenomena. We observe a significant initial voltage loss in continuous cell operation, which could be associated with gas bubble accumulation, transport layer or flow field passivation, and changes in the catalyst oxidation state. Further, we investigate the effects of materials and operational configurations, including the membrane type and thickness, and the electrolyte flow rate, including KOH being fed to both electrodes as well as to the anode only. Furthermore, the effect of membrane drying temperature on ex situ as well as in situ electrochemical performance is evaluated. Finally, we discuss 700 h of AEMWE operation at 1 A/cm 2 , highlighting the underlying degradation phenomena.

25 ENERGY STORAGE↗
DOE OSTI2025

Modulation of thermal conductivity of iron-doped ß-Ga2O3 by helium-ion irradiation

This study examines the impact of helium-ion irradiation on the thermal conductivity of ß-Ga2O3. A laser-based spatial domain thermoreflectance technique is used to investigate thermal conductivity map for both un-irradiated and irradiated ß-Ga2O3, which are then validated against simulation results derived from density functional theory-based phonon transport simulations. Since helium bubble evolution was ob- served at the nanoscale using transmission electron microscopy, the simulation study was carried out on eight distinct helium-induced sites in ß-Ga2O3. Our findings indicate a reduction in thermal conductivity for the irradiated samples. Experimental results show a significant reduction in thermal conductivity in irradiated samples, with de- creases of approximately 25% along the [100] direction and 40% along [001] directions. Phonon transport simulations closely replicate these findings, particularly when helium occupying interstitial sites, predicting reductions of ˜53% along [100] and ˜50% along [001] directions. This work underscores the role of irradiation-induced microstructural changes in the heat transport properties of ß-Ga2O3 which is crucial for its application in sensor devices in extreme environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2020

Aqueous Ni-rich-cathode dispersions processed with phosphoric acid for lithium-ion batteries with ultra-thick electrodes

Lithium-ion battery (LIB) production can benefit both economically and environmentally from aqueous processing. Although these electrodes have the potential to surpass electrodes conventionally processed with N -methyl-2-pyrrolidone (NMP) in terms of performance, significant issues still exist with respect to ultra-thick cathodes ($\gg$4 mAh/cm 2 areal capacities). A major concern for these types of electrodes with high-nickel active material stems from lithium leaching from active material, which drives the pH of the dispersion in excess of 12 and subsequently corrodes the current collector interface. As this corrosion reaction proceeds, hydrogen generation at the interface creates bubbles which cause severe cracking in the dried electrode surface. When areal loadings are increased, this effect becomes more pronounced and is detrimental to both mechanical and electrochemical properties of these electrodes. In this work, a technique for mitigating corrosion at the current collector by adjusting the pH of the dispersion with the addition of phosphoric acid is investigated. Phosphoric acid was added in 0.5 wt% increments between 0.0 and 1.5 wt%, and effects on rheology, adhesion, corrosion, and electrochemical performance were investigated. A technique is reported for producing aqueous processed cathodes with areal loadings of 6–8 mAh/cm 2 with reduced surface cracking and superior high-rate discharge capacity (i.e. high-power performance) for this class of cathode loadings.

25 ENERGY STORAGE↗
DOE OSTI2024

Enhanced rare earth element recovery with cross-linked glutaraldehyde-lanthanide binding peptides in foam-based separations

Lanthanide Binding Tag (LBT) peptides that coordinate selectively with lanthanide ions can be used to replace the energy intensive processes used for the separation of rare earth elements (REEs). These surface-active biomolecules, once selectively complexed with the trivalent REE cations, can adsorb to air/aqueous interfaces of bubbles for foam-based REEs recovery. Glutaraldehyde, an organic compound that is a homobifunctional crosslinker for proteins and peptides, can be used to enhance the adsorption and interfacial stabilization of lanthanide-bound peptides films. The stability of the interfacial cross-linked films was tested by measuring their dilational and shear surface rheological properties. Surface activity of the adsorbed species was analyzed using pendant drop tensiometry, while surface density and molecular arrangement were determined using x-ray reflectivity and x-ray fluorescence near total reflection. Glutaraldehyde cross-linked REE-peptide complexes enhance the adsorption of lanthanides to air-water interfaces, resulting in thicker interfacial structures. Subsequently, these thicker layers enhance the dilational and shear interfacial rheological properties. The interfacial film stabilization and REEs extraction promoted by the cross-linker presented in this work provides an approach to integrate glutaraldehyde as a substitute of common foam stabilizers such as polymers, surfactants, and particles to optimize the recovery of REEs when using biomolecules as extractants.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗
DOE OSTI2025

Conversion of CO 2 from power plant into CaCO 3 nanoparticles

Carbon dioxide (CO 2 ), a main composition of flue gas, represents a significant and largely untapped carbon resource. Herein, mediated by glycine (Gly), we captured and converted CO 2 into CaCO 3 nanoparticles using real flue gas from a power plant, demonstrating for the first time the feasibility of using amino acid to convert CO 2 from power plant flue gasses. The method did not require extraneous energy and CaCO 3 nanoparticles with a size of ∼25 nm were obtained. Moreover, the potential toxicity of CO 2 -converted nanoparticles was investigated. It appeared that both the initial CO 2 loading and the carbamate percentage significantly influence the shape and size of the CaCO 3 particles. Our method was also proven effective for flue gas with varying CO 2 concentrations (4 %, 12 %, and 20 %). By tuning flue gas bubbling time and flow rate to achieve consistent CO 2 loading and carbamate levels, we produced CaCO 3 nanoparticles with similar shapes and sizes across all CO 2 concentrations studied. In addition, our data indicated that although real flue gas contains small amounts of gases like oxygen and CO, they insignificantly influence the shape and size of our nanoparticles but did impact the phase component of CaCO 3 . In conclusion, the toxicity experiments found that CaCO 3 nanoparticles produced from both real flue gas and simulated flue gas exhibited concentration- and time-dependent effects on cell viability.

36 MATERIALS SCIENCE↗