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.     

Stability of convective motion in a vertical layer in the presence of longitudinal vibrations

The influence of vibrations of a cavity containing a fluid on the convective stability of the equilibrium has been investigated on a number of occasions [1]. The stability of convective flows in a modulated gravity field has not hitherto been studied systematically. There is only the paper of Baxi, Arpaci, and Vest [2], which contains fragmentary data corresponding to various values of the determining parameters of the problem. The present paper investigates the linear stability of convective flow in a vertical plane layer with walls at different temperatures in the presence of longitudinal harmonic vibrations of the cavity containing the fluid. It is assumed that the frequency of the vibrations is fairly high; the motion is described by the equations of the averaged convective motion. The stability boundaries of the flow with respect to monotonic perturbations in the region of Prandtl numbers 0 ? P ? 10 are determined. It is found that high-frequency vibrations have a destabilizing influence on the convective motion. At sufficiently large values of the vibration parameter, the flow becomes unstable at arbitrarily small values of the Grashof number, this being due to the mechanism of vibrational convection, which leads to instability even under conditions of weightlessness, when the main flow is absent [3, 4].     

Stability analysis of soret-driven double-diffusive convection of Maxwell fluid in a porous medium

Stability analysis of double-diffusive convection for viscoelastic fluid with Soret effect in a porous medium is investigated using a modified-Maxwell-Darcy model. We use the linear stability analysis to investigate how the Soret parameter and the relaxation time of viscoelastic fluid effect the onset of convection and the selection of an unstable wavenumber. It is found that the Soret effect is to destabilize the system for oscillatory convection. The relaxation time also enhances the instability of the system. The effects of Soret coefficient and relaxation time on the heat transfer rate in a porous medium are studied using the nonlinear stability analysis, the variation of the Nusselt number with respect to the Rayleigh number is derived for stationary and oscillatory convection modes. Some previous results can be reduced as the special cases of the present paper.     

Effects of Brownian Motion and Thermophoresis on the Mixed Convection of Nanofluid in a Porous Channel Including Flow Reversal

This article is devoted to combined convection heat transfer of nanofluids through a vertical channel filled with a homogeneous and isotropic porous medium. The flow is assumed to be fully developed and the “Brinkman extended Darcy” model is used for the flow in the porous media and “clear compatible” viscous dissipation model is considered. Also the model utilized for the nanofluid incorporates the effects of Brownian motion and thermophoresis. The governing momentum, energy, and nanopartices volume fraction equations are solved both analytically and numerically. The effects of the influential dimensionless parameters such as Brownian and thermophoresis parameters, mixed convection parameter (Gr/Re), Brinkman, Darcy and Lewis numbers on dimensionless velocity and temperature distributions and pressure drop are studied. Also, the results of the Nusselt number for the both left and right walls are presented and discussed.     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.     

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.     

多孔介质中热对流的分叉机理研究

本文利用解析分析方法研究了数值模拟发现的多孔介质层中出现的对流分叉机理，指出控制方程中的Ｒａｙｌｅｉｇｈ数，是决定流动的特征参数。当Ｒａｙｌｅｉｇｈ数小于临界数值时，多孔介质内流动处于静止传热状态，并且这种状态是稳定的。如果Ｒａｙｌｅｉｇｈ数大于临界数值，非线性方程出现分叉解，文中指出，存在多个使平凡解失稳而分叉的临界Ｒａｙｌｅｉｇｈ数，当Ｒａｙｌｅｉｇｈ数由小到大经历这些临界数值时，其由平凡解发展起来的分叉解的流态，依次由单回流区转变为双回流区及三回流区。理论分析给出了分叉解和分叉解的振幅方程，阐明了分叉的机理，其结论和数值结果定性一致．     

Vibrational convection in a saturated porous medium in the absence of thermal equilibrium between the solid and liquid phases

Thermal vibrational convection in a saturated porous medium is theoretically studied on the basis of a thermal nonequilibrium model, in accordance with which the temperatures of the porous medium and the saturating liquid can be different. In the high-frequency vibration approximation the averaged equations of convection are derived. The dependence of the vibration force direction on the interphase heat transfer coefficient and the vibration frequency is established. Vibrational convection in a cylindrical layer is studied. It is shown that, depending on the interphase heat transfer coefficient, the flows of two types differing in the liquid circulation direction can exist.     

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.     

Slip effects on a mixed convection flow of a third-grade fluid near the orthogonal stagnation point on a vertical surface

A mixed convection flow of a third-grade fluid near the orthogonal stagnation point on a vertical surface with slip and viscous dissipation effects is investigated. The governing partial differential equations for the third-grade fluid are converted into a system of nonlinear ordinary differential equations by using a similarity transformation. The effects of various parameters, including the Weissenberg number, third-grade parameter, local Reynolds number, Prandtl number, Eckert number, mixed convection parameter, velocity slip, and thermal slip on the velocity and temperature profiles, local skin friction coefficient, and local Nusselt number are discussed.     

