Assessment of dual plane PIV measurements in wall turbulence using DNS data

Experimental dual plane particle image velocimetry (PIV) data are assessed using direct numerical simulation (DNS) data of a similar flow with the aim of studying the effect of averaging within the interrogation window. The primary reason for the use of dual plane PIV is that the entire velocity gradient tensor and hence the full vorticity vector can be obtained. One limitation of PIV is the limit on dynamic range, while DNS is typically limited by the Reynolds number of the flow. In this study, the DNS data are resolved more finely than the PIV data, and an averaging scheme is implemented on the DNS data of similar Reynolds number to compare the effects of averaging inherent to the present PIV technique. The effects of averaging on the RMS values of the velocity and vorticity are analyzed in order to estimate the percentage of turbulence intensity and enstrophy captured for a given PIV resolution in turbulent boundary layers. The focus is also to identify vortex core angle distributions, for which the two-dimensional and three-dimensional swirl strengths are used. The studies are performed in the logarithmic region of a turbulent boundary layer at z ⁺ = 110 from the wall. The dual plane PIV data are measured in a zero pressure gradient flow over a flat plate at Re _τ = 1,160, while the DNS data are extracted from a channel flow at Re _τ = 934. Representative plots at various wall-normal locations for the RMS values of velocity and vorticity indicate the attenuation of the variance with increasing filter size. Further, the effect of averaging on the vortex core angle statistics is negligible when compared with the raw DNS data. These results indicate that the present PIV technique is an accurate and reliable method for the purposes of statistical analysis and identification of vortex structures.     

Asymptotic Analysis for a Mixed Boundary-Value Contact Problem

The variational solution of the nonlinear Signorini contact problem determines also the active contact zone Γ ^c . If the latter is known, then the elastic field is a solution of a linear mixed boundary value problem in which on Γ ^c the normal displacement and tangential traction are given, while on the non-contact part the total traction is zero. Such mixed boundary conditions in general generate singularities of the solution's stress field at the points P ^{( k )} where the boundary conditions change. For smooth data, however, the variational solution of the Signorini contact problem actually belongs to H ²(Ω)², which implies the disappearance of these singularities, i.e., that the corresponding stress intensity factors vanish. This paper is devoted to the characterization of the active contact zone Γ ^c by the vanishing stress intensity factors including their sensitivity with respect to varying Γ ^c for two-dimensional problems provided that Γ ^c consists of a finite number of intervals. We use the method of asymptotic expansions and derive an explicit formula for the sensitivity, which is rigorously justified by employing weighted Sobolev spaces with detached asymptotics. These results can be used to determine the points P ^{( k )} with a corresponding Newton iteration. Accepted July 6, 2000?Published online January 22, 2001     

Near-wall flow in a three-dimensional boundary layer on the endwall of a 30° bend

 Turbulence measurements are reported on the three-dimensional turbulent boundary layer along the centerline of the flat endwall in a 30° bend. Profiles of mean velocities and Reynolds stresses were obtained down to y ⁺≈2 for the mean flow and y ⁺≈8 for the turbulent stresses. Mean velocity data collapsed well on a simple law-of-the-wall based on the magnitude of the resultant velocity. The turbulence intensity and turbulent shear stress magnitude both increased with increased three-dimensionality. The ratio of these two quantities, the a ₁ structure parameter, decreased in the central regions of the boundary layer and showed profile similarity for y ⁺<50. The shear stress vector angle lagged behind the velocity gradient vector angle in the outer region of the boundary layer, however there was an indication that the shear stress vector tends to lead the velocity gradient vector close to the wall. Received: 16 July 1996/Accepted: 14 July 1997     

Steady mixed convection boundary layer flow over a vertical flat plate in a porous medium filled with water at 4°C: case of variable wall temperature

The problem of steady mixed convection boundary layer flow over a vertical impermeable flat plate in a porous medium saturated with water at 4°C (maximum density) when the temperature of the plate varies as x ^m and the velocity outside boundary layer varies as x ^{2 m} , where x measures the distance from the leading edge of the plate and m is a constant is studied. Both cases of the assisting and the opposing flows are considered. The plate is aligned parallel to a free stream velocity U _∞ oriented in the upward or downward direction, while the ambient temperature is T _∞ = T _m (temperature at maximum density). The mathematical models for this problem are formulated, analyzed and simplified, and further transformed into non-dimensional form using non-dimensional variables. Next, the system of governing partial differential equations is transformed into a system of ordinary differential equations using the similarity variables. The resulting system of ordinary differential equations is solved numerically using a finite-difference method known as the Keller-box scheme. Numerical results for the non-dimensional skin friction or shear stress, wall heat transfer, as well as the temperature profiles are obtained and discussed for different values of the mixed convection parameter λ and the power index m. All the numerical solutions are presented in the form of tables and figures. The results show that solutions are possible for large values of λ and m for the case of assisting flow. Dual solutions occurred for the case of opposing flow with limited admissible values of λ and m. In addition, separation of boundary layers occurred for opposing flow, and separation is delayed for the case of water at 4°C (maximum density) compared to water at normal temperature.     

