The structure of parallel layers in steady two-dimensional magnetohydrodynamic flows in sudden duct expansions and contractions

The structure of non-linear, steady, two-dimensional parallel layers at high values of the Hartmann number, Ha, the Reynolds number, Re, and the interaction parameter, N, for duct expansions and contractions has been investigated. The magnetic field is transverse to the flow. For the flow regime, in which the electromagnetic force balances the inertial force in the layer, a viscous sublayer at the solid wall parallel to the magnetic field has been obtained, in which the flow is driven by the pressure gradient induced in the outer, inviscid layer.     

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.      

竖置管流中液固两相脉动特性和颗粒浓度分布

利用激光多普勒分相测量技术,考察了液固两相自下而上通过竖置矩形管时,固、液两相的时均速度、流向及横向的脉动强度和颗粒相的相对浓度分布,证实了颗粒浓度的横向分布主要取决于颗粒的横向脉动强度分布（即npvp′2^-=常数）的分析结果。     

Periodical streaming potential and electro-viscous effects in microchannel flow

This paper presents an analytical solution to periodical streaming potential, flow-induced electric field and velocity of periodical pressure-driven flows in twodimensional uniform microchannel based on the Poisson-Boltzmann equations for electric double layer and Navier-Stokes equation for liquid flow. Dimensional analysis indicates that electric-viscous force depends on three factors： （1） Electric-viscous number representing a ratio between maximum of electric-viscous force and pressure gradient in a steady state, （2） profile function describing the distribution profile of electro-viscous force in channel section, and （3） coupling coefficient reflecting behavior of arnplitude damping and phase offset of electro-viscous force. Analytical results indicate that flow-induced electric field and flow velocity depend on frequency Reynolds number （Re = wh^2/v）. Flow-induced electric field varies very slowly with Re when Re 〈 1, and rapidly decreases when Re 〉 1. Electro-viscous effect on flow-induced electric field and flow velocity are very significant when the rate of the channel width to the thickness of electric double layer is small.     

Periodical streaming potential and electro-viscous effects in microchannel flow

This paper presents an analytical solution to periodical streaming potential,flow-induced electric field and velocity of periodical pressure-driven flows in twodimensional uniform microchannel based on the Poisson.Boltzmann equations for electric double layer and Navier-Stokes equation for liquid flow.Dimensional analysis indicates that electric-Viscous force depends on three factors:(1)Electric-viscous number representing a ratio between maximum of electric-viscous force and pressure gradient in a steady state,(2)profile function describing the distribution profile of electro-viscous forcein channel section,and(3)coupling coefficient reflecting behavior of amplitude damping and phase Offset of electro-viscous force.Analytical results indicate that flow-induced electric field and flow velocity depend on frequency Reynolds number(Re=wh2/v).Flow-induced electric field varies very slowly with Re when Re＜1.and rapidly decreases when Re＞1.Electro-viscous effect on flow-induced electric field and flow velocity are very significant when the rate of the channel width to the thickness of electric double layer is small.     

On the effect of flow direction on mixed convection from a horizontal cylinder

The influence of free stream direction on mixed (natural and forced) convective heat transfer from a circular cylinder is investigated. The cylinder, which has an isothermal surface, is placed with its axis horizontal and normal to the oncoming flow. The free stream direction varies between the vertically upward (parallel flow) and the vertically downward (contraflow) directions. The investigation is based on the time integration of the unsteady, two-dimensional equations of motion and energy until reaching steady conditions. The study is limited to Reynolds numbers up to Re = 40 and Grashoff numbers of Gr = Re². The results are compared with the available experimental data and the agreement is satisfactory.     

Peristaltic transport and heat transfer of a MHD Newtonian fluid with variable viscosity

The influence of temperature‐dependent viscosity and magnetic field on the peristaltic flow of an incompressible, viscous Newtonian fluid is investigated. The governing equations are derived under the assumptions of long wavelength approximation. A regular perturbation expansion method is used to obtain the analytical solutions for the velocity and temperature fields. The expressions for the pressure rise, friction force and the relation between the flow rate and pressure gradient are obtain. In addition to analytical solutions, numerical results are also computed and compared with the analytical results with good agreement. The results are plotted for different values of variable viscosity parameter β, Hartmann number M, and amplitude ratio ?. It is found that the pressure rise decreases as the viscosity parameter β increases and it increases as the Hartmann number M increases. Finally, the maximum pressure rise (σ=0) increases as M increases and β decreases. Copyright © 2009 John Wiley & Sons, Ltd.     

