Gravity Affected Lateral Dispersion and Diffusion in a Stationary Horizontal Porous Medium Shear Flow

Transport of dissolved species by a carrier fluid in a porous medium comprises advection and diffusion/dispersion processes. Hydrodynamic dispersion is commonly characterized by an empirical relationship, in which the dispersion mechanism is described by contributions of molecular diffusion and mechanical dispersion expressed as a function of the molecular Peclét number. Mathematically these two phenomena are modeled by a constant diffusion coefficient and by velocity dependent dispersion coefficients, respectively. Here, the commonly utilized Bear--Scheidegger dispersion model of linear proportionality between mechanical dispersion and velocity, and the more complicated Bear--Bachmat model derived on a streamtube array model porous medium and better describing observed dispersion coefficients in the moderate molecular Peclét number range, will be considered. Analyzing the mixing flow of two parallelly flowing confluent fluids with different concentrations of a dissolved species within the frames of boundary layer theory one has to deal with transverse mixing only. With the Boussinesq approximation being adopted approximate analytical solutions of the corresponding boundary layer system of equations show that there is no effect of density coupling on concentration distributions across the mixing layer in the pure molecular diffusion regime case. With the Peclét number of the oncoming flow growing beyond unity, density coupling has an increasing influence on the mixing zone. When the Peclét number grows further this influence is successively reduced until its disappearance in the pure mechanical dispersion regime.     

Mixed Convection Boundary-Layer Flow in a Porous Medium Filled with Water Close to its Maximum Density

The steady mixed convection boundary-layer flow over a vertical impermeable surface in a porous medium saturated with water close to its maximum density is considered for uniform wall temperature and outer flow. The problem can be reduced to similarity form and the resulting equations are examined in terms of a mixed convection parameter λ and a parameter δ which measures the difference between the ambient temperature and the temperature at the maximum density. Both assisting (λ > 0) and opposing flows (λ < 0) are considered. A value δ₀ is found for which there are dual solutions for a range λ_c < λ < 0 of λ (the value of λ_c dependent on δ) and single solutions for all λ ≥ 0. Another value of δ₁ of δ, with δ₁ > δ₀, is found for which there are dual solutions for a range 0 < λ < λ_c of positive values of λ, with solutions for all λ≤ 0. There is also a range δ₀ <  δ < δ₁ where there are solutions only for a finite range of λ, with critical points at both positive and negative values of λ, thus putting a finite limit on the range of existence of solutions.     

The Convective Boundary-Layer Flow on a Reacting Surface in a Porous Medium

The convective boundary-layer flow on an impermeable vertical surface in a fluid-saturated porous medium is considered where the flow results from the heat released by an exothermic catalytic reaction on the surface converting a reactive component within the convective fluid to an inert product. The reaction is modelled by first-order kinetics with an Arrhenius temperature dependence. Numerical solutions of the governing equations are obtained for a range of parameter values. These show, for large activation energies, that localized rapid changes in wall temperature and localized high reaction rates occur a little way from the leading edge. Asymptotic expansions, valid at large distances from the leading edge, are derived, the form that these expansions take is qualitatively different depending on whether or not reactant consumption is included in the model.     

Effect of Viscous Dissipation on a Quasi-Parallel Free Convection Boundary-Layer Flow over a Vertical Flat Plate in a Porous Medium

Transport in Porous Media -     

Mixed Convection Boundary-Layer Flow Near the Stagnation Point on a Vertical Surface in a Porous Medium: Brinkman Model with Slip

The steady boundary-layer flow near the stagnation point on an impermeable vertical surface with slip that is embedded in a fluid-saturated porous medium is investigated. Using appropriate similarity variables, the governing system of partial differential equations is transformed into a system of ordinary differential equations. This system is then solved numerically. The features of the flow and the heat transfer characteristics for different values of the governing parameters, namely, the Darcy–Brinkman, Γ, mixed convection, λ, and slip, γ, parameters, are analysed and discussed in detail for the cases of assisting and opposing flows. It is found that dual solutions exist for assisting flows, as well as those usually reported in the literature for opposing flows. A stability analysis of the steady flow solutions encountered for different values of the mixed convection parameter λ is performed using a linear temporal stability analysis. This analysis reveals that for γ  =  0 (slip absent) and Γ  =  1 the lower solution branch is unstable while the upper solution branch is stable.     

Mixed Convection on a Horizontal Surface Embedded in a Porous Medium: the Structure of a Singularity

The mixed convection boundary-layer flow on a horizontal impermeable surface embedded in a saturated porous medium and driven by a local heat source is considered. Similarity solutions are obtained for specific outer flow variations and these are shown to have a solution only for parameter values greater than some critical value. When this is not the case the solution develops a singularity at a finite distance from the leading edge. The nature of this singularity is also discussed.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

Modeling Variable Density Flow and Solute Transport in Porous Medium: 1. Numerical Model and Verification

