Effect of rotation on the onset of double diffusive convection in a Darcy porous medium saturated with a couple stress fluid

The effect of rotation on the onset of double diffusive convection in a horizontal couple stress fluid-saturated porous layer, which is heated and salted from below, is studied analytically using both linear and weak nonlinear stability analyses. The extended Darcy model, which includes the time derivative and Coriolis terms, has been employed in the momentum equation. The onset criterion for stationary, oscillatory and finite amplitude convection is derived analytically. The effect of Taylor number, couple stress parameter, solute Rayleigh number, Lewis number, Darcy–Prandtl number, and normalized porosity on the stationary, oscillatory, and finite amplitude convection is shown graphically. It is found that the rotation, couple stress parameter and solute Rayleigh number have stabilizing effect on the stationary, oscillatory, and finite amplitude convection. The Lewis number has a stabilizing effect in the case of stationary and finite amplitude modes, with a destabilizing effect in the case of oscillatory convection. The Darcy–Prandtl number and normalized porosity advances the onset of oscillatory convection. 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. The transient behavior of the Nusselt number and Sherwood number is investigated by solving the finite amplitude equations using Runge–Kutta method.     

An analytical study of linear and nonlinear double diffusive convection in a fluid saturated anisotropic porous layer with Soret effect

The double diffusive convection in a horizontal anisotropic porous layer saturated with a Boussinesq fluid, which is heated and salted from below in the presence of Soret coefficient is studied analytically using both linear and nonlinear stability analyses. The normal mode technique is used in the linear stability analysis while a weak nonlinear analysis based on a minimal representation of double Fourier series method is used in the nonlinear analysis. The generalized Darcy model including the time derivative term is employed for the momentum equation. The critical Rayleigh number, wavenumber for stationary and oscillatory modes and frequency of oscillations are obtained analytically using linear theory. The effect of anisotropy parameters, solute Rayleigh number, Soret parameter and Lewis number on the stationary, oscillatory, finite amplitude convection and heat and mass transfer are shown graphically.     

Lattice Boltzmann simulation of convective heat transfer of non-Newtonian fluids in impeller stirred tank

Lattice Boltzmann simulation of convective heat transfer of non-Newtonian fluids in an impeller stirred tank is performed. The curved and moving boundary methods combined with the unknown-index algorithm are used to solve the flow and thermal fields induced in a cold tank by an oscillating hot impeller. For a given maximum radius of the blades, the simulation results show that a rectangular impeller of large aspect ratio induces stronger heat transfer effect on the tank walls than the small aspect ratio. This is because the latter would cause worse field synergy than the former, i.e. the induced local velocities of fluid are mostly perpendicular to the temperature gradients. The convection effects on the tank walls are also improved as the oscillation amplitude of impeller increases until the swept areas of impeller are close to whole azimuth of the tank, i.e., oscillation amplitude of 90°. The maximum Nusselt number on the tank walls for power-law fluid flows of n = 0.7, 1 and 1.5 occurs at oscillation amplitude of 75°. Finally, it is found that the heat transfer effect on the tank walls is reduced as the power-law index of fluid increases.     

Theoretical model of laminar flow in tubes in rolling motion

The theoretical model of laminar flow in tubes in rolling motion is established. The velocity and temperature correlations are derived, and the frictional resistance coefficient and Nusselt number are also obtained. The oscillation of parameters is induced by the tangential force due to rolling motion. The effect of centrifugal and Coriolis forces on the flow is negligible. The tangential force does not effect on the average parameters. The oscillating amplitude of Nusselt number increase with the Prandtl number increasing. Both the oscillating amplitudes of frictional resistance coefficient and Nusselt number increase with the rolling frequency increasing.     

基于Monte-Carlo随机有限元方法的随机边界条件下自然对流换热不确定性研究

