Moderately Strong Vorticity in a Bathtub-Type Flow

Theoretical and numerical analysis is performed for an inviscid axisymmetric vortical bathtub-type flow. The level of vorticity is kept high so that the image of the flow on the radial–axial plane (r–z plane)is not potential. The most significant findings are: (1) the region of validity of the strong vortex approximation is separated from the drain by a buffer region, (2) the power-law asymptote of the stream function, specified by Δψ∼r ^4/3Δz, appears near the axis when vorticity in the flow is sufficiently strong and (3) the local Rossby number in the region of the power-law is not very sensitive to the changes of the initial vorticity level in the flow and the global Rossby number. Received 3 April 2000 and accepted 29 September 2000     

Microstructural effects on the dynamic rheology of a discotic mesophase pitch

The dynamic rheology and accompanying microstructure of a synthetic mesophase pitch (AR-HP) is reported. The loss modulus (G″) was found to be higher than the storage modulus (G′) at all frequencies (∼0.1 to ∼100 rad/s) and temperatures (280–305 °C) tested. The slope of the low-frequency terminal zone for G′ was found to be approximately 0.8, much lower than a value of 2 that is observed for flexible chain polymers. Cross-polarized optical microscopy with full-wave retardation plate confirmed the presence of different textures of mesophase pitch processed under different conditions. While loss moduli remained fairly unchanged, finer textures led to significantly lower storage moduli. Consistent with this trend, coarsening of the microstructure during textural relaxation led to an increase in storage moduli. Therefore, for the discotic mesophase pitch, the viscous component was found to remain unaffected by the microstructure, but the elastic modulus was dependent on the orientation of layer-planes and size of the texture.     

The rheological characterisation of bread dough using capillary rheometry

An Australian hard wheat flour–water dough has been characterised using parallel plate and capillary rheometers over an extensive range of apparent shear rates (10^− 3–10³ s^− 1) relevant to process conditions. Torsional measurements showed that the shear viscosity of the dough increased with strain to a maximum value and then decreased, suggesting a breakdown of the dough structure. Both torsional and capillary experiments revealed the shear-thinning behaviour of the dough. The wall slip phenomenon in capillary rheometry was investigated and found to be diameter dependent and occurred at a critical shear stress of approximately 5–10 kPa. A two-regime power law behaviour was observed, with the power law index approximately 0.3 in the low shear rate range increasing to 0.67 in the high shear rate range. Pressure fluctuation was observed in the capillary data and increased with shear rate, in particular, at shear rates approaching 10⁴ s^− 1. The results demonstrate that capillary rheometry is a viable means of rheologically testing dough at high shear rates provided pressure fluctuation is carefully monitored and capillary rheometry corrections, including wall slip, are accounted for.     

Three views of viscoelasticity for Cox–Merz materials

A slight rearrangement of the classical Cox and Merz rule suggests that the shear stress value of steady shear flow, , and complex modulus value of small amplitude oscillatory shear, G^∗(ω) = (G^′2 + G^″2)^1/2, are equivalent in many respects. Small changes of material structure, which express themselves most sensitively in the steady shear stress, τ, show equally pronounced in linear viscoelastic data when plotting these with G^∗ as one of the variables. An example is given to demonstrate this phenomenon: viscosity data that cover about three decades in frequency get stretched out over about nine decades in G^∗ while maintaining steep gradients in a transition region. This suggests a more effective way of exploiting the Cox–Merz rule when it is valid and exploring reasons for lack of validity when it is not. The τ −G^∗ equivalence could also further the understanding of the steady shear normal stress function as proposed by Laun.     

Rheo-PIV of a shear-banding wormlike micellar solution under large amplitude oscillatory shear

