Characterizing developing adverse pressure gradient flows subject to surface roughness

An experimental study was conducted to examine the effects of surface roughness and adverse pressure gradient (APG) on the development of a turbulent boundary layer. Hot-wire anemometry measurements were carried out using single and X-wire probes in all regions of a developing APG flow in an open return wind tunnel test section. The same experimental conditions (i.e., T _∞, U _ref, and C _p) were maintained for smooth, k ⁺ = 0, and rough, k ⁺ = 41–60, surfaces with Reynolds number based on momentum thickness, 3,000 < Re _θ < 40,000. The experiment was carefully designed such that the x-dependence in the flow field was known. Despite this fact, only a very small region of the boundary layer showed a balance of the various terms in the integrated boundary layer equation. The skin friction computed from this technique showed up to a 58% increase due to the surface roughness. Various equilibrium parameters were studied and the effect of roughness was investigated. The generated flow was not in equilibrium according to the Clauser (J Aero Sci 21:91–108, 1954) definition due to its developing nature. After a development region, the flow reached the equilibrium condition as defined by Castillo and George (2001), where Λ = const, is the pressure gradient parameter. Moreover, it was found that this equilibrium condition can be used to classify developing APG flows. Furthermore, the Zagarola and Smits (J Fluid Mech 373:33–79, 1998a) scaling of the mean velocity deficit, U _∞δ*/δ, can also be used as a criteria to classify developing APG flows which supports the equilibrium condition of Castillo and George (2001). With this information a ‘full APG region’ was defined.     

Radial Solutions and Phase Separation in a System of Two Coupled Schrödinger Equations

We consider the nonlinear elliptic system

Very high Reynolds number boundary layers over 3D sparse roughness and obstacles: the mean flow

Hot-wire and oil-film interferometry measurements are taken for 3D rough wall boundary layers at very high Reynolds numbers (61,000 < Re θ < 120,000) with low blockage ratios, 10 < δ/H < 135, and high roughness, 100 < H ⁺ < 4,900. The results cover flows over both rough walls and over obstacles and are compared with and provide extension to recent lower Reynolds number results. The validity of the Townsend ‘wall similarity hypothesis’ in relation to consistently increasing 3D roughness is interrogated. In agreement with recent work, Schultz and Flack (J Fluid Mech 580:381–405, 2007) and Castro (J Fluid Mech 585:469–485, 2007) found that, for relatively low roughness, Townsend’s hypothesis holds for the mean velocity field. With increasing roughness, the equilibrium layer diminishes and gradually vanishes. The viscous component of the wall shear stress decreases, while the turbulent component increases as the roughness effects extend across the boundary layer.     

Exact analytic solutions for free convection boundary layers on a heated vertical plate with lateral mass flux embedded in a saturated porous medium

 The problem of the self-similar boundary flow of a “Darcy-Boussinesq fluid” on a vertical plate with temperature distribution T _w(x) = T _∞+A·x ^λ and lateral mass flux v _w(x) = a·x ^(λ−1)/2, embedded in a saturated porous medium is revisited. For the parameter values λ = 1,−1/3 and −1/2 exact analytic solutions are written down and the characteristics of the corresponding boundary layers are discussed as functions of the suction/ injection parameter in detail. The results are compared with the numerical findings of previous authors. Received on 8 March 1999     

On near-wall hot-wire measurements

A specially constructed hot-wire probe was used to obtain very near-wall velocity measurements in both a fully developed turbulent channel flow and flat plate boundary layer flow. The near-wall hot-wire probe, having been calibrated in a specially constructed laminar flow calibration rig, was used to measure the mean streamwise velocity profile, distributions of streamwise and spanwise intensities of turbulence and turbulence kinetic energy k in the viscous sublayer and beyond; these distributions compare very favorably with available DNS results obtained for channel flow. While low Reynolds number effects were clearly evident for the channel flow, these effects are much less distinct for the boundary layer flow. By assuming the dissipating range of eddy sizes to be statistically isotropic and the validity of Taylor's hypothesis, the dissipation rate ɛ _iso in the very near-wall viscous sublayer region and beyond was determined for both the channel and boundary layer flows. It was found that if the convective velocity U _c in Taylor's hypothesis was assumed to be equal to the mean velocity Uˉ at the point of measurement, the value of (ɛ⁺ _iso)₁ thus obtained agrees well with that of (ɛ ⁺)_DNS for y ⁺ ≥ 80 for channel flow; this suggests the validity of assuming U _c=Uˉ and local isotropy for large values of y ⁺. However, if U _c was assumed to be 10.6u _τ , the value of (ɛ⁺ _iso)₂ thus obtained was found to compare reasonably well with the distribution of (ɛ⁺ _iso)_DNS for y ⁺≤ 15. Received: 31 May 1999/Accepted: 20 December 1999     

