Non-Darcian and Anisotropic Effects on Free Convection in a Porous Enclosure

The free convective flow and heat transfer, within the framework of Boussinesq approximation, in an anisotropic fluid filled porous rectangular enclosure subjected to end-to-end temperature difference have been investigated using Brinkman extended non-Darcy flow model. The studies involve simultaneous consideration of hydrodynamic and thermal anisotropy. The flow and temperature fields in general are governed by, Ra, the Rayleigh number, AR, the aspect ratio of the slab, K*, the permeability ratio and k*, the thermal conductivity ratio, and Da, Darcy number. Numerical solutions employing the successive accelerated replacement (SAR) scheme have been obtained for 100 ≤ Ra ≤ 1000, 0.5 ≤ AR ≤ 5, 0.5 ≤ K* ≤ 5, 0.5 ≤ k* ≤ 5, and 0 ≤ Da ≤ 0.1. It has been found that [`(Nu)]{\overline {Nu}}, average Nusselt number increases with increase in K* and decreases as k* increases. However, the magnitude of the change in [`(Nu)]{\overline {Nu}} depends on the parameter Da, characterizing the Brinkman extended non-Darcy flow.     

Secondary flow patterns and mixing in laminar pulsating flow through a curved pipe

Mixing by secondary flow is studied by particle image velocimetry (PIV) in a developing laminar pulsating flow through a circular curved pipe. The pipe curvature ratio is η = r ₀/r _c  = 0.09, and the curvature angle is 90°. Different secondary flow patterns are formed during an oscillation period due to competition among the centrifugal, inertial, and viscous forces. These different secondary-flow structures lead to different transverse-mixing schemes in the flow. Here, transverse mixing enhancement is investigated by imposing different pulsating conditions (Dean number, velocity ratio, and frequency parameter); favorable pulsating conditions for mixing are introduced. To obviate light-refraction effects during PIV measurements, a T-shaped structure is installed downstream of the curved pipe. Experiments are carried out for the Reynolds numbers range 420 ≤ Re_st ≤ 1,000 (Dean numbers 126.6 ≤ Dn ≤ 301.5) corresponding to non-oscillating flow, velocity component ratios 1 ≤ (β = U _max,osc/U _m,st) ≤ 4 (the ratio of velocity amplitude of oscillations to the mean velocity without oscillations), and frequency parameters 8.37 < (α = r ₀(ω/ν)^0.5) < 24.5, where α² is the ratio of viscous diffusion time over the pipe radius to the characteristic oscillation time. The variations in cross-sectional average values of absolute axial vorticity (|ζ|) and transverse strain rate (|ε|) are analyzed in order to quantify mixing. The effects of each parameter (Re_st, β, and α) on transverse mixing are discussed by comparing the dimensionless vorticities (|ζ _P |/|ζ _S |) and dimensionless transverse strain rates (|ε _P |/|ε _S |) during a complete oscillation period.     

Potential Flow Past a Slightly Deformed Porous Circular Cylinder

A uniform potential flow past a porous circular cylinder with a core of different permeability is discussed. The porous circular cylinder is slightly deformed whose radius is r=r₁(1+ecosm q){r=r_1(1+\epsilon \cos m \theta)} , where | e | << 1{\mid\epsilon\mid\ll 1} and m is a positive integer. Here r, θ are the polar coordinates and r ₁ is the characteristic radius of the cylinder. The drag force exerted by the exterior flow on the surface of the cylinder is calculated and it depends on the thickness of the porous material and on the permeabilities of the two porous regions. As special cases, porous cylinder with hollow core, rigid core, and deformed cylinder is discussed.     

Natural Convection Boundary-Layer Flow in a Porous Medium with Temperature-Dependent Boundary Conditions

