Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flow”

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 505 records · Page 28

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗

High-resolution time-resolved experiments on mixing and entrainment of buoyant jets in stratified environments

Fluid jets interacting with a stratified layer play an important role in the safety of several reactor designs. In the containment of nuclear power plants, fluid jets dominate the transport and mixing of gaseous species and consequent hydrogen distribution in case of a severe accident. The mixing phenomena in the containment are driven by buoyant high-momentum injections (jets) and low momentum injection plumes. Mixing near the postulated break is initially dominated by high flow velocities. Plumes with moderate flow velocities are instead relevant in the break compartment during the long-term pressurization phase, or in any of the apertures between two connected compartments if the mass flows are sufficiently high and the density differences between efflux and ambient are sufficiently low. Phenomena of interest include free plumes (as produced by the efflux from the break compartment in a larger room or directly from a break flow), wall plumes (such those produced by low mass flows through inter-compartment apertures), and propagating stratification fronts in the ambient (for any stably stratified conditions). These phenomena have been highly ranked about nuclear reactor design, especially regarding of safety protocols. During a Pressurized Thermal Shock (PTS) scenario, the interaction between the cold ECCS injection plume and the stratified fluid present in the cold (or hot) leg is important in order to determine the temperature at the time-dependent temperature at the inlet of the reactor pressure vessel (RPV) and the potential to cause a thermal shock on the RPV wall. In sodium-cooled fast reactors (SFRs), core channels are typically hydro-dynamically isolated so that there exists a considerable temperature variation at the exit of adjacent fuel assemblies. All the above phenomena are characterized by the interaction of buoyant jets with the stratified flow. In stratified layers baroclinic forces create significant redistribution of turbulent kinetic energy and scales, which leads to anisotropy. This important physical phenomenon is highly three dimensional and is challenging to capture even with high-fidelity CFD simulations, due in part to lack of sufficiently resolved validation data. Furthermore, the experimental data available in the open literature do not feature the level of fidelity needed for an extensive validation of turbulence models in lower order CFD. To shed new lights into the crucial phenomena object of the present research project, it was proposed to conduct coordinated experiments and simulations at the University of Michigan and the George Washington University. The project has resulted in an experimental database of high-resolution time-resolved measurements of jets in uniform and stratified environments. The novel experimental data will be used to validate computational fluid dynamic (CFD) codes, including both Large Eddy Simulations (LES) and unsteady Reynolds-averaged Navier-Stokes equations (URANS) methodologies. In the Experimental and Multiphase flow (ECMF) laboratory at Univerisity of Michigan, we built two experimental facilities to investigate also the effect of scaling. The first facility, DESTROJER (DEnsity Stratified Turbulent ROund free Jet ExpeRiment), featuring a contoured jet nozzle with a diameter of D=12.7mm and a 1m×1m×1m cubic tank, which is made of acrylic glass for optical access. The ratio between the tank width and the nozzle diameter is equal to 78, which ensures that there is no direct interaction between the jet and the side walls. A second, modular experimental facility, features three different tank sizes of size 10×10×30, 20×20×30, 30×30×30 cm3 respectively (all tanks have the same height of 30 cm), and a jet diameter of 2mm. For the modular facility, tank-to-nozzle ratio of 50, 100, and 150 are obtained respectively. Experiments with different-density jet impactions and sharp interface with a density difference of 3.16% In the Laboratory at George Washington University, complementary experiments have been carried out at a facility featuring the same dimensions as DESTROYER. Experiments with different jet diameters have been carried out as well. High-resolution time-resolved measurements have been performed at all facilities using high-speed particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF). In this way, simultaneous measurements of velocity and density fields have been carried out. At the University of Michigan, wire-mesh sensors have been designed and built to gather high-resolution, time-dependent data of the stratified front in the tank. These measurements are important to provide accurate initial conditions for the CFD simulations. An important outcome of the project has been the development of novel techniques to achieve refractive index matching (RIM) of acqueous solutions with high density differences. RIM is crucial if optical techniques such as PIV are used to measure velocity fields. Previous to this project, experiments up to 3% density difference had been reported in the literature. With the methodology developed in this project, we have been able to push the envelop to density differences up to 9%. The experiments have been used to validate Reynolds Averaged Navier-Stokes (RANS) turbulence models and Large Eddy Simulations (LES) models. Both NEK5000 and the commercial CFD code STAR-CCM+ have been used. The experimental data will be made available to the community through the NE-KAMS databasehave been successfully conducted in the scaled facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimal Wetting Angles in Lattice Boltzmann Simulations of Viscous Fingering

