Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pattern formation”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

An Efficient Numerical Algorithm for Solving Coupled Time-Dependent Ginzburg-Landau Equation for Superconductivity and Elasticity

A decoupled finite element algorithm is developed for simulating the vortex dynamics on an elastic superconductor which couples the time-dependent Ginzburg- Landau equation with the complex-valued superconducting order parameter and the vector-valued magnetic potential, and the elasticity equation. We present an iterative algorithm for the decoupled system arising from the time and spatial discretization using a combination of preconditioner, algebraic multigrid method (AMG) and preconditioned conjugate gradient method (PCG). The iterative algorithm allows us to perform large-scale three-dimensional simulations of mesoscale pattern formation during superconducting phase transitions with arbitrary elastic boundary conditions. Here, the performance and efficiency of the algorithm are numerically verified by several benchmark problems, exhibiting up to two orders of magnitude improvement depending on the scale of discrete system compared to the exact solver.

Efficiency↗
DOE OSTI2021

Novel Mechanisms for the Formation of Dislocation Cell Patterns in BCC Metal

In this study, we present the first simulation results of the formation of dislocation cell wall microstructures in tantalum subjected to shock loading. Dislocation patterns and cell wall formation are important to understanding the mechanical properties of the materials in which they spontaneously arise, and yet the processing and self-assembly mechanisms leading to their formation are poorly understood. By employing transmission electron microscopy and discrete dislocation dynamics, we propose a new mechanism involving coplanar dislocations and pseudo-dipole mixed dislocation arrays that is essential to the pattern formation process. Our large-scale 3D DDD simulations demonstrate the self-organization of dislocation networks into cell walls in deformed BCC metal (tantalum) persisting at strain ε = 20%. The simulation analysis captures several crucial aspects of how the dislocation cell pattern affects metal plasticity, as observed in experiments. Although experimental evidence is inconclusive regarding whether cell wall formation takes place at the shock front, after the shock, during release, or when the sample has had enough time to relax post-recovery, our simulations indicate cell wall formation occurs after the shock and before release. The extended Taylor hardening composite model effectively considers the non-uniform dislocation density when cell walls form and accurately describes the corresponding flow stress.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nonequilibrium Physics of Multiphase Flow in Porous Media: Wettability and Disorder

Carbon capture and geologic storage, dissociation of methane hydrates in permafrost, enhanced oil recovery, and water dropout in low-temperature fuel cells, all have something in common: two or more fluids flow simultaneously through a porous medium; and the displacement of one fluid by another is often unstable (either due to gravity or viscous forces). Yet, our ability to model multiphase flow in porous media has remained a challenge. The traditional equations are unable to predict, explain, or even reproduce, the formation of the complex patterns observed in experiments. The overarching goal of the project is to develop new physical understanding of the role of wettability and disorder in multiphase flow through permeable media, and to develop new mathematical and computational models at the pore scale and at the continuum scale.

58 GEOSCIENCES↗
DOE OSTI2023

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurement of Ehrlich-Schwoebel barrier contribution to the self-organized formation of ordered surface patterns on Ge(001)

Normal-incidence 1-keV Ar + ion bombardment leads to amorphization and ultrasmoothing of Ge at room temperature, but at elevated temperatures the Ge surface remains crystalline and is unstable to the formation of self-organized nanoscale patterns of ordered pyramid-shaped pits. The physical phenomenon distinguishing the high-temperature patterning from room-temperature ultrasmoothing is believed to be a surface instability due to the Ehrlich-Schwoebel barrier for diffusing vacancies and adatoms, which is not present on the amorphous material. This real-time grazing-incidence small-angle x-ray scattering study compares smoothing of a prepatterned Ge sample at room temperature with patterning of an initially flat Ge sample at an elevated temperature. In both experiments, when the nanoscale structures are relatively small in height, the average kinetics can be explained by a linear theory. Furthermore, the linear theory coefficients, indicating surface stability or instability, were extracted for both experiments. A comparison between the two measurements allows estimation of the contribution of the Ehrlich-Schwoebel barrier to the self-organized formation of ordered nanoscale patterns on crystalline Ge surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonequilibrium Phenomena in Plasmas in Contact with Liquids

Plasma-liquid interaction is a low-temperature plasma science frontier involving fundamental issues in plasma-materials interaction, multiphase kinetics, and collective effects. Plasmas in contact with liquids present two distinctive types of nonequilibrium: kinetic - across the interface - involving complex physical and chemical kinetics; and dissipative - along the interface - comprising superficial transport that is often conducive to instabilities and the formation of dissipative structures such as spot patterns. The project investigated dissipative and kinetic nonequilibrium phenomena concurrently in plasmas in contact with liquid water electrodes. The project’s objective was to determine the properties of dissipative (macroscopic) and kinetic (microscopic) nonequilibrium phenomena, and the interrelation between them, in low-temperature atmospheric-pressure plasmas in contact with liquid electrodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Rapid Patterning and Advanced Device Structures for Low Cost Manufacturable Crystalline Si IBC Cells

