均匀自由流动的非牛顿流体中连续表面上的磁流体动力学流动和热传递

分析在平行自由流动的非牛顿黏弹性导电流体中,连续平展表面移动时的稳态流和热传递特性,该流动处于横向均匀磁场作用下.以二阶流体构建它的本构方程,得到了速度分布和温度断面图的数值结果.讨论了诸如黏弹性参数、磁场参数和Prandtl数等不同物理参数对诸种动量和热传递特性的影响,并给出相关图示.     

Convective heat and mass transfer in three-dimensional mixed convection flow of viscoelastic fluid in presence of chemical reaction and heat source/sink

Heat and mass transfer effects in the three-dimensional mixed convection flow of a viscoelastic fluid with internal heat source/sink and chemical reaction have been investigated in the present work. The flow generation is because of an exponentially stretching surface. Magnetic field normal to the direction of flow is considered. Convective conditions at the surface are also encountered. Appropriate similarity transformations are utilized to reduce the boundary layer partial differential equations into the ordinary differential equations. The homotopy analysis method is used to develop the solution expressions. Impacts of different controlling parameters such as ratio parameter, Hartman number, internal heat source/sink, chemical reaction, mixed convection, concentration buoyancy parameter and Biot numbers on the velocity, temperature and concentration profiles are analyzed. The local Nusselt and Sherwood numbers are sketched and examined.     

Two-dimensional flow of chemically reactive viscoelastic fluids with or without the influence of thermal convection

Shear and temperature induced dynamics of a viscoelastic liquid is investigated with direct numerical simulations. A numerical algorithm based on the finite difference method is implemented in time and space with the Oldroyd-B constitutive model for the viscoelastic liquids. The temperature effects are modeled by a combination of Arrhenius kinetics and thermal convection. The standard approach which applies periodic boundary conditions is replaced by one in which the unprescribed boundary values are regenerated from the main flow. The algorithm is validated against well-known Newtonian natural convection flow results. The results show that progressive viscoelasticity leads to lower temperature increases within the fluid system, and that the highest temperature increases are recorded for corresponding Newtonian fluids.     

The effect of variable viscosity on MHD viscoelastic fluid flow and heat transfer over a stretching sheet

An analysis has been carried out to study the momentum and heat transfer characteristics in an incompressible electrically conducting non-Newtonian boundary layer flow of a viscoelastic fluid over a stretching sheet. The partial differential equations governing the flow and heat transfer characteristics are converted into highly non-linear coupled ordinary differential equations by similarity transformations. The effect of variable fluid viscosity, Magnetic parameter, Prandtl number, variable thermal conductivity, heat source/sink parameter and thermal radiation parameter are analyzed for velocity, temperature fields, and wall temperature gradient. The resultant coupled highly non-linear ordinary differential equations are solved numerically by employing a shooting technique with fourth order Runge–Kutta integration scheme. The fluid viscosity and thermal conductivity, respectively, assumed to vary as an inverse and linear function of temperature. The analysis reveals that the wall temperature profile decreases significantly due to increase in magnetic field parameter. Further, it is noticed that the skin friction of the sheet decreases due to increase in the Magnetic parameter of the flow characteristics.     

Effects of thermal radiation and magnetic field on unsteady mixed convection flow and heat transfer over an exponentially stretching surface with suction in the presence of internal heat generation/absorption

In this paper, the problem of unsteady laminar two-dimensional boundary layer flow and heat transfer of an incompressible viscous fluid in the presence of thermal radiation, internal heat generation or absorption, and magnetic field over an exponentially stretching surface subjected to suction with an exponential temperature distribution is discussed numerically. The governing boundary layer equations are reduced to a system of ordinary differential equations. New numerical method using Mathematica has been used to solve such system after obtaining the missed initial conditions. Comparison of obtained numerical results is made with previously published results in some special cases, and found to be in a good agreement.     

Analysis of viscoelastic fluid flow with temperature dependent properties in plane Couette flow and thin annuli

The steady state flow in very thin annuli has been studied analytically for the case where the annular gap is much smaller than the radius of the inner cylinder and for the outer cylinder rotating at constant angular speed and the inner cylinder at rest. The cylinders were subjected to two different thermal boundary conditions. The exponential effect of temperature on the relaxation time and the viscosity coefficient was accounted into the governing differential equations using Nahme’s law. Effects of viscous dissipation as well as εDe² (viscoelastic index for SPTT constitutive equation) on the dimensionless velocity and temperature profiles have been investigated. Results show that while the properties of the fluid depend on temperature, the velocity and temperature profiles are different compared to those obtained with constant physical properties. The Nahme–Griffith number increases whereas εDe² as a viscoelastic index decreases when temperature dependent physical properties are considered. In addition, the results indicate that the viscous dissipation has a sensible effect on heat transfer and the Nusselt number decreases with an increase in the Nahme–Griffith number.     