Mixed convection flow over a vertical wedge embedded in a highly porous medium

 The steady mixed convection flow over a vertical wedge with a magnetic field embedded in a porous medium has been investigated. The effects of the permeability of the medium, surface mass transfer and viscous dissipation on the flow and temperature fields have been included in the analysis. The coupled nonlinear partial differential equations governing the flow field have been solved numerically using the Keller box method. The skin friction and heat transfer are found to increase with the parameters characterizing the permeability of the medium, buoyancy force, magnetic field and pressure gradient. However the effect of the permeability and magnetic field on the heat transfer is very small. The heat transfer increases with the Prandtl number, but the skin friction decreases. The buoyancy force which assists the forced convection flow causes an overshoot in the velocity profiles. Both the skin friction and heat transfer increase with suction and the effect of injection is just the reverse. Received on 21 May 1999     

A numerical study of pulsatile laminar flows in a pipe with a ring-type constriction

Numerical simulations have been carried out to study pulsatile laminar flows in a pipe with an axisymmetric ringtype constriction. Three types of pulsatile flows were investigated, namely a physiological flow, a pure sinusoidal flow and a non-zero mean velocity sinusoidal flow. The laminar flow governing equations were solved by the SIMPLE algorithm on a non-staggered grid and a modified Crank-Nicolson approximation was used to discretrize the momentum equations with respect to time. The maximum flow Reynolds numer (Re) is 100. The Womersley number (N_w) ranges from 0 to 50, with the corresponding Strouhal number (St) ranging from 0 to 3·98. The constriction opening ratio (d/D) and thickness ratio (h/D) are fixed at 0·5 and 0·1 respectively. Within the time period investigated, all these pulsatile flows include both forward and backward flows. The unsteady recirculation region and the recirculation points change in size and location with time. For N_w ≤ 1 and St≤ 1·56 x 10^?3 the three pulsatile flows have the same simple relation between the instantaneous flow rate and pressure loss (Δp) across the constriction and the pressure gradient in the axial direction (dp/dz) in the fully developed flow region. The phase angles between the flow rate and pressure loss and the pressure gradient are equal to zero. With increasing N_w and St, the phase angle between the flow rate and the dp/dz becomes larger and has its maximum value of 90° at N_w = 50 and St = 3·98. The three pulsatile flows also show different relations between the flow rate and the pressure gradient. The pure sinusoidal flow has the largest maximum pressure gradient and the non-zero mean velocity sinusoidal flow has the smallest. For larger N_w and St the fully developed velocity profiles in the fully developed flow region have a smaller velocity gradient along the radial direction in the central region. The maximum recirculation length increases for N_w ranging from 0 to 4·2, while this length becomes very small at N_w = 50 and St = 3·98. The deceleration tends to enlarge the recirculation region and this effect appears for N_w ≥ 3 and St ≥ 1·43×10^?2. Linear relations exist between the flow rate and the instantaneous maximum values of velocity, vorticity and shear stress.     

Viscous and Inviscid Instabilities of Flow Along a Streamwise Corner

Here we consider the stability of flow along a streamwise corner formed by the intersection of two large flat plates held perpendicular to each other. Self-similar solutions for the steady laminar mean flow in the corner region have been obtained by solving the boundary layer equations for zero and nonzero streamwise pressure gradients. The stability of the mean flow is investigated using linear stability analysis. An eigensolver has been developed to solve the resulting linear eigenvalue problem either in a global mode to obtain an approximation to all the dominant eigenmodes or in a local mode to refine a particular eigenmode. The stability results indicate that the entire spectrum of two-dimensional and oblique viscous modes of a two-dimensional Blasius boundary layer is active in the case of a corner layer as well, but away from the cornerline. In a corner region of finite spanwise extent, the continuous spectrum of oblique modes degenerates to a discrete spectrum of modes of increasing spanwise wave number. The effect of the corner on the two-dimensional viscous instability is small and decreases the growth rate. The growth rate of outgoing oblique disturbances is observed to decrease, while the growth rate of incoming oblique disturbances is enhanced by the corner. This asymmetry between the outgoing and incoming viscous modes increases with increasing obliqueness of the disturbance. The instability of a zero pressure gradient corner layer is dominated by the viscous modes; however, an inviscid corner mode is also observed. The critical Reynolds number of the inviscid mode rapidly decreases with even a small adverse streamwise pressure gradient and the inviscid mode becomes the dominant one. Received 17 March 1998 and accepted 28 April 1999     