Recent record efficiency Si cells have had an interdigitated back contact (IBC) design where both positive and negative contacts are fabricated on the back side of the Si wafer. Many industry groups have investigated IBC cells but used photolithography to pattern the rear contacts although it is considered impractical for low cost, high volume manufacturing. Alternative patterning methods using lasers and mechanical masking have been utilized in fabrication of patterned regions in IBC Si solar cells to replace photolithography. Another exciting advance in Si solar cell technology has been replacing the high temperature diffusion of doped regions inside the crystalline Si wafer with low temperature deposition of very thin layers of doped amorphous hydrogenated Si (a-Si:H) on the Si wafer. This is called a Si heterojunction (HJ) device. The a-Si:H provides excellent passivation of surface defects and has produced the highest open circuit voltages (VOC) of any Si solar cell device structure. The Institute of Energy Conversion (IEC) demonstrated the first IBC-HJ solar cell combining these two strategies in 2007 and has fabricated IBC-HJ cells with 20% efficiency using three photolithography steps. Despite the demonstrated efficiency potential of this device by industry groups (>26%) there is no commercial production because of the challenges of patterning and processing the structure in an industrial environment. Our work sought to address that challenge. The objective of this 3 year project was to develop the processing for the IBC-HJ Si solar cell using lasers for patterning and contact formation instead of photolithography. Laser patterning enables rapid, contactless manufacturing of patterned regions. We intensively studied laser fired contacts, laser patterning of the a-Si multi-layer stacks and metal layers and the use of plasma shadow masks. After an exhaustive focus in the first year on the laser fired emitter (LFE) and contact (LFC), we were unable to achieve VOC greater than 660 mV, compared to the 720 mV required to meet our milestones. We discovered that an additional issue with our original IBC structure was the presence of an inversion layer at the back surface connecting the p and n regions which reduced the VOC and fill factor (FF). We developed innovative methods for characterizing the inversion layer. These two limitations in the original design lead to development of a new IBC-HJ process sequence and device structure which we called Plasma Masked Laser Processed (PMLP). It retained the original high efficiency features including manufacturability. The LFC was replaced with a standard n-type a-Si HJ contact. The inversion layer was eliminated by replacing the previous p-type stack in the gap with an n-type or SiN stack. The PMLP used laser ablation of a multi-layer a-Si stack followed by chemical etching to open the n-contact. It required deposition of a patterned stack through a mask in the plasma deposition chamber. This turned out to be a source of significant problems due to inevitable unwanted deposition ‘leakage’ under the mask. This formed a blocking contact on the emitter which significantly reduced Voc and FF. Several iterations in PMLP device structure and chemical etching steps resulted in a large improvements, e.g. the efficiency increasing from 3 to 15% and Voc from 450 to 660 mV but they were unable to completely eliminate it. The best IBC-HJ device we fabricated had only 15% efficiency. For perspective, our standard FHJ devices had 20% efficiency.

14 SOLAR ENERGY↗
DOE OSTI2026

QPatLib v1.0 — Measurement-based quantum simulation Pauli string unitary pattern collections

This Zenodo record accompanies the paper “Scalable Measurement-Based Quantum Simulation Patterns for Benchmarking” arXiv.2605.12502 and provides QPatLib v1.0 measurement-pattern datasets in human-readable JSONL together with a ZIP archive of OpenQASM 3.0 circuits used for validation and reproducibility. The patterns and circuits implement Pauli string unitaries for benchmark cases. Cases include all possible string combinations for less than 6 qubits and strings used in Hamiltonians for certain diatomic molecules for 6 or more qubits. Format: Each pattern_*.jsonl file is containins measurement patterns for all subsets for a given model/instance and subset strategy: it begins with a preamble containing model metadata, subset definitions, provenance, and (when feasible) full-pattern test results, followed by one pattern entry per subset. Each subset entry includes a required pattern_ascii field storing the measurement pattern in the measurement-calculus/Graphix standard with signal shifting, written left-to-right in the canonical order nodes → edges → measurements (with signal dependencies) → byproduct corrections (X/Z). The circuits are included as circuit_files.zip. Patterns in this record were validated against the corresponding circuits and checked for causal flow. Codes for generating these patterns can be found at QPatLib repository on Github

Graphix↗
DOE OSTI2021

Transverse barrier formation by electrical triggering of a metal-to-insulator transition