Influence of chemical reaction on heat and mass transfer of non-Newtonian fluid with yield stress by free convection from vertical surface in porous medium considering Soret effect

The effect of chemical reaction on free convection heat and mass transfer for a non-Newtonian power law fluid over a vertical flat plate embedded in a fluid-saturated porous medium has been studied in the presence of the yield stress and the Soret effect. The governing boundary layer equations and boundary conditions are cast into a dimen- sionless form by similarity transformations, and the resulting system of equations is solved by a finite difference method. The results are preSented and discussed for concentration profiles, as well as the Nusselt number and the Sherwood number for various values of the parameters, which govern the problem. The results obtained show that the flow field is influenced appreciably by the presence of the chemical reaction parameter γ the order of.the chemical reaction parameter m, the Soret number St, the buoyancy ratio N, the Lewis number Le, and the dimensionless rheological parameter Ω.     

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.     

Mixed Convection Boundary Layer Flow Near the Stagnation Point on a Vertical Surface Embedded in a Porous Medium with Anisotropy Effect

The steady boundary-layer flow near the stagnation point on a vertical flat plate embedded in a fluid-saturated porous medium characterized by an anisotropic permeability is investigated. Using appropriate similarity transformation, 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 modified mixed convection parameter Λ, and the anisotropy parameter A are analyzed and discussed. It is found that dual solutions exist for both assisting and opposing flows. Moreover, the range of Λ for which the solution exists increases with A.     

Effect of weak rotation on natural convection in a horizontal porous cylinder

Natural convection in a fluid saturated porous medium confined in a horizontal circular cylinder and rotating about its axis, with isothermal boundary conditions and uniform internal heat sink, is studied by both numerical and perturbation methods. No symmetry with respect to the vertical diameter is expected for the flow and temperature fields and the whole region must be involved in the computation. Only the weak rotation regime, for which the centrifugal force is negligible compared to gravity, is considered. Governing equations for the two-dimensional flow field are solved in both rotating and non-rotating coordinate systems. Results indicate that rotation significantly decreases the radial amplitude of fluid particle trajectories in the radial direction and thus reduces the overall heat transfer.     

Effects of chemical reaction on mixed convection flow of a polar fluid through a porous medium in the presence of internal heat generation

This paper reports a detailed numerical investigation on mixed convection flow of a polar fluid through a porous medium due to the combined effects of thermal and mass diffusion. The energy equation accounts for heat generation or absorption, while the nth order homogeneous chemical reaction between the fluid and the diffusing species is included in the mass diffusion equation. The governing equations of the linear momentum, angular momentum, energy and concentration are obtained in a non-similar form by introducing a suitable group of transformations. The final set of non-similar coupled non-linear partial differential equations is solved using an implicit finite-difference scheme in combination with quasi-linearization technique. The effects of various parameters on the velocity, angular velocity, temperature and concentration fields are investigated. Numerical results for the skin friction coefficient, wall stress of angular velocity, Nusselt number and Sherwood number are also presented.     

基于POD-Galerkin降维方法的热毛细对流分岔分析

作为流动与传热相互耦合的非线性过程, 热毛细对流有着复杂的转捩过程, 探究流场和温度场随参数变化而发生的分岔现象, 是热毛细对流研究的一个重要课题. 基于本征正交分解的POD-Galerkin降维方法可以通过提取特征模态, 构建低维模型, 实现流场的快速计算. 数值分岔方法可以通过求解含参数动力系统的分岔方程, 直接计算稳定解和分岔点. 探究了将直接数值模拟方法、POD-Galerkin降维方法、数值分岔方法的优势结合, 以提高热毛细对流转捩过程分析效率的可行性. 利用直接数值模拟得到的流场和温度场数据, 构建了不同体积比下, 二维有限长液层热毛细对流的POD-Galerkin低维模型, 在低维模型上采用数值积分及数值分岔方法计算了分岔点, 得到了低维方程的分岔图. 在一定参数范围内, 在低维模型上模拟热毛细对流, 对雷诺数和体积比进行参数外推, 通过与直接数值模拟的结果对比, 验证了低维模型的准确性与鲁棒性. 说明了低维方程可以定性反映原高维系统的流动特性, 而定量方面, 由低维模型和直接数值模拟计算得到的周期解频率的相对误差大约为5%. 验证了利用POD-Galerkin降维方法研究热毛细对流的可行性.      

Mixed Convection Boundary Layer Flow Past a Horizontal Circular Cylinder Embedded in a Bidisperse Porous Medium