The natural convection boundary-layer flow on a surface embedded in a fluid-saturated porous medium is discussed in the case when the wall heat flux is related to the wall temperature through a power-law variation. The flow within the porous medium is assumed to be described by Darcy’s law and the Boussinesq approximation is assumed for the density variations. Two cases are discussed, (i) stagnation-point flow and (ii) flow along a vertical surface. The possible steady states are considered first with the governing partial equations reduced to ordinary differential equations by similarity transformations and these latter equations further transformed to previously studied free-convection problems. This identifies values of the exponent N in the power-law wall temperature variation, N = 3/2 for stagnation-point flows and 3/2 ≤ N ≤ 3 for the vertical surface, where similarity solutions do not exist. Time development for stagnation-point flows is seen to depend on N, for N <  3/2 the steady state is approached at large times, for N ≥ 3/2 a singularity develops at finite time leading to thermal runaway. Numerical solutions for the vertical surface, where the temperature-dependent boundary condition becomes more significant as the solution develops, show that, for N < 3/2, the corresponding similarity solution is approached, whereas for N >  3/2 the solution breaks down at a finite distance along the surface.     

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.     

Convection in a Lid-Driven Heat-Generating Porous Cavity with Alternative Thermal Boundary Conditions

Mixed convection flow in a two-sided lid-driven cavity filled with heat-generating porous medium is numerically investigated. The top and bottom walls are moving in opposite directions at different temperatures, while the side vertical walls are considered adiabatic. The governing equations are solved using the finite-volume method with the SIMPLE algorithm. The numerical procedure adopted in this study yields a consistent performance over a wide range of parameters that were 10⁻⁴ ≤ Da ≤ 10⁻¹ and 0 ≤ Ra _I ≤ 10⁴. The effects of the parameters involved on the heat transfer characteristics are studied in detail. It is found that the variation of the average Nusselt number is non-linear for increasing values of the Darcy number with uniform or non-uniform heating condition.     

Insights into the Relationships Among Capillary Pressure,Saturation, Interfacial Area and Relative Permeability Using Pore-Network Modeling

To gain insight in relationships among capillary pressure, interfacial area, saturation, and relative permeability in two-phase flow in porous media, we have developed two types of pore-network models. The first one, called tube model, has only one element type, namely pore throats. The second one is a sphere-and-tube model with both pore bodies and pore throats. We have shown that the two models produce distinctly different curves for capillary pressure and relative permeability. In particular, we find that the tube model cannot reproduce hysteresis. We have investigated some basic issues such as effect of network size, network dimension, and different trapping assumptions in the two networks. We have also obtained curves of fluid–fluid interfacial area versus saturation. We show that the trend of relationship between interfacial area and saturation is largely influenced by trapping assumptions. Through simulating primary and scanning drainage and imbibition cycles, we have generated two surfaces fitted to capillary pressure, saturation, and interfacial area (P _c –S _w –a _nw ) points as well as to relative permeability, saturation, and interfacial area (k _r –S _w –a _nw ) points. The two fitted three-dimensional surfaces show very good correlation with the data points. We have fitted two different surfaces to P _c –S _w –a _nw points for drainage and imbibition separately. The two surfaces do not completely coincide. But, their mean absolute difference decreases with increasing overlap in the statistical distributions of pore bodies and pore throats. We have shown that interfacial area can be considered as an essential variable for diminishing or eliminating the hysteresis observed in capillary pressure–saturation (P _c –S _w ) and the relative permeability–saturation (k _r –S _w ) curves.     

The onset of buoyancy-driven convection in a ferromagnetic fluid saturated porous medium

The onset of buoyancy-driven convection in an initially quiescent ferrofluid saturated horizontal porous layer in the presence of a uniform vertical magnetic field is investigated. The Brinkman-Lapwood extended Darcy equation with fluid viscosity different from effective viscosity is used to describe the flow in the porous medium. The lower boundary of the porous layer is assumed to be rigid-paramagnetic, while the upper paramagnetic boundary is considered to be either rigid or stress-free. The thermal conditions include fixed heat flux at the lower boundary, and a general convective–radiative exchange at the upper boundary, which encompasses fixed temperature and fixed heat flux as particular cases. The resulting eigenvalue problem is solved numerically using the Galerkin technique. It is found that increase in the Biot number Bi, porous parameter σ, viscosity ratio Λ, magnetic susceptibility χ, and decrease in the magnetic number M ₁ and non-linearity of magnetization M ₃ is to delay the onset of ferroconvection in a porous medium. Further, increase in M ₁, M ₃, and decrease in χ, Λ, σ and Bi is to decrease the size of convection cells.     

Friction and forced convective heat transfer in a sintered porous channel with obstacle blocks