为分析边界条件不确定性对方腔内自然对流换热的影响,发展了一种求解随机边界条件下自然对流换热不确定性传播的Monte-Carlo随机有限元方法.通过对输入参数场随机边界条件进行Karhunen-Loeve展开及基于Latin(拉丁)抽样法生成边界条件随机样本,数值计算了不同边界条件随机样本下方腔内自然对流换热流场与温度场,并用采样统计方法计算了随机输出场的平均值与标准偏差.根据计算框架编写了求解随机边界条件下方腔内自然对流换热不确定性的MATLAB随机有限元程序,分析了随机边界条件相关长度与方差对自然对流不确定性的影响.结果表明:平均温度场及流场与确定性温度场及流场分布基本相同;随机边界条件下Nu数概率分布基本呈现正态分布,平均Nu数随着相关长度和方差增加而增大;方差对自然对流换热的影响强于相关长度的影响.     

Mixed convection of heat generating nanofluid in a lid-driven cavity with uniform and non-uniform heating of bottom wall

Steady state two-dimensional mixed convection in a lid-driven square cavity filled with Cu  –water nanofluid is investigated numerically in the presence of internal heat generation. In the present investigation, bottom wall is uniformly and non-uniformly heated while two vertical walls are fixed and they are thermally insulated. The top wall is moving from left to right at a constant speed. The governing equations are normalized and solved numerically with boundary conditions by finite volume approach using third order accurate upwinding scheme (deferred QUICK). Effects of the pertinent physical parameters are investigated in terms of the flow and temperature fields, as well as Nusselt number distributions. The presented results show that the solid volume fraction plays a significant role on the flow and thermal fields and the Nusselt number distributions for different flow configurations. It is found that Richardson number strongly affect the fluid flow and heat transfer in the cavity. For Ri<1$\mathit{Ri}<1$, the forced convection becomes dominant in the entire cavity, the natural convection relatively weak.     

Buoyancy-assisted flow of yield stress fluids past a cylinder: Effect of shape and channel confinement

The aiding-buoyancy mixed convection heat transfer in Bingham plastic fluids from an isothermal cylinder of elliptical and circular shape in a vertical adiabatic channel is numerically investigated. For a fixed shape of the elliptical cylinder E = 2 (ratio of major to minor axes), the effect of confinement is studied for three values of blockage ratio, B, defined as the ratio of the channel width to the circumference of the cylinder/π, as 6.5, 2.17 and 1.3. In order to delineate the role of cross-section of the cylinder, results are also presented here for a circular cylinder of the same heat transfer area as the elliptical cylinder. The results presented herein span the range of conditions as: Bingham number, 0 ≤ Bn ≤ 100, Reynolds number, 1 ≤ Re ≤ 40, and Prandtl number, 1 ≤ Pr ≤ 100 over the range of Richardson number Ri = 0 (pure forced convection) to Ri = 10. Extensive results on drag coefficient, local and surface averaged values of the Nusselt number and yield surfaces are presented herein to elucidate the combined effects of buoyancy, blockage ratio and fluid yield stress. The morphology of the yield surfaces shows that the unyielded plug regions formed upstream and downstream of the cylinder grow faster at low Reynolds numbers with the increasing yield stress effects under the weak buoyancy forces, i.e., small values of Grashof or Richardson number. The heat transfer enhancement is observed with the increasing channel-confinement due to the sharpening of the temperature gradients near the surface of the cylinder. The average Nusselt number shows a positive dependence on the Reynolds number, Prandtl number and Richardson number irrespective of the shape of the cylinder or the type of fluid. By employing the modified definitions of the dimensionless parameters (based on the two choices of the overall effective fluid velocity), predictive correlations have been established for estimating the value of the average Nusselt number in a new application.     

A lattice Boltzmann analysis of the conjugate natural convection in a square enclosure with a circular cylinder

In this paper, an improved lattice Boltzmann model that incorporates the effect of heat capacity is adopted to simulate the conjugate natural convection in a square enclosure with a circular cylinder, and the square cavity considered is bounded by four finite thickness and conductive walls. Influences of Rayleigh number, wall thermophysical properties, and wall thickness on streamlines, isotherms and Nusselt number are analyzed. It is found that the heat transfer rate obtained for different wall thicknesses is usually smaller than that for the zero wall thickness, while this difference is not apparent for a larger thermal conductivity ration. Additionally, the study also shown that the increase of the outer wall thickness tends to decrease the heat transfer rate. The effects of the geometric parameter, including aspect ratio and vertical cylinder location, on heat transfer are also studied.     

