Period-doubling route to chaos for a swirling flow in an open cylindrical container with a rotating disk

Conclusion The period-doubling route to chaos for a swirling vortex flow in an open cylindrical container, at an aspect ratio of 2, driven by a rotating bottom disk was recognized by using laser-Doppler velocimetry. The onset of periodic motion for the flow was found when Re was in the range between 1850 and 1900. The flow was subharmonically bifurcated into a double-period motion when Re was about 2150. When the Reynolds number was in the range from 2300 to 2400, the flow bifurcated again through the period-doubling mode. When the Reynolds number was further increased, the flow eventually showed chaotic motion. The existence of a free surface promotes the onset of periodicity, and the difference of the critical Reynolds number with and without a free surface was estimated to be about 600.This work was supported by the National Science Council of the Republic of China under grant No. NSC-82-0410-E-002-191     

Influence of elongational properties on the contraction flow of polyisobutylene in a mixed solvent

Over the last decade several international programmes have been developed around different standard fluids, one of which is the so-called S1 fluid. This is a solution of polyisobutylene in a mixed solvent and the aim of the programme has been to study the rheology of polymer solutions from the dilute solution to the melt. The focus of this paper will be on the flow visualisation of contraction flows of S1 through orifice dies and on the estimation of some of its extensional properties. The contraction ratios range from 24.4:1 and 124.3:1. The measured entry pressure drops will be correlated with contraction ratio and apparent wall shear rate. Experimental evidence will show that, when regarded as a function of wall shear rate, the entry pressure drops are independent of the contraction ratios. The flow fields for different contraction ratios, at any constant apparent wall shear rate, however, differ substantially. The evolution of the flow fields is monitored and it is shown that an initial increase in vortex size is followed by a slower decrease, this happening at a constant Weissenberg number. At the same Weissenberg number, small scale instabilities start occurring near the office. As the shear rate is increased further, these instabilities grow in size until, eventually, the flow structure is destroyed. An extensional viscosity is evaluated using a modified form of the Binding analysis for contraction flows and we show that the results are not only in qualitative agreement with those from other groups, but also that the analysis is able to predict exactly the onset of the aforementioned flow instabilities. Received: 20 March 1997 Accepted: 18 September 1997     

Die-swell,splashing drop and a numerical technique for solving the Oldroyd B model for axisymmetric free surface flows

This work deals with the development of a numerical method for simulating viscoelastic axisymmetric free surface flow of an Oldroyd B fluid. A novel formulation is developed for the computation of the non-Newtonian extra-stress components on rigid boundaries and on the symmetry axis. The full free surface stress conditions are employed. The resulting governing equations are solved by finite differences on a Marker-and-cell (MAC) type grid. Validation is provided by simulating a pipe flow problem. The classical die-swell problem is solved and swelling ratios are provided. The height of the splash caused by a falling liquid drop for various Reynolds and Weissenberg numbers is then studied, and the height of the splash is shown to diminish with increasing viscoelasticity.     

An adaptive viscoelastic stress splitting scheme and its applications: AVSS/SI and AVSS/SUPG

We report an adaptive viscoelastic stress splitting (AVSS) scheme, which was found to be extremely robust in the numerical simulation of viscoelastic flow involving steep stress boundary layers. The scheme is different from the elastic viscous split stress (EVSS) in that the local Newtonian component is allowed to depend adaptively on the magnitude of the local elastic stress. Two decoupled versions of the scheme were implemented for the Upper Convected Maxwell (UCM) model: the streamline integration (AVSS/SI), and the streamline upwind Petrov-Galerkin (AVSS/SUPG) methods of integrating the stress. The implementations were benchmarked against the known analytic Poiseuille solution, and no upper limiting Weissenberg number was found (the computation was stopped at Weissenberg number of O(10⁴)). The flow past a sphere in a tube was solved next, giving convergent results up to a Weissenberg number of 3.2 with the AVSS/SI and 1.55 with the AVSS/SUPG (both were decoupled schemes; without the adaptive scheme, the limiting Weissenberg number for the decoupled streamline integration method was about 0.3). These results show that the decoupled AVSS is more stable than the corresponding EVSS, and the SI is more robust than SUPG in solving the constitutive equation of hyperbolic type. In addition, we found a very long stress wake behind the sphere, and a weak vortex in the rear stagnation region at a Weissenberg number above W_i of about 1.6, which does not seem to increase in size or strength with increasing W_i.     