This work experimentally studies the flow characteristics and forced convective heat transfer in a sintered porous channel that filled with sintered copper beads of three average diameters ( 0.830, and 1.163 mm). The pressure drop and the local temperature measurements can be applied to figure out the distributions of the friction coefficient and the heat transfer coefficient. Three sintered porous channels differ in the arrangement of obstacle blocks. Model A has no obstacle. Models B and C have five obstacle blocks facing down and up, respectively, in a sintered porous channel. The range of experimental parameters, porosity, heat flux, and effect of forced convection are 0.370 ≤ ɛ ≤ 0.385, q=0.228, 0.872, 1.862 W/cm², and 200 ≤ Re _d ≤ 800. The permeability and inertia coefficient of each of the three sintered porous channels are analyzed. The results for Model A agree with those obtained by previous investigations in C _f distribution. The heat transfer of Model C exceeds that of Model A by approximately 20%. Finally, a series of empirical correlation equations were obtained for practical applications and engineering problems.     

A numerical study of non-isothermal turbulent coaxial jets

In this work, we propose to study non isothermal air–air coaxial jets with two different approaches: parabolic and elliptic approaches. The standard k−ε model and the RSM model were applied in this study. The numerical resolution of the equations governing this flow type was carried out for: the parabolic approach, by a “home-made” CFD code based on a finite difference method, and the elliptic approach by an industrial code (FLUENT) based on a finite volume method. In forced convection mode (Fr = ∞), the two turbulence models are valid for the prediction of the mean flow. But for turbulent sizes, k−ε model gives results closer to those achieved in experiments compared to RSM Model. Concerning the limit of validity of the parabolic and elliptic approaches, we showed that for velocities ratio r lower than 1, the results of the two approaches were satisfactory. On the other hand, for r > 1, the difference between the results became increasingly significant. In mixed convection mode (Fr ≅ 20), the results obtained by the two turbulence models for the mean axial velocity were very different even in the plume region. For the temperature and the turbulent sizes the two models give satisfactory results which agree well with the correlations suggested by the experimenters for X ≥ 20. Thus, the second order model with σ _t = 0.85 is more effective for a coaxial jet study in a mixed convection mode.     

Finite Element Modelling of Flow Through a Porous Medium Between Two Parallel Plates Using The Brinkman Equation

Despite the widespread use of the Darcy equation to model porous flow, it is well known that this equation is inconsistent with commonly prescribed no slip conditions at flow domain walls or interfaces between different sections. Therefore, in cases where the wall effects on the flow regime are expected to be significant, the Darcy equation which is only consistent with perfect slip at solid boundaries, cannot predict velocity and pressure profiles properly and alternative models such as the Brinkman equation need to be considered. This paper is devoted to the study of the flow of a Newtonian fluid in a porous medium between two impermeable parallel walls at different Darcy parameters (Da). The flow regime is considered to be isothermal and steady. Three different flow regimes can be considered using the Brinkman equation: free flow (Da > 1), porous flow (high permeability, 1 > Da > 10⁻⁶) and porous flow (low permeability Da < 10⁻⁶). In the present work the described bench mark problem is used to study the effects of solid walls for a range of low to high Darcy parameters. Both no-slip and slip conditions are considered and the results of these two cases are compared. The range of the applicability of the Brinkman equation and simulated results for different cases are shown.     

Velocity-defect scaling for turbulent boundary layers with a range of relative roughness

Velocity profile measurements in zero pressure gradient, turbulent boundary layer flow were made on a smooth wall and on two types of rough walls with a wide range of roughness heights. The ratio of the boundary layer thickness (δ) to the roughness height (k) was 16≤δ/k≤110 in the present study, while the ratio of δ to the equivalent sand roughness height (k _s) ranged from 6≤δ/k _s≤91. The results show that the mean velocity profiles for all the test surfaces agree within experimental uncertainty in velocity-defect form in the overlap and outer layer when normalized by the friction velocity obtained using two different methods. The velocity-defect profiles also agree when normalized with the velocity scale proposed by Zagarola and Smits (J Fluid Mech 373:33–70, 1998). The results provide evidence that roughness effects on the mean flow are confined to the inner layer, and outer layer similarity of the mean velocity profile applies even for relatively large roughness.     

Forced convection flow over a microsphere