A new numerical model for the resolution of density coupled flow and transport in porous media is presented. The model is based on the mixed hybrid finite elements (MHFE) and discontinuous finite elements (DFE) methods. MHFE is used to solve the flow equation and the dispersive part of the transport equation. This method is more accurate in the calculation of velocities and ensures continuity of fluxes from one element to the adjacent one. DFE is used to solve the convective part of the transport equation. Combined with a slope limiting procedure, it avoids numerical instabilities and creates a very limited numerical dispersion, even for high grid Peclet number.Flow and transport equations are coupled by a standard iterative scheme. Residual based criterion is used to stop the iterations. Simulations of an unstable equilibrium show the effects of the criteria used to stop the iterations and the stopping criterion in the solver. The effects are more important for finer grids than for coarser grids.The numerical model is verified by the simulation of standard benchmarks: the Henry and the Elder test cases. A good agreement is found between the revised semianalytical Henry solution and the numerical solution. The Elder test case was also studied. The simulations were similar to those presented in previous works but with significantly less unknowns (i.e. coarser grids). These results show the efficiency of the used numerical schemes.     

An Experimental Study of Miscible Displacements in Porous Media with Variation of Fluid Density and Viscosity

The effect of fluid density and viscosity on dispersion in miscible displacements in porous media is examined. Miscible displacement experiments with fluid pairs having density and viscosity differences are carried out in a Plexiglas column containing a homogeneous and isotropic sand pack. Tracer tests and tests of both stable and unstable miscible displacement are conducted using NaCl and glycerine solutions. Concentration breakthrough curves are measured through an electrical monitoring technique. Following the conventional convection–dispersion formulation, the dispersion coefficient is determined by performing a least squares fit to the measured concentration breakthrough curves. It is found that for stable displacements dispersion coefficients drop continually when density differences increase or when viscosity ratios of the displaced and displacing fluid decrease. In the unstable case the dispersion coefficients increase with both density and viscosity differences.     

Chaos in the Convective Flow of a Fluid Mixture in a Porous Medium

The results of the study of the global behaviour of the convective flow of a binary mixture in a porous medium are presented. Bifurcation diagram, fixed points, periodic, chaotic solutions, stable and unstable manifolds, and basins of attraction have been calculated. Different behaviours (chaos, undecidable behaviour, etc.) have been found.     

The Effect of Magnetic Field Dependent Viscosity on Thermosolutal Convection in a Ferromagnetic Fluid Saturating a Porous Medium

The effect of magnetic field dependent viscosity on thermosolutal convection in a ferromagnetic fluid saturating a porous medium is considered for a fluid layer heated and soluted from below in the presence of uniform magnetic field. Using linearized stability theory and normal mode analysis, an exact solution is obtained for the case of two free boundaries. For case of stationary convection, medium permeability has a destabilizing effect, whereas a stable solute gradient and magnetic field dependent viscosity have a stabilizing effect on the system. In the absence of magnetic field dependent viscosity, the destabilizing effect of non-buoyancy magnetization is depicted but in the presence of magnetic field dependent viscosity non-buoyancy magnetization may have a destabilizing or stabilizing effect on the onset of instability. The critical wave number and the critical magnetic thermal Rayleigh number for the onset of instability are also determined numerically for sufficiently large values of buoyancy magnetization parameter M₁ and the results are depicted graphically. The principle of exchange of stabilities is found to hold true for the ferromagnetic fluid saturating a porous medium heated from below in the absence of stable solute gradient. The oscillatory modes are introduced due to the presence of the stable solute gradient, which were non-existent in its absence. A sufficient condition for the non-existence of overstability is also obtained. The paper also reaffirms the qualitative findings of earlier investigations which are, in fact, limiting cases of the present study.     

Effects of Temperature-Dependent Viscosity with Soret and Dufour Numbers on Non-Darcy MHD Free Convective Heat and Mass Transfer Past a Vertical Surface Embedded in a Porous Medium

An analysis is presented to investigate the effects of temperature-dependent viscosity, thermal dispersion, Soret number and Dufour number on non-Darcy MHD free convective heat and mass transfer of a viscous, incompressible and electrically conducting fluid past a vertical isothermal surface embedded in a saturated porous medium. The governing partial differential equations are transferred into a system of ordinary differential equations, which are solved numerically using a fourth order Runge–Kutta scheme with the shooting method. Comparisons with previously published work by Hong and Tien [Hong, J. T. and Tien, C. L.: 1987, Int. J. Heat Mass Transfer 30, 143–150] and Sparrow et al. [Sparrow, E. M. et al.: 1964, AIAA J. 2 652–659] are performed and good agreement is obtained. Numerical results of the skin friction coefficient, the local Nusselt number and the local Sherwood number as well as the velocity, temperature and concentration profiles are presented for different physical parameters.     

Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid Past a Vertical Plate Immersed in a Porous Medium

The effects of viscous dissipation on unsteady free convection from an isothermal vertical flat plate in a fluid saturated porous medium are examined numerically. The Darcy–Brinkman–Forchheimer model is employed to describe the flow field. A new model of viscous dissipation is used for the Darcy–Brinkman–Forchheimer model of porous media. The simultaneous development of the momentum and thermal boundary layers are obtained by using a finite difference method. Boundary layer and Boussinesq approximation have been incorporated. Numerical calculations are carried out for various parameters entering into the problem. Velocity and temperature profiles as well as local friction factor and local Nusselt number are shown graphically. It is found that as time approaches infinity, the values of friction factor and heat transfer coefficient approach steady state.     

