Numerical simulation of Taylor Couette flow of Bingham fluids

The Bingham fluid flow between two concentric cylinders is studied using numerical simulation. The cylinders are assumed to rotate independently, and with an imposed axial sliding. The flow field is decomposed with linearity arguments of the base circular Couette shear flow and corresponding deviation field. The numerical methods are based on the expression of the deviation field in terms of complete sets of orthogonal functions and Chebyshev series. The Galerkin projection method is used with the pressure term being eliminated. The Adams Bashforth scheme is adopted for time marching. The results show that the vortices are squeezed toward the inner cylinder due to the effect of yield stress. When the outer cylinder is held stationary, the yield stress plays a role in weakening the vortex flow. However, for the co-rotation situation, the vortex flow is initially strengthened with an increase of yield stress, and then weakened as the yield stress is raised large enough. The annular unyielded regions emerge and stick to the outer cylinder. In case of Taylor Couette flow with an imposed axial sliding, a spiral vortex flow is visible with spiral unyielded region being obtained.     

正方形排列四圆柱体绕流特性及互扰效应研究

基于计算流体动力学理论,运用大涡模拟方法对雷诺数Re=3900三维正方形排列四圆柱体结构群的绕流问题进行数值计算,主要分析来流攻角与间距比两个参数对四圆柱体结构群流体参数及流场模态的影响。结果表明:来流攻角与间距比均对四圆柱体结构群绕流特性有较强的影响;来流攻角θ=0°、22.5°、45°下,临界间距比分别为3.5、4.0、3.0;间距比的变化会导致下游圆柱表面压力系数分布发生改变;另一方面,间距比较小时,四圆柱体结构之间的互扰作用均以临近效应为主;随间距比增大,上游圆柱尾流对下游圆柱有显著影响,其互扰作用会转变尾激效应。     

Hysteresis phenomenon at interaction of shock waves generated by a cylinder array

This paper describes investigations of the interaction between bow shock waves generated by cylindrical bodies in a supersonic flow. Numerical simulations are performed using the inviscid Euler equations for cylinders whose axes are parallel to each other and normal to the flow direction. Mostly an infinite periodical cylinder array is considered, but the case of two cylinders is also briefly discussed. Three different regimes of the shock wave interaction, a regular interaction, a Mach interaction, and a choked flow, have been observed for the flow through the periodical cylinder array. In the case of the flow around two bodies, the choked flow is replaced by a regime with a collective bow shock. The transition between different flow regimes is studied by varying the inflow Mach number or the distance between the cylinders. A hysteretic behavior at the transition between the regular and Mach interactions has been observed. The transition is governed by the theoretical detachment and von Neumann criteria based on the local shock wave inclination at the interaction point.     

Hysteresis phenomenon in supersonic flow past a system of cylinders

The interaction between the bow shocks ahead of a system of bodies in a supersonic flow is numerically investigated. Flow past cylinders with parallel axes lying in the plane perpendicular to the flow direction is calculated. Three different shock interaction patterns are obtained in modeling flow past an infinite periodic lattice, namely, the regular and Mach-type regimes and a regime with flow choking. In the case of two bodies a collective bow shock is formed. Transition between the flow patterns with variation in the freestream Mach number and the spacing between the cylinders is studied. Regular-to-Mach-configuration transition and vice versa occurs at particular angles of inclination of the bow shock at the interaction point determined by the criterion of a maximum flow deflection angle and the von Neumann criterion, respectively.     

An experimental study of drop migration in shear flow between concentric cylinders

The phenomenon of migration of liquid drops in Couette flow between concentric cylinders due to non-Newtonian fluid properties and shape deformation has been studied experimentally. The results agree very well with the theory of Chan and Leal, which included the effect of hydrodynamic interaction with the bounding walls, and that of velocity profile curvature in a Couette device. Significant observations that were not reported in previous studies include the migration of a deformable Newtonian drop to an equilibrium position between the centerline and the inner rotor, and the competition between normal stresses and shape deformation effects for the case of a Newtonian drop in a non-Newtonian fluid.     

Instability of Bingham fluids in Taylor-Dean flow between two concentric cylinders at arbitrary gap spacings

Stability of Bingham fluids is investigated numerically in azimuthal pressure-driven flow between two infinitely long concentric cylinders. An infinitesimal perturbation is introduced to the basic flow and its time evolution is monitored using normal mode linear stability analysis. An eigenvalue problem is obtained which is solved numerically using pseudo-spectral collocation method. Numerical results are obtained for two different cases: (i) the inner cylinder is rotating at constant velocity while the outer cylinder is fixed (i.e., the Taylor-Dean flow) and (ii) both cylinders are fixed (i.e., the Dean flow). The results show that the yield stress always has a stabilizing effect on the Taylor-Dean flow. But, for the Dean flow the effect of the yield stress is predicted to be stabilizing or destabilizing depending on the magnitude of the Bingham number and also the gap size.     

