Onset of Double-Diffusive Reaction–Convection in an Anisotropic Rotating Porous Layer

The linear and non-linear stability of a rotating double-diffusive reaction–convection in a horizontal anisotropic porous layer subjected to chemical equilibrium on the boundaries is investigated considering a Darcy model that includes the Coriolis term. The effect of Taylor number, mechanical, and thermal anisotropy parameters, reaction rate, solute Rayleigh number, Lewis number, and normalized porosity on the stability of the system is investigated. We find that the Taylor number has a stabilizing effect, chemical reaction may be stabilizing or destabilizing and that the anisotropic parameters have significant influence on the stability criterion. The effect of various parameters on the stationary, oscillatory, and finite-amplitude convection is shown graphically. A weak nonlinear theory based on the truncated representation of Fourier series method is used to find the finite amplitude Rayleigh number and heat and mass transfer.     

Thermal Convection in a Cattaneo–Fox Porous Material with Guyer–Krumhansl Effects

The Guyer–Krumhansl equation is coupled with the Cattaneo–Fox law for the temperature and heat flux fields to study thermal convection in a fluid-saturated Darcy porous material. In particular the effects of the Guyer–Krumhansl terms on oscillatory convection are studied. It is found that for a certain range of the Guyer–Krumhansl coefficient stationary convection occurs while changing the range results in oscillatory convection. Numerical results quantify this effect.     

Darcy–Carreau Model and Nonlinear Natural Convection for Pseudoplastic and Dilatant Fluids in Porous Media

Transport in Porous Media - The linear and weakly nonlinear stability analyses are carried out to study instabilities in Darcy–Bénard convection for non-Newtonian inelastic fluids. The...     

Thermal Dispersion–Radiation Effects on Non-Darcy Natural Convection in a Fluid Saturated Porous Medium

The effects of thermal dispersion and thermal radiation on the non-Darcy natural convection over a vertical flat plate in a fluid saturated porous medium are studied. Forchheimer extension is considered in the flow equations. The coefficient of thermal diffusivity has been assumed to be the sum of molecular diffusivity and the dispersion thermal diffusivity due to mechanical dispersion. Rosseland approximation is used to describe the radiative heat flux in the energy equation. Similarity solution for the transformed governing equations is obtained. Numerical results for the details of the velocity and temperature profiles which are shown on graphs have been presented. The combined effect of thermal dispersion and thermal radiation, for the two cases Darcy and non-Darcy porous medium, on the heat transfer rate which are entered in tables is discussed.     

Entropy Generation in Double-Diffusive Convection in a Square Porous Cavity using Darcy–Brinkman Formulation

The article reports a numerical study of entropy generation in double-diffusive convection through a square porous cavity saturated with a binary perfect gas mixture and submitted to horizontal thermal and concentration gradients. The analysis is performed using Darcy–Brinkman formulation with the Boussinesq approximation. The set of coupled equations of mass, momentum, energy and species conservation are solved using the control volume finite-element method. Effects of the Darcy number, the porosity and the thermal porous Rayleigh number on entropy generation are studied. It was found that entropy generation considerably depends on the Darcy number. Porosity induces the increase of entropy generation, especially at higher values of thermal porous Rayleigh number.     

Influence of thermal radiation and Joule heating in the Eyring–Powell fluid flow with the Soret and Dufour effects

A two-dimensional magnetohydrodynamic boundary layer flow of the Eyring–Powell fluid on a stretching surface in the presence of thermal radiation and Joule heating is analyzed. The Soret and Dufour effects are taken into account. Partial differential equations are reduced to a system of ordinary differential equations, and series solutions of the resulting system are derived. Velocity, temperature, and concentration profiles are obtained. The skin friction coefficient and the local Nusselt and Sherwood numbers are computed and analyzed.     

MHD mixed convection–radiation interaction along a permeable surface immersed in a porous medium in the presence of Soret and Dufour’s Effects