Application of an electric stimulus to a material with a metal-insulator transition can trigger a large resistance change. Resistive switching from an insulating into a metallic phase, which typically occurs by the formation of a conducting filament parallel to the current flow, is a highly active research topic. Using the magneto-optical Kerr imaging, we found that the opposite type of resistive switching, from a metal into an insulator, occurs in a reciprocal characteristic spatial pattern: the formation of an insulating barrier perpendicular to the driving current. This barrier formation leads to an unusual N-type negative differential resistance in the current-voltage characteristics. We further demonstrate that electrically inducing a transverse barrier enables a unique approach to voltage-controlled magnetism. By triggering the metal-to-insulator resistive switching in a magnetic material, local on/off control of ferromagnetism is achieved using a global voltage bias applied to the whole device.

36 MATERIALS SCIENCE↗
DOE OSTI2022

helixCAM: A platform for programmable cellular assembly in bacteria and human cells

Interactions between cells are indispensable for signaling and creating structure. The ability to direct precise cell-cell interactions would be powerful for engineering tissues, understanding signaling pathways, and directing immune cell targeting. In humans, intercellular interactions are mediated by cell adhesion molecules (CAMs). However, endogenous CAMs are natively expressed by many cells and tend to have cross-reactivity, making them unsuitable for programming specific interactions. Here, we showcase “helixCAM,” a platform for engineering synthetic CAMs by presenting coiled-coil peptides on the cell surface. helixCAMs were able to create specific cell-cell interactions and direct patterned aggregate formation in bacteria and human cells. Based on coiled-coil interaction principles, we built a set of rationally designed helixCAM libraries, which led to the discovery of additional high-performance helixCAM pairs. Finally, we applied this helixCAM toolkit for various multicellular engineering applications, such as spherical layering, adherent cell targeting, and surface patterning.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Wetting transition and fluid trapping in a microfluidic fracture

Immiscible fluid–fluid displacement in confined geometries is a fundamental process occurring in many natural phenomena and technological applications, from geological CO 2 sequestration to microfluidics. Due to the interactions between the fluids and the solid walls, fluid invasion undergoes a wetting transition from complete displacement at low displacement rates to leaving a film of the defending fluid on the confining surfaces at high displacement rates. While most real surfaces are rough, fundamental questions remain about the type of fluid–fluid displacement that can emerge in a confined, rough geometry. Here, we study immiscible displacement in a microfluidic device with a precisely controlled structured surface as an analogue for a rough fracture. We analyze the influence of the degree of surface roughness on the wetting transition and the formation of thin films of the defending liquid. We show experimentally, and rationalize theoretically, that roughness affects both the stability and dewetting dynamics of thin films, leading to distinct late-time morphologies of the undisplaced (trapped) fluid. Finally, we discuss the implications of our observations for geologic and technological applications.

42 ENGINEERING↗
DOE OSTI2020

Anisotropic Dynamics and Mechanics of Macromolecular Crystals Containing Lattice-Patterned Polymer Networks

The mechanical and functional properties of many crystalline materials depend on cooperative changes in lattice arrangements in response to external perturbations. However, the flexibility and adaptiveness of crystalline materials are limited. Additionally, the bottom-up, molecular-level design of crystals with desired dynamic and mechanical properties at the macroscopic level remains a considerable challenge. To address these challenges, we had previously integrated mesoporous, cubic ferritin crystals with hydrogel networks, resulting in hybrid materials (polymer integrated crystals or PIX) which could undergo dramatic structural changes while maintaining crystalline periodicity and display efficient self-healing. Here, the dynamics and mechanics of these ferritin-PIX were devoid of directionality, which is an important attribute of many molecular and macroscopic materials/devices. In this study, we report that such directionality can be achieved through the use of ferritin crystals with anisotropic symmetries (rhombohedral or trigonal), which enable the templated formation of patterned hydrogel networks in crystallo. The resulting PIX expand and contract anisotropically without losing crystallinity, undergo prompt bending motions in response to stimuli, and self-heal efficiently, capturing some of the essential features of sophisticated biological devices like skeletal muscles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Thermal and Catalytic Decomposition of 2-Hydroxyethylhydrazine and 2-Hydroxyethylhydrazinium Nitrate Ionic Liquid

To develop chemical kinetics models for the combustion of ionic liquid-based monopropellants, identification of the elementary steps in the thermal and catalytic decomposition of components such as 2-hydroxyethylhydrazinium nitrate (HEHN) is needed but is currently not well understood. The first decomposition step in protic ionic liquids such as HEHN is typically the proton transfer from the cation to the anion, resulting in the formation of 2-hydroxyethylhydrazine (HEH) and HNO 3 . In the first part of this investigation, the high-temperature thermal decomposition of HEH is probed with flash pyrolysis (<1400 K) and vacuum ultraviolet (10.45 eV) photoionization time-of-flight mass spectrometry (VUV-PI-TOFMS). Next, the investigation into the thermal and catalytic decomposition of HEHN includes two mass spectrometric techniques: (1) tunable VUV-PI-TOFMS (7.4-15 eV) and (2) ambient ionization mass spectrometry utilizing both plasma and laser ionization techniques whereby HEHN is introduced onto a heated inert or iridium catalytic surface and the products are probed. The products can be identified by their masses, their ionization energies, and their collision-induced fragmentation patterns. Further, formation of product species indicates that catalytic surface recombination is an important reaction process in the decomposition mechanism of HEHN. The products and their possible elementary reaction mechanisms are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Aging phenomena in the two-dimensional complex Ginzburg-Landau equation