We explore the behavior of a wormlike micellar solution under both steady and large amplitude oscillatory shear (LAOS) in a cone–plate geometry through simultaneous bulk rheometry and localized velocimetric measurements. First, particle image velocimetry is used to show that the shear-banded profiles observed in steady shear are in qualitative agreement with previous results for flow in the cone–plate geometry. Then under LAOS, we observe the onset of shear-banded flow in the fluid as it is progressively deformed into the non-linear regime—this onset closely coincides with the appearance of higher harmonics in the periodic stress signal measured by the rheometer. These harmonics are quantified using the higher-order elastic and viscous Chebyshev coefficients e _n and v _n , which are shown to grow as the banding behavior becomes more pronounced. The high resolution of the velocimetric imaging system enables spatiotemporal variations in the structure of the banded flow to be observed in great detail. Specifically, we observe that at large strain amplitudes (γ ₀ ≥ 1), the fluid exhibits a three-banded velocity profile with a high shear rate band located in-between two lower shear rate bands adjacent to each wall. This band persists over the full cycle of the oscillation, resulting in no phase lag being observed between the appearance of the band and the driving strain amplitude. In addition to the kinematic measurements of shear banding, the methods used to prevent wall slip and edge irregularities are discussed in detail, and these methods are shown to have a measurable effect on the stability boundaries of the shear-banded flow.     

On Saint-Venant's Principle for a Homogeneous Elastic Arch-Like Region

In this paper, we consider a two-dimensional homogeneous isotropic elastic material state in the arch-like region a ≤ r ≤ b, 0 ≤ θ ≤ α, where (r, θ) denote plane polar coordinates. We assume that three of the edges r = a, r = b, θ = α are traction-free, while the edge θ = 0 is subjected to an (in plane) self-equilibrated load. We define an appropriate measure for the Airy stress function φ and then we establish a clear relationship with the Saint-Venant's principle on such regions. We introduce a cross-sectional integral function I(θ) which is shown to be a convex function and satisfies a second-order differential inequality. Consequently, we establish a version of the Saint-Venant principle for such a curvilinear strip, without requiring of any condition upon the dimensions of the arch-like region.     

Conoscopic observations of shear-induced rotations in nematic liquid crystals

Orientational changes in monodomains of flow-aligning liquid crystals, 4-n-pentyl-4′-cyanobiphenyl and N-(4-methoxybenzylidene)-4-butylaniline, were studied during shear and recovery in a linear shearing device fitted to an optical microscope. Planar alignment (director in the shear plane) allows the study of twist effects and was generated by strong planar anchoring at the walls with orientations in a range of 0–90° with the shear direction. While being held back by the anchoring walls, shear caused the bulk director to rotate towards a steady-state alignment angle in the shear direction (Leslie angle θ_L). The transient director rotation was observed with conoscopy. It was found that increasing the initial alignment towards the vorticity direction increased the measured θ_L. Upon stopping the flow, the bulk director returned to its initial state. With initial alignment orientation changing from parallel to perpendicular to the flow direction, the rate of the twist-driven recovery process increases. This rate increase is not seen in the splay-driven recovery which is constant and consistently faster than twist-driven recovery at all orientations studied. Received: 10 December 1998/Accepted: 7 June 1999     

Viscoelastic properties of ultra-high viscosity alginates

Ultra-high viscosity alginates were extracted from the brown seaweeds Lessonia nigrescens (UHV_N, containing 61% mannuronate (M) and 2% guluronate (G)) and Lessonia trabeculata (UHV_T, containing 22% M and 78% G). The viscoelastic behavior of the aqueous solutions of these alginates was determined in shear flow in terms of the shear stress σ ₂₁, the first normal stress difference N ₁, and the shear viscosity η in isotonic NaCl solutions (0.154 mol/L) at T = 298 K in dependence of the shear rate [(g)\dot]\dot{\gamma} for solutions of varying concentrations and molar masses (3–10 × 10⁵ g/mol, homologous series was prepared by ultrasonic degradation). Data obtained in small-amplitude oscillatory shear (SAOS) experiments obey the Cox–Merz rule. For comparison, a commercial alginate with intermediate chemical composition was additionally characterized. Particulate substances which are omnipresent in most alginates influenced the determination of the material functions at low shear rates. We have calculated structure–property relationships for the prediction of the viscosity yield, e.g., η–M _w–c–[(g)\dot]\dot{\gamma} for the Newtonian and non-Newtonian region. For the highest molar masses and concentrations, the elasticity yield in terms of N ₁ could be determined. In addition, the extensional flow behavior of the alginates was measured using capillary breakup extensional rheometry. The results demonstrate that even samples with the same average molar mass but different molar mass distributions can be differentiated in contrast to shear flow or SAOS experiments.     