This work is focused on the numerical modeling of steady, laminar, heat and mass transfer by MHD mixed convection from a semi-infinite, isothermal, vertical and permeable surface immersed in a uniform porous medium in the presence of thermal radiation and Dufour and Soret effects. A mixed convection parameter for the entire range of free-forced-mixed convection is employed and the governing equations are transformed into non-similar equations. These equations are solved numerically by an efficient, implicit, iterative, finite-difference scheme. The obtained results are checked against previously published work on special cases of the problem and are found to be in excellent agreement. A parametric study illustrating the influence of the thermal radiation coefficient, magnetic field, porous medium inertia parameter, concentration to thermal buoyancy ratio, and the Dufour and Soret numbers on the fluid velocity, temperature and concentration as well as the local Nusselt and the Sherwood numbers is conducted. The obtained results are shown graphically and the physical aspects of the problem are discussed.     

Beyond Kozeny–Carman: Predicting the Permeability in Porous Media

Transport in Porous Media - Various processes such as heterogeneous reactions or biofilm growth alter a porous medium’s underlying geometric structure. This significantly affects its...     

The Onset of Prandtl–Darcy–Prats Convection in a Horizontal Porous Layer

We consider the effect of finite Prandtl–Darcy numbers of the onset of convection in a porous layer heated isothermally from below and which is subject to a horizontal pressure gradient. A dispersion relation is found which relates the critical Darcy–Rayleigh number and the induced phase speed of the cells to the wavenumber and the imposed Péclet and Prandtl–Darcy numbers. Exact numerical solutions are given and these are supplemented by asymptotic solutions for both large and small values of the governing nondimensional parameters. The classical value of the critical Darcy–Rayleigh number is $4\pi ^2$ 4 π 2 , and we show that this value increases whenever the Péclet number is nonzero and the Prandtl–Darcy number is finite simultaneously. The corresponding wavenumber is always less than $\pi $ π and the phase speed of the convection cells is always smaller than the background flux velocity.     

Diffusion–Convection in a Porous Medium with Impervious Inclusions at Low Flow Rates

The aim of this paper is to develop a macroscopic model for the transport of a passive solute, by diffusion and convection, in a heterogeneous medium consisting of impervious solids periodically distributed in a porous matrix. In the porous part, the flow is described by Darcy's law. Attempt is made to derive the macroscopic equation governing the average concentration field in the equivalent macroscopic medium and the macroscopic transport parameters. The analysis is conducted in the case when convection and diffusion are of the same order of magnitude at the macroscopic level, that is, when the Péclet number is of order 1. The proposed macroscopic model is obtained using the homogenization method for periodic structures with a double scale asymptotic expansion, in which the small parameter is the ratio between the two characteristic lengths l (the period scale of the impervious bodies distribution) and L (the scale of the macroscopic sample). The macroscopic parameters, which characterize the multiporous medium, depend solely on the transport parameters in the porous matrix and on the geometry of the impervious inclusions without any phenomenological assumption. Numerical computations are performed using a finite element method for several geometries of the solid inclusions, in two- and three-dimensional cases.     

Impact of the Buoyancy–Viscous Force Balance on Two-Phase Flow in Layered Porous Media

Motivated by geological carbon storage and hydrocarbon recovery, the effect of buoyancy and viscous forces on the displacement of one fluid by a second immiscible fluid, along parallel and dipping layers of contrasting permeability, is characterized using five independent dimensionless numbers and a dimensionless storage or recovery efficiency. Application of simple dimensionless models shows that increased longitudinal buoyancy effects increase storage efficiency by reducing the distance between the leading edges of the injected phase in each layer and decreasing the residual displaced phase saturation behind the leading edge of the displacing phase. Increased transverse buoyancy crossflow increases storage efficiency if it competes with permeability layering effects, but reduces storage efficiency otherwise. When both longitudinal and transverse buoyancy effects are varied simultaneously, a purely geometrical dip angle group defines whether changes in storage efficiency are dominated by changes in the longitudinal or transverse buoyancy effects. In the limit of buoyancy-segregated flow, we report an equivalent, unidimensional flow model which allows rapid prediction of storage efficiency. The model presented accounts for both dip and layering, thereby generalizing earlier work which accounted for each of these but not both together. We suggest that the predicted storage efficiency can be used to compare and rank geostatistical realizations, and complements earlier heterogeneity measures which are applicable in the viscous limit.     