A two-dimensional model of the flow of ordered suspensions of rods

Analysis of creeping flow over an array of freely-rotating cylinders sandwiched between two sliding parallel plates is studied using a finite-difference and a least-squares numerical technique. The flow pattern was found to be very much influenced by the cylinder-to-cylinder spacing and by the gap width of the parallel plates. The shear stress on the cylinder surface and on the parallel plates was found to be a strong function of position. The viscosity of a suspension composed of an array of freely-rotating cylinders was deduced from the applied shear rate and the evaluated shear stress on the parallel plates. Experimental results confirm the numerical findings.     

Drift due to two obstacles in different arrangements

We study the drift induced by the passage of two cylinders through an unbounded extent of inviscid incompressible fluid under the assumption that the flow is two dimensional and steady in the moving frame of reference. The goal is to assess how the resulting total particle drift depends on the parameters of the geometric configuration, namely the distance between the cylinders and their angle with respect to the direction of translation. This problem is studied by numerically computing, for different cylinder configurations, the trajectories of particles starting at various initial locations. The velocity field used in these computations is expressed in closed form using methods of the complex function theory, and the accuracy of calculations is carefully verified. We identify cylinder configurations which result in increased and decreased drift with respect to the reference case when the two cylinders are separated by an infinite distance. Particle trajectories shed additional light on the hydrodynamic interactions between the cylinders in configurations resulting in different drift values. This ensemble of results provides insights about the accuracy of models used to study biogenic transport.     

Flow of a viscoelastic fluid between eccentric cylinders: impact on flow stability

The flow of an elastic fluid between eccentric cylinders was experimentally and theoretically examined. A purely elastic flow instability in the eccentric cylinder flow geometry has been discovered 1, 2. The focus of this paper is an investigation of how the characteristics of the base flow ultimately influence the flow stability. The velocity profiles of a Boger fluid between eccentric cylinders was measured by Laser Doppler velocimetry and compared to the theoretical predictions from a lubrication analysis. The base flow stress profiles between eccentric cylinders were calculated by integrating the polymeric stresses along streamlines. It is shown that the convection of polymeric stresses alter the hoop stresses in the flow. Implications of this observation to the flow stability are discussed.     

Taylor–Dean instability of yield-stress fluids at large gaps

Centrifugal instability of Bingham fluids is investigated in Taylor–Dean flow when the gap size is large compared to the cylinders radii. To determine conditions for the onset of instability, an infinitesimal axisymmetric disturbance was introduced to the basic flow and its evolution in time was monitored using a normal-mode linear stability analysis. To avoid the problem with the stress discontinuity at the location of the yield surface(s), use was made of the Papanastasiou’s regularized variation of the Bingham model in order to obtain the basic flow velocity profiles. An eigenvalue problem was obtained for the exact Bingham fluid which was solved numerically using the collocation method. A plot of the neutral instability curve at different Bingham numbers suggests that the yield stress can have a stabilizing or destabilizing effect on Taylor–Dean flow depending on the sign and magnitude of the pressure gradient, and also on the sense of rotation of the two cylinders with respect to each other.     

On the flow behaviour of confined finite-length wavy cylinders

Confined aspect-ratio of 6 wavy cylinders with a mean blockage-ratio of 0.5 were studied using time-resolved particle-image velocimetry at a sub-critical Reynolds number of 2700. Wavelengths and wave amplitudes of 2–4 and 0.1–0.3 mean diameters respectively were investigated. Results show that vortices are generally shed from the wavy cylinder and channel walls regularly, reminiscent of the unsteady symmetric flow configuration in confined non-wavy cylinders. Furthermore, vortex formation lengths for confined wavy cylinders are generally shorter than their unconfined counterparts, though their variations with respect to geometrical changes remain consistent with unconfined flow conditions. Gross cross-stream flow behaviour does not differ significantly between confined and unconfined wavy cylinders, indicating that finite-length effects are independent of the present confinement. Confined wavy cylinder wake regions are more sensitive towards geometrical changes and a combination of small wavelength and large wave amplitude leads to significant suppression of coherent cylinder and wall vortex-shedding. This is supported by phase-averaged flow reconstructions derived from Proper Orthogonal Decomposition analysis. Lastly, larger wave amplitudes lead to redistributions of dominant flow energy further downstream and to higher mode numbers, which suggests a causal link to the formation of stronger and more coherent streamwise vortices.     

