Boundary layer flow of Maxwell fluid due to torsional motion of cylinder: modeling and simulation

This paper investigates the boundary layer flow of the Maxwell fluid around a stretchable horizontal rotating cylinder under the influence of a transverse magnetic field. The constitutive flow equations for the current physical problem are modeled and analyzed for the first time in the literature. The torsional motion of the cylinder is considered with the constant azimuthal velocity E. The partial differential equations (PDEs) governing the torsional motion of the Maxwell fluid together with energy transport are simplified with the boundary layer concept. The current analysis is valid only for a certain range of the positive Reynolds numbers. However, for very large Reynolds numbers, the flow becomes turbulent. Thus, the governing similarity equations are simplified through suitable transformations for the analysis of the large Reynolds numbers. The numerical simulations for the flow, heat, and mass transport phenomena are carried out in view of the bvp4c scheme in MATLAB. The outcomes reveal that the velocity decays exponentially faster and reduces for higher values of the Reynolds numbers and the flow penetrates shallower into the free stream fluid. It is also noted that the phenomenon of stress relaxation, described by the Deborah number, causes to decline the flow fields and enhance the thermal and solutal energy transport during the fluid motion. The penetration depth decreases for the transport of heat and mass in the fluid with the higher Reynolds numbers. An excellent validation of the numerical results is assured through tabular data with the existing literature.     

垂直湍流液-固流中大颗粒的相对速度

通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式．用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数．实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。     

GPU accelerated simulations of three-dimensional flow of power-law fluids in a driven cube

Newtonian fluid flow in two- and three-dimensional cavities with a moving wall has been studied extensively in a number of previous works. However, relatively a fewer number of studies have considered the motion of non-Newtonian fluids such as shear thinning and shear thickening power law fluids. In this paper, we have simulated the three-dimensional, non-Newtonian flow of a power law fluid in a cubic cavity driven by shear from the top wall. We have used an in-house developed fractional step code, implemented on a Graphics Processor Unit. Three Reynolds numbers have been studied with power law index set to 0.5, 1.0 and 1.5. The flow patterns, viscosity distributions and velocity profiles are presented for Reynolds numbers of 100, 400 and 1000. All three Reynolds numbers are found to yield steady state flows. Tabulated values of velocity are given for the nine cases studied, including the Newtonian cases.     

A correlation for yield stress fluid flow through packed beds

Relatively few correlations are available for non-Newtonian fluid flows through packed beds, even though such fluids are frequently used in industry. In this paper, a correlation is presented for yield stress fluid flow through packed beds. The correlation is developed by introducing the yield stress model in place of the Newtonian model used in deriving Erguns equation. The resulting model has three parameters that are functions of the geometry and roughness of the particle surfaces. Two of the parameters can be deduced in the limit as the yield stress becomes negligible and the model reduces to Erguns equation for Newtonian fluids. The third model parameter is determined from experimental data. The correlation relates a defined friction factor to the dimensionless Reynolds and Hedstrom numbers and can be used to predict pressure drop for flow of a yield stress fluid through a packed bed of spherical particles. Conditions for flow or no-flow are also determined in the correlation. Comparison of model calculations, between a Newtonian and a yield stress fluid for flow penetration into a packed bed of spheres, shows the yield stress fluid initially performs similar to the Newtonian fluid, at large Reynolds numbers. At lower Reynolds numbers the yield stress effect becomes important and the flow rate significantly decreases when compared to the Newtonian fluid.     

Finite element simulations of incompressible flow past a heated/cooled sphere

This paper reports numerical simulation of the flow past a heated/cooled sphere. A Galerkin finite element method is used to solve the 3D incompressible Boussinesq equations in primitive variable form. Numerical simulations of flow around the sphere for a range of Grashof numbers and moderate Reynolds numbers, were conducted. The drag coefficient for adiabatic flow shows good agreement with standard correlations over the range of the Reynolds numbers investigated. It is shown that the drag can vary considerably with heating of the sphere and that computational fluid dynamics methods can be used to derive constitutive laws for macroscopic momentum and heat exchange in multiphase flow. © 1998 John Wiley & Sons, Ltd.     

Stability of Poiseuille flow in the presence of a longitudinal magnetic field

The stability of the plane flow of an electrically conducting fluid with respect to small perturbations was studied at large Reynolds numbers in the presence of a longitudinal magnetic field. The dependence of the critical Reynolds number on the electrical conductivity is investigated. At large Reynolds numbers, a new branch of instability and a sudden change in the critical Reynolds numbers is found. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 45–53, May–June, 2008.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

