MHD free-convective flow of micropolar and Newtonian fluids through porous medium in a vertical channel

We have studied the fully-developed free-convective flow of an electrically conducting fluid in a vertical channel occupied by porous medium under the influence of transverse magnetic field. The internal prefecture of the channel is divided into two regions; one region filled with micropolar fluid and the other region with a Newtonian fluid or both the regions filled by Newtonian fluids. Analytical solutions of the governing equations of fluid flow are found to be in excellent agreement with analytical prediction. Analytical results for the details of the velocity, micro-rotation velocity and temperature fields are shown through graphs for various values of physical parameters. It is noticed that Newtonian fluids prop up the linear velocity of the fluid in contrast to micropolar fluid. Also the skin friction coefficient at both the walls is derived and its numerical values are offered through tables.     

Two fluid flow and heat transfer in an inclined channel containing porous and fluid layer

Convective flow and heat transfer in an inclined channel bounded by two rigid plates held at constant different temperatures with one region filled with porous matrix saturated with a viscous fluid and another region with a clear viscous fluid different from the fluid in first region is studied analytically. The coupled nonlinear governing equations are solved using regular perturbation method. It is found that the presence of porous matrix in one of the region reduces the velocity and temperature. Results have been presented for a wide range of governing parameters such as Grashof number, porous parameter, angle of inclination, ratio of heights of the two layers and also the ratio of viscosities.     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

纳米尺度通道内的界面流体粘性研究

含液微纳米孔隙在自然界中普遍存在,在发展 日趋精密化、微型化的工业中也有着广泛的应用,深刻理解流体在微纳米通道内的物性变化对于相关自然现象以及工业应用具有重要的指导意义.本文基于分子动力学方法,建立了由金属铂板构成的二维纳米尺度通道分子模型,分别考察了受限Lennard-Jones流体和水的物性变化.根据密度、剪切应力...     

Forced convection — axial conduction between parallel walls with unequal heat fluxes

Laminar flow heat transfer is computed for a situation in which a fluid moves along a parallel plate channel with unequal wall heat fluxes (one wall is insulated). The fluid enters the heating section through an upstream region which is perfectly insulated. This situation serves to describe an upper bound for the commonly encountered case of double pipe heat exchangers with identical thermal conditions. A control volume approach has been employed for the numerical work enabling a fast calculation for the thermally developing regime in the parallel plate channel. The merits of the adopted procedure are assessed by comparison with other results available in the literature for the one-region and for the two-region problem.     

Flow in rotating radial channels at small Rossby and Ekman numbers

The problem of flow of a viscous fluid in a rotating channel is considered in the region of very small Rossby and Ekman numbers and moderately large Reynolds numbers. Asymptotic expressions with respect to the Ekman number are found for the velocity components and the longitudinal pressure gradient by solving a system of linear differential equations using Fourier series. The stability limits of such flow are predicted. Attention is drawn to a similarity between the velocity profiles of these flows and flows of a magnetic fluid and a fluid executing longitudinal oscillations in a fixed channel.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 11–15, January–February, 1984.     

On the orientation of ellipsoidal particles in a turbulent shear flow

Direct numerical simulation (DNS) of small prolate ellipsoidal particles suspended in a turbulent channel flow is reported. The coupling between the fluid and the particles is one-way. The particles are subjected to the hydrodynamic drag force and torque valid for creeping flow conditions. Six different particle cases with varying particle aspect ratios and equivalent response times are investigated. Results show that, in the near-wall region, ellipsoidal particles tend to align with the mean flow direction, and the alignment increases with increasing particle aspect ratio. When the particle inertia increases, the particles are less oriented in the spanwise direction and more oriented in the wall-normal direction. In the core region of the channel, the orientation becomes isotropic.     

Development of flow in the entrance region of a converging channel

An approximate solution for the laminar flow of an incompressible viscous fluid in the entrance region of a converging channel is obtained. The radial velocity distribution at the entrance of the channel is taken to be a symmetric but otherwise arbitrary function of the angular position. Expressions for the velocity components and pressure are given. The case of the uniformly distributed entrance velocity is considered as an example.     

VIBRATION AND STABILITY OF VERTICAL UPWARD-FLUID-CONVEYING PIPE IMMERSED IN RIGID CYLINDRICAL CHANNEL

A theoretical model is developed for the vibration and stability of a vertical pipe subjected concurrently to two dependent axial flows. The external fluid, after exiting the outer annular region between the pipe and a rigid cylindrical channel, is conveyed upwards inside the pipe. This configuration thus resembles of a pipe that aspirating fluid. The equation of planar mo- tion is solved by means of the differential quadrature method （DQM）. Calculations are conducted for a slender drill-string-like and a bench-top-size system, for different confinement conditions of the outer annular channel. It is shown that the vibrations of these two systems are closely related to the degree of confinement of the outer annular channel. For a drill-string-like system with narrow annuli, buckling instability may occur in the second and third modes. For a bench-top-size system, however, both buckling and flutter may occur in the lowest three modes. The form of instability depends on the annuli size.     

