Microscopic-based fluid flow simulation of invasion on a two-dimensional lattice. II. Mobilization and cohesion

An algorithm for modeling secondary invasion processes in porous media is presented. Mobilization of trapped defender fluid is accomplished through interfacial interaction rules. Cohesive forces are also included within the defender phase. A series of simulation runs are performed using two-dimensional lattices and examined to determine optimal conditions for secondary invasions that sweep the trapped defender phase from the porous medium.     

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.     

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.     

On hydraulic permeability of random packs of monodisperse spheres: Direct flow simulations versus correlations

Hydraulic permeability is studied in porous media consisting of randomly distributed monodisperse spheres by means of computational fluid dynamics (CFD) simulations. The packing of spheres is generated by inserting a certain number of nonoverlapping spherical particles inside a cubic box at both low and high packing fractions using proper algorithms. Fluid flow simulations are performed within the interparticulate porous space by solving Navier-Stokes equations in a low-Reynolds laminar flow regime. The hydraulic permeability is calculated from the Darcy equation once the mean values of velocity and pressure gradient are calculated across the particle packing. The simulation results for the pressure drop across the packing are verified by the Ergun equation for the lower range of porosities (<0.75), and the Stokes equation for higher porosities (∼1). Using the results of simulations, the effects of porosity and particle diameters on the hydraulic permeability are investigated. Simulations precisely specified the range of applicability of empirical or semi-empirical correlations for hydraulic permeability, namely the Carman-Kozeny, Rumpf-Gupte, and Howells-Hinch formulas. The number of spheres in the model is gradually decreased from 2000 to 20 to discover the finite-size effect of pores on the hydraulic permeability of spherical packing, which has not been clearly addressed in the literature. In addition, the scale dependence of hydraulic permeability is studied via simulations of the packing of spheres shrunk to lower scales. The results of this work not only reveal the validity range of the aforementioned correlations, but also show the finite-size effect of pores and the scale-independence of direct CFD simulations for hydraulic permeability.     

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.     

Rate processes on fractals: Theory,simulations, and experiments

Heterogeneous kinetics are shown to differ drastically from homogeneous kinetics. For the elementary reaction A + A products we show that the diffusion-limited reaction rate is proportional tot ^– h[A]² or to [A]^x, whereh=1- d _s/2, X=1+2/d _s =(h-2)(h-1), andd _s is the effective spectral dimension. We note that ford = d _s =1, h =1/2 andX = 3, for percolating clustersd _s = 4/3,h = 1/3 andX = 5/2, while for dust ds <1, 1 >h > 1/2 and >X > 3. Scaling arguments, supercomputer simulations and experiments give a consistent picture. The interplay of energetic and geometric heterogeneity results in fractal-like kinetics and is relevant to excitation fusion experiments in porous membranes, films, and polymeric glasses. However, in isotopic mixed crystals, the geometric fractal nature (percolation clusters) dominates.     

Numerical simulation of heat transfer and fluid flow of Water-CuO Nanofluid in a sinusoidal channel with a porous medium

This study has compared the convection heat transfer of Water-based fluid flow with that of Water-Copper oxide (CuO) nanofluid in a sinusoidal channel with a porous medium. The heat flux in the lower and upper walls has been assumed constant, and the flow has been assumed to be two-dimensional, steady, laminar, and incompressible. The governing equations include equations of continuity, momentum, and energy. The assumption of thermal equilibrium has been considered between the porous medium and the fluid. The effects of the parameters, Reynolds number and Darcy number on the thermal performance of the channel, have been investigated. The results of this study show that the presence of a porous medium in a channel, as well as adding nanoparticles to the base fluid, increases the Nusselt number and the convection heat transfer coefficient. Also the results show that As the Reynolds number increases, the temperature gradient increases. In addition, changes in this parameter are greater in the throat of the flow than in convex regions due to changes in the channel geometry. In addition, porous regions reduce the temperature difference, which in turn increases the convective heat transfer coefficient.     

Tunable-slip boundaries for coarse-grained simulations of fluid flow