Velocity-defect scaling for turbulent boundary layers with a range of relative roughness

Velocity profile measurements in zero pressure gradient, turbulent boundary layer flow were made on a smooth wall and on two types of rough walls with a wide range of roughness heights. The ratio of the boundary layer thickness (δ) to the roughness height (k) was 16≤δ/k≤110 in the present study, while the ratio of δ to the equivalent sand roughness height (k _s) ranged from 6≤δ/k _s≤91. The results show that the mean velocity profiles for all the test surfaces agree within experimental uncertainty in velocity-defect form in the overlap and outer layer when normalized by the friction velocity obtained using two different methods. The velocity-defect profiles also agree when normalized with the velocity scale proposed by Zagarola and Smits (J Fluid Mech 373:33–70, 1998). The results provide evidence that roughness effects on the mean flow are confined to the inner layer, and outer layer similarity of the mean velocity profile applies even for relatively large roughness.     

Dependence of the zero shear-rate viscosity and the viscosity function of linear high-density polyethylenes on the mass-average molar mass and polydispersity

Linear high-density polyethylenes with molar masses M _w between 240 and 1,000,000 g/mol, obtained by metallocene catalysts, were characterized in shear using oscillatory and creep tests. The polydispersities of the molar mass distributions (MMDs) lay between 1 and 16. The resulting zero shear-rate viscosities η₀ covered a range from 2.5×10⁻³ to around 10⁸ Pas. Above a critical molar mass of M _c≈2,900 g/mol, the experimental results can be described by the relation η₀ ∼ M _w^3.6, independently of the MMD. The oscillatory data were fitted with a Carreau–Yasuda equation. The resulting parameters were correlated to molecular structure. The parameter a, being a quantity for the width of the transition between the Newtonian and the non-Newtonian regime, showed a dependence on the molar mass M _w but not on M _w/M _n. The parameter λ of the Carreau-Yasuda equation was found to be the reciprocal crossover frequency for all samples with a log-Gaussian MMD. λ depends on the molar mass M _w and also on M _w/M _n.

Viscous dissipation effects on thermal entrance region heat transfer in pipes with uniform wall heat flux

The analytical solution to Graetz problem with uniform wall heat flux is extended by including the viscous dissipation effect in the analysis. The analytical solution obtained reduces to that of Siegel, Sparrow and Hallman neglecting viscous dissipation as a limiting case. The sample developing temperature profiles, wall and bulk temperature distributions and the local Nusselt number variations are presented to illustrate the viscous dissipation effects. It is found that the role of viscous dissipation on thermal entrance region heat transfer is completely different for heating and cooling at wall. In the case of cooling at wall, a critical value of Brinkman number, Br _c=−11/24, exists beyond which (−11/24<Br<0) the fluid bulk temperature will always be less than the uniform entrance temperature indicating the predominance of cooling effect over the viscous heating effect. On the other hand, with Br < Br _c the bulk temperature T _b will approach the wall temperature T _w at some downstream position and from there onward the bulk temperature T _b becomes less than the wall temperature T _w with T _w > B _b > T ₀ indicating overall heating effect for the fluid. The numerical results for the case of cooling at wall Br < 0 are believed to be of some interest in the design of the proposed artctic oil pipeline.     

Rheo-dielectric behavior of low molecular weight liquid crystal 2. Behavior of 8CB in nematic and smectic states

Rheo-dielectric behavior was examined for 4−4^′−n-octyl-cyanobiphenyl (8CB) having large dipoles parallel to its principal axis (in the direction of the C≡N bond). In the quiescent state at all temperatures (T) examined, orientational fluctuation of the 8CB molecules was observed as dielectric dispersions at characteristic frequencies ω_c>10⁶ s⁻¹. In the isotropic state at high T, no detectable changes of the complex dielectric constant ɛ^*(ω) were found under slow flow at shear rates ˙γ≫ω_c. In the nematic state at intermediate T, the terminal relaxation intensity of ɛ^*(ω) was decreased under such slow flow. In the smectic state at lower T, the flow effect became much less significant. These results were related to the flow-induced changes of the liquid crystalline textures in the nematic and smectic states, and the differences of the rheo-dielectric behavior in these states are discussed in relation to a difference of the symmetry of molecular arrangements in the nematic and smectic textures. Received: 1 October 1998 Accepted: 13 January 1999     

Stability of an Initially, Stably Stratified Fluid Subjected to a Step Change in Temperature