A continuum mechanical theory for turbulence: a generalized Navier–Stokes-<Emphasis Type="Italic">α</Emphasis> equation with boundary conditions

We develop a continuum-mechanical formulation and generalization of the Navier–Stokes-α equation based on a recently developed framework for fluid-dynamical theories involving higher-order gradient dependencies. Our flow equation involves two length scales α and β. The first of these enters the theory through the specific free-energy α ²|D|², where D is the symmetric part of the gradient of the filtered velocity, and contributes a dispersive term to the flow equation. The remaining scale is associated with a dissipative hyperstress which depends linearly on the gradient of the filtered vorticity and which contributes a viscous term, with coefficient proportional to β ², to the flow equation. In contrast to Lagrangian averaging, our formulation delivers boundary conditions and a complete structure based on thermodynamics applied to an isothermal system. For a fixed surface without slip, the standard no-slip condition is augmented by a wall-eddy condition involving another length scale ℓ characteristic of eddies shed at the boundary and referred to as the wall-eddy length. As an application, we consider the classical problem of turbulent flow in a plane, rectangular channel of gap 2h with fixed, impermeable, slip-free walls and make comparisons with results obtained from direct numerical simulations. We find that α/β ~ Re ^0.470 and ℓ/h ~ Re ^−0.772, where Re is the Reynolds number. The first result, which arises as a consequence of identifying the specific free-energy with the specific turbulent kinetic energy, indicates that the choice β = α required to reduce our flow equation to the Navier–Stokes-α equation is likely to be problematic. The second result evinces the classical scaling relation η/L ~ Re ^−3/4 for the ratio of the Kolmogorov microscale η to the integral length scale L.      

Variation of the stress intensity factor along the front of a 3-D rectangular crack subjected to mixed-mode load

Summary  The singular integral equation method is applied to the calculation of the stress intensity factor at the front of a rectangular crack subjected to mixed-mode load. The stress field induced by a body force doublet is used as a fundamental solution. The problem is formulated as a system of integral equations with r ⁻³-singularities. In solving the integral equations, unknown functions of body-force densities are approximated by the product of polynomial and fundamental densities. The fundamental densities are chosen to express two-dimensional cracks in an infinite body for the limiting cases of the aspect ratio of the rectangle. The present method yields rapidly converging numerical results and satisfies boundary conditions all over the crack boundary. A smooth distribution of the stress intensity factor along the crack front is presented for various crack shapes and different Poisson's ratio. Received 5 March 2002; accepted for publication 2 July 2002     

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     

Free convection from a vertical permeable circular cone with non-uniform surface heat flux

 In this paper, the problem of laminar free convection from a vertical permeable circular cone maintained with non-uniform surface heat flux is considered. The governing boundary layer equations are reduced non-similar boundary layer equations with surface heat flux proportional to x ⁿ (where x is the distance measured from the leading edge). The solutions of the reduced equations are obtained by using three distinct solution methodologies; namely, (i) perturbation solution for small transpiration parameter, ξ, (ii) asymptotic solution for large ξ, and (iii) the finite difference solutions for all ξ. The solutions are presented in terms of local skin-friction and local Nusselt number for smaller values of Prandtl number and heat flux gradient and are displayed in tabular form as well as graphically. Effects of pertinent parameters on velocity and temperature profiles are also shown graphically. Solutions obtained by finite difference method are also compared with the perturbation solutions for small and large ξ and found to be in excellent agreement. Received on 1 October 1999     

Rheology of lyotropic and thermotropic lamellar phases

We investigate the steady-state rheological behaviour of the lamellar phase of a lyotropic system (CpCl, hexanol, brine) and of a thermotropic system (8CB). Power laws characterize the behaviour of the imposed stress as a function of the measured shear rate and similarities are observed for both systems; the same regime γ˙∼σ ^m with m≈1.7 is obtained at low shear stresses corresponding to a texture of oily streaks oriented in the direction of the flow, as shown by microscopic observations. The “onion state” only exists in the case of dilute samples of the lyotropic lamellar phase; the stress then varies as γ˙∼σ ^m with m≈4.8, as already observed by other groups with different systems. Rheological measurements at different temperatures allow determination of different activation energies relating to the still badly understood processes involved in the different rheological regimes. We propose a model which reproduces the experimental power laws and which is based on an analogy with the theory of high-temperature creep in metals and alloys. Received: 19 October 1999/Accepted: 1 November 1999     