Simultaneous determination of shear and extensional viscosities and wall slip by on-line rheometry with parameter fitting: an application to EPDM rubber

Shear and extensional viscosities and wall slip are determined simultaneously under extrusion processing conditions using an on-line rheometer. Because it is not possible to independently control flow rate and temperature, classical methods for interpretation of capillary data cannot be used with on-line rheometry. This limitation is overcome using computational optimization to fit parameters in a flow model. This consists of three parts, representing shear viscosity, extensional viscosity, and wall slip. Three-parameter, power law forms, based on local instantaneous deformation rates and including temperature dependence, are used for each, and analytic solutions applied for entry flow and flow in the capillary. For entry flow, the Cogswell–Binding approach is used, and for developed flow in the capillary a solution incorporating wall slip is derived. The rheometer, with interchangeable capillaries, is mounted in place of the die on a rubber profile extrusion line. Pressure drops and temperatures for extrusion of an EPDM rubber through 2 mm diameter capillaries of length 0, 2, 3, 4, and 5 mm are logged and flow rates determined for a range of extruder speeds (5 to 20 rpm). Pressures ranged from 60 to 75 bar and temperatures from 86 to 116 °C. Mean flow velocity in the capillaries was between 5 × 10⁻³ and 5 × 10⁻¹ m s⁻¹. The nine material parameters are optimized for best fit of the analytic pressure drops to experimental data, using about 100 data points, with the Levenberg–Marquardt method. It is concluded that flow is dominated by extension and wall slip. Shear flow appears to play little part. The slip model indicates that slip velocity increases much more rapidly than the wall shear stress (in the range 0.5–1 MPa) and decreases with temperature for a given stress level. Results for the (uniaxial) extensional viscosity represent an engineering approximation to this complex phenomenon at the high strains (approximately 200) and high extension rates (up to 800 s⁻¹) applying in the extrusion. Results indicate a slight extension hardening and a decrease with temperature. Results are put into the context of the available studies in the literature, which, particularly with regard to wall-slip and extensional flow, consider conditions far removed from those applying in industrial extrusion. The present methods provide a powerful means for flow characterization under processing conditions, providing data suitable for use in computer simulations of extrusion and optimization of die design.     

Numerical simulation of fluid flow and heat transfer characteristics in channel with V corrugated plates

A detailed numerical study is carried out to investigate fluid flow and heat transfer characteristics in a channel with heated V corrugated upper and lower plates. The parameters studied include the Reynolds number (Re = 2,000–5,500), angles of V corrugated plates (θ = 20°, 40°, 60°), and constant heat fluxs (q″ = 580, 830, 1,090 W/m²). Numerical results have been validated using the experimented data reported by Naphon, and a good agreement has been found. The angles of V corrugated plates (θ) and the Reynolds number are demonstrated to significantly affect the fluid flow and the heat transfer rate. Increasing the angles of V corrugated plates can make the heat transfer performance become better. The increasing Reynolds number leads to a more complex fluid flow and heat transfer rate. The numerical calculations with a non-equilibrium wall function have a better accuracy than with a standard wall function for solving high Reynolds numbers or complex flow problems.     

Near-ground tornado-like vortex structure resolved by particle image velocimetry (PIV)

The near-ground flow structure of tornadoes is of utmost interest because it determines how and to what extent civil structures could get damaged in tornado events. We simulated tornado-like vortex flow at the swirl ratios of S = 0.03–0.3 (vane angle θ_v = 15°–60°), using a laboratory tornado simulator and investigated the near-ground-vortex structure by particle imaging velocimetry. Complicated near-ground flow was measured in two orthogonal views: horizontal planes at various elevations (z = 11, 26 and 53 mm above the ground) and the meridian plane. We observed two distinct vortex structures: a single-celled vortex at the lowest swirl ratio (S = 0.03, θ_v = 15°) and multiple suction vortices rotating around the primary vortex (two-celled vortex) at higher swirl ratios (S = 0.1–0.3, θ_v = 30°–60°). We quantified the effects of vortex wandering on the mean flow and found that vortex wandering was important and should be taken into account in the low swirl ratio case. The tangential velocity, as the dominant velocity component, has the peak value about three times that of the maximum radial velocity regardless of the swirl ratio. The maximum velocity variance is about twice at the high swirl ratio (θ_v = 45°) that at the low swirl ratio (θ_v = 15°), which is contributed significantly by the multiple small-scale secondary vortices. Here, the results show that not only the intensified mean flow but greatly enhanced turbulence occurs near the surface in the tornado-like vortex flow. The intensified mean flow and enhanced turbulence at the ground level, correlated with the ground-vortex interaction, may cause dramatic damage of the civil structures in tornadoes. This work provides detailed characterization of the tornado-like vortex structure, which has not been fully revealed in previous field studies and laboratory simulations. It would be helpful in improving the understanding of the interaction between the tornado-like vortex structure and the ground surface, ultimately leading to better predictions of tornado-induced wind loads on civil structures.     

