Determination of flow activation energy at viscosity maximum for spherical and wormlike micelles of different lengths and flexibility

Steady shear rheological measurements were carried out on aqueous solutions containing 15 mM cetyltrimethylammonium bromide (CTABr) and a constant value of [MX] and temperature for MX = 2,3-; 2,4-; 2,5-; 2,6-; 3,4-; and 3,5-Cl₂BzNa with Bz^?representing C₆H₃CO_2?. Plots of zero shear viscosity (η ₀) vs. [MX] at 35 °C and 15 mM CTABr show the presence of single maximum and double maxima for MX = 2,3- and 3,5-Cl₂BzNa, respectively. Turbidity data (absorbance at 600 nm vs. [MX]) coupled with η ₀vs. [MX] data at 35 °C reveal indirectly the presence of vesicles along with wormlike micelles (WM) at MX / CTABr > 0. 7 for MX = 3,5-Cl₂BzNa. Temperature dependence of η ₀in the vicinity of the viscosity maximum shows nonlinear and linear Arrhenius behavior, within the temperature range of 20–55 °C, for MX = 2,3-; 2,4-; 2,5-; 3,4-; and 3,5-Cl₂BzNa, respectively. The values of η ₀, $\dot {\gamma }_{\text {cr}} $ (critical shear rate), and flow activation energy correlate with CTABr micellar binding constants of counterions.     

Influence of single rectangular groove on the flow past a circular cylinder

In the present study, flow control mechanism of single groove on a circular cylinder surface is presented experimentally using Particle image velocimetry (PIV). A square shaped groove is patterned longitudinally on the surface of the cylinder with a diameter of 50 mm. The flow characteristics are studied as a function of angular position of the groove from the forward stagnation point of the cylinder within 0° ≤ θ ≤ 150°. In the current work, instantaneous and time-averaged flow data such as vorticity, ω streamline, Ψ streamwise, u/U_o and transverse, v/U_o velocity components, turbulent kinetic energy, TKE and RMS of streamwise, u_rms and transverse, v_rms velocity components are utilized in order to present the results of quantitative analyses. Furthermore, Strouhal numbers are calculated using Karman vortex shedding frequency, f_k obtained from single point spectral analysis. It is concluded that a critical angular position of the groove, θ = 80° is observed. The flow separation is controlled within 0° ≤ θ < 80°. At θ = 80°, the flow separation starts to occur in the upstream direction. The instability within the shear layer is also induced on grooved side of the cylinder with frequencies different than Karman vortex shedding frequency, f_k.     

An eddy viscosity model with near‐wall modifications

An extended version of the isotropic k–ε model is proposed that accounts for the distinct effects of low‐Reynolds number (LRN) and wall proximity. It incorporates a near‐wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

Co-current flow effects on a rising Taylor bubble

The effects of co-current flows on a rising Taylor bubble are systematically investigated by a front tracking method coupled with a finite difference scheme based on a projection approach. Both the upward (the co-current flows the same direction as the buoyancy force) and the downward (the co-current moves in the opposite direction of the buoyancy force) co-currents are examined. It is found that the upward co-current tends to elongate the bubble, while the downward co-current makes the bubble fatter and shorter. For large N_f (the inverse viscosity number), the upward co-current also elongates the skirted tail and makes the tail oscillate, while the downward co-current shortens the tail and even changes a dimpled bottom to a round shape. The upward co-current promotes the separation at the tail, while the downward co-current suppresses the separation. The terminal velocity of the Taylor bubble rising in a moving flow is a linear combination of the mean velocity (U_C) of the co-current and the terminal velocity (U₀) of the bubble rising in the stagnant liquid, and the constant is around 2 which agrees with the literature. The wake length is linearly proportional to the velocity ratio (U_C/U₀). The co-currents affect the distribution of the wall shear stresses near the bubble, but not the maximum.     

DNS Scrutiny of the <Emphasis Type="Italic">ζ</Emphasis>-<Emphasis Type="Italic">f</Emphasis> Elliptic-Relaxation Eddy-Viscosity Model in Channel Flows with a Moving Wall