Heat transfer in MHD viscoelastic boundary layer flow over a stretching sheet with thermal radiation and non-uniform heat source/sink

An analysis has been carried out to study the flow and heat transfer characteristics for MHD viscoelastic boundary layer flow over an impermeable stretching sheet with space and temperature dependent internal heat generation/absorption (non-uniform heat source/sink), viscous dissipation, thermal radiation and magnetic field due to frictional heating. The flow is generated due to linear stretching of the sheet and influenced by uniform magnetic field, which is applied vertically in the flow region. The governing partial differential equations for the flow and heat transfer are transformed into ordinary differential equations by a suitable similarity transformation. The governing equations with the appropriate conditions are solved exactly. The effects of viscoelastic parameter and magnetic parameter on skin friction and the effects of viscous dissipation, non-uniform heat source/sink and the thermal radiation on heat transfer characteristics for two general cases namely, the prescribed surface temperature (PST) case and the prescribed wall heat flux (PHF) case are presented graphically and discussed. The numerical results for the wall temperature gradient (the Nusselt number) are presented in tables and are discussed.     

Nonaxisymmetric Free Convection Flow over a Rotating Disk in a Viscoelastic fluid with Magnetic Field

The non axisymmetric motion produced by a buoyancy-induced secondary flow of a viscoelastic fluid over an infinite rotating disk in a verticalplane with a magnetic field applied normal to the disk has been studied.The governing Navier Stokes equations and the energy equation admit a self similar solution. The system of ordinary differential equations has been solved numerically using Runge-Kutta Gill subroutine.The turning moment for the viscoelastic fluid is found to be less than that of the Newtonian fluid but the turning moment is increased due to the magnetic parameter. The resultant force due to the buoyancy-induced secondary flow increases with the magnetic parameter but reduces as the viscoelastic parameter increases. The quantity of fluid, which is pumped outwards due to the centrifuging action of the disk, for the viscoelastic fluid is more than that of the Newtonian fluid. The buoyancy-induced secondary flow boundary layer is much thicker than the primary boundary layer thickness. The thermal boundary layer due to the primary flow increases with the magnetic parameter decreases as the viscoelastic parameter increases. The heat transfer increases with the viscoelastic parameter but decreases as the magnetic parameter increases. The effect of the viscoelastic parameter is more pronounced on the secondary flow than on the primary flow.     

Effect of Hall current on transient hydromagnetic Couette–Poiseuille flow of a viscoelastic fluid with heat transfer

The unsteady Couette–Poiseuille flow of an electrically conducting incompressible non-Newtonian viscoelastic fluid between two parallel horizontal non-conducting porous plates is studied with heat transfer considering the Hall effect. A sudden uniform and constant pressure gradient, an external uniform magnetic field that is perpendicular to the plates and uniform suction and injection through the surface of the plates are applied. The two plates are kept at different but constant temperatures while the Joule and viscous dissipations are taken into consideration. Numerical solutions for the governing momentum and energy equations are obtained using finite difference approximations. The effect of the Hall term, the parameter describing the non-Newtonian behavior, and the velocity of suction and injection on both the velocity and temperature distributions is examined.     

Effect of variable density on hydromagnetic mixed convection flow of a non-Newtonian fluid past a moving vertical plate

The effect of temperature-dependent density on MHD mixed convection flow of power-law fluid past a moving semi-infinite vertical plate for high temperature differences between the plate and the ambient fluid is studied. The fluid density is assumed to decrease exponentially with temperature. The usual Boussinesq approximations are not considered due to the large temperature differences. The surface temperature of the moving plate was assumed to vary according to a power-law form, that is, T_w(x) = T_∞ + Ax^γ. The fluid is permeated by a uniform magnetic field imposed perpendicularly to the plate on the assumption of small magnetic Reynolds number. A numerical shooting algorithm for two unknown initial conditions with fourth-order Runge–Kutta integration scheme has been used to solve the coupled non-linear boundary value problem. The effects of various parameters on the velocity and temperature profiles as well as the local skin-friction coefficient and the local Nusselt number are presented graphically and in the tabular form. The results show that application of Boussinesq approximations in a non-Newtonian fluid subjected to high temperature differences gives a significant error in the values of the skin-friction coefficient and the application of an external magnetic field reduces this error markedly in the case of shear-thickening fluid.     