An analysis of turbulent film condensation on a sphere with variable wall temperature under the effect of local shear stress

This paper presents an investigation into turbulent film condensation on a sphere with variable wall temperature. Under the wide range of vapor velocity, the wall temperature and the local film shear stress were considered. The result shows that under the high velocity vapor, the increase of the temperature amplitude will bring out a larger Nusselt number, and the increase is about 2.7–5.6%. Besides, under the effect of the local film shear stress, the mean Nusselt number will decrease about 0.65–0.8%. Furthermore, the paper then discusses the influence of shears and temperature amplitudes on the local dimensionless film thickness and heat transfer characteristics. Finally, the results developed in the current study are compared with those generated by previous theoretical results.     

Stationary statistical properties of Rayleigh‐Bénard convection at large Prandtl number

This is the third in a series of our study of Rayleigh‐Bénard convection at large Prandtl number. Here we investigate whether stationary statistical properties of the Boussinesq system for Rayleigh‐Bénard convection at large Prandtl number are related to those of the infinite Prandtl number model for convection that is formally derived from the Boussinesq system via setting the Prandtl number to infinity. We study asymptotic behavior of stationary statistical solutions, or invariant measures, to the Boussinesq system for Rayleigh‐Bénard convection at large Prandtl number. In particular, we show that the invariant measures of the Boussinesq system for Rayleigh‐Bénard convection converge to those of the infinite Prandtl number model for convection as the Prandtl number approaches infinity. We also show that the Nusselt number for the Boussinesq system (a specific statistical property of the system) is asymptotically bounded by the Nusselt number of the infinite Prandtl number model for convection at large Prandtl number. We discover that the Nusselt numbers are saturated by ergodic invariant measures. Moreover, we derive a new upper bound on the Nusselt number for the Boussinesq system at large Prandtl number of the form which asymptotically agrees with the (optimal) upper bound on Nusselt number for the infinite Prandtl number model for convection. © 2007 Wiley Periodicals, Inc.     

柔性圆柱形微管道内的电动流动及传热研究

研究了在纯压力驱动下,流体通过壁面带有某种电荷的聚电解质层（PEL）的微管道，即柔性微管道的电动流动和热传输特性.基于先前得到的电势和速度的解析解以及流向势的数值解,在热充分发展的情况下, 假设壁面热流恒定,利用有限差分法求解了包括黏性耗散和Joule(焦耳)热影响下的能量方程，获得了无量纲温度数值解.通过数值计算，给出了相关的无量纲参数对速度、温度以及Nusselt(努赛尔)数的影响.研究表明,当其他参数固定时,无量纲速度和温度随着无量纲聚电解质层厚度d的增大而减小,随着聚电解质层中等效双电层厚度与双电层厚度之比K_λ的增大而增大；Nusselt数随着Joule热系数S的增大而减小,随无量纲聚电解质层厚度d的增大而减小,随着K_λ的增大而增大.     

Viscous dissipation and Joule heating effects on MHD-free convection from a vertical plate with power-law variation in surface temperature in the presence of Hall and ion-slip currents

An analysis is presented to study the effects of viscous dissipation and Joule heating on MHD-free convection flow past a semi-infinite vertical flat plate in the presence of the combined effect of Hall and ion-slip currents for the case of power-law variation of the wall temperature. The fluid is permeated by a strong transverse magnetic field imposed perpendicularly to the plate on the assumption of a small magnetic Reynolds number. The governing differential equations are transformed by introducing proper non-similarity variables and solved numerically. The effects of various parameters on the velocity and temperature profiles as well as the local wall shear stresses and the local Nusselt number are presented graphically and in tabular form. It is found that the magnetic field acts as a retarding force on the tangential flow but has a propelling effect on the induced lateral flow. The skin-friction factor for the tangential flow and the local Nusselt number decrease but the skin-friction factor for the lateral flow increases as the magnetic field increases. The skin-friction factor for the tangential and lateral flows are increased while the local Nusselt number is decreased if the effect of viscous dissipation, Joule heating and heat generation are considered. The opposite trend was observed as the temperature power coefficient n is increased. Also, the skin-friction factor for the tangential flow and the local Nusselt number are increased due to the Hall and ion-slip currents, whereas the skin-friction factor for the tangential flow increases when Hall values increase to one and decreases for values of Hall greater than one, but reduces by rising ion-slip values.     

