Slow visco-elastic flow past a cylinder in a rectangular channel

This paper is concerned with the flow of a visco-elastic liquid through a rectangular channel containing a cylindrical obstruction placed either in a symmetric or asymmetric position with respect to the centre of the channel. Numerical predictions of the flow are obtained using a well established finite element Galerkin mixed formulation. The influence of elasticity on the streamline pattern is found to be negligible, and one only observes changes due to different geometries, i.e. relative positions of cylinder and channel. However, both elasticity and a variable viscosity are found to have significant effects on the forces exerted on the cylinder.     

A new constitutive model for fibre suspensions: flow past a sphere

A new phenomenological constitutive equation for homogeneous suspensions of macrosized fibres is proposed. In the model, the local averaged orientation of the fibres is represented by a director field, which evolves in time in a manner similar to the rotation of a prolate spheroid. The stress is linear in the strain rate, but the viscosity is a fourth-order tensor that is directly related to the director field. In the limit of low-volume fractions of fibres, the model reduces properly to the leading terms of the constitutive equation for dilute suspensions of spheroids. The model has three parameters: the aspect ratio R of the fibres, the volume fraction , and A, which plays the role of the maximum-volume fraction of the fibres. Experimental shear data are used to estimate the parameter A, and the resulting model is used in a boundary-element program to study the flow past a sphere placed at the centre line of a cylinder for the whole range of volume fractions from 0.01 to near maximum volume fraction. The agreement with experimental data from Milliken et al. [1] is good.     

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.     

Conjugate heat transfer in the unbounded flow of a viscoelastic fluid past a sphere

This work addresses the conjugate heat transfer of a simplified PTT fluid flowing past an unbounded sphere in the Stokes regime (Re = 0.01). The problem is numerically solved with the finite-volume method assuming axisymmetry, absence of natural convection and constant physical properties. The sphere generates heat at a constant and uniform rate, and the analysis is conducted for a range of Deborah (0 ≤ De ≤ 100), Prandtl (10⁰ ≤ Pr ≤ 10⁵) and Brinkman (0 ≤ Br ≤ 100) numbers, in the presence or absence of thermal contact resistance at the solid–fluid interface and for different conductivity ratios (0.1 ≤ κ ≤ 10). The drag coefficient shows a monotonic decrease with De, whereas the normalized stresses on the sphere surface and in the wake first increase and then decrease with De. A negative wake was observed for the two solvent viscosity ratios tested (β = 0.1 and 0.5), being more intense for the more elastic fluid. In the absence of viscous dissipation, the average Nusselt number starts to decrease with De after an initial increase. Heat transfer enhancement relative to an equivalent Newtonian fluid was observed for the whole range of conditions tested. The dimensionless temperature of the sphere decreases and becomes more homogeneous when its thermal conductivity increases in relation to the conductivity of the fluid, although small changes are observed in the Nusselt number. The thermal contact resistance at the interface increases the average temperature of the sphere, without affecting significantly the shape of the temperature profiles inside the sphere. When viscous dissipation is considered, significant changes are observed in the heat transfer process as Br increases. Overall, a simplified PTT fluid can moderately enhance heat transfer compared to a Newtonian fluid, but increasing De does not necessarily improve heat exchange.     

A power-law fluid flow past a porous sphere

Summary A model has been developed for the flow of a non-Newtonian fluid past a porous sphere. The drag force exerted on a porous sphere moving in a power-law fluid is obtained by an approximate solution of equations of motion in the creeping flow regime. It is predicted that the effect of the pseudoplastic anomaly on the drag force is more pronounced at large porosity parameters.

---

Zusammenfassung Es wird ein Modell für die Strömung einer nichtnewtonschen Flüssigkeit längs einer porösen Kugel entwickelt. Die auf die in einer Ostwald-DeWaele-Flüssigkeit bewegte Kugel ausgeübte Reibungskraft wird durch eine Näherungslösung der Bewegungsgleichungen für schleichende Strömung gewonnen. Man findet, daß der Einfluß der Abweichung vom newtonschen Verhalten um so ausgeprägter wird, je größer die Porosität ist.