In this paper we present a Direct Numerical Simulations (DNS) of channel flow with stationary and moving walls. Three cases, Poiseuille-type with U_W/U_b = 0.75, intermediate-type with U_W/U_b = 1.215, and Couette-type with U_W/U_b = 1.5 (U_W and U_b are the wall and the bulk velocity), were compared with the pure Poiseuille U_W/U_b = 0, at a bulk Reynolds number equal to 4,800 corresponding to Re _\uptau = 288_{\uptau} =288. The DNS results were used to scrutinize the capabilities of ζ-f eddy viscosity model (based on the elliptic relaxation concept) in reproducing the near-wall turbulence in non conventional flows where the shear stress structures are strongly different with respect to the cases used for models calibration. The ζ-f model (also in its basic formulation) demonstrated to have good prospects to reproduce the main phenomenology of such class of flows due to its built-in capabilities to account separately for the different (and opposite) near wall effects on turbulence: the damping due to viscosity and pressure reflection. The results of the computations demonstrated that standard ζ-f model can reasonably reproduce the phenomenology of these flows in terms of velocity and turbulent kinetic energy profiles and budgets.     

Calibration of a four-hole pyramid probe and area traverse measurements in a short-duration transonic turbine cascade tunnel

A four-hole pyramid probe has been calibrated for use in a short-duration transonic turbine cascade tunnel. The probe is used to create area traverse maps of total and static pressure, and pitch and yaw angles of the flow downstream of a transonic annular cascade. This data is unusual in that it was acquired in a short-duration (5 s of run time) annular cascade blowdown tunnel. A four-hole pyramid probe was used which has a 2.5 mm section head, and has the side faces inclined at 60° to the flow to improve transonic performance. The probe was calibrated in an ejector driven, perforated wall transonic tunnel over the Mach number range 0.5–1.2, with pitch angles from -20° to + 20° and yaw angles from-23° to +23°. A computer driven automatic traversing mechanism and data collection system was used to acquire a large probe calibration matrix (～ 10,000 readings) of non dimensional pitch, yaw, Mach number, and total pressure calibration coefficients. A novel method was used to transform the probe calibration matrix of the raw coefficients into a probe application matrix of the physical flow variables (pitch, yaw, Mach number etc.). The probe application matrix is then used as a fast look-up table to process probe results. With negligible loss of accuracy, this method is faster by two orders of magnitude than the alternative of global interpolation on the raw probe calibration matrix. The blowdown tunnel (mean nozzle guide vane blade ring diameter 1.1 m) creates engine representative Reynolds numbers, transonic Mach numbers and high levels (≈ 13%) of inlet turbulence intensity. Contours of experimental measurements at three different engine relevant conditions and two axial positions have been obtained. An analysis of the data is presented which includes a necessary correction for the finite velocity of the probe. Such a correction is non trivial for the case of fast moving probes in compressible flow.     

The Effect of Sweep-Angle Variation on the Turbulence Structure in a Separated,Three-Dimensional Flow

A three-dimensional separated flow behind a swept, backward-facing step is investigated by means of DNS for Re _H = C _∞ H/ν = 3000 with the purpose to identify changes in the statistical turbulence structure due to a variation of the sweep angle α from 0° up to 60°. With increasing sweep angle, the near-wall turbulence structure inside the separation bubble and downstream of reattachment changes due to the presence of an edge-parallel mean flow component W. Turbulence production due to the spanwise shear ∂W/∂y at the wall becomes significant and competes with the processes caused by impingement of the separated shear-layer. Changes due to a sweep angle variation can be interpreted in terms of two competing velocity scales which control the global budget of turbulent kinetic energy: the step-normal component U _∞ = C _∞cosα throughout the separated flow region and the velocity difference C _∞ across the entire shear-layer downstream of reattachment. As a consequence, the significance of history effects for the development into a two-dimensional boundary layer decreases with increasing sweep angle. For α ≥50°, near-wall streaks tend to form inside the separated flow region. Received 7 November 2000 and accepted 9 July 2002 Published online 3 December 2002 RID="*" ID="*" Part of this work was funded by the Deutsche Forschungsgemeinschaft within Sfb 557. Computer time was provided by the Konrad-Zuse Zentrum (ZIB), Berlin. Communicated by R.D. Moser     

Scaling of flow separation on a pitching low aspect ratio plate

We use two different dye injection approaches, in two different water tunnels, to visualize the formation and subsequent evolution of leading-edge vortices and related separated structures, for a pitching low aspect ratio plate. The motion is a smoothed linear pitch ramp from 0° to 40° incidence, brief hold, and return to 0°, executed at reduced pitch rates ranging from 0.1 to 0.35 and about various pivot locations. All cases evince a leading edge vortex with pronounced axial flow, which leads to formation of large-scale, three-dimensional flow structures, culminating in a large vortical structure centered at the wing symmetry plane. Pitch is also compared to plunge, where the plunge-induced angle of attack is taken as the geometric pitch incidence angle, ignoring pitch-rate effects. At successively increasing values of convective time C/U, the three-dimensional patterns of the flow structure are remarkably similar for the pitching and plunging motions. The similarity of these patterns persists, though they are shifted in time, for variation of either the location of the pitching axis or the dimensionless pitch rate.     