Hyperbolicity and change of type in the flow of viscoelastic fluids through channels

We consider steady flow of an upper convected Maxwell fluid through a channel with wavy walls. The vorticity of this flow will change type when the velocity in the center of the channel is larger than a critical value defined by the propagation of shear waves. There is then a central region of the channel in which the vorticity equation is hyperbolic and a low speed region near the walls where the vorticity equation is elliptic. We linearize the problem for small amplitude waviness and the linearized problem is solved in detail. The characterstic nets depend on the viscoelastic “Mach” number which is the ratio (M = U/c) of the unperturbed maximum velocity U to the speed of shear waves c into the fluid at rest and the elasticity number E. There is a supercritical (hyperbolic) region around the center of the channel when M > 1. When M ? 1, the thickness of this hyperbolic region is small when E is large, and large when E is small. Regions of positive and negative vorticity are swept out along forward facing characteristics in the hyperbolic region. There is rapid damping of vorticity in the hyperbolic region away from the boundary when M ? 1 and the Weissenberg number

. (The Weissenberg number is proportional to the relaxation time of the fluid.)The rate of damping of vorticity decreases as W is increased. Flows with high M appear to be more “elastic” when W is large in the sense that the damping is suppressed as the relaxation time of the fluid is increased.     

Simulations of magnetorheological suspensions in Poiseuille flow

Particle-level simulations are conducted to study magnetorheological fluids in plane Poiseuille flow. The importance of the boundary conditions for the particles at the channel walls is examined by considering two extreme cases: no friction and infinite coefficient of friction. The inclusion of friction produces Bingham fluid behavior, as commonly observed experimentally for MR suspensions. Lamellar structures, similar to those reported for electrorheological fluids in shear flow, are observed in the post-yield region for both particle boundary conditions. The formation of these lamellae is accompanied by an increase in the bulk fluid velocity. The slip boundary condition produces higher fluid velocities and thicker lamellar structures.     

Dynamics of a long tubular cantilever conveying fluid downwards,which then flows upwards around the cantilever as a confined annular flow

A theoretical model is developed for the dynamics of a hanging tubular cantilever conveying fluid downwards; the fluid, after exiting from the free end, is pushed upwards in the outer annular region contained by the cantilever and a rigid cylindrical channel. This configuration thus resembles that of a drill-string with a floating fluid-powered drill-bit. The linear equation of motion is solved by means of a hybrid Galerkin–Fourier method, as well as by a conventional Galerkin method. Calculations are conducted for a very slender system with parameters appropriate for a drill-string, for different degrees of confinement of the outer annular channel; and also for another, bench-top-size experiment. For wide annuli, the dynamics is dominated by the internal flow and, for low flow velocities, the flow increases the damping associated with the presence of the annular fluid. For narrow annuli, however, the annular flow is dominant, tending to destabilize the system, giving rise to flutter at remarkably low flow velocities. The mechanisms underlying the dynamics are also considered, in terms of energy transfer from the fluid to the cantilever and vice versa, as are possible applications of this work.     

Motion of a circular cylinder in a viscous liquid between parallel plates

The two-dimensional motion of a cylinder in a viscous fluid between two parallel walls of a vertical channel is studied. It is found that when the cylinder moves very closely along one of the channel walls, it always rotates in the direction opposite to that of contact rolling along the nearest wall. When the cylinder is away from the walls, its rotation depends on the Reynolds number of the flow. In this study two numerical methods were used. One is for the unsteady motion of a sedimenting cylinder initially released from a position close to one of the channel walls, where the Navier-Stokes equations are solved for the fluid and Newton's equations of motion are solved for the rigid cylinder. The other method is for the steady flow in which a cylinder is fixed in a uniform flow field where the channel walls are sliding past the cylinder at the speed of the approaching flow, or equivalently a cylinder is moving with a constant velocity in a quiescent fluid. The flow field, the drag, the side force (lift), and the torque experienced by the cylinder are studied in detail. The effects of the cylinder location in the channel, the size of the channel relative to the cylinder diameter, and the Reynolds number of the flow are examined. In the limit when the cylinder is translating very closely along one of the walls, the flow in the gap between the cylinder and the wall is solved analytically using lubrication theory, and the numerical solution in the other region is used to piece together the whole flow field.This research was supported by NSF DMR91-20668 through the Laboratory for Research on the Structure of Matter at the University of Pennsylvania and from the Research Foundation of the University of Pennsylvania.     

Finite Element Analysis of Convection Flow Through a Porous Medium in a Horizontal Channel

We analyse the convection flow of a viscous fluid through a horizontal channel enclosing a fully saturated porous medium. The Galerkin finite element analysis is used to discuss the flow and heat transfer through the porous medium using serendipity elements. The velocity, the temperature distributions and the rate of heat transfer are analysed for variations in the governing parameters. The profiles at different vertical levels are asymmetric curves, exhibiting reversal flow everywhere except on the midplane. In a given porous medium, for fixed G or N, the temperature in the fluid region at any position in fluids with a higher Prandtl number, is much higher than in fluids with a lower Prandtl number. Likewise, other parameters being fixed, lesser the permeability of the medium, lower the temperature in the flow field. Nu reduces across the flow at all axial positions, while it enhances along the axial direction of the channel. Nu reduces with decrease in the Darcy parameter D, and thus lesser the permeability of the medium, lesser the rate of heat transfer across the boundary at any axial position of the channel.     