放/吸热流体在两个充满多孔材料的竖直平行板之间作不稳定的自然对流

在一个由两块无限竖直平行板组成的管道中,充满着多孔的介质材料,使用Darcy模型(Brinkman模型的推广)的动量方程,连同能量方程,计算不可压缩、粘性、放/吸热流体在该管道中的不稳定自然对流,即Couette流动.流动是由于边界平板有不对称的加热,以及作加速运动所引起.选用合理的无量纲参数,对控制方程进行简化,通过Laplace变换进行解析求解,得到闭式的速度和温度分布曲线解,随后导出表面摩擦力和传热率.发现在竖直管道中的不同剖面,流体的流动及温度分布曲线随着时间而增加,且在运动平板附近更高.特别是,流体的速度和温度随着平板间距的增加而增加,但是,表面摩擦力和热传导率随着平板间距的增加而减小.     

Analytical solution of stagnation-point flow of a viscoelastic fluid towards a stretching surface

In this paper the problem of stagnation-point flow of a viscoelastic fluid towards a stretching surface [T.R. Mahapatra, A.S. Gupta, Stagnation-point flow of a viscoelastic fluid towards a stretching surface, Int. J. Non-Linear Mech. 39 (2004) 811] is solved analytically by using the homotopy analysis method (HAM). The results for velocity and temperature profiles are obtained. It is noted that the behavior of the HAM solution for velocity and temperature profiles is in good agreement with the numerical solution given in reference [T.R. Mahapatra, A.S. Gupta, Stagnation-point flow of a viscoelastic fluid towards a stretching surface, Int. J. Non-Linear Mech. 39 (2004) 811].     

Effects of partial slip,viscous dissipation and Joule heating on Von Kármán flow and heat transfer of an electrically conducting non-Newtonian fluid

The steady Von Kármán flow and heat transfer of an electrically conducting non-Newtonian fluid is extended to the case where the disk surface admits partial slip. The fluid is subjected to an external uniform magnetic field perpendicular to the plane of the disk. The constitutive equation of the non-Newtonian fluid is modeled by that for a Reiner–Rivlin fluid. The momentum equations give rise to highly non-linear boundary value problem. Numerical solutions for the governing non-linear equations are obtained over the entire range of the physical parameters. The effects of slip, magnetic parameter and non-Newtonian fluid characteristics on the velocity and temperature fields are discussed in detail and shown graphically. Emphasis has been laid to study the effects of viscous dissipation and Joule heating on the thermal boundary layer. It is interesting to find that the non-Newtonian cross-viscous parameter has an opposite effect to that of the slip and the magnetic parameter on the velocity and the temperature fields.     

Numerical Studies of Viscoelastic Flow Using the Software OpenFOAM

Viscoelastic fluids are a special class of non-Newtonian fluids. There are several types of viscoelastic fluid models, and all of them have a complex rheological response in comparison to Newtonian fluids. This response can be viewed as a combination of viscous and elastic effects and non-linear phenomena. This complex physics makes a numerical simulation a rather challenging task, even in simple test-cases. Studies presented in this paper are numerical studies of the viscoelastic fluid flow in several test cases. These studies have been done in OpenFOAM, an open-source CFD package. Implementation of viscoelastic models and a solver is only available in a community driven version of software (OpenFOAM-ext). One of the goals of research in this paper was to test the solver and models on some simple test cases. We considered start-up and pulsating flows of viscoelastic fluid in a channel and a circular pipe. The important thing is that an analytical solution can be found in these cases, making in possible to test all aspects of numerical simulation in OpenFOAM. Obtained results showed an excellent agreement with the analytical solution for both velocity and stress components. These results encouraged authors' motivation and a choice to use OpenFOAM for simulation of viscoelastic flows. We hope that our research will make a contribution to the OpenFOAM community. Our plan for the further research is a simulation of blood flow in arteries with the viscoelastic solver. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

微极液体在两个多孔圆盘间的MHD流动及其热传导

两个平行的无限大多孔圆盘,圆盘表面有均匀注入时,数值地研究圆盘间不可压缩导电微极流体,在横向外加磁场作用下的轴对称稳定层流.运用von Krmn的相似变换,将非线性运动的控制方程转化为无量纲形式.使用基于有限差分格式的算法,在相应的边界条件下,求解简化后耦合的常微分方程组.讨论Reynolds数、磁场参数、微极参数和Prandtl数,对流动速度和温度分布的影响.在特殊情况下,所得结果与已有文献的工作有着很好的一致性.研究表明,圆盘表面的传热率随着Rynolds数、磁场参数和Prandtl数的增加而增加;剪切应力随着注入的增加而减少,但它随着外部磁场的加强而增加.和Newton流体相比较,微极流体的剪切应力因素较弱,有利于聚合体加工过程中流动和温度的控制.     