Forced convection heat transfer characteristics around a microsphere subjected to uniform heat flux boundary condition is numerically investigated in this study. Moderate to high values of Reynolds number and a wide range of Prandtl number are considered. The analysis assumes that the continuity assumption is valid and hence the Navier–Stokes equations are solved for the range of Knudsen number of 0.001 ≤ Kn ≤ 0.1. The appropriate boundary conditions at the surface of the microsphere; the velocity slip and temperature jump are applied. The effect of the flow parameters: Re, Pr and Kn on the velocity and temperature distribution is presented and hence a better control on the boundary layer thickness can be achieved in the microscale level. Furthermore, the effect of the controlling parameters on the delay of flow separation, reduced shear stress, drag coefficient and on the Nusselt number profiles is also presented in the results.     

A New Algorithm for Estimating Three-Phase Relative Permeability from Unsteady-State Core Experiments

We present a new history matching method based on a Genetic Algorithm to estimate three-phase k _r (relative permeability) from unsteady-state coreflood experiments. In this method, relative permeabilities (k _r) are represented by quadratic B-Spline functions. Adjustable coefficients in k _r functions are changed in an iterative process to minimize an objective function. The objective function is defined as the difference between the measures and simulated values of the pressure drop across the core and fluids recovery during the experiment. One of the main features of this approach is that water and gas relative permeabilities (k _rw and k _rg) are assumed to be functions of two independent saturations as opposed to most of the existing empirical k _r models in which k _rw and k _rg are assumed to be only dependent of their own saturations. Another important aspect of this algorithm is that it considers inequality constrains to ensure that physically acceptable k _r curves are maintained throughout the iterative optimization process. A three-phase coreflood simulator has been developed based on this methodology that generates best k _r values by matching experimental data. The integrity of the developed software was first successfully verified by using two sets of experimental three-phase k _r data published in the literature. Then, the results of some three-phase coreflood experiments carried out in our laboratory were used to obtain three-phase k _r curves by this approach.     

Mixed Convection Heat Transfer in the Annulus Between Two Concentric Vertical Cylinders Using Porous Layers

A numerical investigation of the steady-state, laminar, axi-symmetric, mixed convection heat transfer in the annulus between two concentric vertical cylinders using porous inserts is carried out. The inner cylinder is subjected to constant heat flux and the outer cylinder is insulated. A finite volume code is used to numerically solve the sets of governing equations. The Darcy–Brinkman–Forchheimer model along with Boussinesq approximation is used to solve the flow in the porous region. The Navier–Stokes equation is used to describe the flow in the clear flow region. The dependence of the average Nusselt number on several flow and geometric parameters is investigated. These include: convective parameter, λ, Darcy number, Da, thermal conductivity ratio, K _r, and porous-insert thickness to gap ratio (H/D). It is found that, in general, the heat transfer enhances by the presence of porous layers of high thermal conductivity ratios. It is also found that there is a critical thermal conductivity ratio on which if the values of Kr are higher than the critical value the average Nusselt number starts to decrease. Also, it found that at low thermal conductivity ratio (K _r ≈ 1) and for all values of λ the porous material acts as thermal insulation.     

The optimal dimensions of convective-radiating circular fins

The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r _o/r) ⁿ ] is studied, while variation of thermal conductivity is of the form k=k _o[1+ɛ((T−T _∞)/ (T _b−T _∞)) ^m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r _o)/(r _e−r _o)]^λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999     

Non-Darcy buoyancy driven flows in a fluid saturated porous medium: the use of asymptotic computational fluid dynamics (ACFD) approach

Natural convection in a fluid saturated porous medium has been numerically investigated using a generalized non-Darcy approach. The governing equations are solved by using Finite Volume approach. First order upwind scheme is employed for convective formulation and SIMPLE algorithm for pressure velocity coupling. Numerical results are presented to study the influence of parameters such as Rayleigh number (10⁶ ≤Ra ≤10⁸), Darcy number (10⁻⁵ ≤ Da ≤ 10⁻²), porosity (0.4 ≤ ɛ ≤ 0.9) and Prandtl number (0.01 ≤ Pr ≤ 10) on the flow behavior and heat transfer. By combining the method of matched asymptotic expansions with computational fluid dynamics (CFD), so called asymptotic computational fluid dynamics (ACFD) technique has been employed to generate correlation for average Nusselt number. The technique is found to be an attractive option for generating correlation and also in the analysis of natural convection in porous medium over a fairly wide range of parameters with fewer simulations for numerical solutions.     