Direct numerical simulation of particle-laden plane turbulent wall jet and the influence of Stokes number

The distribution and motion of inertial particles in plane turbulent wall jet are investigated using direct numerical simulation, under the assumption of one-way coupling. To our knowledge, this appears to be the first direct numerical simulation of a particle-laden plane turbulent wall jet. It is shown that, in outer part of the wall jet, the behaviour of particles closely resembles that of a free plane jet. Due to the streamwise decay of particle Stokes number, the particle streaks formed in the near wall region of the wall jet are characterized by their intensity variation, which differs significantly from those in the channel flow. The streamwise growth of the particle velocity half-width is approximately equal to that of the fluid velocity half-width and the maximum velocity of particles decays slower than that of fluid due to inertia. The outer scaling can collapse the mean particle velocity in both the inner and outer region for heavier particles. In the buffer region, the particle–fluid velocity difference can be negative or positive depending on the Stokes number since there are two competing effects, namely the memory effect and turbophoresis. In the viscous region, the larger particles are on average faster than fluid and the velocity difference is found to be self-similar depending on outer Stokes number. The near-wall distribution of velocity difference is significantly correlated with the presence of high-momentum particles which are entrained by vortical structures generated in the outer region of the wall jet. These results are useful for environmental and engineering applications.     

Viscoplastic flow in slightly varying channels with wall slip pertaining to a magnetorheological (MR) polishing process

In this paper, the first analytical endeavor into the fluid dynamic modeling of an MR polishing process is reported. The velocity and shear stress fields of an MR fluid running through a thin slippery channel with a slightly varying height are analytically solved using a bi-viscosity constitutive model and a Navier slip model. Estimations of the mechanical power density and the total power per unit depth applied onto the channel surfaces are also presented. Analytical solutions for the Couette–Poiseuille flow behavior of a bi-viscous fluid flowing through either parallel or non-uniform channels are obtained, and the associated necessary and sufficient conditions characterizing a total of 5 types of flow are derived.The behaviors of the fluid are examined through the use of a parametric diagram of Bingham number (Bn) and Couette number (Co), i.e., Bn–Co or 1/Bn–1/Co diagram, by changing the geometric and operating conditions. Using these diagrams, variations in the rheological characteristics of the flow are investigated in great detail, with a special focus on the movement of the pseudo-core region. Finally, the mechanical power density field obtained for the flow in a converging–diverging channel is used to explain the wear mechanism in the MR polishing process. The effects of the power density field and the total power on the material removal rate (MRR) and the within-workpiece nonuniformity (WIWNU) with respect to various geometric and operating conditions are evaluated.     

A Direct theory of viscous fluid flow in channels

This paper deals with an Eulerian formulation of the theory of directed fluid sheets appropriate for incompressible, linear viscous fluid flow in channels with arbitrary shapes for their major boundaries which may be moving or fixed. Special cases of the theory are applied to a number of two-dimensional fluid flow problems and these solutions are in general discussed for unsteady flow. Specific applications include fluid flow in a channel whose boundaries are symmetric with respect to a middle plane in the channel, subjected to time-dependent pressure gradient at one end; and to lubrication problems in a general shaped channel when one of the channel walls is a fixed plane while the other is moving with a constant velocity. Flow of a viscous fluid with a free surface over a fixed boundary is also discussed.Dedicated to J. L. Ericksen on the occasion of his Sixtieth Birthday     

Mathematical model of micropolar fluid in two-phase immiscible fluid flow through porous channel

This paper is concerned with the flow of two immiscible fluids through a porous horizontal channel. The fluid in the upper region is the micropolar fluid/the Eringen fluid, and the fluid in the lower region is the Newtonian viscous fluid. The flow is driven by a constant pressure gradient. The presence of micropolar fluids introduces additional rotational parameters. Also, the porous material considered in both regions has two different permeabilities. A direct method is used to obtain the analytical solution of the concerned problem. In the present problem, the effects of the couple stress, the micropolarity parameter, the viscosity ratio, and the permeability on the velocity profile and the microrotational velocity are discussed. It is found that all the physical parameters play an important role in controlling the translational velocity profile and the microrotational velocity. In addition, numerical values of the different flow parameters are computed. The effects of the different flow parameters on the flow rate and the wall shear stress are also discussed graphically.     

Streamline upwind finite element method for conjugate heat transfer problems

This paper presents a combined finite element method for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow. The streamline upwind finite element method is used for the analysis of thermal viscous flow in the fluid region, whereas the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the proposed method is to consistently couple heat transfer along the fluid-solid interface. Three test cases, i.e. conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the efficiency of the present method. The English text was polished byYunming Chen.