We conduct pore-scale simulations of two-phase flow using the 2D Rothman–Keller colour gradient lattice Boltzmann method to study the effect of wettability on saturation at breakthrough (sweep) when the injected fluid first passes through the right boundary of the model. We performed a suite of 189 simulations in which a “red” fluid is injected at the left side of a 2D porous model that is initially saturated with a “blue” fluid spanning viscosity ratios M=ν r /ν b ∈[0.001,100] and wetting angles θ w ∈[0°,180°]. As expected, at low-viscosity ratios M=ν r /ν b $\ll$1 we observe viscous fingering in which narrow tendrils of the red fluid span the model, and for high-viscosity ratios M$\gg$1, we observe stable displacement. The viscous finger morphology is affected by the wetting angle with a tendency for more rounded fingers when the injected fluid is wetting. However, rather than the expected result of increased saturation with increasing wettability, we observe a complex saturation landscape at breakthrough as a function of viscosity ratio and wetting angle that contains hills and valleys with specific wetting angles at given viscosity ratios that maximize sweep. This unexpected result that sweep does not necessarily increase with wettability has major implications to enhanced oil recovery and suggests that the dynamics of multiphase flow in porous media has a complex relationship with the geometry of the medium and the hydrodynamical parameters.

Engineering↗

Evaluating the performance of parallel subsurface simulators: An illustrative example with PFLOTRAN

To better inform the subsurface scientist on the expected performance of parallel simulators, this work investigates performance of the reactive multiphase flow and multicomponent biogeochemical transport code PFLOTRAN as it is applied to several realistic modeling scenarios run on the Jaguar supercomputer. After a brief introduction to the code's parallel layout and code design, PFLOTRAN's parallel performance (measured through strong and weak scalability analyses) is evaluated in the context of conceptual model layout, software and algorithmic design, and known hardware limitations. PFLOTRAN scales well (with regard to strong scaling) for three realistic problem scenarios: (1) in situ leaching of copper from a mineral ore deposit within a 5-spot flow regime, (2) transient flow and solute transport within a regional doublet, and (3) a real-world problem involving uranium surface complexation within a heterogeneous and extremely dynamic variably saturated flow field. Weak scalability is discussed in detail for the regional doublet problem, and several difficulties with its interpretation are noted.

54 ENVIRONMENTAL SCIENCES↗

Transient seafloor venting on continental slopes from warming‐induced methane hydrate dissociation

Methane held in frozen hydrate cages within marine sediment comprises one of the largest carbon reservoirs on the planet. Recent submarine observations of widespread methane seepage may record hydrate dissociation due to oceanic warming, which consequently may further amplify climate change. Here, in this work, we simulate the effect of seafloor warming on marine hydrate deposits using a multiphase flow model. We show that hydrate dissociation, gas migration, and subsequent hydrate formation can generate temporary methane venting into the ocean through the hydrate stability zone. Methane seeps venting through the hydrate stability zone on the eastern Atlantic margin may record this process due to warming begun thousands of years ago. Our results contrast with the traditional view that venting occurs only updip of the hydrate stability zone.

hydrate↗

Experimental analysis of spatial correlation effects on capillary trapping of supercritical CO 2 at the intermediate laboratory scale in heterogeneous porous media