Instability of the Benard–Marangoni Convection in a Porous Layer Affected by a Non-Vertical Magnetic Field

The onset of the Benard–Marangoni convection in a horizontal porous layer permeated by a magnetohydrodynamic fluid with a nonlinear magnetic permeability is examined. The porous layer is assumed to be governed by the Brinkman model; it is bounded by a rigid surface from below and by a non-deformable free surface from above and subjected to a non-vertical magnetic field. The critical effective Marangoni number and the critical Rayleigh number are obtained for different values of the effective Darcy number, Biot number, Chandrasekhar number, nonlinear magnetic parameter, and angle from the vertical axis for the cases of stationary convection and overstability. The related eigenvalue problem is solved by using the first-order Chebyshev polynomial method.     

The Falkner–Skan Wedge Flows of Power-Law Fluids Embedded in a Porous Medium

The non-Darcy flow characteristics of power-law non-Newtonian fluids past a wedge embedded in a porous medium have been studied. The governing equations are converted to a system of first-order ordinary differential equations by means of a local similarity transformation and have been solved numerically, for a number of parameter combinations of wedge angle parameter m, power-law index of the non-Newtonian fluids n, first-order resistance A and second-order resistance B, using a fourth-order Runge–Kutta integration scheme with the Newton–Raphson shooting method. Velocity and shear stress at the body surface are presented for a range of the above parameters. These results are also compared with the corresponding flow problems for a Newtonian fluid. Numerical results show that for the case of the constant wedge angle and material parameter A, the local skin friction coefficient is lower for a dilatant fluid as compared with the pseudo-plastic or Newtonian fluids.     

A Three-Dimensional Model of Two-Phase Flows in a Porous Medium Accounting for Motion of the Liquid–Liquid Interface

A new three-dimensional hydrodynamic model for unsteady two-phase flows in a porous medium, accounting for the motion of the interface between the flowing liquids, is developed. In a minimum number of interpretable geometrical assumptions, a complete system of macroscale flow equations is derived by averaging the microscale equations for viscous flow. The macroscale flow velocities of the phases may be non-parallel, while the interface between them is, on average, inclined to the directions of the phase velocities, as well as to the direction of the saturation gradient. The last gradient plays a specific role in the determination of the flow geometry. The resulting system of flow equations is a far generalization of the classical Buckley–Leverett model, explicitly describing the motion of the interface and velocity of the liquid close to it. Apart from propagation of the two liquid volumes, their expansion or contraction is also described, while rotation has been proven negligible. A detailed comparison with the previous studies for the two-phase flows accounting for propagation of the interface on micro- and macroscale has been carried out. A numerical algorithm has been developed allowing for solution of the system of flow equations in multiple dimensions. Sample computations demonstrate that the new model results in sharpening the displacement front and a more piston-like character of displacement. It is also demonstrated that the velocities of the flowing phases may indeed be non-collinear, especially at the zone of intersection of the displacement front and a zone of sharp permeability variation.     

Influence of Viscous Fingering on Dynamic Saturation–Pressure Curves in Porous Media

We report on results from primary drainage experiments on quasi-two-dimensional porous models. The models are transparent, allowing the displacement process and structure to be monitored in space and time during primary drainage experiments carried out at various speeds. By combining detailed information on the displacement structure with global measurements of pressure, saturation and the capillary number Ca, we obtain a scaling relation relating pressure, saturation, system size and capillary number. This scaling relation allows pressure–saturation curves for a wide range of capillary numbers to be collapsed on the same master curve. We also show that in the case of primary drainage, the dynamic effect in the capillary pressure–saturation relationship observed on partially water saturated soil samples might be explained by the combined effect of capillary pressure along the invasion front of the gaseous phase, and pressure changes caused by viscous effects in the wetting fluid phase.     