---

A, B, C, D a, b, c, d coefficients in eqs. [10] and [18] - F _D drag force - K consistency index in power-law model - k ₁ ,k ₂ coefficients defined by eq. [18] - m porosity parameter - n flow index in power-law model - P pressure - P ^* dimensionless pressure defined by eq. [4] - P pressure difference - R radius of porous sphere - r radial distance from the center of the sphere - U velocity of uniform stream - u _i velocity component - u _i ^* dimensionless velocity component defined by eq. [4] - Y drag force correction factor defined by eq. [27] - _ij rate of deformation tensor - _ij ^* dimensionless rate of deformation tensor defined by eq. [4] - , spherical coordinates - dimensionless radial distance defined by eq. [4] - second invariant of rate of deformation tensor - ^* dimensionless second invariant of rate of deformation tensor defined by eq. [4] - _ij stress tensor - _ij ^* dimensionless stress tensor defined by eq. [4] - stream function - ^* dimensionless stream function defined by eq. [4] - i inside the surface of the sphere - o outside the surface of the sphere With 1 figure and 1 table     

O(2)‐symmetry breaking bifurcation: with application to the flow past a sphere in a pipe

The steady, axisymmetric laminar flow of a Newtonian fluid past a centrally‐located sphere in a pipe first loses stability with increasing flow rate at a steady O(2)‐symmetry breaking bifurcation point. Using group theoretic results, a number of authors have suggested techniques for locating singularities in branches of solutions that are invariant with respect to the symmetries of an arbitrary group. These arguments are presented for the O(2)‐symmetry encountered here and their implementation for O(2)‐symmetric problems is discussed. In particular, how a bifurcation point may first be detected and then accurately located using an ‘extended system’ is described. Also shown is how to decide numerically if the bifurcating branch is subcritical or supercritical. The numerical solutions were obtained using the finite element code ENTWIFE. This has enabled the computation of the symmetry breaking bifurcation point for a range of sphere‐to‐pipe diameter ratios. A wire along the centerline of the pipe downstream of the sphere is also introduced, and its effect on the critical Reynolds number is shown to be small. Copyright © 2000 John Wiley & Sons, Ltd.     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.     

Numerical investigation of flow structures around a sphere

A numerical investigation of flow around a sphere is performed and compared with previous studies. Here, a second-order accurate, finite volume method is used in order to predict the instantaneous and time-averaged flow characteristics using large eddy simulation (LES) on the multi-block grid system. Namely, the objectives of this article are: (i) the presentation of flow structures in the wake region downstream of the sphere with a wide variety of flow properties such as the distribution of velocity vectors, patterns of streamlines, Reynolds stress correlations, root mean square of velocity components and other time-averaged flow data in order to reveal the vortical flow structures in detail and (ii) to demonstrate the abilities of computational methods in simulation of vortical flow data. Finally, it has been concluded that there are good agreements between the experimental results and numerical predictions.     

The flow past a sphere in a cylindrical tube: effects of intertia,shear-thinning and elasticity

Numerical solutions are presented for the flow past a sphere placed at the centreline of a cylindrical tube for Reynolds numbers ranging from 0 to 150, using a boundary element method. Fluids are modelled by a variety of constitutive equations including the Newtonian, the Carreau and the Phan-Thien-Tanner models. The influence of inertia, shear-thinning and fluid elasticity on the flow field, drag and the pressure drop force-drag ratio is examined. Some results are compared with available experimental data.     

圆柱管道内运动的小球受到的黏滞阻力的数值分析