Some new aspects of the slow flow of a viscous fluid through an axisymmetric duct expansion or contraction. II — Experimental part

The flow of a newtonian viscous fluid through a cylindrical duct in the vicinity of a section discontinuity is studied using a visualization technique. The evolution with the Reynolds number of the features of the stationary vortex cell is given, in particular the reversibility of the flow is verified for the very small Reynolds numbers. In the creeping regime a detailed analysis of the velocity and strain-rate fields is performed.     

Diffusion flux to a distorted gas bubble at large Reynolds numbers

The diffusion flux to a distorted gas bubble situated in a uniform viscous incompressible fluid flow is determined for large Reynolds and Péclet numbers and finite Weber numbers. The bubble has the shape of an ellipsoid of revolution, oblate in the flow direction, making it possible to use the flow field derived by Moore [1] in the form of a two-term expansion with respect to the flow parameter =R^–1/2 (R is the Reynolds number; the zeroth term of the expansion corresponds to potential flow). The dependence of the diffusion flux onto the bubble surface on the Weber and Reynolds numbers is determined. The results of Winnikow [2] and Sy and Lightfoot [3] are thus generalized to the case of finite Weber numbers and a broader range of Reynolds numbers.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 70–76, July–August, 1976.     

Numerical simulation of confined laminar jet flows

Two‐dimensional incompressible jet development inside a duct has been studied in the laminar flow regime, for cases with and without entrainment of ambient fluid. Results have been obtained for the flow structure and critical Reynolds number values for steady asymmetric jet development and for the onset of temporal oscillations, at various values of the duct‐to‐jet width ratio (aspect ratio). It is found that at low aspect ratios and Reynolds numbers, jet development inside the duct is symmetric. For larger aspect ratios and Reynolds numbers, the jet flow at steady state becomes asymmetric with respect to the midplane, and for still higher values, it becomes oscillatory with respect to time. When entrainment is present, the instabilities of asymmetric development and temporal oscillations occur at a much higher critical Reynolds number for a given aspect ratio, indicating that the stability of the jet flow is higher with entrainment. Copyright © 2000 John Wiley & Sons, Ltd.     

Unsteady Laminar to Turbulent Flow in a Spacer-Filled Channel

A combined numerical and experimental investigation has been carried out to study the flow behaviour in a spacer-filled channel, representative of those used in spiral-wound membrane modules. Direct numerical simulation and particle image velocimetry were used to investigate the fluid flow characteristics inside a 2 × 2 cell at Reynolds numbers that range between 100 and 1000. It was found that the flow in this geometry moves parallel to and also rotates between the spacer filaments and that the rate of rotation increases with Reynolds number. The flow mechanisms, transition process and onset of turbulence in a spacer-filled channel are investigated including the use of the velocity spectra at different Reynolds numbers. It is found that the flow is steady for Re < 200 and oscillatory at Re ～ 250 and increasingly unsteady with further increases in Re before the onset of turbulent flow at Re ～ 1000.     

Extension of the Darcy�CForchheimer Law for Shear-Thinning Fluids and Validation via Pore-Scale Flow Simulations

Flow of non-Newtonian fluids through porous media at high Reynolds numbers is often encountered in chemical, pharmaceutical and food, as well as petroleum and groundwater engineering, and in many other industrial applications. Under the majority of operating conditions typically explored, the dependence of pressure drops on flow rate is non-linear and the development of models capable of describing accurately this dependence, in conjunction with non-trivial rheological behaviors, is of paramount importance. In this work, pore-scale single-phase flow simulations conducted on synthetic two-dimensional porous media are performed via computational fluid dynamics for both Newtonian and non-Newtonian fluids and the results are used for the extension and validation of the Darcy?CForchheimer law, herein proposed for shear-thinning fluid models of Cross, Ellis and Carreau. The inertial parameter ?? is demonstrated to be independent of the viscous properties of the fluids. The results of flow simulations show the superposition of two contributions to pressure drops: one, strictly related to the non-Newtonian properties of the fluid, dominates at low Reynolds numbers, while a quadratic one, arising at higher Reynolds numbers, is dependent on the porous medium properties. The use of pore-scale flow simulations on limited portions of the porous medium is here proposed for the determination of the macroscale-averaged parameters (permeability K, inertial coefficient ?? and shift factor ??), which are required for the estimation of pressure drops via the extended Darcy?CForchheimer law. The method can be applied for those fluids which would lead to critical conditions (high pressures for low permeability media and/or high flow rates) in laboratory tests.     