Hall effects on the magnetohydrodynamic shear flow past an infinite porous flat plate subjected to uniform suction or blowing

An analysis is made of Hall effects on the steady shear flow of a viscous incompressible electrically conducting fluid past an infinite porous plate in the presence of a uniform transverse magnetic field. It is shown that for suction at the plate, steady shear flow solution exists only when S²<Q, where S and Q are the suction and magnetic parameters, respectively. The primary flow velocity decreases with increase in Hall parameter m. But the cross-flow velocity first increases and then decreases with increase in m. Similar results are obtained for variation of the induced magnetic field with m. It is further found that for blowing at the plate, steady shear flow solution exists only when , where S₁ is the blowing parameter.     

Local properties of vertical mercury-argon two-phase flow in a circular tube under transverse magnetic field

Experimental data are presented in this paper on the profiles of local void fraction, bubble impaction rate, bubble velocity and its spectrum, and also bubble length and its spectrum, of mercury-argon two-phase slug flow flowing upwards in a vertical circular tube in the presence of a transverse magnetic field. Decrease in void fraction and increase in bubble velocity are significant when the magnetic flux density is larger than $\text{0.3～0.4}$$\text{T}\text{(Ha ? 100)}$. This effect is discussed by analyzing the bubble size distribution. Recovery of local void fraction profile in the downstream of an obstacle and diffusion of void injected from only one nozzle in the presence of magnetic field are also discussed.     

Laminar flow through channels with porous walls and with an applied transverse magnetic field

Summary The problem of two-dimensional steady laminar flow of a viscous incompressible and electrically conducting fluid through a channel with two equally porous walls in the presence of a transverse magnetic field has been extended to include all values of Hartmann number and small suction velocity at the walls. Expressions for the velocity components, the pressure and the wall friction in terms of the Hartmann number and the suction Reynolds number are given. It is found that the pressure drop in the major flow direction and the wall friction decrease with the increase in suction and increase with the increase in the strength of the magnetic field.     

Determination of gradients of flow variables behind three-dimensional stationary and pseudo stationary curved shock waves in conducting gases

Summary In this paper we have obtained the gradients of magnetic field, velocity, pressure and density behind a shock wave in three dimensional steady motion of a conducting gas. For the shock configuration, we take a continuous differentiable function of coordinates and it is assumed that the components of the magnetic field H ⁱ , velocity components u ⁱ , pressure p and density behind the shock-surface are differentiable functions. Moreover we take H ⁱ , u ⁱ , p and in front of the shock-wave as constant quantities. In § 4 we have obtained the gradients of flow and field quantities behind the pseudostationary shockwave. § 5 is devoted to the calculation of gradients of flow and field quantities in cases where the normal component of the magnetic field is zero on both sides of the shock wave. In § 6 the relation between the curvature k of the shock-surface and the curvature K of the stream line just behind the shock surface in two dimensional steady motion has been derived. § 7 deals with the determination of the ratio K/k for an attached shock in the case of a wedge.     

Pulsatile magneto-biofluid flow and mass transfer in a non-Darcian porous medium channel

A numerical study of pulsatile flow and mass transfer of an electrically conducting Newtonian biofluid via a channel containing porous medium is considered. The conservation equations are transformed and solved under boundary conditions prescribed at both walls of the channel, using a finite element method with two-noded line elements. The influence of magnetic field on the flow is studied using the dimensionless hydromagnetic number, Nm, which defines the ratio of magnetic (Lorentz) retarding force to the viscous hydrodynamic force. A Darcian linear impedance for low Reynolds numbers is incorporated in the transformed momentum equation and a second order drag force term for inertial (Forchheimer) effects. Velocity and concentration profiles across the channel width are plotted for various values of the Reynolds number (Re), Darcy parameter (λ), Forchheimer parameter (Nf), hydro-magnetic number (Nm), Schmidt number (Sc) and also with dimensionless time (T). Profiles of velocity varying in space and time are also provided. The conduit considered is rigid with a pulsatile pressure applied via an appropriate pressure gradient term. Increasing the hydromagnetic number (Nm) from 1 to 15 considerably depresses biofluid velocity (U) indicating that a magnetic field can be used as a flow control mechanism in, for example, medical applications. A rise in Nf from 1 to 20 strongly retards the flow development and decreases the velocity, U, across the width of the channel. The effects of other parameters on the flowfield are also discussed at length. The flow model also has applications in the analysis of electrically conducting haemotological fluids flowing through filtration media, diffusion of drug species in pharmaceutical hydromechanics, and also in general fluid dynamics of pulsatile systems.     