The Unsteady Wake of a Circular Cylinder near a Free Surface

The behaviour of the wake Strouhal number for flow past a cylinder close to a free surface at both low and moderate Froude numbers is investigated numerically. For the low Froude number case (i.e., gravity-dominated), the results obtained are similar to those for flow past a cylinder close to an adjacent no-slip boundary. As the distance between the wall and the cylinder is reduced, the Strouhal number, as measured from the time varying lift, increases to a maximum at a gap ratio of 0.70. Further gap reduction leads to a rapid decrease in the Strouhal number, with shedding finally ceasing altogether at gap ratios below 0.16. The agreement between the results for a free surface and a no-slip boundary suggests that the mechanism behind the suppression of vortex shedding is common. For flow at a fixed gap ratio and a moderate Froude number, two distinctly different wake states are observed with the flow passing over the cylinder tending to switch from a state of attachment to the free surface, to one of separation from it, and then back again in a pseudo-periodic fashion. Even though there is a significant difference in Reynolds number, the predicted numerical two-dimensional behaviour is found to compare favourably with the experimental observations at higher Reynolds number.     

Flow above the free end of a surface-mounted finite-height circular cylinder: A review

The wake of a surface-mounted finite-height circular cylinder and the associated vortex patterns are strongly dependent on the cylinder aspect ratio and the thickness of the boundary layer on the ground plane relative to the dimensions of the cylinder. Above a critical aspect ratio, the mean wake is characterized by streamwise tip vortex structures and Kármán vortex shedding from the sides of the cylinder. Below a critical aspect ratio, a unique mean wake structure is observed. Recent experimental studies in the literature that used phase-averaged techniques, as well as recent numerical simulations, have led to an improved physical understanding of the near-wake vortex flow patterns. However, the flow above the free end of the finite circular cylinder, and its relationship to the near wake, has not been systematically studied. The effects of aspect ratio and boundary layer thickness on the free-end flow field are also not completely understood, nor has the influence of Reynolds number on the free-end flow field been fully explored. Common features associated with the free end include separation from the leading edge, a mean recirculation zone containing a prominent cross-stream arch (or mushroom) vortex, and reattachment onto the free-surface. Other flow features that remain to be clarified include a separation bubble near the leading edge, one or two cross-stream vortices within this separation bubble, the origins of the streamwise tip or trailing vortices, and various critical points in the near-surface flow topology. This paper reviews the current understanding of the flow above the free end of a surface-mounted finite-height circular cylinder, with a focus on models of the flow field, surface oil flow visualization studies, pressure and heat flux distributions on the free-end surface, measurements of the local velocity field, and numerical simulations, found in the literature.     

Drop oscillations under simple shear in a highly viscoelastic matrix

Recent two-dimensional numerical simulations and experiments have shown that, when a drop undergoes shear in a viscoelastic matrix liquid, the deformation can undergo an overshoot. I implement a volume-of-fluid algorithm with a paraboloid reconstruction of the interface for the calculation of the surface tension force for three-dimensional direct numerical simulations for a Newtonian drop in an Oldroyd-B liquid near criticalities. Weissenberg numbers up to 1 at viscosity ratio 1 and retardation parameter 0.5 are examined. Critical capillary numbers rise with the Weissenberg number. Just below criticality, drop deformation begins to undergo an overshoot when the Weissenberg number is sufficiently high. The overshoot becomes more pronounced, and at higher matrix Weissenberg numbers, such as 0.8, drop deformation undergoes novel oscillations before settling to a stationary shape. Breakup simulations are also described.     

Self-sustained global oscillations in a jet in crossflow

A jet in crossflow with an inflow ratio of 3, based on the maximum velocity of the parabolic jet profile, is studied numerically. The jet is modeled as an inhomogeneous boundary condition at the crossflow wall. We find two fundamental frequencies, pertaining to self-sustained oscillations in the flow, using full nonlinear direct numerical simulation (DNS) as well as a modal decomposition into global linear eigenmodes and proper orthogonal decomposition (POD) modes; a high frequency which is characteristic for the shear-layer vortices and the upright vortices in the jet wake, and a low frequency which is dominant in the region downstream of the jet orifice. Both frequencies can be related to a region of reversed flow downstream of the jet orifice. This region is observed to oscillate predominantly in the wall-normal direction with the high frequency, and in the spanwise direction with the low frequency. Moreover, the steady-state solution of the governing Navier?CStokes equations clearly shows the horseshoe vortices and the corresponding wall vortices further downstream, and the emergence of a distinct counter-rotating vortex pair high in the free stream. It is thus found that neither the inclusion of the jet pipe nor unsteadiness is necessary to generate the characteristic counter-rotating vortex pair.     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermocapillary convection in an annular cylindrical container