Flow-induced crystallization of high-density polyethylene: the effects of shear and uniaxial extension

The effects of shear, uniaxial extension and temperature on the flow-induced crystallization of two different types of high-density polyethylene (a metallocene and a ZN-HDPE) are examined using rheometry. Shear and uniaxial extension experiments were performed at temperatures below and well above the peak melting point of the polyethylenes in order to characterize their flow-induced crystallization behavior at rates relevant to processing (elongational rates up to 30 s^− 1 and shear rates 1 to 1,000 s^− 1 depending on the application). Generally, strain and strain rate found to enhance crystallization in both shear and elongation. In particular, extensional flow was found to be a much stronger stimulus for polymer crystallization compared to shear. At temperatures well above the melting peak point (up to 25°C), polymer crystallized under elongational flow, while there was no sign of crystallization under simple shear. A modified Kolmogorov crystallization model (Kolmogorov, Bull Akad Sci USSR, Class Sci, Math Nat 1:355–359, 1937) proposed by Tanner and Qi (Chem Eng Sci 64:4576–4579, 2009) was used to describe the crystallization kinetics under both shear and elongational flow at different temperatures.     

Measurements and scaling of wall shear stress fluctuations

Measurements of velocity and wall shear stress fluctuations were made in an external turbulent boundary layer developed over a towed surface-piercing flat plate. An array of eight flush-mounted wall shear stress sensors was used to compute the space-time correlation function. A methodology for in situ calibration of the sensors for ship hydrodynamic applications is presented. The intensity of the wall shear stress fluctuations, τ _rms/τ _avg was measured as 0.25 and 0.36 for R _θ =3,150 and 2,160 respectively. The probability density is shown to exhibit positive skewness, and lack of flow reversals at the wall. Correlations between velocity and wall shear stress fluctuations are shown to collapse with outer boundary layer length and velocity scales, verifying the existence of large-scale coherent structures which convect and decay along the wall at an angle of inclination varying from 10 to 13° over the range of Reynolds numbers investigated. The wall shear stress convection velocity determined from narrow band correlation measurements is shown to scale with outer variables. The space-time correlation of the wall shear is shown to exhibit a well-defined convective ridge, and to decay 80% over approximately 3δ for R _θ =3,150. Published online: 7 November 2002     

Three-dimensional elastic behavior of a nonhomogeneous layer

Summary  This paper deals with the theoretical treatment of a three-dimensional elastic problem governed by a cylindrical coordinate system (r,θ,z) for a medium with nonhomogeneous material property. This property is defined by the relation G(z)=G ₀(1+z/a) ^m where G ₀,a and m are constants, i.e., shear modulus of elasticity G varies arbitrarily with the axial coordinate z by the power product form. We propose a fundamental equation system for such nonhomogeneous medium by using three kinds of displacement functions and, as an illustrative example, we apply them to an nonhomogeneous thick plate (layer) subjected to an arbitrarily distributed load (not necessarily axisymmetric) on its surfaces. Numerical calculations are carried out for several cases, taking into account the variation of the nonhomogeneous parameter m. The numerical results for displacement and stress components are shown graphically. Received 10 May 1999; accepted for publication 15 August 1999     

Dynamics of hairpin vortices generated by a mixing tab in a channel flow