Two-Phase Flow Properties Prediction from Small-Scale Data Using Pore-Network Modeling

The objective of this work is to evaluate the prediction accuracy of network modeling to calculate transport properties of porous media based on the interpretation of mercury invasion capillary pressure curves only. A pore-scale modeling approach is used to model the multi-phase flow and calculate gas/oil relative permeability curves. The characteristics of the 3-D pore-network are defined with the requirement that the network model satisfactorily reproduces the capillary pressure curve (P_c curve), the porosity and the permeability. A sensitivity study on the effect of the input parameters on the prediction of capillary pressure and gas/oil relative permeability curves is presented. The simulations show that different input parameters can lead to similarly good reproductions of the experimental P_c, although the predicted relative permeabilities K_r are somewhat widespread. This means that the information derived from a mercury invasion P_c curve is not sufficient to characterize transport properties of a porous medium. The simulations indicate that more quantitative information on the wall roughness and the node/bond aspect ratio would be necessary to better constrain the problem. There is also evidence that in narrow pore size distributions pore body volume and pore throat radius are correlated while in broad pore size distributions they would be uncorrelated.     

Two-Phase Flow in Porous Media: Predicting Its Dependence on Capillary Number and Viscosity Ratio

Motivated by the need to determine the dependencies of two-phase flow in a wide range of applications from carbon dioxide sequestration to enhanced oil recovery, we have developed a standard two-dimensional, pore-level model of immiscible drainage, incorporating viscous and capillary effects. This model has been validated through comparison with several experiments. For a range of stable viscosity ratios (M = μ _injected,nwf/μ _{defending, wf} ≥ 1), we had increased the capillary number, N _c and studied the way in which the flows deviate from fractal capillary fingering at a characteristic time and become compact for realistic capillary numbers. This crossover has enabled predictions for the dependence of the flow behavior upon capillary number and viscosity ratio. Our results for the crossover agreed with earlier theoretical predictions, including the universality of the leading power-law indicating its independence of details of the porous medium structure. In this article, we have observed a similar crossover from initial fractal viscous fingering (FVF) to compact flow, for large capillary numbers and unstable viscosity ratios M < 1. In this case, we increased the viscosity ratio from infinitesimal values, and studied the way in which the flows deviate from FVF at a characteristic time and become compact for non-zero viscosity ratios. This crossover has been studied using both our pore-level model and micro-fluidic flow-cell experiments. The same characteristic time, τ = 1/M ^0.7, satisfactorily describes both the pore-level results for a range of large capillary numbers and the micro-fluidic flow cell results. This crossover should lead to predictions similar to those mentioned above.     

Relative permeability and capillary pressure functions of porous media as related to the displacement growth pattern

Visualization experiments of the unsteady immiscible displacement of a fluid by another are performed on glass-etched pore networks of well-controlled morphology by varying the fluid system and flow conditions. The measured transient responses of the fluid saturation and pressure drop across the porous medium are introduced into numerical solvers of the macroscopic two-phase flow equations to estimate the non-wetting phase, k_rnw, and wetting phase, k_rw, relative permeability curves and capillary pressure, P_c, curve. The correlation of k_rnw, k_rw, and P_c with the displacement growth pattern is investigated. Except for the capillary number, wettability, and viscosity ratio, the immiscible displacement growth pattern in a porous medium may be governed by the shear-thinning rheology of the injected or displaced fluid, and the porous sample length as compared to the thickness of the frontal region. The imbibition k_rnw increases as the flow pattern changes from compact displacement to viscous fingering or from viscous to capillary fingering. The imbibition k_rw increases as the flow pattern changes from compact displacement or capillary fingering to viscous fingering. As the shear-thinning behaviour of the NWP strengthens and/or the contact angle decreases, then the flow pattern is gradually dominated by irregular interfacial configurations, and the imbibition k_rnw increases. The imbibition P_c is a decreasing function of the capillary number or increasing function of the injected phase viscosity in agreement with the linear thermodynamic theory.