Rheology of SiO 2 /(acrylic polymer/epoxy) suspensions. III. Uniaxial elongational viscosity

Uniaxial elongational viscosity of SiO₂/(acrylic polymer/epoxy (AP/EP)) suspensions with various SiO₂ volume fractions (?) in a blend of acrylic polymer and epoxy was investigated at various temperatures (T). The matrix polymer ((AP/EP) blend) contained 70?vol.% of EP. At ?????35?vol.% at T????80°C, where the suspensions were in sol state, strain-hardening behavior was observed. This strain hardening of the suspensions is attributable to the elongational flow properties of (AP/EP) medium. At critical gel state (??=?35?vol.% and T?=?100°C) and in gel state (?????40?vol.%), the elongational viscosity exhibited the strain-softening behavior. These results strongly suggest that the strain softening results from the strain-induced disruption of the network structure of the SiO₂ particles therein.     

Temperature dependence of the viscosity of highly starch-filled poly(hydroxy ester ether) biodegradable composites

 The temperature dependence of the viscosity of starch-filled poly(hydroxy ester ether) (PHEE) biodegradable composites was analyzed using Arrhenius and WLF equations. Corn starch/PHEE materials were extruded using a twin screw extruder with starch volume fractions from 0.27 to 0.66. Dynamic strain sweep measurements were carried out at 10 rad/s at six different temperatures from 100 °C to 150 °C. Both Arrhenius and WLF equations model equally well the temperature effect on viscosity of PHEE and starch/PHEE composites with starch volume fractions up to 0.36. Arrhenius equation with stress correction describes the stress dependence of viscosity of starch/PHEE composites with higher starch volume fractions. The activation energy using both Arrhenius equation and Arrhenius equation with stress correction is 62.7 kJ/mol for pure PHEE and starch/PHEE composites. Received: 10 September 1999 Accepted: 27 March 2000     

The analysis of the channel flow sensitivity to the parameters of the k–ε method

This paper deals with the problem of using sensitivity analysis for fluid mechanics solutions to the constants of the standard k–ε method for 2D, incompressible and steady flows. The problem is described and analysed on the basis of a channel flow. Sensitivity coefficients of the following properties were determined: a pressure, two components of a velocity, a turbulence kinetic energy, a dissipation rate of turbulence kinetic energy and a turbulence dynamic viscosity. The calculated property values depend on five model constants that are parameters of the sensitivity analysis in this paper. Sensitivity coefficients are derivatives of the above properties, for individual parameters. In this paper these coefficients are determined using a finite difference approximation to the sensitivities coefficients. The author of this paper compares three models of the boundary layer with regard to the sensitivity of properties to the parameters. Irrespective of the boundary layer model used here, the analysis of sensitivity coefficients for the channel flow properties shows that the most sensitive property is the turbulence dissipation rate. Next properties of consequence, although of significantly smaller values of sensitivity coefficients, are the turbulence viscosity and the turbulence kinetic energy. All flow properties are mostly sensitive to the C_µ parameter. One of the final conclusions in this paper is that the analysis of sensitivity coefficient fields allows the reliable checking of results and indicates those areas most prone to calculation difficulties. Copyright © 2008 John Wiley & Sons, Ltd.     

Flow around a cylinder surrounded by a permeable cylinder in shallow water

The change in flow characteristics downstream of a circular cylinder (inner cylinder) surrounded by an outer permeable cylinder was investigated in shallow water using particle image velocimetry technique. The diameter of the inner cylinder and the water height were kept constant during the experiments as d?=?50?mm and h _w ?=?25?mm, respectively. The depth-averaged free-stream velocity was also kept constant as U?=?170?mm/s which corresponded to a Reynolds number of Re_d?=?8,500 based on the inner cylinder diameter. In order to examine the effect of diameter and porosity of the outer cylinder on flow characteristics of the inner cylinder, five different outer cylinder diameters (D?=?60, 70, 80, 90 and 100?mm) and four different porosities (???=?0.4, 0.5, 0.6 and 0.7) were used. It was shown that both porosity and outer cylinder diameter had a substantial effect on the flow characteristics downstream of the circular cylinder. Turbulent statistics clearly demonstrated that in comparison with the bare cylinder (natural case), turbulent kinetic energy and Reynolds stresses decreased remarkably when an outer cylinder was placed around the inner cylinder. Thereby, the interaction of shear layers of the inner cylinder has been successfully prevented by the presence of outer cylinder. It was suggested by referring to the results that the outer cylinder having 1.6????D/d????2.0 and 0.4????D/d????0.6 should be preferred to have a better flow control in the near wake since the peak magnitude of turbulent kinetic energy was considerably low in comparison with the natural case and it was nearly constant for these mentioned porosities ??, and outer cylinder to inner cylinder diameter ratios D/d.     