Numerical simulation of viscous oscillatory flow past four cylinders in square arrangement

The results of a numerical study of the viscous oscillating flow past four circular cylinders, for a constant frequency parameter equal to 50 and KC ranging between 0.2 and 10, are presented. The cylinders were placed on the vertices of a square, two sides of which were perpendicular and two parallel to the oncoming flow, for pitch ratios, P/D, ranging between 2 and 5. The finite-element method was employed for the solution of the Navier-Stokes equations, in the formulation where the stream function and the vorticity are the field variables, whereas the pressure distribution throughout the computational domain was obtained from the solution of Poisson’s equation. When the Keulegan-Carpenter number is lower than 4, the flow remains symmetrical with respect to the horizontal axis of symmetry of the solution domain and periodic at consecutive cycles. As KC increases to 4, the flow becomes aperiodic in different cycles, although symmetry with respect to the horizontal central line of the domain is preserved. For KC equal to 5, asymmetries appear intermittently in the flow, which are eventually amplified as KC increases still further. These asymmetries, in association with the aperiodicity of flow in different cycles, lead to an almost chaotic configuration as KC grows larger. For characteristic cases the flow pattern and the time histories of the in-line and transverse forces exerted on the cylinders are presented. The mean transverse forces acting on the cylinders, the r.m.s. values of the in-line and transverse forces, together with the drag and inertia coefficients of the in-line force, were evaluated for each pitch ratio in the range of Keulegan-Carpenter numbers examined and are presented in diagrams.     

基于中间构形的大变形弹塑性模型

在大变形弹塑性本构理论中,一个基本的问题是弹性变形和塑性变形的分解.通常采用两种分解方式,一是将变形率(或应变率)加法分解为弹性和塑性两部分,其中,弹性变形率与Kirchhoff应力的客观率通过弹性张量联系起来构成所谓的次弹性模型,而塑性变形率与Kirchhoff应力使用流动法则建立联系;另一种是基于中间构形将变形梯度进行乘法分解,它假定通过虚拟的卸载过程得到一个无应力的中间构形,建立所谓超弹性-塑性模型.研究了基于变形梯度乘法分解并且基于中间构形的大变形弹塑性模型所具有的若干性质,包括:在不同的构形上,塑性旋率的存在性、背应力的对称性、塑性变形率与屈服面的正交性以及它们之间的关系.首先,使用张量函数表示理论,建立了各向同性函数的若干特殊性质,并导出了张量的张量值函数在中间构形到当前构形之间进行前推后拉的简单关系式.然后,基于这些特殊性质和关系式,从热力学定律出发,建立模型在不同构形上的数学表达,包括客观率表示的率形式和连续切向刚度等,从而获得模型所具有的若干性质.最后,将模型与4种其他模型进行了比较分析.      

Numerical and experimental study of two turbulent opposed plane jets

The turbulent interaction between two opposed plane jets separated by a distance H is experimentally studied by using a PIV (Particle Image Velocimetry) method and numerically investigated by means of a finite volume code. Two turbulence models have been tested: the standard k-ɛ model and a second-order model. The validation of the numerical study was performed by comparing the results with experimental data obtained for the case of two interacting opposed jets at ambient temperature (isothermal case). The effect of the angle of inclination of the jets is studied. Conclusions of the validation are then used to study the interaction between two jets, one being maintained at ambient temperature whereas the other is heated. Results show that the stagnation point moves towards the heated jet. It is shown that the heating induces a stabilizing effect on the flow.     

Perturbation theory for viscoelastic fluids between eccentric rotating cylinders

The circumferential and radial profiles of velocity, pressure and stress are derived for the flow of model viscoelastic liquids between two slightly eccentric cylinders with the inner one rotating. Singular perturbation methods are used to derive expansions valid for small gaps between the cylinders, but for all Deborah numbers. Results for Newtonian, second-order, Criminale-Ericksen-Filbey, upper-convected Maxwell, and White-Metzner constitutive equation separate the effects of elasticity, memory, and shear thinning on the development of the large stress gradients that hinder numerical solutions with these models in more complicated geometries. The effect of the constitutive equation on the critical Deborah number for flow separation is presented.     

A numerical study on the flow patterns of two oscillating cylinders

