A cumulant analysis for non-Gaussian displacement distributions in Newtonian and non-Newtonian flows through porous media

We use displacement encoding pulsed field gradient (PFG) nuclear magnetic resonance to measure Fourier components S(q) of flow displacement distributions P(zeta) with mean displacement (zeta) for Newtonian and non-Newtonian flows through rocks and bead packs. Displacement distributions are non-Gaussian; hence, there are finite terms above second order in the cumulant expansion of ln(S(q)). We describe an algorithm for an optimal self-consistent cumulant analysis of data, which can be used to obtain the first three (central) moments of a non-Gaussian P(zeta), with error bars. The analysis is applied to Newtonian and non-Newtonian flows in rocks and beads. Flow with shear-thinning xanthan solution produces a 15.6+/-2.3% enhancement of the variance sigma(2) of displacement distributions when compared to flow experiments with water.     

Direct pore-level modeling of incompressible fluid flow in porous media

We present a dynamic particle-based model for direct pore-level modeling of incompressible viscous fluid flow in disordered porous media. The model is capable of simulating flow directly in three-dimensional high-resolution micro-CT images of rock samples. It is based on moving particle semi-implicit (MPS) method. We modify this technique in order to improve its stability for flow in porous media problems. Using the micro-CT image of a rock sample, the entire medium, i.e., solid and fluid, is discretized into particles. The incompressible Navier–Stokes equations are then solved for each particle using the MPS summations. The model handles highly irregular fluid–solid boundaries effectively. An algorithm to split and merge fluid particles is also introduced. To handle the computational load, we present a parallel version of the model that runs on distributed memory computer clusters. The accuracy of the model is validated against the analytical, numerical, and experimental data available in the literature. The validated model is then used to simulate both unsteady- and steady-state flow of an incompressible fluid directly in a representative elementary volume (REV) size micro-CT image of a naturally-occurring sandstone with 3.398 μm resolution. We analyze the quality and consistency of the predicted flow behavior and calculate absolute permeability using the steady-state flow rate.     

Quantum effects on the Rayleigh-Taylor instability of stratified fluid/plasma through porous media

Quantum effect on Rayleigh-Taylor instability of stratified plasma layer through a porous medium are investigated. The linear growth rate is obtained analytically and is analyzed. In the presence of quantum effect, both the porosity of porous medium and the medium permeability has different influence on the coup point () for stability, but they do not have influence on the critical point () for stability. The quantum effect plays the principal role of the complete stability case for the system considered.     

Elliptic random-walk equation for suspension and tracer transport in porous media

We propose a new approach to transport of the suspensions and tracers in porous media. The approach is based on a modified version of the continuous time random walk (CTRW) theory. In the framework of this theory we derive an elliptic transport equation. The new equation contains the time and the mixed dispersion terms expressing the dispersion of the particle time steps. The properties of the new equation are studied and the fundamental analytical solutions are obtained. The solution of the pulse injection problem describing a common tracer injection experiment is studied in greater detail. The new theory predicts delay of the maximum of the tracer, compared to the velocity of the flow, while its forward “tail” contains much more particles than in the solution of the classical parabolic (advection-dispersion) equation. This is in agreement with the experimental observations and predictions of the CTRW theory.     

Density profile and flow of miscible fluid with dissimilar constituent masses

A computer simulation model is used to study the density profile and flow of a miscible gaseous fluid mixture consisting of differing constituent masses (m_A=m_B/3) through an open matrix. The density profile is found to decay with the height ∝exp (−m_A(B)h), consistent with the barometric height law. The flux density shows a power-law increase ∝(p_c−p)^μ with μ2.3 at the porosity 1−p above the pore percolation threshold 1−p_c.     

Conditional methods for continuum reacting flows in porous media

A special version of conditional moment closure—PCMC—is suggested for modeling reacting flows in porous media. The model involves conditioning on a special tracer scalar, which is introduced to characterize scalar transport in the gaseous phase. (i.e., for the flow in the interparticle space or in the pores). The model accounts for interparticle variations of species concentrations and emulates diffusion in the interparticle space. Special boundary conditions that are consistent with conventional conditions at the phase interface are obtained for the PCMC model. The model is tested against complete direct simulation of a reacting flow in porous media with favourable results.     

An efficient S-DDM iterative approach for compressible contamination fluid flows in porous media

In this paper, we develop an efficient splitting domain decomposition method (S-DDM) for compressible contamination fluid flows in porous media over multiple block-divided sub-domains by combining the non-overlapping domain decomposition, splitting, linearization and extrapolation techniques. The proposed S-DDM iterative approach divides the large domain into multiple block sub-domains. In each time interval, the S-DDM scheme is applied to solve the water head equation, in which an efficient local multilevel scheme is used for computing the values of water head on the interfaces of sub-domains, and the splitting implicit scheme is used for computing the interior values of water head in sub-domains; and the S-DDM scheme is then proposed to solve the concentration equation by combining the upstream volume technique. Numerical experiments are performed and analyzed to illustrate the efficiency of the S-DDM iterative approach for simulating compressible contamination fluid flows in porous media. The developed method takes the excellent attractive advantages of both the non-overlapping domain decomposition and the splitting technique, and reduces computational complexities, large memory requirements and long computational durations.     

Transition in the propagation mechanism during flame acceleration in porous media

The combustion of stoichiometric hydrogen-air at various initial pressures was investigated in a 7.62 cm square cross-section channel filled with 1.27 cm diameter beads. The flame time-of-arrival and pressure time history along the channel were obtained by ionization probes and piezoelectric pressure transducers. Flame acceleration was found to be very rapid, e.g. at an initial pressure of 45 kPa the flame achieves a velocity of over 600 m/s in roughly 0.3 m. It was determined that at this high speed a well defined planar shock wave precedes a thick reaction zone. It was also shown that there is a transition in the flame propagation mechanism, similar to that observed in an obstacle laden channel [G. Ciccarelli and C. Johansen, The role of shock-flame interactions on flame acceleration in an obstacle laden channel, Proc. 22nd International Colloquium on the Dynamics of Explosions and Reactive Systems, Minsk, 2009]. By varying the initial pressure of the mixture, changes in the axial location of the transition between combustion propagation regimes was also observed. A soot foil technique was used to identify the transition in the propagation mechanism, as well as to provide information concerning the local flow field around the beads and the overall average flow direction.     

A numerical simulation method for dissolution in porous and fractured media

We describe an algorithm for simulating reactive flows in porous media, in which the pore space is mapped explicitly. Chemical reactions at the solid–fluid boundaries lead to dissolution (or precipitation), which makes it necessary to track the movement of the solid–fluid interface during the course of the simulation. We have developed a robust algorithm for constructing a piecewise continuous (C₁) surface, which enables a rapid remapping of the surface to the grid lines. The key components of the physics are the Navier–Stokes equations for fluid flow in the pore space, the convection–diffusion equation to describe the transport of chemical species, and rate equations to model the chemical kinetics at the solid surfaces. A lattice-Boltzmann model was used to simulate fluid flow in the pore space, with linear interpolation at the solid boundaries. A finite-difference scheme for the concentration field was developed, taking derivatives along the direction of the local fluid velocity. When the flow is not aligned with the grid this leads to much more accurate convective fluxes and surface concentrations than a standard Cartesian template. A robust algorithm for the surface reaction rates has been implemented, avoiding instabilities when the surface is close to a grid point. We report numerical tests of different aspects of the algorithm and assess the overall convergence of the method.     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.