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.     

Convection currents in a vertical layer with a space-modulated temperature distribution on the boundaries

Steady convective motions in a plane vertical fluid layer are investigated. The temperature along the boundaries of the layer varies harmonically and has different average values on each of the boundaries. Thus space-period modulation of the temperature of the walls is assigned along with average lateral heating of the layer. The form of the plane steady motions and regions of existence of through currents and currents of cellular structure are found for various values of the parameters of the problem by the finite difference grid-point method. The dependence of the main characteristics of fluid motion on the Grashof number is determined. The results presented in the article pertain to the case when the period of modulation of the temperature of the boundaries coincides with the wavelength of the critical mode of a plane-parallel current. A numerical investigation of supercritical motions in a vertical layer with plane isothermal boundaries heated to a different temperature was carried out in [1–3]. The effect of a space-periodic inhomogeneity due to curvature of walls on the form and stability of convective motions in a vertical layer with lateral heating was examined in [4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 20–25, September–October, 1978.The author thanks E. M. Zhukhovitskii for formulating the problem and supervising the work and G. Z. Gershuni for discussions and useful comments.     

The Onset and Nonlinear Regimes of Convection in a Two-Layer System of Fluid and Porous Medium Saturated by the Fluid

The onset of convection and its nonlinear regimes in a heated from below two-layer system consisting of a horizontal pure fluid layer and porous medium saturated by the same fluid is studied under the conditions of static gravitational field. The problem is solved numerically by the finite-difference method. The competition between the long-wave and short-wave convective modes at various ratios of the porous layer to the fluid layer thicknesses is analyzed. The data on the nature of convective motion excitation and flow structure transformation are obtained for the range of the Rayleigh numbers up to quintuple supercriticality. It has been found that in the case of a thick porous layer the steady-state convective regime occurring after the establishment of the mechanical equilibrium becomes unstable and gives way to the oscillatory regime at some value of the Rayleigh number. As the Rayleigh number grows further the oscillatory regime of convection is again replaced by the steady-state convective regime.     

Analysis of Flow and Heat Transfer in a Vertical Rectangular Duct Using a Non-Darcy Model

Numerical investigation of steady natural convection flow through a fluid-saturated porous medium in a vertical rectangular duct is investigated. The Darcy-Forchheimer-Brinkman model is used to represent the fluid transport within the porous medium. One of the vertical walls of the duct is cooled to a constant temperature, while the other wall is heated to constant but different temperature. The other two sides of the duct are insulated. The finite difference method of second-order accuracy is used to solve the non-dimensional governing equations. The results are presented graphically to show the effects of the Darcy number, inertial parameter, Grashof number, Brinkman number, aspect ratio, and viscosity ratio. It is found that an increase in the Darcy number and inertial parameter reduces the flow intensity whereas an increase in the Grashof number, Brinkman number, aspect ratio, and viscosity ratio increases the flow intensity.     

Asymptotic form of the solution to the problem of multicomponent mixture flow through a porous medium with a boundary layer

In problems of two-phase mixture flow through a porous medium in a subterranean stratum a boundary layer phenomenon arises caused by the fact that relative phase motion exists in the system, and so having no analogy with the single-phase case. The physical nature of boundary layer phenomena is explained, and an asymptotic solution is constructed for the self-similar problem with an arbitrary number of components in the system, by using the method of matched asymptotic forms. The conditions are established for the motions of a multicomponent and a binary mixture to be equivalent, and a study is made of the role of convective factors in the formation of averaged working indices for the stratum.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 94–100, July–August, 1985.     

Fully Developed Mixed Convection Flow in a Vertical Channel Containing Porous and Fluid Layer with Isothermal or Isoflux Boundaries

An analysis of fully developed combined free and forced convective flow in a fluid saturated porous medium channel bounded by two vertical parallel plates is presented. The flow is modeled using Brinkman equation model. The viscous and Darcy dissipation terms are also included in the energy equation. Three types of thermal boundary conditions such as isothermal–isothermal, isoflux–isothermal, and isothermal–isoflux for the left–right walls of the channel are considered. Analytical solutions for the governing ordinary differential equations are obtained by perturbation series method. In addition, closed form expressions for the Nusselt number at both the left and right channel walls are derived. Results have been presented for a wide range of governing parameters such as porous parameter, ratio of Grashof number and Reynolds number, viscosity ratio, width ratio, and conductivity ratio on velocity, and temperature fields. It is found that the presence of porous matrix in one of the region reduces the velocity and temperature.     

Convective cooling/heating induced thermal stresses in a fluid saturated porous medium undergoing local thermal non-equilibrium

Local thermal non-equilibrium (LTNE) may have profound effects on the pore pressure and thermal stresses in fluid saturated porous media under transient thermal loads. This work investigates the temperature, pore pressure, and thermal stress distributions in a porous medium subjected to convective cooling/heating on its boundary. The LTNE thermo-poroelasticity equations are solved by means of Laplace transform for two fundamental problems in petroleum engineering and nuclear waste storage applications, i.e., an infinite porous medium containing a cylindrical hole or a spherical cavity subjected to symmetrical thermo-mechanical loads on the cavity boundary. Numerical examples are presented to examine the effects of LTNE under convective cooling/heating conditions on the temperature, pore pressure and thermal stresses around the cavities. The results show that the LTNE effects become more pronounced when the convective heat transfer boundary conditions are employed. For the cylindrical hole problem of a sandstone formation, the thermally induced pore pressure and the magnitude of thermal stresses are significantly higher than the corresponding values in the classical poroelasticity, which is particularly true under convective cooling with moderate Biot numbers. For the spherical cavity problem of a clay medium, the LTNE effect may become significant depending on the boundary conditions employed in the classical theory.     

Different types of nonlinear convective oscillation in a multilayer system

The nonlinear development of oscillatory instability under the joint action of buoyant and thermocapillary effects in a multilayer system is investigated. The nonlinear two-dimensional convective regimes are studied by the finite-difference method. Rigid heat-insulated lateral walls are considered. Different types of nonlinear flow, symmetric and asymmetric oscillations, have been found. It is shown that the oscillatory motion takes place in an interval of the Grashof number values bounded both from below (by the mechanical equilibrium) and from above (by the steady state). Cavities with different lengths are considered.     

Non-darcy natural convection on a vertical cylinder in a saturated porous medium

The steady free convection boundary layer flow of non-Darcy fluid along an isothermal vertical cylinder embedded in a saturated porous medium using the Ergun model has been studied. The partial differential equations governing the flow have been solved numerically using an implicit finite-difference scheme developed by Keller. It is found that the heat transfer is strongly affected by the modified Grashof number which characterizes the non-Darcy fluid, and the curvature parameter. Also the heat transfer is found to be more than that of the flat plate.     

Numerical investigation of convective motion in a closed cavity

The investigation of thermal convection in a closed cavity is of considerable interest in connection with the problem of heat transfer. The problem may be solved comparatively simply in the case of small characteristic temperature difference with heating from the side, when equilibrium is not possible and when slow movement is initiated for an arbitrarily small horizontal temperature gradient. In this case the motion may be studied using the small parameter method, based on expanding the velocity, temperature, and pressure in series in powers of the Grashof number—the dimensionless parameter which characterizes the intensity of the convection [1–4]. In the problems considered it has been possible to find only two or three terms of these series. The solutions obtained in this approximation describe only weak nonlinear effects and the region of their applicability is limited, naturally, to small values of the Grashof number (no larger than 10³).With increase of the temperature difference the nature of the motion gradually changes—at the boundaries of the cavity a convective boundary layer is formed, in which the primary temperature and velocity gradients are concentrated; the remaining portion of the liquid forms the flow core. On the basis of an analysis of the equations of motion for the plane case, Batchelor [4] suggested that the core is isothermal and rotates with constant and uniform vorticity. The value of the vorticity in the core must be determined as the eigenvalue of the problem of a closed boundary layer. A closed convective boundary layer in a horizontal cylinder and in a plane vertical stratum was considered in [5, 6] using the Batchelor scheme. The boundary layer parameters and the vorticity in the core were determined with the aid of an integral method. An attempt to solve the boundary layer equations analytically for a horizontal cylinder using the Oseen linearization method was made in [7].However, the results of experiments in which a study was made of the structure of the convective motion of various liquids and gases in closed cavities of different shapes [8–13] definitely contradict the Batchelor hypothesis. The measurements show that the core is not isothermal; on the contrary, there is a constant vertical temperature gradient directed upward in the core. Further, the core is practically motionless. In the core there are found retrograde motions with velocities much smaller than the velocities in the boundary layer.The use of numerical methods may be of assistance in clarifying the laws governing the convective motion in a closed cavity with large temperature differences. In [14] the two-dimensional problem of steady air convection in a square cavity was solved by expansion in orthogonal polynomials. The author was able to progress in the calculation only to a value of the Grashof numberG=10⁴. At these values of the Grashof numberG the formation of the boundary layer and the core has really only started, therefore the author's conclusion on the agreement of the numerical results with the Batchelor hypothesis is not justified. In addition, the bifurcation of the central isotherm (Fig. 3 of [14]), on the basis of which the conclusion was drawn concerning the formation of the isothermal core, is apparently the result of a misunderstanding, since an isotherm of this form obviously contradicts the symmetry of the solution.In [5] the method of finite differences is used to obtain the solution of the problem of strong convection of a gas in a horizontal cylinder whose lateral sides have different temperatures. According to the results of the calculation and in accordance with the experimental data [9], in the cavity there is a practically stationary core. However, since the authors started from the convection equations in the boundary layer approximation they did not obtain any detailed information on the core structure, in particular on the distribution of the temperature in the core.In the following we present the results of a finite difference solution of the complete nonlinear problem of plane convective motion in a square cavity. The vertical boundaries of the cavity are held at constant temperatures; the temperature varies linearly on the horizontal boundaries. The velocity and temperature distributions are obtained for values of the Grashof number in the range 0<G4·10⁵ and for a value of the Prandtl number P=1. The results of the calculation permit following the formation of the closed boundary layer and the very slowly moving core with a constant vertical temperature gradient. The heat flux through the cavity is found as a function of the Grashof number.     

Investigation of the Relative Motion of a Viscous Fluid and a Porous Medium Using the Brinkman Equations

Exact solutions are obtained for the following three problems in which the Brinkman filtration equations are used: laminar fluid flow between parallel plane walls, one of which is rigid while the other is a plane layer of saturated porous medium, motion of a plane porous layer between parallel layers of viscous fluid, and laminar fluid flow in a cylindrical channel bounded by an annular porous layer.     

Radiation effects on combined convection over a vertical flat plate embedded in a porous medium of variable porosity

The present paper is concerned with the study of radiation effects on the combined (forced-free) convection flow of an optically dense viscous incompressible fluid over a vertical surface embedded in a fluid saturated porous medium of variable porosity with heat generation or absorption. The effects of radiation heat transfer from a porous wall on convection flow are very important in high temperature processes. The inclusion of radiation effects in the energy equation leads to a highly non-linear partial differential equations which are transformed to a system of ordinary differential equations using non-similarity transformation. These equations are then solved numerically using implicit finite-difference method subject to appropriate boundary and matching conditions. A parametric study of the physical parameters such as the particle diameter-based Reynolds number, the flow based Reynolds number, the Grashof number, the heat generation or absorption co-efficient and radiation parameter is conducted on temperature distribution. The effects of radiation and other physical parameters on the local skin friction and on local Nusselt number are shown graphically. It is interesting to observe that the momentum and thermal boundary layer thickness increases with the radiation and decrease with increase in the Prandtl number.     

Three dimensional unsteady convection and mass transfer flow through porous medium

The problem of three dimensional unsteady convection flow through a porous medium, with effect of mass transfer bounded by an infinite vertical porous plate is discussed, when the suction at the plate is transverse sinusoidal and the plate temperature oscillates in time about a constant mean. Assuming the free stream velocity to be uniform, approximate solutions are obtained for the flow field, the temperature field, the skin-friction and the rate of heat transfer. The dependence of solution on Pr (Prandtl number), Gr (Grashof number based on temperature), Gc (modified Grashof number based on concentration difference), Sc (Schimdt number), the frequency and the permeability parameter is also investigated.     

Throughflow Effects on Penetrative Convection in Superposed Fluid and Porous Layers

The effect of vertical throughflow on the onset of penetrative convection simulated via internal heating in a two-layer system in which a layer of fluid overlies and saturates a layer of porous medium is studied. Flow in the porous medium is governed by Forchheimer-extended Darcy equation, and Beavers?CJoseph slip condition is applied at the interface between the fluid and the porous layers. The boundaries are considered to be rigid, however permeable, and insulated to temperature perturbations. The eigenvalue problem is solved using a regular perturbation technique with wave number as a perturbation parameter. The ratio of fluid layer thickness to porous layer thickness, ??, the direction of throughflow, and the presence of volumetric internal heat source in fluid and/or porous layer play a decisive role on the stability characteristics of the system. In addition, the influence of Prandtl number arising due to throughflow is also emphasized on the stability of the system. It is observed that both stabilizing and destabilizing factors can be enhanced because of the simultaneous presence of a volumetric heat source and vertical throughflow so that a more precise control (suppress or augment) of thermal convective instability in a layer of fluid or porous medium is possible.     

Effect of Thermal Modulation on the Onset of Convection in Walters B Viscoelastic Fluid-Saturated Porous Medium

The linear stability of Walters B viscoelastic fluid-saturated horizontal porous layer is examined theoretically when the walls of the porous layer are subjected to time-periodic temperature modulation. Three types of boundary temperature modulations are considered namely, symmetric, asymmetric, and only the lower wall temperature is modulated while the upper wall is held at constant temperature. A regular perturbation method based on small amplitude of applied temperature field is used to compute the critical values of Rayleigh number and the corresponding wave number. The shift in critical Rayleigh number is calculated as a function of modulation frequency, viscoelastic parameter, and Prandtl number. The effect of all three types of modulations is found to be destabilizing as compared to the unmodulated system. This result is in contrast to the system with other types of fluids. Besides, the influence of physical parameters on the control of convective instability of the system is discussed.     

Unsteady Flow Through a Highly Porous Medium in the Presence of Radiation

The unsteady natural convection flow of a viscous and incompressible fluid through a porous medium with high porosity bounded by a vertical infinite stationary plate in the presence of radiation has been investigated. The fluid is assumed to be a gray, emitting and absorbing radiation, but non-scattering medium. The effects of the radiation parameter, Grashof number and permeability parameter of the medium on the velocity field as well as the effects of the radiation parameter and Prandtl number on the temperature field have been included in the analysis.     

Stability of steady nonisothermal flow in a vertical layer with permeable boundaries in the microconvection model

The problem of the stability of steady convective viscous incompressible fluid flow in a vertical layer with boundaries at different temperatures is considered in the presence of transverse blowing through the layer. The complete spectral problem is solved for a silicon melt. The neutral curve is constructed and the critical Grashof number is found. The numerical calculations show that the presence of transverse blowing significantly affects the flow stability. As compared with the Oberbeck-Boussinesq model, in the microconvection model the instability develops at lower wavenumbers.     

Similarity solutions for buoyancy induced flows in a saturated porous medium adjacent to impermeable horizontal surfaces

In the present paper similarity solutions for the convective flow induced by buoyancy in a saturated porous medium adjacent to horizontal impermeable surfaces are obtained. The analysis incorporates the variation of permeability from the wall and expressions for boundary layer thickness, local and overall surface heat-flux are obtained. Applications of the results to convective flows in a geothermal reservoir are discussed.     

The Onset of Convection in an Anisotropic Porous Layer Using a Thermal Non-Equilibrium Model

The stability of a horizontal fluid saturated anisotropic porous layer heated from below and cooled from above is examined analytically when the solid and fluid phases are not in local thermal equilibrium. Darcy model with anisotropic permeability is employed to describe the flow and a two-field model is used for energy equation each representing the solid and fluid phases separately. The linear stability theory is implemented to compute the critical Rayleigh number and the corresponding wavenumber for the onset of convective motion. The effect of thermal non-equilibrium and anisotropy in both mechanical and thermal properties of the porous medium on the onset of convection is discussed. Besides, asymptotic analysis for both very small and large values of the interphase heat transfer coefficient is also presented. An excellent agreement is found between the exact and asymptotic solutions. Some known results, which correspond to thermal equilibrium and isotropic porous medium, are recovered in limiting cases.     

Double diffusive convection of anomalous density fluids in a porous cavity

A numerical study has been performed to analyze the combined effect of temperature and species gradients on the buoyancy-driven natural convection flow of cold water near its density extremum contained in a porous cavity. The governing equations are descretized using the finite volume method. The results of the investigation are presented in the form of steady-state streamlines, velocity vectors, isotherms, and isoconcentrationlines. The results are discussed for different porosities, Darcy numbers, and Grashof numbers. The heat and mass transfer rates calculated are found to behave nonlinearly with hot wall temperature. The heat and mass transfer are increased with increasing Darcy number and porosity. It is found that the convective heat and mass transfer rate are greatly affected by the presence of density maximum.