为了更好地运用落球法测量研究流体的黏滞系数,研究小球在黏性流体中下落的受力情况,本文对小球在充满黏性流体的圆柱管道的下落过程进行分析。利用COMSOL4.4仿真模拟,建立了合理的仿真模型,并分析了小球受到的黏滞阻力与小球的大小、下落位置的关系。结果表明:选择速度项二阶近似、压强项一阶近似的离散化方法,可以得到和理论值非常相符的仿真结果;当下落过程中小球球心始终在圆柱轴线上时,小球受到的黏滞阻力相对于Stokes力的修正系数,是小球半径与圆柱管道半径的比例函数,本文得到了更大范围的符合理论解的修正系数;当下落过程中小球的球心偏离圆柱轴线时,对于同样大小的小球,黏滞阻力、压强力、黏性力均随着球心到轴线的距离先减小后增大,且具有不同的极小值点。     

Rolling resistance of a rigid sphere with viscoelastic coatings

We present a novel three-dimensional boundary-element formulation that fully characterizes the mechanical behavior of the external boundary of a multi-layered viscoelastic coating attached to a hard rotating spherical core. The proposed formulation incorporates both, the viscoelastic, and the inertial effects of the steady-state rolling motion of the sphere, including the Coriolis effect. The proposed formulation is based on Fourier-domain expressions of all mechanical governing equations. It relates two-dimensional Fourier series expansions of surface displacements and stresses, which results in the formation of a compliance matrix for the outer boundary of the deformable coating, discretized into nodes. The computational cost of building such a compliance matrix is optimized, based on configurational similarities and symmetry. The proposed formulation is applied, in combination with a rolling contact solving strategy, to evaluate the viscoelastic rolling friction of a coated sphere on a rigid plane. Steady-state results generated by the proposed model are verified by comparison to those obtained from running dynamic simulations on a three-dimensional finite element model, beyond the transient. A detailed application example includes a verification of convergence and illustrates the dependence of rolling resistance on the applied load, the thickness of the coating, and the rolling velocity.     

并列双椭圆柱绕流的格子Boltzmann-虚拟区域方法的模拟研究

采用格子Boltzmann-虚拟区域方法对并列双椭圆柱绕流进行了模拟研究。首先,通过与并列双圆柱的结果进行对比,验证了数值方法的有效性。其次对雷诺数为200时两种间距(g=0.5a和2.0a,g为柱体表面间距,a为椭圆长轴)的情况进行了研究,考察了椭圆长轴与短轴之比,即α值对柱体升阻力系数及涡结构的影响。研究发现,与圆柱相比,对于g=0.5a椭圆柱的升阻力系数可能出现两种变化,一是升阻力随时间演化较规则,接近周期性;二是流场可能长时间偏向于其中一个椭圆柱,这些变化与α的值有关。对于g=2.0a,两个椭圆柱后的某一区域内会出现四列涡街,经过一段时间,四列涡街又会演化成两列向两侧扩张的涡街。     

Aerodynamic coefficients of a spinning sphere in a rarefied-gas flow

A three-dimensional rarefied-gas flow past a spinning sphere in the transitional and near-continuum flow regimes is studied numerically. The rarefaction and compressibility effects on the lateral (Magnus) force and the aerodynamic torque exerted on the sphere are investigated for the first time. The coefficients of the drag force, the Magnus force, and the aerodynamic torque are found for Mach numbers ranging from 0.1 to 2 and Knudsen numbers ranging from 0.05 to 20. In the transitional regime, at a certain Knudsen number depending on the Mach number the Magnus force direction changes. This change is attributable to the increase in the role of normal stresses and the decrease in the contribution of the shear stresses to the Magnus force with decrease in the Knudsen number. A semi-empirical formula for the calculation of the Magnus force coefficient in the transitional flow regime is proposed.     

Isothermal flow past a blowing sphere