Impact comminution of solids due to local kinetic energy of high shear strain rate: I. Continuum theory and turbulence analogy

The modeling of high velocity impact into brittle or quasibrittle solids is hampered by the unavailability of a constitutive model capturing the effects of material comminution into very fine particles. The present objective is to develop such a model, usable in finite element programs. The comminution at very high strain rates can dissipate a large portion of the kinetic energy of an impacting missile. The spatial derivative of the energy dissipated by comminution gives a force resisting the penetration, which is superposed on the nodal forces obtained from the static constitutive model in a finite element program. The present theory is inspired partly by Grady's model for expansive comminution due to explosion inside a hollow sphere, and partly by analogy with turbulence. In high velocity turbulent flow, the energy dissipation rate gets enhanced by the formation of micro-vortices (eddies) which dissipate energy by viscous shear stress. Similarly, here it is assumed that the energy dissipation at fast deformation of a confined solid gets enhanced by the release of kinetic energy of the motion associated with a high-rate shear strain of forming particles. For simplicity, the shape of these particles in the plane of maximum shear rate is considered to be regular hexagons. The particle sizes are assumed to be distributed according to the Schuhmann power law. The condition that the rate of release of the local kinetic energy must be equal to the interface fracture energy yields a relation between the particle size, the shear strain rate, the fracture energy and the mass density. As one experimental justification, the present theory agrees with Grady's empirical observation that, in impact events, the average particle size is proportional to the (−2/3) power of the shear strain rate. The main characteristic of the comminution process is a dimensionless number B_a (Eq. (37)) representing the ratio of the local kinetic energy of shear strain rate to the maximum possible strain energy that can be stored in the same volume of material. It is shown that the kinetic energy release is proportional to the (2/3)-power of the shear strain rate, and that the dynamic comminution creates an apparent material viscosity inversely proportional to the (1/3)-power of that rate. After comminution, the interface fracture energy takes the role of interface friction, and it is pointed out that if the friction depends on the slip rate the aforementioned exponents would change. The effect of dynamic comminution can simply be taken into account by introducing the apparent viscosity into the material constitutive model, which is what is implemented in the paper that follows.     

A note on capillary viscometer theory: Start-up flow contribution

In gravity flow type capillary viscometer, the non-linear relationship between kinematic viscosity and efflux time can be due to two factors — the well known kinetic energy correction and the start-up flow during which the flow within the capillary attains its Poiseuillean character. The start-up contribution, however, appears to be an order of magnitude lower than that introduced by the kinetic energy correction for commonly used low-viscosity liquids encountered in practice.     

Entropy generation for pipe flow of a third grade fluid with Vogel model viscosity

The flow of fluid-solid mixtures in a pipe can be treated as non-Newtonian fluids of third grade. Depending upon the fluid viscosity, entropy generation in the flow system varies. In the present study, flow of third grade fluid in a pipe is considered. The Vogel model is introduced to account for the temperature-dependent viscosity. Entropy generation due to fluid friction and heat transfer in the flow system is formulated. The influence of viscosity parameters A and B on the entropy generation number is investigated. It is found that increasing viscosity parameter A reduces the entropy generation number and opposite is true for increasing viscosity parameter B.     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

Flow structures behind a vibrissa-shaped cylinder at different angles of attack: Complication on vortex-induced vibration