The complex Ginzburg-Landau equation with additive noise is a stochastic partial differential equation that describes a remarkably wide range of physical systems which include coupled non-linear oscillators subject to external noise near a Hopf bifurcation instability and spontaneous structure formation in non-equilibrium systems, e.g., in cyclically competing populations or oscillatory chemical reactions. Here, we employ a finite-difference method to numerically solve the noisy complex Ginzburg-Landau equation on a two-dimensional domain with the goal to investigate its non-equilibrium dynamics when the system is quenched into the “defocusing spiral quadrant”. We observe slow coarsening dynamics as oppositely charged topological defects annihilate each other, and characterize the ensuing aging scaling behavior. We conclude that the physical aging features in this system are governed by non-universal aging scaling exponents. We also investigate systems with control parameters residing in the “focusing quadrant”, and identify slow aging kinetics in that regime as well. Finally, we provide heuristic criteria for the existence of slow coarsening dynamics and physical aging behavior in the complex Ginzburg-Landau equation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

From classical to quantum dynamics of atomic and ionic species interacting with graphene and its analogue

Graphene and its analogues offer a broad range of application opportunities for (opto)-electronics, sensing, catalysis, phase separation, energy conversion and storage, etc. Engineering graphene properties often relies on its controllable functionalization, defect formation and patterning, and reactive gas etching. In this chapter, we survey the dynamics of graphene using classical and quantum-classical dynamics methods. We discuss the reactivity, scattering, and transmission of atomic and ionic species including Ar cluster ion, H/D, and H+/D+ on graphene flakes of various sizes, focusing on the atomic-scale motion and energy dissipation pathways involved in forming and breaking covalent bonding. Discussions on the nuclear quantum effects of light species, the effects of isotopic substitution, and the methodologies for such modeling are also included.

Garashchuk, Sofiya V.↗
DOE OSTI2024

On principles of emergent organization

After more than a century of concerted effort, physics still lacks basic principles of spontaneous organization. To appreciate why, we first state the problem, outline historical approaches, and survey the present state of the physics of self-organization. This frames the particular challenges arising from mathematical intractability and the resulting need for computational approaches, as well as those arising from a chronic failure to define structure. Then an overview of two modern mathematical formulations of organization—intrinsic computation and evolution operators—lays out a way to overcome these challenges. Together, the vantage point they afford shows how to account for the emergence of structured states via a statistical mechanics of systems arbitrarily far from equilibrium. The result is a constructive path forward to principles of organization that builds on mathematical identification of structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Impacts of Mixed-Wettability on Brine Drainage and Supercritical CO 2 Storage Efficiency in a 2.5-D Heterogeneous Micromodel

Geological carbon storage (GCS) involves unstable drainage processes, the formation of patterns in a morphologically unstable interface between two fluids in a porous medium during drainage. The unstable drainage processes affect CO 2 storage efficiency and plume distribution and can be greatly complicated by the mixed-wet nature of rock surfaces common in hydrocarbon reservoirs where supercritical CO 2 (scCO 2 ) is used in enhanced oil recovery. Here, we performed scCO 2 injection (brine drainage) experiments at 8.5 MPa and 45°C in heterogeneous micromodels, two mixed-wet with varying water- and intermediate-wet patches, and one water-wet. The flow regime changes from capillary fingering through crossover to viscous fingering in the micromodels of the same pore geometry but different wetting surfaces at displacement rates with logCa (capillary number) increasing from –8.1 to –4.4. While the mixed-wet micromodel with uniformly distributed intermediate-wet patches yields ~0.15 scCO 2 saturation increase at both capillary fingering and crossover flow regimes (–8.1 ≤ logCa ≤ –6.1), the one heterogeneous wetting to scCO 2 results in ~0.09 saturation increase only at the crossover flow regime (–7.1 ≤ logCa ≤ –6.1). The interconnected flow paths in the former are quantified and compared to the channelized scCO 2 flow through intermediate-wet patches in the latter by topological analysis. At logCa > –6.1 (near well), the effects of wettability and pore geometry are suppressed by strong viscous force. Both scCO 2 saturation and distribution suggest the importance of wettability on CO 2 storage efficiency and plume shape in reservoirs and capillary leakage through caprock at GCS conditions.

58 GEOSCIENCES↗