The onset of convective instability in an initially quiescent, stably stratified fluid layer between two horizontal plates is analyzed with linear theory. The bottom boundary is heated suddenly from below, subjected to a step change in surface temperature. The critical time t _c to mark the onset of Rayleigh-Bénard convection is predicted by propagation theory. This theory uses the length scaled by , where α denotes thermal diffusivity. Under the normal mode analysis the dimensionless disturbance equations are obtained as a function of τ(=αt/d ²) and ζ(=Z/), where d is the fluid layer depth and Z is the vertical distance. The resulting equations are transformed to self-similar ones by using scaling and finally fixing τ as τ_c under the frame of coordinates τ and ζ. For a given γ, Pr and τ_c, the minimum value of Ra is obtained from the marginal stability curve. Here γ denotes the temperature ratio to represent the degree of stabilizing effect, Pr is the Prandtl number and Ra is the Rayleigh number. With γ=0, the minimum Ra value approaches the well-known value of 1708 as τ_c increases. However, it is inversely proportional to τ_c ^3/2 as τ_c decreases. With increasing γ, the system becomes more stable. It is interesting that in the present system, propagation theory produces the stability criteria to bound the available experimental data over the whole domain of time. Received 5 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 RID="*" ID="*" This work has been supported by both SK Chemicals Co. Ltd. and LG Chemical Ltd., Seoul under the Brain Korea 21 Project of the Ministry of Education. Communicated by H.J.S. Fernando     

Global Semigroup of Conservative Solutions of the Nonlinear Variational Wave Equation

We prove the existence of a global semigroup for conservative solutions of the nonlinear variational wave equation u _tt − c(u)(c(u)u _x ) _x  = 0. We allow for initial data u| _{t = 0} and u _t | _t=0 that contain measures. We assume that 0 < k^-1 \leqq c(u) \leqq k{0 < \kappa^{-1} \leqq c(u) \leqq \kappa}. Solutions of this equation may experience concentration of the energy density (u_t²+c(u)²u_x²)dx{(u_t^2+c(u)^2u_x^2){\rm d}x} into sets of measure zero. The solution is constructed by introducing new variables related to the characteristics, whereby singularities in the energy density become manageable. Furthermore, we prove that the energy may focus only on a set of times of zero measure or at points where c′(u) vanishes. A new numerical method for constructing conservative solutions is provided and illustrated with examples.     

A Beale–Kato–Majda Criterion for Three Dimensional Compressible Viscous Heat-Conductive Flows

We prove a blow-up criterion in terms of the upper bound of (ρ, ρ ⁻¹, θ) for a strong solution to three dimensional compressible viscous heat-conductive flows. The main ingredient of the proof is an a priori estimate for a quantity independently introduced in Haspot (Regularity of weak solutions of the compressible isentropic Navier–Stokes equation, arXiv:1001.1581, 2010) and Sun et al. (J Math Pure Appl 95:36–47, 2011), whose divergence can be viewed as the effective viscous flux.     

Computation of conjugate heat transfer in the turbulent mixed convection regime in a vertical channel with multiple heat sources

This paper presents the results of a comprehensive numerical study to analyze conjugate, turbulent mixed convection heat transfer from a vertical channel with four heat sources, uniformly flush-mounted to one of the channel walls. The results are presented to study the effect of various parameters like thermal conductivity of wall material (k _s), thermal conductivity of flush-mounted discrete heat source (k _c), Reynolds number of fluid flow (Re _s), modified Richardson number (Ri ⁺) and aspect ratio (AR) of the channel. The standard k-ε turbulence model, modified by including buoyancy effects with physical boundary conditions, i.e. without wall functions, has been used for the analysis. Semi-staggered, non-uniform grids are used to discretise the two dimensional governing equations, using finite volume method. A correlation, encompassing a wide range of parameters, is developed for the non-dimensional maximum temperature (T ^*) using the asymptotic computational fluid dynamics (ACFD) technique.     

Control of Flow and Heat Transfer in a Porous Enclosure due to an Adiabatic Thin fin on the Hot Wall