Numerical investigation of natural convection in a rectangular enclosure due to partial heating and cooling at vertical walls

A comprehensive numerical investigation on the natural convection in a rectangular enclosure is presented. The flow is induced due to the constant partial heating at lower half of the left vertical wall and partial cooling at upper half of the right vertical wall along with rest walls are adiabatic. In this investigation the Special attention is given to understand the effect of aspect ratio and heat source intensity i.e. Rayleigh number, Ra, on the fluid flow configuration as well as on the local and average heat transfer rates. The range of Rayleigh (Ra) and aspect ratio (A) is taken [10³, 10⁶] and [0.5, 4] respectively. The results are presented in terms of stream function (ψ), temperature (θ) and heat transfer rates (local Nusselt numbers Nu_L, and average Nusselt numbers Nu). The numerical experiments show that increasing of Ra implies the enhancement of thermal buoyancy force, which in turn increases the thermal convection in the cavity. As a result, the local as well as average heat transfer rate is expected to increase. The local transfer rate (Nu_L) is increases in the small region near the left vertical wall of the left wall of the cavity and after that it is decreases in the middle portion of heated region. And, it start to increase near to the middle point of left wall. It is also observed that the local heat transfer is increases as increases the aspect ratio. The average heat transfer rate (Nu) is increases as the aspect ratio A increases from 0.5 to 1 and beyond that it is decreases smoothly. It is also found that the heat transfer rate attains its maximum value at aspect ratio one.     

Effects of variable viscosity on non-Darcy MHD free convection along a non-isothermal vertical surface in a thermally stratified porous medium

An analysis is performed for non-Darcy free convection flow of an electrically conducting fluid over an impermeable vertical plate embedded in a thermally stratified, fluid saturated porous medium for the case of power-law surface temperature. The present work examines the effects of non-Darcian flow phenomena, variable viscosity, Hartmann–Darcy number and thermal stratification on free convective transport and demonstrates the variation in heat transfer prediction based on three different flow models. The wall effect on porosity variation is approximated by an exponential function. The effects of thermal dispersion and variable stagnant thermal conductivity are taken into consideration in the energy equation. The resulting non-similar system of equations is solved using a finite difference method. Results are presented for velocity, temperature profiles and local Nusselt number for representative values of different controlling parameters.     

Effect of magnetic field dependent viscosity on ferroconvection in the presence of dust particles

This paper deals with the theoretical investigation of the effect of magnetic field dependent (MFD) viscosity on the thermal convection in a ferromagnetic fluid in the presence of dust particles. For a flat ferromagnetic fluid layer contained between two free boundaries, the exact solution is obtained using a linear stability analysis and a normal mode analysis method. For the case of stationary convection, dust particles always have a destabilizing effect, whereas the MFD viscosity has a stabilizing effect on the onset of convection. In the absence of MFD viscosity, the destabilizing effect of magnetization is depicted but in the presence of MFD viscosity, non-buoyancy magnetization may have a destabilizing or a stabilizing effect on the onset of convection. The critical wave number and critical magnetic thermal Rayleigh number for the onset of stationary convection are also determined numerically for sufficiently large values of buoyancy magnetization parameter M ₁. Graphs have been plotted by giving numerical values to the parameters to depict the stability characteristics. It is observed that the critical magnetic thermal Rayleigh number is reduced solely because the heat capacity of clean fluid is supplemented by that of the dust particles. The principle of exchange of stabilities is found to hold true for the ferromagnetic fluid heated from below in the absence of dust particles. The oscillatory modes are introduced due to the presence of the dust particles, which were non-existent in their absence. A sufficient condition for the non-existence of overstability is also obtained.     

Mixed convection of a viscous dissipating fluid about a vertical flat plate