Steady, axisymmetric, isothermal, incompressible flow past a sphere with uniform blowing out of the surface is investigated for Reynolds numbers in the range 1 to 100 and surface velocities up to 10 times the free stream value. A stream-function-velocity formulation of the flow equations in spherical polar co-ordinates is used and the equations are solved by a Galerkin finite-element method. Reductions in the drag coefficients arising from blowing are computed and the effects on the viscous and pressure contributions to the drag considered. Changes in the surface pressure, surface vorticity and flow patterns for two values of the Reynolds number (1 and 40) are examined in greater detail. Particular attention is paid to the perturbation to the flow field far from the sphere.     

Uniform flow of viscoelastic fluids past a confined falling cylinder

Uniform steady flow of viscoelastic fluids past a cylinder placed between two moving parallel plates is investigated numerically with a finite-volume method. This configuration is equivalent to the steady settling of a cylinder in a viscoelastic fluid, and here, a 50% blockage ratio is considered. Five constitutive models are employed (UCM, Oldroyd-B, FENE-CR, PTT and Giesekus) to assess the effect of rheological properties on the flow kinematics and wake patterns. Simulations were carried out under creeping flow conditions, using very fine meshes, especially in the wake of the cylinder where large normal stresses are observed at high Deborah numbers. Some of the results are compared with numerical data from the literature, mainly in terms of a drag coefficient, and significant discrepancies are found, especially for the constant-viscosity constitutive models. Accurate solutions could be obtained up to maximum Deborah numbers clearly in excess of those reported in the literature, especially with the PTT and FENE-CR models. The existence or not of a negative wake is identified for each set of model parameters.     

An analytic approximation of the drag coefficient for the viscous flow past a sphere

We give an analytic solution at the 10th order of approximation for the steady-state laminar viscous flows past a sphere in a uniform stream governed by the exact, fully non-linear Navier-Stokes equations. A new kind of analytic technique, namely the homotopy analysis method, is applied, by means of which Whitehead's paradox can be easily avoided and reasonably explained. Different from all previous perturbation approximations, our analytic approximations are valid in the whole field of flow, because we use the same approximations to express the flows near and far from the sphere. Our drag coefficient formula at the 10th order of approximation agrees better with experimental data in a region of Reynolds number R_d<30, which is considerably larger than that (R_d<5) of all previous theoretical ones.     

Oscillatory flow of second grade fluid in cylindrical tube

The unsteady oscillatory flow of an incompressible second grade fluid in a cylindrical tube with large wall suction is studied analytically. Flow in the tube is due to uniform suction at the permeable walls, and the oscillations in the velocity field are due to small amplitude time harmonic pressure waves. The physical quantities of interest are the velocity field, the amplitude of oscillation, and the penetration depth of the oscillatory wave. The analytical solution of the governing boundary value problem is obtained, and the effects of second grade fluid parameters are analyzed and discussed.     

Stokes flow of a rotating sphere around the axis of a circular orifice or a circular disk

This paper presents an infinite series solution to the creeping flow equations for the axisymmetric motion of a sphere of arbitrary size rotating in a quiescent fluid around the axis of a circular orifice or a circular disk whose diameters are either larger or smaller than that of the sphere. Numerical tests of the convergence are passed and the comparison with the exact solution and other computational results shows an agreement to five significant figures for the torque coefficients in both cases. The torque coefficients are obtained for the sphere located up to a position tangent to the wall plane containing either the orifice or the disk. It is concluded that the torque coefficients of the sphere and the disk are monotonically increasing with the decrease of the distance from the disk or the orifice plane in both cases.     

Numerical studies of slow viscous rotating flow past a sphere. II

The Navier–Stokes equations, which are the governing equations for a steady, viscous, incompressible fluid rotating about the z-axis with angular velocity ω, are linearized using the Oseen approximation. Two parameters, namely the Reynolds number Re = Ua/v and Re_ω = 2ωa²/v (the Reynolds number w.r.t. rotation), enter the linearized equations. These equations are solved by the Peaceman–Rachford ADI method and the resulting algebraic equations are solved by the SOR method. Streamlines are plotted and compared with the Oseen solution for the non-rotating case. The magnitude of the vorticity vector with increasing θ is also plotted.