Several numerical studies have demonstrated that the heterogeneous nature of typical sedimentary formations can favorably dampen the accumulation of mobile CO 2 phase underneath the caprock. Core flooding experiments have also shown that contrasts in capillary entry pressure can lead to buildup of nonwetting fluid phase (NWP) at interfaces between facies. Explicit representation of geological heterogeneity at the intermediate (cm-to-m) scale is a powerful approach to identify the key mechanisms that control multiphase flow dynamics in porous media. The ability to carefully control flow regime and permeability contrast at a scale that is relevant to CO 2 plume dynamics in saline formations offers valuable information to understand immiscible displacement processes and provides a benchmark for mathematical models. To provide insight into the impact of capillary heterogeneity on flow dynamics and trapping efficiency of supercritical CO 2 under successive drainage and imbibition conditions, we present an experimental investigation conducted in a synthetic sand reservoir. By mimicking the interplay of governing forces at reservoir conditions via application of surrogate fluids, we performed three immiscible displacement experiments to observe the entrapment of NWP in heterogeneous porous media. Capillary trapping performance is evaluated for each scenario through spatial and temporal variations of NWP saturation; for this reason we adopted X-ray attenuation to precisely measure phase saturation throughout the flow domain and apply spatial moment analysis. Furthermore, the sweeping performance of two different permeability fields with comparable variance but distinct spatial correlation was compared against a homogeneous base case with equivalent mean permeability by means of spatial moment analysis.

58 GEOSCIENCES↗

An extended JBN method of determining unsteady-state two-phase relative permeability

Relative permeability is the reduction of permeability of porous media when subjected to multiphase flow and a key parameter in subsurface hydrology. The JBN method is a well-known method of obtaining relative permeability, which measures the overall pressure drop and the effluent phase ratio versus time during two-phase displacements. By assuming no capillary pressure or gravity, the JBN method obtains the relative permeabilities to both phases at the core outlet. Since data across a range of saturations are acquired in a relatively short time, this method is widely used. This work extends the JBN method by having (1) section-wise pressure drop measurements between the core inlet, four pressure taps on the core and the outlet, (2) local saturation measurements, and (3) local phase fluxes. With these data, the extended JBN method can determine relative permeabilities to both phases at each pressure tap of the core (not just at the core outlet). The JBN extension is shown using a data set where CO 2 invades a brine-filled core. From this it is found that the advantages of the extended JBN method over the regular JBN method are: (1) four times more data are obtained, and (2) data are more accurate because the capillary end effect is experimentally avoided. Avoiding the end effect results in tripling the saturation range, and obtaining relative permeabilities that are consistent with steady-state measurements and roughly 40% higher than those from the regular JBN method.

54 ENVIRONMENTAL SCIENCES↗

Length scale selects directionality of droplets on vibrating pillar ratchet

Directional control of droplet motion at room temperature is of interest for applications such as microfluidic devices, self-cleaning coatings, and directional adhesives. Here, arrays of tilted pillars ranging in height from the nanoscale to the microscale are used as structural ratchets to directionally transport water at room temperature. Water droplets deposited on vibrating chips with a nanostructured ratchet move preferentially in the direction of the feature tilt while the opposite directionality is observed in the case of microstructured ratchets. This remarkable switch in directionality is consistent with changes in the contact angle hysteresis. To glean further insights into the length scale dependent asymmetric contact angle hysteresis, the contact lines formed by a nonvolatile room temperature ionic liquid placed onto the tilted pillar arrays were visualized and analyzed in situ in a scanning electron microscope. As a result, the ability to tune droplet directionality by merely changing the length scale of surface features all etched at the same tilt angle would be a versatile tool for manipulating multiphase flows and for selecting droplet directionality in other lap-on-chip applications.

42 ENGINEERING↗

Experimental investigation on multiscale hydrodynamics in a novel gas–Liquid–Solid three phase jet-Loop reactor