In this study, the effect of the viscous dissipation in steady, laminar mixed convection heat transfer from a heated/cooled vertical flat plate is investigated in both aiding and opposing buoyancy situations. The external flow field is assumed to be uniform. The governing systems of partial differential equations are solved numerically using the finite difference method. A parametric study is performed in order to illustrate the interactive influences of the governing parameters, mainly, the Richardson number, Ri (also known as the mixed convection parameter) and the Eckert number, Ec on the velocity and temperature profiles as well as the friction and heat transfer coefficients. Based on the facts the free stream is either in parallel or reverse to the gravity direction and the plate is heated or cooled, different flow situations are identified. The influence of the viscous dissipation on the heat transfer varied according to the situation. For some limiting cases, the obtained results are validated by comparing with those available from the existing literature. An expression correlating Nu in terms of Pr, Ri and Ec is developed.     

MHD-mixed convection from a vertical slender cylinder

The problem of steady laminar magnetohydrodynamic (MHD) mixed convection heat transfer about a vertical slender cylinder is studied numerically. A uniform magnetic field is applied perpendicular to the cylinder. The resulting governing equations are transformed into the non-similar boundary layer equations and solved using the Keller box method. The velocity and temperature profiles as well as the local skin friction and the local heat transfer parameters are determined for different values of the governing parameters, mainly the transverse curvature parameter, the magnetic parameter, the electric field parameter and the Richardson number. For some specific values of the governing parameters, the results agree very well with those available in the literature. Generally, it is determined that the local skin friction coefficient and the local heat transfer coefficient increase, increasing the Richardson number, Ri (i.e. the mixed convection parameter), electric field parameter E₁ and magnetic parameter Mn.     

Effect of dust particles on a layer of micropolar ferromagnetic fluid heated from below saturating a porous medium

The problem of the effect of dust particles on the thermal convection in micropolar ferromagnetic fluid saturating a porous medium subject to a transverse uniform magnetic field has been investigated theoretically. Linear stability analysis and normal mode analysis methods are used to find an exact solution for a flat micropolar ferromagnetic fluid layer contained between two free boundaries. In case of stationary convection, the effect of various parameters like medium permeability, dust particles, non-buoyancy magnetization, coupling parameter, spin-diffusion parameter and micropolar heat conduction parameter are analyzed. For sufficiently large values of magnetic parameter M₁, the critical magnetic thermal Rayleigh number for the onset of instability is determined numerically and results are depicted graphically. It is also observed that the critical magnetic thermal Rayleigh number is reduced solely because the heat capacity of clean fluid is supplemented by that of the dust particles. The principle of exchange of stabilities is found to hold true for the micropolar ferromagnetic fluid saturating a porous medium heated from below in the absence of micropolar viscous effect, microinertia and dust particles.     

Influence of variable heat flux on natural convection along a corrugated wall in porous media

In this work the coupled non-linear partial differential equations, governing the free convection from a wavy vertical wall under a power law heat flux condition, are solved numerically. For both Darcy and Forchheimer extended non-Darcy models, a wavy to flat surface transformation is applied and the governing equations are reduced to boundary layer equations. A finite difference scheme based on the Keller Box approach has been used in conjunction with a block tri-diagonal solver for obtaining the solution. Detailed simulations are carried out to investigate the effect of varying parameters such as power law heat flux exponent m, wavelength–amplitude ratio a and the transformed Grashof number Gr′. Both surface undulations and inertial forces increase the temperature of the vertical surface while increasing m reduces it. The wavy pattern observed in surface temperature plots, become more prominent with increasing m or a but reduces as Gr′ increases.     

Effects of temperature-dependent viscosity and variable thermal conductivity on MHD non-Darcy mixed convective diffusion of species over a stretching sheet

A numerical model is presented to study the effects of temperature-dependent viscosity and variable thermal conductivity on mixed convection problem. Two important types of wall heating conditions namely, prescribed surface temperature and prescribed wall heat flux which arise in polymer industries are considered. The problem is solved numerically by using the fifth-order Runge–Kutta Fehlberg method with shooting technique. It is found that the Prandtl number is to decrease the skin friction coefficient, local Nusselt number and local Sherwood number. The effects of non-uniform heat source/sink and porous parameter are analyzed on velocity, temperature, skin friction co-efficient, Nusselt and Sherwood numbers.