Cattaneo–Christov heat flux model for three-dimensional magnetohydrodynamic flow of an Eyring Powell fluid over an exponentially stretching surface with convective boundary condition

The present model concentrates on three-dimensional steady incompressible flow of an Eyring-Powell nanofluid past an exponentially stretching sheet with magnetic field. The Cattaneo–Christov heat flux with convective boundary condition is accounted for. Shooting method is the instrumental for obtaining numerical solution of the transformed-converted system of the mathematical models. Behavior of the determining thermo-physical parameters on the velocity, temperature, skin friction, heat transfer rate, and finally isotherms are considered. The major relevant outcomes of the current investigation are that increment in Eyring-Powell parameter uplifts flow velocity, while that peters out the fluid temperature. Enhanced values of the mixed convection parameter weakened the skin friction coefficient while it slightly strengthened the rate of heat transfer.     

Biomagnetic viscoelastic fluid flow over a stretching sheet

Of concern in this paper is an investigation of biomagnetic flow of a non-Newtonian viscoelastic fluid over a stretching sheet under the influence of an applied magnetic field generated owing to the presence of a magnetic dipole. The viscoelasticity of the fluid is characterised by Walter’s B fluid model. The applied magnetic field has been considered to be sufficiently strong to saturate the ferrofluid. The magnetization of the fluid is considered to vary linearly with temperature as well as the magnetic field intensity. The theoretical treatment of the physical problem consists of reducing it to solving a system of non-linear coupled differential equations that involve six parameters, which are solved by developing a finite difference technique. The velocity profile, the skin-friction, the wall pressure and the rate of heat transfer at the sheet are computed for a specific situation. The study shows that the fluid velocity increases as the rate of heat transfer decreases, while the local skin-friction and the wall pressure increase as the magnetic field strength is increased. It is also revealed that fluid viscoelasticity has an enhancing effect on the local skin-friction. The study will have an important bearing on magnetic drug targeting and separation of red cells as well as on the control of blood flow during surgery.     

Momentum and heat transfer over a continuously moving surface with a parallel free stream in a viscoelastic fluid

The flow and heat transfer characteristics for a continuous moving surface in a viscoelastic fluid are investigated. Constitutive equations of viscoelastic fluid obey the second‐grade model. Analytic expressions to velocity and temperature have been developed by employing homotopy analysis method. The criterion to the convergence of the solution is properly discussed. Furthermore, the values of skin friction coefficient and the local Nusselt number have been computed and discussed. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Falling film flow of a viscoelastic fluid along a wall

In this paper, we study the heat transfer in the fully developed flow of a viscoelastic fluid, a slag layer, down a vertical wall. A new constitutive relation for the stress tensor of this fluid is proposed, where the viscosity depends on the volume fraction, temperature, and shear rate. For the heat flux vector, we assume the Fourier's law of conduction with a constant thermal conductivity. The model is also capable of exhibiting normal stress effects. The governing equations are non‐dimensionalized and numerically solved to study the effects of various dimensionless parameters on the velocity, temperature, and volume fraction. The effect of the exponent in the Reynolds viscosity model is also discussed. The different cases of shear‐thinning and shear‐thickening, cooling and heating, are compared and discussed. The results indicate that the viscous dissipation and radiation (at the free surface) cause the temperature to be higher inside the flow domain. Copyright © 2013 John Wiley & Sons, Ltd.     

Nonlinear evolution of the travelling waves at the surface of a thin viscoelastic falling film

The problem of hydrodynamic instability of a thin condensate viscoelastic liquid film flowing down on the outer surface of an axially moving vertical cylinder is investigated. In order to improve the accuracy of numerical results, the viscoelastic and heat transfer parameters have been included into the governing equations. Also, the analytical solutions are obtained by utilizing the long-wave perturbation method. The influence of some physical parameters is discussed in both linear and nonlinear steps of the problem. It has been revealed that the stability of the film flow is weakened when the radius of cylinder and the temperature difference are reduced. Moreover, it is found that the increment of down-moving motion of the cylinder can enhance the flow stability. Further, the thin film flow can be destabilized by the viscoelastic property. The results show that both supercritical stability and subcritical instability can take place within the film flow system given appropriate conditions. Moreover, the absence of Reynolds number leads to an obvious difference in the behavior of some physical parameters.