The liquid–gas interface of a liquid in an annular container is subjected to a temperature gradient. Since shear stress on the free liquid surface depends on the temperature it shall create a variable shear stress on the surface which in turn yields by viscous traction a thermocapillary convection in the bulk of the liquid. For constant temperature T ₂ > T ₁ on the outer cylindrical wall and T ₁ on the inner wall a steady convection shall emerge. The morphology of the ensuing flow is investigated as a function of the diameter ratio and the liquid height ratio and exhibits the growth of bottom rim vortex rings, which unite to single vortex rings of opposite flow direction to the original single strong vortex ring above it. The results also show how such a system depends on the magnitude of the diameter ratio and the liquid height ratio. The evaluation of the results exhibits the evolution of vortex rings.     

Simultaneous measurements of velocity and deformation in flows through compliant diaphragms

Flow through a circular orifice in a deformable diaphragm mounted in a pipe was studied experimentally as a simple yet suitable case for validating numerical fluid/structure interaction (FSI) codes including structures with significant deformation and strain. The flow was characterized using pressure taps, particle image velocimetry (PIV), and hot-film anemometry while deformation of the compliant diaphragm was determined directly from PIV images. The diaphragm material properties were measured independently by a uniaxial tensile testing machine. The diaphragm material modulus, orifice diameter, and pipe Reynolds number were varied over ranges appropriate for simulations of flows through heart valves. Pipe Reynolds numbers ranged from 600 (laminar upstream condition) to 8800 (turbulent upstream condition). The pressure drop across the diaphragm resulted in a concave deformation for all cases studied. For the range of Reynolds number tested, the Euler number decreased with increasing Reynolds number as a result of orifice expansion. The flow immediately downstream of compliant diaphragms was jet-like with strong inward radial velocity components and vena contracta. Laminar low Reynolds number flow (Re=600) through both rigid and compliant diaphragms yielded early and regular roll up of coherent vortex rings at a fixed frequency in contrast to turbulent higher Reynolds number flow (Re=3900), which yielded a broad range of vortex passage frequencies. Expansion of the compliant orifice for Re=3900 resulted in an initially broader slower jet with delayed shear layer development compared with the equivalent rigid case.     

Observations of vortex breakdown in an open cylindrical container with a rotating bottom

A laser induced fluorescent dye technique was used to visualize the steady-state flow driven by a rotating bottom in an open, cylindrical container. The flow behaviour and the vortex breakdown conditions were studied as a function of the container aspect ratio H/R and the Reynolds number Re = R ²/v. Like in the closed container configuration, previously studied by Vogel (1968) and Escudier (1984), vortex breakdown occurs in a certain parameter range (H/R, Re). However, in the free surface configuration vortex breakdown conditions as well as the forms of the breakdown bubbles differ notably from what is observed in the closed container configuration. In particular, it is found that as Re is increased, the breakdown bubbles get attached to the free surface and grow in diameter.     

Particle removal from a surface by a bounded vortex flow