Permeable plate at angle of attack in steady flow of a viscous fluid

An approximate solution is presented for the problem of the resistance of a permeable plate of widthl at an angle of attack in a steady plane flow of an incompressible viscous fluid for the case of both small and very large Reynolds numbers with different permeability laws. The results obtained in the case of large Reynolds numbers are compared with the corresponding results for flow past plane rod grids.     

Lattice Boltzmann simulations for flow and heat/mass transfer problems in a three‐dimensional porous structure

The lattice Boltzmann method (LBM) for a binary miscible fluid mixture is applied to problems of transport phenomena in a three‐dimensional porous structure. Boundary conditions for the particle distribution function of a diffusing component are described in detail. Flow characteristics and concentration profiles of diffusing species at a pore scale in the structure are obtained at various Reynolds numbers. At high Reynolds numbers, the concentration profiles are highly affected by the flow convection and become completely different from those at low Reynolds numbers. The Sherwood numbers are calculated and compared in good agreement with available experimental data. The results indicate that the present method is useful for the investigation of transport phenomena in porous structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Two-dimensional unsteady laminar flow of a power law fluid across a square cylinder

The two-dimensional and unsteady free stream flow of power law fluids past a long square cylinder has been investigated numerically in the range of conditions 60≤Re≤160 and 0.5≤n≤2.0. Over this range of Reynolds numbers, the flow is periodic in time. A semi-explicit finite volume method has been used on a non-uniform collocated grid arrangement to solve the governing equations. The global quantities such as drag coefficients, Strouhal number and the detailed kinematic variables like stream function, vorticity and so on, have been obtained for the above range of conditions. While, over this range of Reynolds number, the flow is known to be periodic in time for Newtonian fluids, a pseudo-periodic flow regime displaying more than one dominant frequency in the lift is observed for shear-thinning fluids. This seems to occur at Reynolds numbers of 120 and 140 for n=0.5 and 0.6, respectively. Broadly speaking, the smaller the value of the power law index, lower is the Reynolds number of the onset of the pseudo-periodic regime. This work is concerned only with the fully periodic regime and, therefore, the range of Reynolds numbers studied varies with the value of the power law index. Not withstanding this aspect, in particular here, the effects of Reynolds number and of the power law index have been elucidated in the unsteady laminar flow regime. The leading edge separation in shear-thinning fluids produces an increase in drag values with the increasing Reynolds number, while shear-thickening fluid behaviour delays this separation and shows the lowering of the drag coefficient with the Reynolds number. Also, the preliminary results suggest the transition from the steady to unsteady flow conditions to occur at lower Reynolds numbers in shear-thinning fluids than that in Newtonian fluids.     

Vortex motions in stratified wake flows

The vortex wakes of obstacles (circular cylinder and sphere) moving through a linearly stratified fluid have been investigated, by means of flow visualization, at Reynolds numbers smaller than 800 and non-dimensional buoyancy frequencies smaller than 6. Vortex shedding from a horizontally suspended circular cylinder is suppressed when the fluid is stratified. The wake of a sphere is affected by lee waves when the Reynolds number exceeds about 200.     

Velocity field in a cylindrical hydrocyclone

The flow in a cylindrical hydrocyclone is investigated at moderate Reynolds numbers. Such a flow regime is realized in a hydrocyclone when, for example, certain solid particles are separated. The flow is calculated in the frame-work of the theory of an inviscid incompressible fluid, since the influence of turbulent pulsations on the flow structure in the investigated flow regime (Reynolds numbers of the order of a few thousand) is slight. Comparison of the results with experimental data indicates good qualitative agreement.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 12–20, November–December, 1980.I thank Yu. P. Gupalo and Yu. S. Ryazantsev for discussing the work and for valuable comments.     

Polymer “threads” and drag reduction

A number of different polymer fluids were ejected on the centerline of a water pipe-flow facility. Two distinct flow regions were identified: Reynolds numbers above 25000, where centerline injection acted as a rather efficient mixing device for water-soluble polymer — and no drag-reduction resulted from non water-soluble materials; and Reynolds numbers from 10000 to 25000, where strong evidence exists that under certain conditions, a viscoelastic fluid thread can interact with turbulence eddies and reduce the overall flow friction in the pipe.On Sabbatical leave from San Diego State University.