To better understand mixing by hairpin vortices, time-series particle image velocimetry (PIV) was applied to the wake of a trapezoidal-shaped passive mixing tab mounted at the bottom of a square turbulent channel (Re _h =2,080 based on the tab height). Instantaneous velocity/vorticity fields were obtained in sequences of 10 Hz in the tab wake in the center plane (x–y) and in a plane (x–z) parallel to the wall. Periodically-shed hairpin vortices were clearly identified and seen to rise as they advected downstream. Experimental evidence shows that the vortex-induced ejection of the near-wall viscous fluid to the immediate upstream is important to the dynamics of hairpin vortices. It can increase the strength of the hairpin vortices in the near tab region and cause generation of secondary hairpin vortices further downstream when the hairpin heads are farther away from the wall. Measurements also reveal the existence of a type of new secondary vortice with the opposite-sign spanwise vorticity. The distribution of vortex loci in the x–y plane shows that the hairpin vortices and the reverse vortices are spatially segregated in distinct layers. Turbulence statistics, including mean velocity profiles, Reynolds stresses, and turbulent kinetic energy dissipation rate distributions, were obtained from the PIV data. These statistical quantities clearly reveal imprints of the identified vortex structures and provide insight into mixing effectiveness. Received: 24 February 2000/Accepted: 24 October 2000     

Effect of local forcing on a turbulent boundary layer

An experimental study is performed to analyze flow structures behind local suction and blowing in a flat-plate turbulent boundary layer. The local forcing is given to the boundary layer flow by means of a sinusoidally oscillating jet issuing from a thin spanwise slot at the wall. The Reynolds number based on the momentum thickness is about Re _θ =1700. The effects of local forcing are scrutinized by altering the forcing frequency (0.011 ≤ f ⁺≤ 0.044). The forcing amplitude is fixed at A ₀=0.4. It is found that a small local forcing reduces the skin friction and the skin friction reduction increases with the forcing frequency. A phase-averaging technique is employed to capture the large-scale vortex evolution. An organized spanwise vortical structure is generated by the local forcing. The cross-sectional area of vortex and the time fraction of vortex are examined by changing the forcing frequency. An investigation of the random fluctuation components reveals that turbulent energy is concentrated near the center of vortical structures. Received: 17 March 2000/Accepted: 3 April 2001     

Numerical investigation of turbulent flow,heat transfer and entropy generation in a helical coiled tube with larger curvature ratio

Three-dimensional turbulent forced convective heat transfer and flow characteristics, and the non-dimensional entropy generation number in a helical coiled tube subjected to uniform wall temperature are simulated using the k–ε standard turbulence model. A finite volume method is employed to solve the governing equations. The effects of Reynolds number, curvature ratio, and coil pitch on the average friction factor and Nusselt number are discussed. The results presented in this paper cover a Reynolds number range of 2 × 10⁴ to 6 × 10⁴, a pitch range of 0.1–0.2 and a curvature ratio range of 0.1–0.3. The results show that the coil pitch, curvature ratio and Reynolds number have different effects on the average friction factor and Nusselt number at different cross-sections. In addition, the flow and heat transfer characteristics in a helical coiled tube with a larger curvature ratio for turbulent flow are different from that of smaller curvature ratio for laminar and turbulent flow in certain ways. Some new features that are not obtained in previous researches are revealed. Moreover, the effects of Reynolds number, curvature ratio, and coil pitch on the non-dimensional entropy generation number of turbulent forced convection in a helical coiled tube are also discussed.     

Characterization of the flow properties of sodium carboxymethylcellulose via mechanical and optical techniques

Sodium carboxymethylcellulose (NaCMC) in solution represents a complex rheological system, since it forms aggregates and associations and hence higher-level structures and, depending on the synthesis, is only found in a molecularly dispersed form in exceptional cases. Rheo-mechanical investigations of the viscoelasticity showed that the Cox-Merz rule is not fulfilled. The aim was therefore to examine whether rheo-optics could be employed to provide more detailed conclusions about the parameters that influence the flow behavior of NaCMC than has hitherto been available with mechanical methods. The flow birefringence, Δn ^′, rises as the degree of polymerization increases, and exhibits the same dependence on molar mass as does the viscosity: Δn ^′∝M _w ^3.4. As the degree of polymerization increases while the shear rate remains constant, the polymer segments become more distinctly aligned in the direction of shear. Hence increasing the degree of polymerization also affects the solution structure, i.e. the interaction of the molecules with one another. The stress-optical rule only applies to a limited extent for this system. The stress-optical coefficient, C, is almost independent of the shear rate, but is strongly influenced by the concentration and attains a limiting value of 3 × 10⁻⁸ Pa⁻¹. C was determined for a polymer in dilute solution and the curve obtained also enabled transitions in the solution structure to be recognized. Received: 1 May 1998 Accepted: 5 October 1998     