Adopting a two-temperature and two-velocity model, appropriate to a bidisperse porous medium (BDPM) proposed by Nield and Kuznetsov (2008), the classical steady, mixed convection boundary layer flow about a horizontal, isothermal circular cylinder embedded in a porous medium has been theoretically studied in this article. It is shown that the boundary layer analysis leads to expressions for the flow and heat transfer characteristics in terms of an inter-phase momentum parameter, a thermal diffusivity ratio, a thermal conductivity ratio, a permeability ratio, a modified thermal capacity ratio, and a buoyancy or mixed convection parameter. The transformed partial differential equations governing the flow and heat transfer in the f-phase (the macro-pores) and the p-phase (the remainder of the structure) are solved numerically using a very efficient implicit finite-difference technique known as Keller-box method. A good agreement is observed between the present results and those known from the open literature in the special case of a traditional Darcy formulation (monodisperse system).     

On the stability of natural convection in a porous vertical slab saturated with an Oldroyd-B fluid

The stability of the conduction regime of natural convection in a porous vertical slab saturated with an Oldroyd-B fluid has been studied. A modified Darcy’s law is utilized to describe the flow in a porous medium. The eigenvalue problem is solved using Chebyshev collocation method and the critical Darcy–Rayleigh number with respect to the wave number is extracted for different values of physical parameters. Despite the basic state being the same for Newtonian and Oldroyd-B fluids, it is observed that the basic flow is unstable for viscoelastic fluids—a result of contrast compared to Newtonian as well as for power-law fluids. It is found that the viscoelasticity parameters exhibit both stabilizing and destabilizing influence on the system. Increase in the value of strain retardation parameter \(\Lambda _2 \) portrays stabilizing influence on the system while increasing stress relaxation parameter \(\Lambda _1\) displays an opposite trend. Also, the effect of increasing ratio of heat capacities is to delay the onset of instability. The results for Maxwell fluid obtained as a particular case from the present study indicate that the system is more unstable compared to Oldroyd-B fluid.     

Free Convection Heat Transfer From A Sphere In A Porous Medium Using A Thermal Non-equilibrium Model

This work studies the free convection heat transfer from a sphere with constant wall temperature embedded in a fluid-saturated porous medium using a thermal non-equilibrium model. The governing equations are transformed into boundary-layer partial differential equations by the coordinate transform, and the obtained governing equations are then solved by the cubic spline collocation method. The temperature distributions for fluid and solid phases are shown for different values of the porosity scaled thermal conductivity ratio, the interphase heat transfer parameter, and the streamwise coordinate. The effects of the porosity scaled thermal conductivity ratio and the interphase heat transfer parameter between solid and fluid phases on the local Nusselt numbers for fluid and solid phases are examined. Results show the local Nusset number for the porous medium can be increased by increasing the porosity scaled thermal conductivity ratio. Moreover, the thermal non-equilibrium effect is more significant for low values of the porosity scaled thermal conductivity ratio or the interphase heat transfer parameter.     

The Onset of Double-Diffusive Convection in a Superposed Fluid and Porous Layer Under High-Frequency and Small-Amplitude Vibrations

We numerically simulate the initiation of an average convective flow in a system composed of a horizontal binary fluid layer overlying a homogeneous porous layer saturated with the same fluid under gravitational field and vibration. In the layers, fixed equilibrium temperature and concentration gradients are set. The layers execute high-frequency oscillations in the vertical direction. The vibration period is small compared with characteristic timescales of the problem. The averaging method is applied to obtain vibrational convection equations. Using for computation the shooting method, a numerical investigation is carried out for an aqueous ammonium chloride solution and packed glass spheres saturated with the solution. The instability threshold is determined under two heating conditions—on heating from below and from above. When the solution is heated from below, the instability character changes abruptly with increasing solutal Rayleigh number, i.e., there is a jump-wise transition from the most dangerous shortwave perturbations localized in the fluid layer to the long-wave perturbations covering both layers. The perturbation wavelength increases by almost 10 times. Vibrations significantly stabilize the fluid equilibrium state and lead to an increase in the wavelength of its perturbations. When the fluid with the stabilizing concentration gradient is heated from below, convection can occur not only in a monotonous manner but also in an oscillatory manner. The frequency of critical oscillatory perturbations decreases by 10 times, when the long-wave instability replaces the shortwave instability. When the fluid is heated from above, only stationary convection is excited over the entire range of the examined parameters. A lower monotonic instability level is associated with the development of perturbations with longer wavelength even at a relatively large fluid layer thickness. Vibrations speed up the stationary convection onset and lead to a decrease in the wavelength of most dangerous perturbations of the motionless equilibrium state. In this case, high enough amplitudes of vibration are needed for a remarkable change in the stability threshold. The results of numerical simulation show good agreement with the data of earlier works in the limiting case of zero fluid layer thickness.