Multiphase flow hydrodynamics in a novel gas–liquid–solid jet-loop reactor (JLR) were experimentally investigated at the macroscales and mesoscales. The chord length distribution was measured by an optical fiber probe and transformed for bubble size distribution through the maximum entropy method. The impacts of key operating conditions (superficial gas and liquid velocity, solid loading) on hydrodynamics at different axial and radial locations were comprehensively investigated. JLR was found to have good solid suspension ability owing to the internal circulation of bubbles and liquid flow. The gas holdup, axial liquid velocity, and bubble velocity increase with gas velocity, while liquid velocity has little influence on them. Compared with the gas–liquid JLRs, solids decrease the gas holdup and liquid circulation, reduces the bubble velocity and delays the flow development due to the enhanced interaction between bubbles and particles (Stokes number >1). This work also provides a benchmark data for computational fluid dynamics (CFD) model validation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of drilling fluid rheology and wellbore pressure on rock cuttings removal performance: numerical investigation

Drilling fluids exhibit non-Newtonian behavior which is generally described in the literature by two rheological models: Power Law and Herschel-Bulkley models. We have used Eulerian- Eulerian multiphase flow model to evaluate the drilling fluid flow patterns and filter cake formation patterns for both models for ultra-deep drilling conditions. Simulated drilling fluid flow patterns clearly show that the Herschel - Bulkley model exhibits fewer vortices around the drilling pipe, indicating better rock cuttings removal performance. In addition to mud rheology, we have analyzed the pressure drop variation with the filter cake formation on the wellbore during ultradeep drilling process. We have shown that filter cake thickness at the wall linearly increases with differential pressure between the well and formation. This implies that the differential pressure needs to be optimized. Very thin filter cake may lead to excessive fluid seepage into the permeable subterranean vertical zone resulting in a loss of fluid circulation and hence in a loss of rock cuttings removal. Yet, very thick filter cake reduces the effective diameter of the hole which may lead to a poor rock cuttings removal performance and stuck pipe. The mathematical model presented here needs to be experimentally validated when filter cake thickness data in deep drilling conditions become available. The model could be useful to provide drilling engineers insight on the filter cake formation, drilling fluid flow pattern, and give guide-line to optimize the filter cake as field experimental data are scarce.

02 PETROLEUM↗

Eigenparticles: characterizing particles using eigenfaces

The shape characteristics of particles have a pinnacle role in microsopic and macroscopic features of a system. Several studies have highlighted the need for considering deviations from a spherical representation of particles for accurate modeling of granular and multiphase flow systems. Using a shape factor, sphericity or roundness parameter alone is proven to be inadequate to capture the physical phenomena. In the present study we propose a novel metric based on the pattern recognition method Eigenfaces, coining the technique ‘Eigenparticles’. Using this technique we create a single statistical distribution of basis shapes to describe the morphological composition. The proposed technique is successfully validated with test shapes and applied to real particles. As a result, when compared with a state-of-the-art Fourier based method, ‘Eigenparticles’ performs favorably, clearly distinguishing the different particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theory and Applications of Macroscale Models in Porous Media

Systems dominated by heterogeneity over a multiplicity of scales, like porous media, still challenge our modeling efforts. The presence of disparate length- and time-scales that control dynamical processes in porous media hinders not only models predictive capabilities, but also their computational efficiency. Macrosopic models, i.e., averaged representations of pore-scale processes, are computationally efficient alternatives to microscale models in the study of transport phenomena in porous media at the system, field or device scale (i.e., at a scale much larger than a characteristic pore size). We present an overview of common upscaling methods used to formally derive macroscale equations from pore-scale (mass, momentum and energy) conservation laws. This review includes the volume averaging method, mixture theory, thermodynamically constrained averaging, homogenization, and renormalization group techniques. We apply these methods to a number of specific problems ranging from food processing to human bronchial system, and from diffusion to multiphase flow, to demonstrate the methods generality and flexibility in handling different applications. The primary intent of such an overview is not to provide a thorough review of all currently available upscaling techniques, nor a complete mathematical treatment of the ones presented, but rather a primer on some of the tools available for upscaling, the basic principles they are based upon, and their specific advantages and drawbacks, so to guide the reader in the choice of the most appropriate method for particular applications and of the most relevant technical literature.

42 ENGINEERING↗

Characterization of immiscible fluid displacement processes with various capillary numbers and viscosity ratios in 3D natural sandstone