A bounded vortex flow consists of an axisymmetric vortex that is confined top and bottom between two plates (the “confinement plate” and “impingement plate”, respectively) and surrounded laterally by a swirling annular slot jet. The bottom of the vortex terminates on the boundary layer along the impingement plate and the top of the vortex is drawn into a suction port positioned at the center of the confinement plate. The circumferential flow within the annular jet is important for supplying circulation to the central wall-normal vortex. This flow field is proposed as a method for mitigation of dust build-up on a surface, where the vortex–jet combination supplements the more traditional vacuum port by enhancing the surface shear stress and related particle transport rate. The paper reports on a computational study of the velocity field and particle transport by a bounded vortex flow. Fluid flow computations are performed using a finite-volume approach for an incompressible fluid and particle transport is simulated using a discrete-element method. Computations are performed for different values of two dimensionless parameters – the ratio of the plate separation distance and the average radial location of the jet inlet (the dimensionless confinement height) and the ratio of flow rate withdrawn at the suction outlet and that injected by the jet (the flow rate ratio). For small values of the flow rate ratio, the impinging jet streamlines pass down to the boundary layer along the bottom surface and then travel up the vortex core. By contrast, for large values of flow rate ratio, the annular jet is quickly entrained into the suction outlet and no wall-normal vortex is formed. Particles are observed to roll along the impingement surface in a direction determined by the fluid shear stress lines. Particles roll outward when they lie beyond a separatrix curve of the surface shear stress lines, where particles within this separatrix curve roll inward, piling up at the center of the flow field. A toroidal vortex ring forms for the small confinement height case with flow rate ratio equal to unity, which yields double separatrix curves in the shear stress lines. The inward rolling particles intermittently lift up due to collision forces and burst away from the impingement surface, eventually to become entrained into the flow out the suction port or resettling back onto the impingement surface.     

Dynamics of annular gas–liquid two-phase swirling jets

The dynamics of annular gas–liquid two-phase swirling jets have been examined by means of direct numerical simulation and proper orthogonal decomposition. An Eulerian approach with mixed-fluid treatment, combined with an adapted volume of fluid and a continuum surface force model, was used to describe the two-phase flow system. The unsteady, compressible, three-dimensional Navier–Stokes equations have been solved by using highly accurate numerical methods. Two computational cases have been performed to examine the effects of liquid-to-gas density ratio on the flow development. It was found that the higher density ratio case is more vortical with larger spatial distribution of the liquid, in agreement with linear theories. Proper orthogonal decomposition analysis revealed that more modes are of importance at the higher density ratio, indicating a more unstable flow field. In the lower density ratio case, both a central and a geometrical recirculation zone are captured while only one central recirculation zone is evident at the higher density ratio. The results also indicate the formation of a precessing vortex core at the high density ratio, indicating that the precessing vortex core development is dependent on the liquid-to-gas density ratio of the two-phase flow, apart from the swirl number alone.     

Application of a vortex method to free surface flows

Vortex methods have found wide applications in various practical problems. The use of vortex methods in free surface flow problems, however, is still very limited. This paper demonstrates a vortex method for practical computation of non-linear free surface flows produced by moving bodies. The method is a potential flow formulation which uses the exact non-linear free surface boundary condition at the exact location of the instantaneous free surface. The position of the free surface, on which vortices are distributed, is updated using a Lagrangian scheme following the fluid particles on the free surface. The vortex densities are updated by the non-linear dynamic boundary condition, derived from the Euler equations, with an iterative Lagrangian numerical scheme. The formulation is tested numerically for a submerged circular cylinder in unsteady translation. The iteration is shown to converge for all cases. The results of the unsteady simulations agree well with classical linearized solutions. The stability of the method is also discussed.     

Interaction of a transverse gas jet with supersonic flow in a dihedral

The results of an experimental investigation of supersonic flow at the Mach number M₁=3 past a transverse gas jet flowing from an orifice in the edge of a dihedral with a linear angle of 90° are presented. The Mach number of the jet was varied from 1 to 3, and the ratio of the total pressure in the jet to the free stream pressure from 90 to 760. Visualization of the flow near the faces of the dihedral revealed the existence of internal lines of flow convergence and divergence in the region of three-dimensional separated flow, which indicates the presence of complex vortex structures. The dependence of the dimensions of the separated flow zone and the characteristic pressures in it on the jet parameters is explored.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 161–165, November–December, 1994.     

A finite volume approach for unsteady viscoelastic fluid flows

A finite volume, time‐marching for solving time‐dependent viscoelastic flow in two space dimensions for Oldroyd‐B and Phan Thien–Tanner fluids, is presented. A non‐uniform staggered grid system is used. The conservation and constitutive equations are solved using the finite volume method with an upwind scheme for the viscoelastic stresses and an hybrid scheme for the velocities. To calculate the pressure field, the semi‐implicit method for the pressure linked equation revised method is used. The discretized equations are solved sequentially, using the tridiagonal matrix algorithm solver with under‐relaxation. In both, the full approximation storage multigrid algorithm is used to speed up the convergence rate. Simulations of viscoelastic flows in four‐to‐one abrupt plane contraction are carried out. We will study the behaviour at the entrance corner of the four‐to‐one planar abrupt contraction. Using this solver, we show convergence up to a Weissenberg number We of 20 for the Oldroyd‐B model. No limiting Weissenberg number is observed even though a Phan Thien–Tanner model is used. Several numerical results are presented. Smooth and stable solutions are obtained for high Weissenberg number. Copyright © 2002 John Wiley & Sons, Ltd.     