Hydromagnetic free convection flow with induced magnetic field effects

An exact solution is presented for the hydromagnetic natural convection boundary layer flow past an infinite vertical flat plate under the influence of a transverse magnetic field with magnetic induction effects included. The transformed ordinary differential equations are solved exactly, under physically appropriate boundary conditions. Closed-form expressions are obtained for the non-dimensional velocity (u), non-dimensional induced magnetic field component (B _x ) and wall frictional shearing stress i.e. skin friction function (τ _x ) as functions of dimensionless transverse coordinate (η), Grashof free convection number (G _r ) and the Hartmann number (M). The bulk temperature in the boundary layer (Θ) is also evaluated and shown to be purely a function of M. The Rayleigh flow distribution (R) is derived and found to be a function of both Hartmann number (M) and the buoyant diffusivity parameter (ϑ ^*). The influence of Grashof number on velocity, induced magnetic field and wall shear stress profiles is computed. The response of Rayleigh flow distribution to Grashof numbers ranging from 2 to 200 is also discussed as is the influence of Hartmann number on the bulk temperature. Rayleigh flow is demonstrated to become stable with respect to the width of the boundary layer region and intensifies with greater magnetic field i.e. larger Hartman number M, for constant buoyant diffusivity parameter ϑ ^*. The induced magnetic field (B _x ), is elevated in the vicinity of the plate surface with a rise in free convection (buoyancy) parameter G _r , but is reduced over the central zone of the boundary layer regime. Applications of the study include laminar magneto-aerodynamics, materials processing and MHD propulsion thermo-fluid dynamics.     

Unsteady free convection on a vertical cylinder with variable heat and mass flux

The unsteady natural convection boundary layer flow over a semi-infinite vertical cylinder is considered with combined buoyancy force effects, for the situation in which the surface temperature T ^′ _w(x) and C ^′ _w(x) are subjected to the power-law surface heat and mass flux as K(T ^′/r) = −ax ⁿ and D(C ^′/r) = −bx ^m . The governing equations are solved by an implicit finite difference scheme of Crank-Nicolson method. Numerical results are obtained for different values of Prandtl number, Schmidt number ‘n’ and ‘m’. The velocity, temperature and concentration profiles, local and average skin-friction, Nusselt and Sherwood numbers are shown graphically. The local Nusselt and Sherwood number of the present study are compared with the available result and a good agreement is found to exist. Received on 7 July 1998     

Oscillatory Couette flow in the presence of an inclined magnetic field

Oscillatory MHD Couette flow of electrically conducting fluid between two parallel plates in a rotating system in the presence of an inclined magnetic field is considered when the upper plate is held at rest and the lower plate oscillates non-torsionally . An exact solution of the governing equations has been obtained by using Laplace transform technique. Asymptotic behavior of the solution is analyzed for M ² ≪1, K ² ≪1 and ω ≪1 and for large M ², K ² and ω. Numerical results of velocities are depicted graphically and the frictional shearing stresses are presented in tables. It is found that a thin boundary layer is formed near the lower plate, for large values of rotation parameter K ², Hartman number M ² and frequency parameter ω. The thickness of this boundary layer increases with increase in inclination of the magnetic field with the axis of rotation.     

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).     

Natural Convection in a Cavity Filled with Porous Layers on the Top and Bottom Walls

Natural convection in a partially filled porous square cavity is numerically investigated using SIMPLEC method. The Brinkman-Forchheimer extended model was used to govern the flow in the porous medium region. At the porous-fluid interface, the flow boundary condition imposed is a shear stress jump, which includes both the viscous and inertial effects, together with a continuity of normal stress. The thermal boundary condition is continuity of temperature and heat flux. The results are presented with flow configurations and isotherms, local and average Nusselt number along the cold wall for different Darcy numbers from 10⁻¹ to 10⁻⁶, porosity values from 0.2 to 0.8, Rayleigh numbers from 10³ to 10⁷, and the ratio of porous layer thickness to cavity height from 0 to 0.50. The flow pattern inside the cavity is affected with these parameters and hence the local and global heat transfer. A modified Darcy–Rayleigh number is proposed for the heat convection intensity in porous/fluid filled domains. When its value is less than unit, global heat transfer keeps unchanged. The interfacial stress jump coefficients β ₁ and β ₂ were varied from  −1 to +1, and their effects on the local and average Nusselt numbers, velocity and temperature profiles in the mid-width of the cavity are investigated.     

On corner flows of Oldroyd-B fluids