Flow around two oscillating cylinders in a side-by-side arrangement at Reynolds number (Re)=185 is simulated using the immersed boundary method. The purpose of this study is to investigate the combined effects of the gap between the two cylinders and their oscillation in the flow. The cylinders oscillate transversely to a uniform cross-flow with a prescribed sinusoidal function in the opposite direction, with the oscillation amplitude equal to 20% of the cylinder diameter. The gap between the two cylinders and the oscillating frequency are chosen as major variables for the parametric study to investigate their influence on the flow pattern. The ratio of mean gap distance between the two oscillating cylinders to the cylinder diameter is chosen to be 0.6, 1.0, 1.4, and 1.8, and the ratio of oscillating frequencies to the natural vortex shedding frequency of a fixed cylinder is 0.8, 1.0, and 1.2. Wake patterns and the drag and lift coefficients are described and compared with those from a single oscillating cylinder and two stationary cylinders. The wake patterns of two oscillating cylinders can be explained by flow mechanisms of two stationary cylinders, a single oscillating cylinder, and their combinations, and are in agreement with classifications of flow over two stationary cylinders presented in previous studies. In the case of two oscillating cylinders, the modulation phenomenon appears from a lower excitation frequency than in a single oscillating cylinder. Generally, oscillating cylinders have higher drag and root-mean-square (r.m.s.) values of drag coefficients than stationary cylinders.     

Weissenberg effect and its dependence upon the experimental geometry

The flow of a second-order fluid with a free surface between two coaxially mounted cylinders of finite length, the inner one of which is rotating, is being studied. In the case of slow flow and small shear rates the flow can be divided into a primary flow in the plane perpendicular to the axis of rotation and a secondary flow in the meridional plane. These flow components are numerically calculated and the results are compared with the analytical results for the semi-infinite cylinder approximation. The influence of the finiteness of the cylinders (end effect) upon the free surface deformation is analysed. The numerical results for the secondary flow are compared with results obtained by flow visualisation.     

Three-dimensional numerical simulations of cross-flow around four cylinders in an in-line square configuration

This paper presents a numerical study of three-dimensional (3-D) laminar flow around four circular cylinders in an in-line square configuration. The investigation focuses on effects of spacing ratio (L/D) and aspect ratio (H/D) on 3-D flow characteristics, and the force and pressure coefficients of the cylinders. Extensive 3-D numerical simulations were performed at Reynolds number of 200 for L/D from 1.6 to 5.0 at H/D=16 and H/D from 6 to 20 at L/D=3.5. The results show that the 3-D numerical simulations have remedied the inadequacy of 2-D simulations and the results are in excellent agreement with the experimental results. The relation between 3-D flow patterns and pressure characteristics around the four cylinders is examined and discussed. The critical spacing ratio for flow pattern transformation was found to be L/D=3.5 for H/D=16, while a bistable wake pattern was observed at L/D=1.6 for the same aspect ratio. Moreover, a transformation of flow pattern from a stable shielding flow pattern to a vortex shedding flow pattern near the middle spanwise positions of the cylinders was observed and was found to be dependent on the aspect ratio, spacing ratio, and end wall conditions. Due to the highly 3-D nature of the flows, different flow patterns coexist over different spanwise positions of the cylinders even for the same aspect ratio. It is concluded that spacing ratio, aspect ratio, and the no-slip end wall condition have important combined effects on free shear layer development of the cylinders and hence have significant effects on the pressure field and force characteristics of the four cylinders with different spacing ratios and aspect ratios.     

Effect of yield stress on the flow of a Casson fluid in a homogeneous porous medium bounded by a circular tube

The effect of yield stress on the flow characteristics of a Casson fluid in a homogeneous porous medium bounded by a circular tube is investigated by employing the Brinkman model to account for the Darcy resistance offered by the porous medium. The non-linear coupled implicit system of differential equations governing the flow is first transformed into suitable integral equations and are solved numerically. Analytical solution is obtained for a Newtonian fluid in the case of constant permeability, and the numerical solution is verified with that of the analytic solution. The effect of yield stress of the fluid and permeability of the porous medium on shear stress and velocity distributions, plug flow radius and flow rate are examined. The minimum pressure gradient required to start the flow is found to be independent of the permeability of the porous medium and is equal to the yield stress of the fluid.     

Marangoni effect in the Couette flow

Viscous heating in an axisymmetric creeping flow of a second-order fluid with free surface between two coaxially mounted cylinders produces a radial temperature gradient in the fluid. The dependence of the surface tension upon temperature is the cause for a secondary flow in the meridional plane of the flow field. This secondary flow (Marangoni effect), and its influence upon the shape of the free surface are studied. The deformation of the free surface caused by the Marangoni effect is compared with the deformation caused by inertia and normal stress differences.