Mixed Convection in a Darcy–Brinkman Porous Medium with a Constant Convective Thermal Boundary Condition

The characteristics of the boundary layer flow past a plane surface adjacent to a saturated Darcy–Brinkman porous medium are investigated in this paper. The flow is driven by an external free stream moving with constant velocity. The surface is heated with a convective boundary condition with constant heat transfer coefficient. The problem is non-similar and is investigated numerically by a finite difference method. The problem is governed by four non-dimensional parameters, that is, the convective Darcy number, the convective Grashof number, the Prandtl number, and the axial distance along the plate. The influence of these parameters on the results is investigated, and the results are presented in tables and figures. The Darcy term and the Grashof term in the momentum equation contradict each other and this contradiction makes the problem complicated. However, the wall shear stress and the wall temperature increase continuously along the plate and the wall temperature always tends to 1.     

Three-Dimensional Lattice Boltzmann Simulations of Single-Phase Permeability in Random Fractal Porous Media with Rough Pore–Solid Interface

Single-phase permeability k has intensively been investigated over the past several decades by means of experiments, theories and simulations. Although the effect of surface roughness on fluid flow and permeability in single pores and fractures as well as networks of fractures was studied previously, its influence on permeability in a random mass fractal porous medium constructed of pores of different sizes remained as an open question. In this study, we, therefore, address the effect of pore–solid interface roughness on single-phase flow in random fractal porous media. For this purpose, we apply a mass fractal model to construct porous media with a priori known mass fractal dimensions \(2.579 \le D_{\mathrm{m}} \le 2.893\) characterizing both solid matrix and pore space. The pore–solid interface of the media is accordingly roughened using the Weierstrass–Mandelbrot approach and two parameters, i.e., surface fractal dimension \(D_{\mathrm{s}}\) and root-mean-square (rms) roughness height. A single-relaxation-time lattice Boltzmann method is applied to simulate single-phase permeability in the corresponding porous media. Results indicate that permeability decreases sharply with increasing \(D_{\mathrm{s}}\) from 1 to 1.1 regardless of \(D_{\mathrm{m}}\) value, while k may slightly increase or decrease, depending on \(D_{\mathrm{m}}\), as \(D_{\mathrm{s}}\) increases from 1.1 to 1.6.     

Pore-Scale Imaging and Modelling of Reactive Flow in Evolving Porous Media: Tracking the Dynamics of the Fluid–Rock Interface

Transport in Porous Media - Fluid–mineral and fluid–rock interfaces are key parameters controlling the reactivity and fate of fluids in reservoir rocks and aquifers. The interface...     

Effect of Conduction in Bottom Wall on Darcy–Bénard Convection in a Porous Enclosure

Conjugate natural convection-conduction heat transfer in a square porous enclosure with a finite-wall thickness is studied numerically in this article. The bottom wall is heated and the upper wall is cooled while the verticals walls are kept adiabatic. The Darcy model is used in the mathematical formulation for the porous layer and the COMSOL Multiphysics software is applied to solve the dimensionless governing equations. The governing parameters considered are the Rayleigh number (100 ≤ Ra ≤ 1000), the wall to porous thermal conductivity ratio (0.44 ≤ K _r ≤ 9.90) and the ratio of wall thickness to its height (0.02 ≤ D ≤ 0.4). The results are presented to show the effect of these parameters on the heat transfer and fluid flow characteristics. It is found that the number of contrarotative cells and the strength circulation of each cell can be controlled by the thickness of the bottom wall, the thermal conductivity ratio and the Rayleigh number. It is also observed that increasing either the Rayleigh number or the thermal conductivity ratio or both, and decreasing the thickness of the bounded wall can increase the average Nusselt number for the porous enclosure.