Exact series solutions for planar creeping flows of Oldroyd-B fluids in the neighbourhood of sharp corners are presented and discussed. Both reentrant and non-reentrant sectors are considered. For reentrant sectors it is shown that more than one type of series solution can exist formally, one type exhibiting Newtonian-like asymptotic behaviour at the corner, away from walls, and another type exhibiting the same kind of asymptotics as an Upper Convected Maxwell (UCM) fluid. The solutions which are Newtonian-like away from walls are shown to develop non-integrable stress singularities at the walls when the no-slip velocity boundary condition is imposed. These mathematical solutions are therefore inadmissible from the physical viewpoint under no-slip conditions. An inadmissible solution, with stress singularities which are not everywhere integrable, is identified among the solutions of UCM-type. For a 270° reentrant sector the radial behaviour of the normal stress is everywhere r^−0.613. In the viscometric region near a wall, the radial normal stress σ^rr behaves like (rε)^−0.613, where ε is the angle made with the wall. In addition σ^rθ is infinite (not integrable) at the wall even when r is non-zero. Another UCM-type solution has a normal stress behaviour away from walls which is r^−0.985 for 270° sector. Again, this solution has a non-integrable stress singularity and is therefore inadmissible. Finally, for non-reentrant sectors it is shown that the flow is always Newtonian-like away from walls.     

A new approach on MHD natural convection boundary layer flow past a flat plate of finite dimensions

A new approach on MHD natural convection boundary layer flow from a finite flat plate of arbitrary inclination in a rotating environment, is presented. This problem plays a significant role on boundary layer flow control. It is shown that taking into account the pressure rise region at the leading edge of the plate leads to avoid separation and the back flow is reduced by the strong magnetic field. It is also shown that the frictional drag at the leading edge of the plate is reduced when the inclination angle α=π/4. In the case of isothermal flat plate, the bulk temperature becomes identical for any value of Gr (Grashof number) when the value of M ² (Hartmann number) and K ² (rotation parameter) are kept fixed.     

Exact solutions for the flow of second grade fluid in annulus between torsionally oscillating cylinders

The velocity field and the associated shear stress corresponding to the torsional oscillatory flow of a second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. At time t = 0, the fluid and both the cylinders are at rest and at t = 0 + , cylinders suddenly begin to oscillate around their common axis in a simple harmonic way having angular frequencies ω 1 and ω 2 . The obtained solutions satisfy the governing differential equation and all imposed initial and boundary conditions. The solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for Newtonian fluid are also obtained as limiting cases of our general solutions.     

Tolerance averaging of heat conduction in transversally graded laminates

It is considered a heat conduction in a layer made of two conductors distributed in the form of laminas with varied thicknesses. Macroscopic (averaged) properties of the layer are continuously “transversally” graded across its thickness (TGL layer), cf. Fig. 1. The aim of the paper is to present and apply an averaged model of the heat conduction, obtained within the tolerance averaging technique, discussed in the book edited by Woźniak et al. (Thermomechanics of microheterogeneous solids and structures. Tolerance averaging approach, Łódź, Wydawnictwo Politechniki Łódzkiej, 2008). It is shown that the proposed model describes the microstructural effect on the heat conduction of the TGL layer. Moreover, results obtained within this model are compared to results by the higher order theory (cf. Aboudi et al., Composites B, 30:777–832, 1999).     

Some Reynolds number effects on two- and three-dimensional turbulent boundary layers

Experimental data for a two-dimensional (2-D) turbulent boundary layer (TBL) flow and a three-dimensional (3-D) pressure-driven TBL flow outside of a wing/body junction were obtained for an approach Reynolds number based on momentum thickness of Re _θ =23,200. The wing shape had a 3:2 elliptical nose, NACA 0020 profiled tail, and was mounted on a flat wall. Some Reynolds number effects are examined using fine spatial resolution (Δy ⁺=1.8) three-velocity-component laser-Doppler velocimeter measurements of mean velocities and Reynolds stresses at nine stations for Re _θ =23,200 and previously reported data for a much thinner boundary layer at Re _θ =5,940 for the same wing shape. In the 3-D boundary layers, while the stress profiles vary considerably along the flow due to deceleration, acceleration, and skewing, profiles of the parameter correlate well and over available Reynolds numbers. The measured static pressure variations on the flat wall are similar for the two Reynolds numbers, so the vorticity flux and the measured mean velocities scaled on wall variables agree closely near the wall. The stresses vary similarly for both cases, but with higher values in the outer region of the higher Re _θ case. The outer layer turbulence in the thicker high Reynolds number case behaves similarly to a rapid distortion of the flow, since stream-wise vortical effects from the wall have not diffused completely through the boundary layer at all measurement stations. Received: 9 June 2000/Accepted: 26 January 2001     

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 × 10³ to 6.5 × 10⁵ to study the effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics. In the subcritical range of Re (<2 × 10⁵), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 × 10⁵), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag. Received: 8 September 1998 / Accepted: 1 January 2000