On the micro- and nanoscale, classical hydrodynamic boundary conditions such as the no-slip condition no longer apply. Instead, the flow profiles exhibit "slip" at the surface, which is characterized by a finite slip length (partial slip). We present a new, systematic way of implementing partial-slip boundary conditions with arbitrary slip length in coarse-grained computer simulations. The main idea is to represent the complex microscopic interface structure by a spatially varying effective viscous force. An analytical equation for the resulting slip length can be derived for planar and for curved surfaces. The comparison with computer simulations of a DPD (dissipative particle dynamics) fluid shows that this expression is valid from full slip to no slip.     

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.     

The effects of porosity and permeability on fluid flow and heat transfer of multi walled carbon nano-tubes suspended in oil (MWCNT/Oil nano-fluid) in a microchannel filled with a porous medium

The forced convection heat transfer and laminar flow in a two-dimensional microchannel filled with a porous medium is numerically investigated. The nano-particles which have been used are multi walled carbon nano-tubes (MWCNT) suspended in oil as the based fluid. The assumption of no-slip condition between the base fluid and nano-particles as well as the thermal equilibrium between them allows us to study the nanofluid in a single phase. The nanofluid flow through the microchannel has been modeled using the Darcy–Forchheimer equation. It is also assumed that there is a thermal equilibrium between the solid phase and the nanofluid for energy transfer. The walls of the microchannel are under the influence of a fluctuating heat flux. Also, the slip velocity boundary condition has been assumed along the walls. The effects of Darcy number, porosity and slip coefficients and Reynolds number on the velocity and temperature profiles and Nusselt number will be studied in this research.     

超声波对多孔介质中两相流动的影响

近年来，超声技术已被应用于采油工程中，在油井解堵，水井增注等方面发挥了重要的作用。     

Lattice-gas simulation of escaping pedestrian flow in corridor

A lattice-gas model of biased-random walkers is used to simulate the escaping pedestrian flow under the open boundary condition in corridor. Given that the total number of people is unchanging, we have studied the evolution of pedestrian flow by varying parameters of system size. Relationships between parameters of system size and the transition time are discussed in this paper. Scaling behaviour is found as follows: the transition time t_c scales as t_c∝W^{-0.85±0.04}, and t_c∝D, where W is the width of corridor, and D is the strength of drift. However, the other parameters are found to have little influence on the transition time.     

Gravity driven instabilities in miscible non-Newtonian fluid displacements in porous media

Gravity driven instabilities in model porous packings of 1 mm diameter spheres are studied by comparing the broadening of the displacement front between fluids of slightly different densities in stable and unstable configurations. Water, water–glycerol and water–polymer solutions are used to vary independently viscosity and molecular diffusion and study the influence of shear-thinning properties. Both injected and displaced solutions are identical but for a different concentration of NaNO₃ salt used as an ionic tracer and to introduce the density contrast. Dispersivity in stable configuration increases with polymer concentration – as already reported for double porosity packings of porous grains. Gravity-induced instabilities are shown to develop below a same threshold Péclet number Pe for water and water–glycerol solutions of different viscosities and result in considerable increases of the dispersivity. Measured threshold Pe values decrease markedly on the contrary with polymer concentration. The quantitative analysis demonstrates that the development of the instabilities is controlled by viscosity through a characteristic gravity number G (ratio between hydrostatic and viscous pressure gradients). A single threshold value of G accounts for results obtained on Newtonian and non-Newtonian solutions.     

多孔介质中卡森流体的分形分析

研究了非牛顿流体中的卡森流体在多孔介质中的流动特性.基于服从分形分布的弯曲毛细管束模型,运用分形几何理论推导出了该流体在多孔介质中流动的流量、流速、启动压力梯度和有效渗透率的分形解析解.模型中的每一个参数都有明确的物理意义,它将卡森流体在多孔介质中的流动特性与多孔介质的微结构参数有机联系起来.文中给出了卡森流体的流速、启动压力梯度和有效渗透率随着各影响因素的变化趋势,并进行了讨论.所得分形模型可以更深刻地理解卡森流体在多孔介质中流动的内在物理机理. 关键词： 多孔介质 卡森流体 分形     

MHD oscillatory flow on free convection–radiation through a porous medium with constant suction velocity