Three-dimensional Turbulent Boundary Layer in a Shrouded Rotating System

A thre-dimensional direct numerical simulation is combined with a laboratory study to describe the turbulent flow in an enclosed annular rotor-stator cavity characterized by a large aspect ratio G = (b − a)/h = 18.32 and a small radius ratio a/b = 0.152, where a and b are the inner and outer radii of the rotating disk and h is the interdisk spacing. The rotation rate Ω considered is equivalent to the rotational Reynolds number Re = Ωb ²/ν= 9 .5 × 10⁴ (ν the kinematic viscosity of water). This corresponds to a value at which experiment has revealed that the stator boundary layer is turbulent, whereas the rotor boundary layer is still laminar. Comparisons of the computed solution with velocity measurements have given good agreement for the mean and turbulent fields. The results enhance evidence of weak turbulence by comparing the turbulence properties with available data in the literature (Lygren and Andersson, J Fluid Mech 426:297–326, 2001). An approximately self-similar boundary layer behavior is observed along the stator. The wall-normal variations of the structural parameter and of characteristic angles confirm that this boundary layer is three-dimensional. A quadrant analysis (Kang et al., Phys Fluids 10:2315–2322, 1998) of conditionally averaged velocities shows that the asymmetries obtained are dominated by Reynolds stress-producing events in the stator boundary layer. Moreover, Case 1 vortices (with a positive wall induced velocity) are found to be the major source of generation of special strong events, in agreement with the conclusions of Lygren and Andersson (J Fluid Mech 426:297–326, 2001).     

Reynolds Number Effects on the Coherent Dynamics of the Turbulent Horseshoe Vortex System

The adverse pressure gradient induced by a surface-mounted obstacle in a turbulent boundary layer causes the approaching flow to separate and form a dynamically rich horseshoe vortex system (HSV) in the junction of the obstacle with the wall. The Reynolds number of the flow (Re) is one of the important parameters that control the rich coherent dynamics of the vortex, which are known to give rise to low-frequency, bimodal fluctuations of the velocity field (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). We carry out detached eddy simulations (DES) of the flow past a circular cylinder mounted on a rectangular channel for Re = 2.0 × 10⁴ and 3.9 × 10⁴ (Dargahi, Exp Fluids 8:1–12, 1989) in order to systematically investigate the effect of the Reynolds number on the HSV dynamics. The computed results are compared with each other and with previous experimental and computational results for a related junction flow at a much higher Reynolds number (Re = 1.15 × 10⁵) (Devenport and Simpson, J Fluid Mech 210:23–55, 1990; Paik et al., Phys Fluids 19:045107, 2007). The computed results reveal significant variations with Re in terms of the mean-flow quantities, turbulence statistics, and the coherent dynamics of the turbulent HSV. For Re = 2.0 × 10⁴ the HSV system consists of a large number of necklace-type vortices that are shed periodically at higher frequencies than those observed in the Re = 3.9 × 10⁴ case. For this latter case the number of large-scale vortical structures that comprise the instantaneous HSV system is reduced significantly and the flow dynamics becomes quasi-periodic. For both cases, we show that the instantaneous flowfields are dominated by eruptions of wall-generated vorticity associated with the growth of hairpin vortices that wrap around and disorganize the primary HSV system. The intensity and frequency of these eruptions, however, appears to diminish rapidly with decreasing Re. In the high Re case the HSV system consists of a single, highly energetic, large-scale necklace vortex that is aperiodically disorganized by the growth of the hairpin mode. Regardless of the Re, we find pockets in the junction region within which the histograms of velocity fluctuations are bimodal as has also been observed in several previous experimental studies.