Natural convection flow in a differentially heated square enclosure filled with porous matrix with a solid adiabatic thin fin attached at the hot left wall is studied numerically. The Brinkman–Forchheimer-extended Darcy model is used to solve the momentum equations, in the porous medium. The numerical investigation is done through streamlines, isotherms, and heat transfer rates. A parametric study is carried out using the following parameters: Darcy number (Da) from 10⁻⁴ to 10⁻², dimensionless thin fin lengths (L _p) 0.3, 0.5, and 0.7, dimensionless positions (S _p) 0.25, 0.5, and 0.75 with Prandtl numbers (Pr) 0.7 and 100 for Ra = 10⁶. For Da = 10⁻³ and Pr = 0.7, it is observed that there is a counter clock-wise secondary flow formation around the tip of the fin for S _p = 0.5 for all lengths of L _p. Moreover when Da = 10⁻² the secondary circulation behavior has been observed for S _p = 0.25 and 0.75 and there is another circulation between the top wall and the fin that is separated from the primary circulation. However, these secondary circulations features are not observed for Pr = 100. It is also found that the average Nusselt number decreases as the length of the fin increases for all locations. However, the rate of decrease of average Nusselt number becomes slower as the location of fin moves from the bottom wall to the top wall. The overall heat transfer rate can be controlled with a suitable selection of the fin location and length.     

Heat transfer characteristics of boundary-layer flows induced by continuous surfaces stretched with prescribed skin friction

A continuous surface stretched with velocity u _w=u _w (x) and having the temperature distribution T _w=T _w (x) interacts with the viscous fluid in which it is immersed both mechanically and thermally. The thermal interaction is characterized by the surface heat flux q _w=q _w (x) and the mechanical one by the skin friction τ _w=τ _w (x). In the whole previous theoretical research concerned with such processes, either (u _w and T _w) or (u _w and q _w) have been prescribed as known boundary conditions. The goal of the present paper is to initiate the investigation of the boundary layer flows induced by stretching processes for which either (τ _w and T _w ) or (τ _w and q _w) are the prescribed quantities. The case of an isothermal surface stretched with constant skin friction, (τ _w=const., T _w=const. ≠ T _∞) is worked out in detail. The corresponding flow and heat transfer characteristics are compared to those obtained for the (well known) case of a uniformly moving isothermal surface (u _w=const., T _w=const. ≠ T _∞).     

Influence of molecular parameters on the stress dependence of viscous and elastic properties of polypropylene melts in shear

The stress dependencies of the steady-state viscosity η and, particularly, that of the steady-state elastic compliance J _e of various linear isotactic polypropylenes (PP) and one long-chain branched PP are investigated using creep-recovery tests. The creep stresses applied range from 2 to 10,000 Pa. In order to discuss the stress-dependent viscosity η and elastic compliance J _e with respect to the influence of the weight average molar mass M _w and the polydispersity factor M _w/M _n the PP are characterized by SEC–MALLS. For the linear PP, linear steady-state elastic compliances J_e⁰J_{\rm e}^0 in the range of 10^− 5–10^− 3 Pa^− 1 are obtained depending on the molar mass distribution. J_e⁰J_{\rm e}^0 of the LCB-PP is distinctly higher and comes to lie at around 10^− 2 Pa^− 1. J_e⁰J_{\rm e}^0 is found to be independent of M _w but strongly dependent on polydispersity. η and J _e decrease with increasing stress. For the linear PP, J _e as a function of the stress τ is temperature independent. The higher M _w/M _n the stronger is the shear thinning of η and the more pronounced is the stress dependence of J _e. For the LCB-PP, the strongest stress dependence of η and J _e is observed. Furthermore, for all PP J _e reacts more sensitively to an increasing stress than η. A qualitative explanation for the stronger stress dependence of J _e compared to η is given by analyzing the contribution of long relaxation times to the viscosity and elasticity.     

Scaling law in saddle-node bifurcations for one-dimensional maps: a complex variable approach

The study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, τ, scales according to an inverse square-root power law, τ∼(μ−μ _c )^−1/2, as the bifurcation parameter μ, is driven further away from its critical value, μ _c . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.     

Linear viscoelastic behavior of perfluorooctyl sulfonate micelles: effects of counter-ions