Unsteady two-dimensional hydromagnetic free convection and thermal radiation flow of an electrically conducting viscous-incompressible fluid, through a highly porous medium bounded by a vertical plane surface of constant temperature are presented. The Rosseland diffusion approximation is used to describe the radiative heat flux in the energy equation. Expressions for the velocity and temperature are obtained. The free-stream velocity of the fluid vibrates about a mean constant value and the surface absorbs the fluid with constant velocity. Effects of varying R (radiative parameter), G (Grashof number), k′ (permeability of the porous medium) and M (magnetic parameter upon the velocity field and the effect of varying R and Pr (Prandtl number) on the temperature are discussed.     

Large-eddy simulations of temporally accelerating turbulent channel flow

The unsteady turbulent channel flow subject to the temporal acceleration is considered in this study. Large-eddy simulations were performed to study the response of the turbulent flow to the temporal acceleration. The simulations were started with the fully developed turbulent channel flow at an initial Reynolds number of Re₀ = 3500 (based on the channel half-height and the bulk-mean velocity), and then a constant temporal acceleration was applied. During the acceleration, the Reynolds number of the channel flow increased linearly from the initial Reynolds number to the final Reynolds number of Re₁ = 22,600. The effect of grid resolution, domain size, time step size on the simulation results was assessed in a preliminary study using simulations of the accelerating turbulent flow as well as simulations of the steady turbulent channel flow at various Reynolds numbers. Simulation parameters were carefully chosen from the preliminary study to ascertain the accuracy of the simulation. From the accelerating turbulent flow simulations, the delays in the response of various flow properties to the temporal acceleration were measured. The distinctive features of the delays responsible for turbulence production, energy redistribution, and radial propagation were identified. Detailed turbulence statistics including the wall shear stress response during the acceleration were examined. The results reveal the changes in the near-wall structures during the acceleration. A self-sustaining mechanism of turbulence is proposed to explain the response of the turbulent flow to the temporal acceleration. Although the overall flow characteristics are similar between the channel and pipe flows, some differences were observed between the two flows.     

Analysis of magnetohydrodynamic flow of a fractional viscous fluid through a porous medium

The aim of this paper is to investigate the exact general solutions of the incompressible viscous fluid flow by using the time-fractional Caputo–Fabrizio derivative. The flow of the fluid is subject to the motion of a plane wall, embedded in a porous medium under the influence of magnetic field. The corresponding non-dimensional governing fractional differential equation with appropriate initial and boundary conditions is solved by means of integral transforms namely, Laplace and Fourier transforms. Solutions are expressed as a sum of steady and transient parts, for the sinusoidal oscillations of the plane wall. The influence of involved physical parameters are discussed graphically. Specifically, it has been observed that the effective permeability K_eff reduces the time taken to reach the steady state.     

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.     

Diffusion of oriented particles in porous media

Diffusion of particles in porous media often shows subdiffusive behavior. Here, we analyze the dynamics of particles exhibiting an orientation. The features we focus on are geometrical restrictions and the dynamical consequences of the interactions between the local surrounding structure and the particle orientation. This interaction can lead to particles getting temporarily stuck in parts of the structure. Modeling this interaction by a particular random walk dynamics on fractal structures we find that the random walk dimension is not affected while the diffusion constant shows a variety of interesting and surprising features.     

Magnetic resonance velocity imaging of liquid and gas two-phase flow in packed beds

Single-phase liquid flow in porous media such as bead packs and model fixed bed reactors has been well studied by MRI. To some extent this early work represents the necessary preliminary research to address the more challenging problem of two-phase flow of gas and liquid within these systems. In this paper, we present images of both the gas and liquid velocities during stable liquid–gas flow of water and SF₆ within a packing of 5 mm spheres contained within columns of diameter 40 and 27 mm; images being acquired using ¹H and ¹⁹F observation for the water and SF₆, respectively. Liquid and gas flow rates calculated from the velocity images are in agreement with macroscopic flow rate measurements to within 7% and 5%, respectively. In addition to the information obtained directly from these images, the ability to measure liquid and gas flow fields within the same sample environment will enable us to explore the validity of assumptions used in numerical modelling of two-phase flows.