A low‐dimensional spectral approach for transient axisymmetric free‐surface flow inside thin cavities of arbitrary shape

A low‐dimensional spectral method is used to solve the transient axisymmetric free surface flow inside thin cavities of arbitrary shape. The flow field is obtained on the basis of the lubrication equations, which are expanded in terms of orthonormal functions over the cavity gap. The formulation accounts for nonlinearities stemming from inertia and front location. The work is of close relevance to the filling stage during die casting, and injection molding, or the flow inside annular (extrusion) dies. Both flows under an imposed flow rate, and an imposed pressure at the cavity entrance are examined. The influence of inertia, aspect ratio, gravity, and wall geometry on the evolution of the front, flow rate, and pressure is assessed particularly in the early stage of flow, when a temporal behavior of the ‘boundary‐layer’ type develops. The multiple‐scale method is applied to obtain an approximate solution at small Reynolds number, Re. Comparison with the exact (numerical) solution indicates a wide range of validity for the multiple‐scale approach, including the moderately small Re range. Copyright © 2002 John Wiley & Sons, Ltd.     

Simulating air entrainment and vortex dynamics in a hydraulic jump

The air entrainment characteristics of three separate Froude number hydraulic jumps are investigated numerically using an unsteady RANS, realizable k–ε turbulence model, with a Volume of Fluid treatment for the free surface. Mean velocity profiles, average void fraction, and Sauter mean diameter compare favorably with experimental data reported in literature. In all simulations, time-averaged void fraction profiles show good agreement with experimental values in the turbulent shear layer and an accurate representation of interfacial aeration at the free surface. Sauter mean diameter is well represented in the shear layer, and free surface entrainment results indicate bubble size remains relatively unchanged throughout the depth of the jump. Several different grid resolutions are tested in the simulations. Significant improvements in void fraction and bubble size comparison are seen when the diameter to grid size ratio of the largest bubbles in the shear layer surpasses eight. A three-dimensional simulation is carried out for one Froude number jump, showing an improvement in the prediction of entrained air and bubble size compared with two-dimensional results at a substantial increase in computation time. An analysis of three-dimensional vorticity shows a complex interaction between spanwise and streamwise vortical structures and entrained air bubbles. The jump is similar to a turbulent mixing layer, constrained by the free surface, with vortex pairing and subsequent fluctuations in free surface elevation. Downstream fluctuations of the toe are associated with a roll up of the primary spanwise vortex, fluctuations of the free surface, and counter-rotating streamwise vortex pairs. The action of these flow structures is likely responsible for the improvement in three-dimensional results.     

Non-linear dynamics of semi-dilute polydisperse polymer solutions in microfluidics: effects of flow geometry

The non-linear dynamics of a semi-dilute (c/c* = 15) polydisperse polyethylene oxide (PEO) solution in microfluidics are studied experimentally using benchmark contraction–expansion flow geometries with three contraction–expansion ratios (4:1:4, 8:1:8 and 16:1:16) and two narrow channel lengths (L _c/D _h = 53 and 5.3, where L _c is the length of the narrow channel and D _h is its hydraulic diameter). Complex flows over a range of elasticity numbers (El), Weissenberg numbers (Wi) and Reynolds numbers (Re) are characterized using micro-particle image velocimetry ( \upmu\upmu-PIV) and pressure drop measurements. The evolution of vortex formation and dynamics has been visualized through a step-flow-rate experiment. Various flow dynamics regimes have been quantified and are presented in a Wi–Re diagram. The experimental results reveal that the contraction ratio can result in qualitatively different vortex dynamics of semi-dilute polymer solutions in microfluidics, whereas the length of the narrow channel merely affects the dynamics at a quantitative level. A single elasticity number, if defined by the size of the narrow channel, is not sufficient to account for the effects of contraction ratio on the non-linear vortex dynamics.