Linear viscoelastic behavior was investigated for aqueous solutions of perfluorooctyl sulfonate (C₈F₁₇SO⁻ ₃; abbreviated as FOS) micelles having a mixture of tetraethylammonium (N⁺(C₂H₅)₄; TEA) and lithium (Li⁺) ions as the counter-ions. The solutions had the same FOS concentration (0.1 mol l⁻¹) and various Li⁺ fractions in the counter-ions, φ_Li = 0−0.6, and the FOS micelles in these solutions formed threads which further organized into dendritic networks. At T ≤ 15 °C, the terminal relaxation time τ and the viscosity η, governed by thermal scission of the networks, increased with increasing φ_Li up to 0.55. A further increase of φ_Li resulted in decreases of τ and η and in broadening of the relaxation mode distribution. These rheological changes are discussed in relation to the role of TEA ions in thermal scission: Previous NMR studies revealed that only a fraction of TEA ions were tightly bound to the FOS micellar surfaces and these bound ions stabilized the thread/network structures. The concentration of non-bound TEA ions, C_TEA ^*, decreased and finally vanished on increasing φ_Li up to φ_Li ^* ≅ 0.6, and the concentration of the bound TEA ions significantly decreased on a further increase of φ_Li. The non-bound TEA ions appeared to catalyze the thermal scission of the FOS threads, and the observed increases of τ and η for φ_Li < 0.55 were attributed to the decrease of C_TEA ^*. On the other hand, the decreases of τ and η as well as the broadening of the mode distribution, found for φ_Li > 0.55 (where C_TEA ^* ≅ 0), were related to destabilization of the FOS threads/networks due to a shortage of the bound TEA ions and to the existence of concentrated Li⁺ ions. Viscoelastic data of pure FOSTEA and FOSTEA/FOSLi/TEACl solutions lent support to these arguments for the role of TEA ions in the relaxation of FOSTEA/FOSLi solutions. Received: 12 October 1999/Accepted: 1 November 1999     

The optimal dimensions of convective-radiating circular fins

The optimal dimensions of convective-radiating circular fins with variable profile, heat-transfer coefficient and thermal conductivity, as well as internal heat generation are obtained. A profile of the form y=(w/2) [1+(r _o/r) ⁿ ] is studied, while variation of thermal conductivity is of the form k=k _o[1+ɛ((T−T _∞)/ (T _b−T _∞)) ^m ]. The heat-transfer coefficient is assumed to vary according to a power law with distance from the bore, expressed as h=K[(r−r _o)/(r _e−r _o)]^λ. The results for λ=0 to λ=1.9, and −0.4≤ɛ≤0.4, have been expressed by suitable dimensionless parameters. A correlation for the optimal dimensions of a constant and variable profile fins is presented in terms of reduced heat-transfer rate. It is found that a (quadratic) hyperbolic circular fin with n=2 gives an optimum performance. The effect of radiation on the fin performance is found to be considerable for fins operating at higher base temperatures, whereas the effect of variable thermal conductivity on the optimal dimensions is negligible for the variable profile fin. It is also observed, in general, that the optimal fin length and the optimal fin base thickness are greater when compared to constant fin thickness. Received on 22 February 1999     

Measurements of fully developed turbulent flow in a trapezoidal duct

The turbulence characteristics of fully developed isothermal air flows through a symmetric trapezoidal duct were examined experimentally using Pitot tube and hot-wire anemometry over a Reynolds number range of 3.7–11.6×10⁴. The measurements included local wall shear stress and the cross-sectional distributions of mean axial velocity, secondary velocities and Reynolds stresses. Four secondary flow cells were detected in a symmetric half of the duct. Although secondary velocity components were typically less than about 1% of the bulk axial velocity, their effect was especially pronounced on the distributions of turbulent kinetic energy and local wall shear stress.List of symbols a, b, c, d trapezoidal duct dimensions (Fig. 1) - A, B coefficients in log law (Table 1) - D _h equivalent hydraulic diameter - f Darcy friction factor, (2D _h /U _b ² ) (dP/dx) - k turbulent kinetic energy per unit mass, - k ⁺ dimensionless turbulent kinetic energy, k/( ^*)² - P static pressure - Re Reynolds number, U _b D _h / - s distance along inclined wall, measured from top corner (Fig. 1) - u, v, w fluctuating components of the velocities in the x, y, z directions - u ⁺, v ⁺, w ⁺ dimensionless turbulence intensities; u ²/ ^*, v ²/ ^*, w ²/ ^* - u ^* local friction velocity, ( _w /)^1/2 - ^* average friction velocity, (¯gt/)^1/2 - axial mean velocity (time-average) - U _b average mean axial velocity - U _sec resultant of ¯V and ¯W, (¯V ²+¯ ²)^1/2 - U ⁺ dimensionless velocity, /u ^* - ¯V, ¯W mean velocities in the y, z directions (secondary velocities) - x axial direction - y, 2 horizontal and vertical directions (Fig. 1) - z ⁺ dimensionless distance from (and normal to) a wall, zu ^*/v - distance from wall (at y=0) to location of maximum axial velocity - laminar dynamic viscosity - v kinematic viscosity - air density - _w local wall shear stress - ¯ _w average of local wall shear stresses over all walls - ¯ average wall shear stress, (dP/dx) (D _h /4) - corner angle of trapezoidal duct (Fig. 1) A version of this paper was presented at the 10th Symposium on Turbulence, University of Missouri-Rolla, Sept. 22–24, 1986