Application of the Dorodnitsyn finite element method to swirling boundary layer flow

The Dorodnitsyn finite element method for turbulent boundary layer flow with surface mass transfer is extended to include axisymmetric swirling internal boundary layer flow. Turbulence effects are represented by the two-layer eddy viscosity model of Cebeci and Smith¹ with extensions to allow for the effect of swirl. The method is applied to duct entry flow and a 10 degree included-angle conical diffuser, and produces results in close agreement with experimental measurements with only 11 grid points across the boundary layer. The introduction of swirl (w_e/u_e = 0.4) is found to have little effect on the axial skin friction in either a slightly favourable or adverse pressure gradient, but does cause an increase in the displacement area for an adverse pressure gradient. Surface mass transfer (blowing or suction) causes a substantial reduction (blowing) in axial skin friction and an increase in the displacement area. Both suction and the adverse pressure gradient have little influence on the circumferential velocity and shear stress components. Consequently in an adverse pressure gradient the flow direction adjacent to the wall is expected to approach the circumferential direction at some downstream location.     

Forced Convection with Laminar Pulsating Counterflow in a Saturated Porous Circular Tube

An analytical solution is obtained for forced convection in a circular tube occupied by a core–sheath-layered saturated porous medium with counterflow produced by pulsating pressure gradients. The case of the constant heat-flux boundary conditions is considered, and the Brinkman model is employed for the porous medium. A perturbation approach is used to obtain analytical expressions for the velocity, temperature distribution, and transient Nusselt number for convection produced by an applied pressure gradient that fluctuates with small amplitude harmonically in time about a non-zero mean. It is shown that the fluctuating part of the Nusselt number alters in magnitude and phase as the dimensionless frequency increases. The magnitude increases from zero, goes through a peak, and then decreases to zero. The height of the peak depends on the values of various parameters. The phase (relative to that of the steady component) decreases from π/2 to − π/2 as the frequency increases.     

Steady mixed convection stagnation-point flow of upper convected Maxwell fluids with magnetic field

The steady MHD mixed convection flow of a viscoelastic fluid in the vicinity of two-dimensional stagnation point with magnetic field has been investigated under the assumption that the fluid obeys the upper-convected Maxwell (UCM) model. Boundary layer theory is used to simplify the equations of motion, induced magnetic field and energy which results in three coupled non-linear ordinary differential equations which are well-posed. These equations have been solved by using finite difference method. The results indicate the reduction in the surface velocity gradient, surface heat transfer and displacement thickness with the increase in the elasticity number. These trends are opposite to those reported in the literature for a second-grade fluid. The surface velocity gradient and heat transfer are enhanced by the magnetic and buoyancy parameters. The surface heat transfer increases with the Prandtl number, but the surface velocity gradient decreases.     

局部狭窄血管中血液振荡流的速度和切应力

本文建立一种分析局部缓慢狭窄血管中血液振荡流的数学模型，给出了血液的轴向流速，径向流速和切应力的包含压力梯度项的解析表达式，并讨论了血管内由局部狭窄引起的压力梯度沿轴向变化的规律。文章以局部余弦狭窄为例进行数值计算，详细讨论上游均匀管段压力梯度的定常部分和不同次谐波对狭窄管段内流速和切应力的影响。数值结果表明，与均匀管情况相比，在狭窄段内，血液振荡流轴向流速无论平均值还是脉动幅值均明显增大，且径向流速不再为零。但径向流速仍远小于轴向流速。同时，切应力也不再仅由轴向流速梯度提供，径向流速梯度也将产生切应力，但是在计算管壁切向上的切应力时，径向流速梯度的贡献仍相当大。与均匀管管壁切应力沿流运方向保持恒定不同。狭窄管管壁切应力（平均值和脉动值）将随着狭窄高度的增大而增大，在狭窄最大高度处达到最大，因而沿流动方向产生了较大的切应力梯度。