csm-logo-headerColorado School of Mines
207 Marquez Hall
Golden, Colorado 80401
Ph: (303) 384-2449


Theoretical models for particulate flows often treat particles and fluid as two interpenetrating continua. In such models, particle-fluid interactions are coarse-grained to the continuum level, giving rise to constitutive models that are, nowadays, often obtained through detailed numerical solutions of the Navier-Stokes equation and particles' equations of motion. Why accurate constitutive models are important? Without them, predictions of continuum-scale models cannot be quantitative. In simulations of proppant transport, the use of an accurrate sedimentation correlation can significantly improve the accuracy of prediction. Using detailed simulations to evaluate / improve constitutive models for particulate flows is a significant research direction of our group. A current focus of our studies is proppant transport.


proppant distribution

Simulated formation of proppant dunes and comparison with experiments (courtesy of Prof. Dongxiao Zhang and Prof. Qingzhi Wen at Peking University).


Many constitutive models operate on the assumption that particles are homogeneously distributed in the fluid, which is usually not true for gas-solid suspensions and certain non-Newtonian suspensions. Take dilute gas-solid flows as an example, particle clusters may arise because 1) particle-particle collisions can be inelastic; 2) viscous dissipation of particles' fluctuating motion increases nonlinearly with solid volume fraction; 3) particle-fluid drag increases nonlinearly with solid volume fraction. Because of computational cost, it is usually not possible to use very small grid size to ensure that the clusters are adequately resolved. Hence, there are many studies that focus on modeling the effect of clusters on sub-grid length scales. Using detailed simulations to extract properties of heterogeneous suspension flows with clusters and voids is another direction of interest to us.


Sedimentation of particles (only centers are shown) through a dense gas in a quasi-2D periodic domain. Simulation is conducted using SUSP3D code. Ar = 71, solid fraction = 0.15, normal restitution = 0.9, particle-fluid density ratio = 100, particle diameter = 6 lattice units.


Heat and mass transfers in particulate-laden flows are important for many of their intended applications. We are developing and using a random-walk based method, coupled with detailed solution of hydrodynamics, to simulate heat and mass transfer in suspensions. The motion of particles in the fluid could agigate the fluid enough to achieve tens to hundreds of fold of enhanced transfer especially when Pe is large. It is therefore important to be able to quantify such effects using detailed simulations and this constitutes the third direction of interest to us in the area of suspension dynamics and transport.




Our studies on porous media focus on pore- and core-scale phenomena. The objective of these studies is to improve existing flow and transport models (e.g. relative permeability) for continuum-scale simulations. In pore-scale studies, we use a range of computational and experimental tools. We conduct direct numerical simulations of fluid flow and transport through porous media using lattice Boltzmann methods and random-walk particle tracking. We build model structures of porous media using a Voronoi-diagram based method. Experimentally, we use transparent microfluidic devices to visualize flow and transport on the scale of micrometers, observations made from which are used to direct or calibrate pore-scale simulations. High-performance computing and microfluidics reaching the size of centimeter give us ability to upscale flow and transport to the core level. These results are compared to real core experiments for full validation.


Movie - Dynamics of water-oil displacement in a microfluidic porous medium without surfactant.


Movie - Dynamics of water-oil displacement in the same microfluidic porous medium with surfactant.


In pores of "conventional" sizes (above O(100) nm), transport of fluids is largely governed by hydrodynamics. In pores of smaller sizes, which are abundant in tight oil and gas reservoirs termed "unconventional", surface effects can moderate Navier-Stokes hydrodynamics and, depending on pore size and fluid-surface interactions, can even become dominant. Transports that cannot be explained and modeled by Navier-Stokes hydrodynamics alone include slippage of gases and chromatographical separation of oils in nanoporous media. These phenomena are critical for reservoir-scale movements of hydrocarbons but little is known on how to model them. Thus, specific rock-fluid interactions and their influences on transport must be measured, the results of which could then lead to predictive models. We approach these problems by using a combination of (mini)core experiments and pore-scale molecular and hydrodynamic simulations.



Filtration of C10/C17 binary mixture through a Niobrara shale sample.


Transport of soluble (i.e. stable dispersion) macromolecules or colloid-size particles, such as hydrated polymer coils, silica-alumina clay particles, and asphaltene aggregates, through porous media is often complicated by interactions between macromolecules / colloids and surfaces of porous media. These interactions are electrostatic in nature and can be, as a leading approximation, modeled as a combination of monopole (i.e. charge) and dipole (i.e. van der Waals) potentials. Particle-tracking method is an efficient way to simulate transport of such "particulates" in porous media. However, interactions between these particles and surfaces of porous media under specific liquid environment (oil or water) must be correctly modeled or measured and then incorporated into a particle-tracking simulator. Research along this direction is being carried out in collaboration with Dr. Ning Wu's group in Chemical Engineering.


breakthrough curve

Breakthrough curve of colloidal nanoparticles through a packed bed that contains a small fraction of collectors - comparison between microfluidic experiment and pore-scale simulation.

Many of our pore-scale simulations use MPI-based parallel computing. The size of computational domain is critical for porous media flow simulations, because porous media often contain pores of different sizes, and simulations must incorporate adequate number of pores of varying sizes. Simulations therefore must admit large domains and parallel computing is indispensable. Our pore-scale simulation programs have all been parallelized.

Liquid-to-gas and gas-to-liquid phase transitions are important dynamics of oil and condensate gas reservoirs. In unconventional reservoirs, capillary pressure between gas and liquid phases can be as high as several atms. As a result, onset of gas phase (bubble point) in an oil reservoir during depressurization may be delayed, and liquid dropout (dew point) in a condensate gas reservoir may be expedited. In modeling, this process is complicated by the physics that both gas and liquid can develop adsorbed layers that attach to the surface of rock, affecting the material balance. We write multiphase equilibrium models to describe the process of depressurization and phase transition in porous media. In addition to modeling, we also conduct core experiments with condensate gas, using a simple oscillation method, to generate data for validation.



Initial Condition - Liquid fully saturates the porous medium.



As piston moves to the right, pressure is reduced and gas is generated in the porous medium.



Critical gas saturation is reached and gas is produced from the porous medium.