Buoyant Flow with Viscous Heating in a Vertical Circular Duct Filled with a Porous Medium

Combined, forced, and free flow in a vertical circular duct filled with a porous medium is investigated according to the Darcy–Boussinesq model. The effect of viscous dissipation is taken into account. It is shown that a thermal boundary condition compatible with fully developed and axisymmetric flow is either a linearly varying wall temperature in the axial direction or, only in the case of uniform velocity profile, an axial linear-exponential wall temperature change. The case of a linearly varying wall temperature corresponds to a uniform wall heat flux and includes the uniform wall temperature as a special case. A general analytical solution procedure is performed, by expressing the seepage velocity profile as a power series with respect to the radial coordinate. It is shown that, for a fixed thermal boundary condition, i.e., for a prescribed slope of the wall temperature, and for a given flow rate, there exist two solutions of the governing balance equations provided that the flow rate is lower than a maximum value. When the maximum value is reached, the dual solutions coincide. When the flow rate is higher than its maximum, no axisymmetric solutions exist. E. Magyari is on leave from the Institute of Building Technology, ETH—Zürich.     

Aiding and Opposing Mixed Convection Flow in Melting From a Vertical Flat Plate Embedded in a Porous Medium

The problem of melting from a vertical flat plate embedded in a porous medium is studied. The main focus is to determine the effect of mixed convection flow in the liquid phase on the melting phenomenon. Both aiding and opposing flows are considered. The conservation equations that govern the problem are reduced to a system of nonlinear ordinary differential equations. The governing equations are solved numerically. Numerical results are obtained for the temperature and flow fields in the melting region. The melting phenomenon decreases the local Nusselt number at the solid–liquid interface.     

Analytic Series Solution for Unsteady Mixed Convection Boundary Layer Flow Near the Stagnation Point on a Vertical Surface in a Porous Medium

In this paper, we solve the unsteady mixed convection flow near the stagnation point on a heated vertical flat plate embedded in a Darcian fluid-saturated porous medium by means of an analytic technique, namely the Homotopy Analysis Method. Different from previous perturbation results, our analytic series solutions are accurate and uniformly valid for all dimensionless times and for all possible values of mixed convection parameter, and besides agree well with numerical results. This provides us with a new analytic approach to investigate related unsteady problems.     

Effects of Viscous Dissipation and Flow Work on Forced Convection in a Channel Filled by a Saturated Porous Medium

Fully developed forced convection in a parallel plate channel filled by a saturated porous medium, with walls held either at uniform temperature or at uniform heat flux, with the effects of viscous dissipation and flow work included, is treated analytically. The Brinkman model is employed. The analysis leads to expressions for the Nusselt number, as a function of the Darcy number and Brinkman number.     

Analytical Study of Fluid Flow and Heat Transfer during Forced Convection in a Composite Channel Partly Filled with a Brinkman–Forchheimer Porous Medium

In this paper, the problem of fully developed forced convection in a parallel-plate channel partly filled with a homogeneous porous material is considered. The porous material is attached to the walls of the channel, while the center of the channel is occupied by clear fluid. The flow in the porous material is described by a nonlinear Brinkman–Forchheimer-extended Darcy equation. Utilizing the boundary-layer approach, analytical solutions for the flow velocity, the temperature distribution, as well as for the Nusselt number are obtained. Dependence of the Nusselt number on several parameters of the problem is extensively investigated.     

Numerical Simulation of Thermal and Concentration Natural Convection in a Porous Cavity in the Presence of an Opposing Flow

Two-dimensional steady-state thermal concentration convection in a rectangular porous cavity is simulated numerically. The temperature and concentration gradients are horizontal and the buoyancy forces act either in the same or in opposite directions. The flow through the porous medium is described by the Darcy-Brinkman or Forchheimer equations. The SIMPLER numerical algorithm based on the finite volume approach is used for solving the problem in the velocity-pressure variables.Numerous series of calculations were carried out over the range Ra _t =3·10⁶ and 3·10⁷, 10^-6 < Da < 1, 1 < N < 20, Le=10 and 100, where Ra, Da, Le, and N are the Rayleigh, Darcy, and Lewis numbers and the buoyancy ratio, respectively. It is shown that the main effect of the presence of the porous medium is to reduce the heat and mass transfer and attenuate the flow field with decrease in permeability. For a certain combination of the Ra, Le, and N numbers the flow has a multicellular structure. The mean Nusselt and Sherwood numbers are presented as functions of the governing parameters.     

Inertial and Second Grade Effects on Macroscopic Equations for the Flow of a Fluid in a Porous Medium

Homogenization techniques are used to upscale from pore to laboratory or field scale viscous and second grade nonNewtonian flow in a porous medium. Nonlinear forms of Darcy's law are obtained and analysed under a series of symmetry properties. The general case of displacement of one of these fluids by another with different properties is considered and a linear stability analysis is performed.