We characterize the influence of reservoir conditions upon multiphase flow, we calculated fluid displacements (drainage processes) in 3D pore spaces of Berea sandstone using two-phase lattice Boltzmann (LB) simulations. The results of simulations under various conditions were used to classify the resulting two-phase flow behavior into three typical fluid displacement patterns on the diagram of capillary number (Ca) and viscosity ratio of the two fluids (M). In addition, the saturation of the nonwetting phase was calculated and mapped on the Ca–M diagram. We then characterized dynamic pore-filling events (i.e., Haines jumps) from the pressure variation of the nonwetting phase, and linked this behavior to the occurrence of capillary fingering. The results revealed the onset of capillary fingering in 3D natural rock at a higher Ca than in 2D homogeneous granular models, with the crossover region between typical displacement patterns broader than in the homogeneous granular model. Furthermore, saturation of the nonwetting phase mapped on the Ca–M diagram significantly depends on the rock models. These important differences between two-phase flow in 3D natural rock and in 2D homogeneous models could be due to the heterogeneity of pore geometry in the natural rock and differences in pore connectivity. By quantifying two-phase fluid behavior in the target reservoir rock under various conditions (e.g., saturation mapping on the Ca–M diagram), our approach could provide useful information for investigating suitable reservoir conditions for geo-fluid management (e.g., high CO 2 saturation in CO 2 storage).

58 GEOSCIENCES↗

CO 2 -brine relative permeability and capillary pressure of Tuscaloosa sandstone: Effect of anisotropy

Relative permeability and capillary pressure are known as essential properties that have substantial impacts on the accuracy of reservoir simulations. The effect of small-scale heterogeneity and lamination in the rock structure is often ignored during the measurement of capillary pressure and relative permeability curves in core samples. This study highlights the remarkable impact of anisotropy on the multiphase flow properties of stratified formations. A series of steady-state CO 2 -brine drainage and imbibition tests are conducted at reservoir conditions in horizontal and vertical core samples of the Tuscaloosa sandstone from the Cranfield CO 2 injection site in Mississippi. The relative permeability curves represent an anisotropic behavior influenced by the heterogeneous and laminated structure of realistic rock samples. The CO 2 saturation profiles during drainage and imbibition cycles indicate that the phase distribution in the pore space is controlled by core-scale heterogeneity in the porosity distribution among the laminations that causes capillary pressure inhomogeneity. Using the saturation profile during the imbibition cycle, the trapping characteristic of the horizontal and vertical rock samples are compared and we found that the capillary trapping is less likely in the vertical direction. Furthermore, the centrifuge-measured capillary pressure demonstrates distinctive characteristics for horizontal and vertical core samples. Since the flooding experiments are performed under capillary controlled flow, the capillary pressure contrast in the laminated structure of the rock strongly affects the relative permeability. The presented results can potentially improve the accuracy of the large-scale simulations for the CO 2 post-injection period, in which the vertical displacement has an important role in the plume migration.

58 GEOSCIENCES↗

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

Mixed imbibition controls the advance of wetting fluid in multiscale geological media

Imbibition plays a central role during multiphase flows in geological media; yet, our understanding of imbibition dynamics across the multiscale features of rocks, including fractures, microcracks, and heterogeneous matrices is limited. Here, we deployed a unique combination of micro- and clinical X-ray CT imaging to probe directly the spatiotemporal evolution of imbibition in a natural, multiscale shale sample. The method enabled unraveling of mixed imbibition behavior, with spontaneous imbibition rates in the matrix pores and microcracks driven by capillarity proceeding as $\propto t^{-1/2}$, whereas pressure-driven flows in the macroscopic fractures yield apparent imbibition rates with $\propto t^{-1}$ scaling. Regardless of the aperture, we showed that spontaneous imbibition proceeds uniformly away from the fractures. The effects of microcracks are further investigated by comparing the dynamics of the imbibition fronts with and without microcracks. Finally, these results highlight the paramount role of microcracks on anomalous front roughening and the need to account for microcracks appropriately in the development of predictive theoretical models.

58 GEOSCIENCES↗