Complementary experimental studies have been conducted with a vibrissa-shaped cylinder at different angles of attack, through vortex-induced vibration (VIV) test in a wind tunnel, along with extensive measurements of wake dynamics in a water channel using time-resolved particle image velocimetry (TR-PIV). The VIV responses of an elastically mounted vibrissa-shaped cylinder are experimentally compared at various angles of attack in the range of θ = 0°–90°. At the reduced velocity of U₀/f₀D_h = 3–10 (f₀ being the system's natural frequency), the cross-flow displacement of the cylinders convincingly demonstrates that the vibrissa-shaped cylinder at a small angle of attack (θ ≤ 30°) is stable, and without appreciable displacement. Beyond θ = 30°, a prominent three-branched VIV response is readily identified, and increasing the angle of attack results in an upward shift of the synchronized region and a considerable intensification of the peak amplitude. Subsequently, TR-PIV measurements are made of the wake flow behind the vibrissa-shaped cylinder, to determine the spatio-temporally varying flow fields in two spanwise planes, i.e., the saddle and the nodal planes. Four systems with different angles of attack are chosen for comparison at Re_D = 1.8 × 10³, i.e., θ = 0°, 30°, 60° and 90°. In the two systems with θ = 0° and 30°, the wake regions feature weak velocity fluctuations in highly limited areas. However, increasing the angles of attack (to θ = 60° and 90°) gives rise to expanded recirculation zones, highly unstable flow reversals immediately behind the cylinder, and strengthened velocity fluctuations in the bulk wake regions. Cross-correlation of the fluctuating longitudinal velocities shows that at θ = 60° and 90° the energetic large-scale vortical structures form earlier, and they exert considerable influence on the near-wake fluid behind the cylinder. Finally, a sophisticated data-driven dynamic mode decomposition (DMD) process is used to extract the dominant unsteady structures in the four systems with different angles of attack. In the system with θ = 0°, two dominant DMD modes at frequencies St= 0.23 and 0.30 are identified in the saddle and the nodal planes, respectively, and those frequencies are St= 0.18 and 0.19 in the system with θ = 30°. The interaction between these dominant events at different frequencies tends to disrupt the formation of a strong vortex-shedding process. Therefore, the hydrodynamic force on the cylinder does not make a concerted contribution to suppressing the VIV behavior along the spanwise direction. In the systems with θ = 60° and 90°, the corresponding DMD modes exhibit much more synchronous, organized characteristics in the saddle and nodal planes, and unsteady events at the same frequencies are detected in both planes, reaching St = 0.14 (for θ = 60°) and 0.12 (for θ = 90°). These effects, along with the intensified vortex-shedding processes in the saddle and nodal planes, exert a concerted hydrodynamic force on the cylinder, causing it to start with an oscillatory state.     

The reduction of noise induced by flow over an open cavity

Large eddy simulations were performed and Lighthill's analogy applied in order to predict the noise induced by the flow of a turbulent boundary layer (Re_D = 12,000 or Re_θ = 300) over an open cavity (length/depth = 1). The effects of three simple geometrical changes, namely lids at the upstream edge (C10 and C20) and a chamfered edge (C23) on the downstream side of the cavity, on the internal circulation flows and shear-layer flow oscillations were determined and their effectiveness in reducing flow-induced noise was assessed. These small modifications were found to weaken the circulating flows inside the cavity and to suppress the pressure fluctuations near its downstream edge. As the lid length increases, the turbulence intensity inside the cavity is suppressed. The acoustic sources are concentrated at the downstream edge. The wall pressure fluctuations at the fundamental frequency are slightly higher for C23. The energy densities at the fundamental frequency are lower by 87.5% (C10), 30% (C20), and 18.7% (C23) than that of C00. The contributions to the acoustic directivity patterns of the dipole and quadrupole sources are oriented toward 135 and 45° respectively. The ratios of decrease of the dipole and quadrupole sources are 93% (C10), 39% (C20), and 23% (C23) with respect to C00.     

Effects of aspect ratio on the drag of a wall-mounted finite-length cylinder in subcritical and critical regimes

The effect of cylinder aspect ratio (??H/d, where H is the cylinder height or length, and d is the cylinder diameter) on the drag of a wall-mounted finite-length circular cylinder in both subcritical and critical regimes is experimentally investigated. Two cases are considered: a smooth cylinder submerged in a turbulent boundary layer and a roughened cylinder immersed in a laminar uniform flow. In the former case, the Reynolds number Re _d (??U _?? d/??, with U _?? being the free-stream velocity and ?? the fluid viscosity) was varied from 2.61?×?10⁴ to 2.87?×?10⁵, and two values of H/d (2.65 and 5) were examined; in the latter case, Re _d ?=?1.24?×?10⁴?C1.73?×?10⁵ and H/d?=?3, 5 and 7. In the subcritical regime, both the drag coefficient C _D and the Strouhal number St are smaller than their counterparts for a two-dimensional cylinder and reduce monotonously with decreasing H/d. The presence of a turbulent boundary layer causes an early transition from the subcritical to critical regime and considerably enlarges the Re _d range of the critical regime. No laminar separation bubble occurs on the finite-length cylinder immersed in the turbulent boundary layer, and consequently, the discontinuity is not observed in the C _D?CRe _d and St?CRe _d curves. In the roughened cylinder case, the Re _d range of the critical regime grows gradually with decreasing H/d, while the C _D crisis becomes less obvious. In both cases, H/d has a negligible effect on the critical value of Re _d at which transition occurs from the subcritical to critical regime.