Dynamic response of micro-pillar sensors measuring fluctuating wall-shear-stress

We present in this paper test results of flexible micro-pillars and pillar arrays for wall shear stress measurements in flows with fluctuating wall shear stress such as unsteady separated flows or turbulent flows. Previous papers reported on the sensing principle and fabrication process. Static calibrations have shown this sensor to have a maximum nonlinearity of 1% over two orders of wall-shear-stress. For measurements in flows with fluctuating wall shear stress the dynamic response has been experimentally verified in an oscillating pipe flow and compared to a calculated response based on Stokes’ and Oseen’s solution for unsteady flow around a cylinder. The results demonstrate good agreement under the given boundary conditions of cylindrical micro-pillars and the limit of viscous Stokes-flow around the pillar. Depending on the fluid and pillar geometry, different response curves result ranging from a flat low-pass filtered response to a strong resonant behavior. Two different methods are developed to detect the frequency content and the directional wall shear stress information from image processing of large sensor films with arrays of micro-pillars of different geometry. Design rules are given to achieve the optimal conditions with respect to signal-to-noise ratio, sensitivity and bandwidth for measurements in turbulent flows.     

Compressibility effects due to turbulent fluctuations

This paper deals with intrinsic effects of compressibility, i.e. with dilatation fluctuations in response to pressure fluctuations. Three different types of turbulent flows are considered in more detail: homogeneous turbulent shear flow, wall-bounded turbulent shear flow and shock/turbulence interaction. A survey of the present knowledge in this field, mainly based on DNS data, is given. Using the linear inviscid perturbation equations a direct link between fluctuations of dilatation and of velocity in the direction of mean shear is presented for homogeneous shear flow. This relation might form the basis for a more universal pressure-dilatation model. It is conjectured that the insignificance of intrinsic compressibility effects in wall-bounded supersonic shear flow is mainly due to the impermeability constraint of the wall. To this end, a linear stability analysis of supersonic channel flow along cooled, but permeable walls has been performed based on Coleman et al.'s [5] mean flow data. It shows an increase in the moduli of eigenfunctions related to compressibility, like pressure, and in moduli of quantities derived from eigenfunctions such as ‘pressure dilatation’ and squared dilatation. Although these results do not prove our hypothesis they provide hints in this direction. Shock/turbulence interaction is viewed as a source of compressibility. Former DNS data of Hannappel and Friedrich [10] for shock/isotropic turbulence interaction showing the effect of compressibility on the amplification of fluctuations are interpreted based on linear perturbation equations.     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part II: influence of operational parameters

The present paper is the subsequent research of the first part (Theor Comput Fluid Dyn, 2009). It investigates the boundary film shear elastic modulus effect in a hydrodynamic contact in different operating conditions. The hydrodynamic contact is one-dimensional, composed of two parallel plane surfaces, which are respectively rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact consists of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. The hydrodynamic contact in a micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity of the film–contact interface large enough. In boundary film area, the viscosity, density, and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling which depends on the film thickness. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. It is found from the simulation results that the boundary film shear elastic modulus influences are normally negligible on the mass flow through the contact, the carried load of the contact and the overall film thickness of the contact, and the boundary film shear elastic modulus would normally influence the local film thickness in an elastic contact when the local film thickness is on the film molecule diameter scale. It is also found that the boundary film shear elastic modulus effect has the tendency of being increased with the reduction of the width of a micro contact. It is increased with the reduction of the boundary film–contact interfacial shear strength or with the increase of the critical boundary film thickness, while it is strongest at certain values of the contact surface roughness, the width ratio of fluid film area to boundary film area, and the lubricant film shear elastic modulus.

---

     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part I: theoretical analysis and typical solution

Boundary film shear elastic modulus effect is analyzed in a hydrodynamic contact. The contact is one-dimensional composed of two parallel plane surfaces, which are, respectively, rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact is consisted of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. Analysis is thus carried out for a micro Raleigh bearing in this contact. The hydrodynamic contact in this micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity large enough at the film–contact interface. In boundary film area, the viscosity, density and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling, which depend on the film thickness. The analytical approach proposed by Zhang (J Mol Liq 128:60–64, 2006) and Zhang et al. (Int J Fluid Mech Res 30:542–557, 2003) is used for boundary film area. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. Conventional hydrodynamic analysis is used for fluid film area. The present paper presents the theoretical analysis and a typical solution. It is found that for the simulated case the boundary film shear elastic modulus effects on the mass flow through the contact, the overall film thickness of the contact and the carried load of the contact are negligible but the boundary film shear elastic modulus effect on the local film thickness of the contact may be significant when the boundary film thickness is on the 1 nm scale and the contact surfaces are elastic. In Part II will be presented detailed results showing boundary film shear elastic modulus effects in different operating conditions.

---

     

Evaluation of three techniques for wall-shear measurements in three-dimensional flows

Recent improvements in three techniques for measuring skin friction in two- and three-dimensional turbulent wall-bounded shear flows are presented. The techniques are: oil-film interferometry, hot wires mounted near the wall, and surface hot-film sensors based on MEMS technology. First, we demonstrate that the oil-film interferometry technique can be used to measure the skin-friction magnitude and its direction in two- and three-dimensional wall-bounded shear flows. Second, a simple method is outlined to measure the skin friction with a wall wire located outside of the viscous sublayer. Finally, a systematic study of the parameters influencing wall-friction measurements with MEMS sensors is presented. The results demonstrate that accurate measurements of the mean skin friction with MEMS sensors are possible in two- and three-dimensional wall flows. Measurements by the three techniques are compared to each other and to past measurements in the same facility.     

A general one-equation turbulence model for free shear and wall-bounded flows

The purpose of thiswork is to introduce a complete and general one-equation model capable of correctly predicting a wide class of fundamental turbulent flows like boundary layer, wake, jet, and vortical flows. The starting point is the mature and validated two-equation k−ω turbulence model of Wilcox. The newly derived one-equation model has several advantages and yields better predictions than the Spalart-Allmaras model for jet and vortical flows while retaining the same efficiency and quality of the results for near-wall turbulent flows without using a wall distance. The derivation and validation of the new model using findings computed by the Spalart-Allmaras and the k−ω models are presented and discussed for several free shear and wall-bounded flows.     

Near-wall scaling for turbulent boundary layers with adverse pressure gradient

A new extended inner scaling is proposed for the wall layer of wall-bounded flows under the influence of both wall shear stress and streamwise pressure gradient. This scaling avoids problems of the classical wall coordinates close to flow separation and reattachment. Based on the proposed extended velocity and length scales a universal nondimensional family of velocity profiles is derived for the viscous region in the vicinity of a wall that depend on wall distance and a parameter α quantifying the importance of the streamwise pressure gradient with respect to the wall shear stress in the momentum balance. The performance of the proposed extended scaling is investigated in two different flow fields, a separating and reattaching turbulent boundary layer and a turbulent flow over a periodic arrangement of smoothly contoured hills. Both flows are results of highly resolved direct numerical simulation (DNS). The results show that the viscous assumptions are valid up to about two extended wall units. If the profiles are scaled by the extended inner coordinates, they seem to behave in a universal way. This gives rise to the hope that a universal behavior of velocity profiles can be found in the proposed extended inner coordinates even beyond the validity of the extended viscous law of the wall.      

Direct numerical simulation of low Reynolds number flows in an open‐channel with sidewalls

A direct numerical simulation of low Reynolds number turbulent flows in an open‐channel with sidewalls is presented. Mean flow and turbulence structures are described and compared with both simulated and measured data available from the literature. The simulation results show that secondary flows are generated near the walls and free surface. In particular, at the upper corner of the channel, a small vortex called inner secondary flows is simulated. The results show that the inner secondary flows, counter‐rotating to outer secondary flows away from the sidewall, increase the shear velocity near the free surface. The secondary flows observed in turbulent open‐channel flows are related to the production of Reynolds shear stress. A quadrant analysis shows that sweeps and ejections are dominant in the regions where secondary flows rush in toward the wall and eject from the wall, respectively. A conditional quadrant analysis also reveals that the production of Reynolds shear stress and the secondary flow patterns are determined by the directional tendency of the dominant coherent structures. Copyright © 2009 John Wiley & Sons, Ltd.     

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

A carbon nanotube sensor for wall shear stress measurement

A novel carbon nanotube (CNT) sensor is being developed to measure the mean and fluctuating wall shear stress (WSS) in a turbulent boundary layer. The CNT WSS sensor is based on the thermal principle and featured by high spatial and temporal resolutions (in the order of nm and kHz, respectively), low power consumption (in the order of μW), and a compact fabrication process compared with traditional WSS measurement sensors. The CNT WSS-sensing element was characterized in detail before its calibration. The CNT sensor was operated under a constant temperature (CT) operation mode and an overheat ratio range of −0.15 to −0.19 and calibrated in a fully developed turbulent channel flow. It has been observed for the first time in a macroscopic flow that the sensor output power is approximately proportional to the 1/3 powered WSS, as expected for a thermal-principle-based WSS sensor, and the wall shear stress measurement is demonstrated for a low Reynolds number flow.     

多孔介质壁面剪切湍流速度时空关联的研究

作为一个基础统计量,时空关联函数在湍流问题的研究中有着广泛的应用,是研究湍流噪声、湍流中物质扩散和大涡模拟亚格子模型等问题的重要参考.本文通过建立三维多孔结构壁面剪切湍流模型,采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对无穷大多孔介质平行板之间壁湍流进行了数值模拟,进而研究其速度脉动时空关联函数的统计特性.一方面,根据计算得到的流场数据,对比分析了常规槽道湍流与多孔介质壁面槽道湍流的时间关联函数.另一方面,计算并讨论了不同孔隙率和渗透率的多孔介质壁面对速度脉动时空关联性的影响.通过研究表明:多孔结构壁面剪切湍流的时空关联函数等值线与椭圆理论相符;在研究参数范围内,多孔介质壁面的速度时空关联系数随着孔隙率增大而增大,随着渗透率增大而减小.同时发现在槽道壁面的近壁区、过渡区、对数律区和中心区等不同位置处,速度时空关联呈现较大差异性:越远离壁面位置(对数律区和中心区),其时空关联函数所呈现的关联等值线椭圆越细长,高值相关等值线越集中.多孔介质主要改变速度时空关联椭圆图像的椭圆率,说明多孔介质壁面主要影响湍流横扫速度.     

Predictive Capabilities of an Improved Cubic k–ε Model for Inert Steady Flows

Through an improved ε transport equation, a major quality enhancement of the cubic k–ε model, earlier developed in[13], is obtained. The ε-equation of [13],yielding good results for wall-bounded and rotating flows, is combined with the one derived by Shih et al. [20], which produces good results for free shear flows (e.g. the plane jet–round jet anomaly is resolved).Results are presented for the following flows: fully developed stationary and rotating channel and pipe, backward-facing step, sudden pipe expansion, smooth channel expansion and contraction, plane and round jet. Heat transfer predictions in turbulent impinging jets are also discussed. Accurate results are obtained for the mean flow quantities for all test cases, without case dependent model tuning. This revised version was published online in July 2006 with corrections to the Cover Date.     

Large eddy simulation of compressible channel flow

Astract The present study is a contribution to the analysis of wall-bounded compressible flows, including a special focus on wall modeling for compressible turbulent boundary layer in a plane channel. large eddy simulation (LES) of fully developed isothermal channel flows at Re = 3,000 and Re = 4,880 with a sufficient mesh refinement at the wall are carried out in the Mach number range 0.3 ≤ M ≤ 3 for two different source term formulations: first the classical extension of the incompressible configuration by Coleman et al. (J. Fluid Mech. 305:159–183, 1995), second a formulation presently derived to model both streamwise pressure drop and streamwise internal energy loss in a spatially developed compressible channel flow. It is shown that the second formulation is consistent with the spatial problem and yields a much stronger cooling effect at the wall than the classical formulation. Based on the present LES data bank, compressibility and low Reynolds number effects are analysed in terms of coherent structure and statistics. A study of the universality of the structure of the turbulence in non-hypersonic compressible boundary layers (M≤5) is performed in reference to Bradshaw (Annu. Rev. Fluid. Mech. 9:33–54, 1977). An improvement of the van Driest transformation is proposed; it accounts for both density and viscosity changes in the wall layer. Consistently, a new integral wall scaling (y ^c+) which accounts for strong temperature gradients at the wall is developed for the present non-adiabatic compressible flow. The modification of the strong Reynolds analogy proposed by Huang et al. (J. Fluid Mech. 305:185–218, 1995) to model the correlation between velocity and temperature for non-adiabatic wall layers is assessed on the basis of a Crocco–Busemann relation specific to channel flow. The key role of the mixing turbulent Prandtl number Pr _m is pointed out. Results show very good agreement for both source formulations although each of them involve a very different amount of energy transfer at the wall. The present work was performed within the framework of the French–German research initiative “large eddy simulation of complex flows’ (UR 507). The computing resources were provided by IDRIS-France. The authors gratefully acknowledge the financial support from the Centre National de la Recherche Scientifique (CNRS), the Centre d’été Mathématique de Recherche Avancée en Calcul Scientifique (CEMRACS) and the Direction Générale de l’Armement (DGA/D4S).     

Hydrodynamics study in a plane ultrafiltration module using an electrochemical method and particle image velocimetry visualization

 The wall shear stress is determined at the surface of a plane sheet of Plexiglas, taking the place of a membrane, using an electrochemical method. Several microelectrodes are mounted flush to this plane plate, and maps of shear stress are determined for two inlet and outlet configurations and three channel heights. The heterogeneity of the wall shear stress is observed for both configurations. Furthermore, the study of the turbulence features of the flow shows a decreasing fluctuating rate of velocity gradient when the channel height is decreased. The wall velocity gradients and turbulent intensity rates analysis are confirmed by a flow visualization using the particle image velocimetry method. Received: 25 September 2000 / Accepted: 23 April 2001     

A vorticity stretching diagnostic for turbulent and transitional flows

Vorticity stretching in wall-bounded turbulent and transitional flows has been investigated by means of a new diagnostic measure, denoted by Γ, designed to pick up regions with large amounts of vorticity stretching. It is based on the maximum vorticity stretching component in every spatial point, thus yielding a three-dimensional scalar field. The measure was applied in four different flows with increasing complexity: (a) the near-wall cycle in an asymptotic suction boundary layer (ASBL), (b) K-type transition in a plane channel flow, (c) fully turbulent channel flow at Re _τ = 180 and (d) a complex turbulent three-dimensional separated flow. Instantaneous data show that the coherent structures associated with intense vorticity stretching in all four cases have the shape of flat ‘pancake’ structures in the vicinity of high-speed streaks, here denoted ‘h-type’ events. The other event found is of ‘l-type’, present on top of an unstable low-speed streak. These events (l-type) are further thought to be associated with the exponential growth of streamwise vorticity in the turbulent near-wall cycle. It was found that the largest occurrence of vorticity stretching in the fully turbulent wall-bounded flows is present at a wall-normal distance of y ⁺?=?6.5, i.e. in the transition between the viscous sublayer and buffer layer. The associated structures have a streamwise length of ~200–300 wall units. In K-type transition, the Γ-measure accurately locates the regions of interest, in particular the formation of high-speed streaks near the wall (h-type) and the appearance of the hairpin vortex (l-type). In the turbulent separated flow, the structures containing large amounts of vorticity stretching increase in size and magnitude in the shear layer upstream of the separation bubble but vanish in the backflow region itself. Overall, the measure proved to be useful in showing growing instabilities before they develop into structures, highlighting the mechanisms creating high shear region on a wall and showing turbulence creation associated with instantaneous separations.     

The surface hot wire as a means of measuring mean and fluctuating wall shear stress

 The paper describes some applications of a wall shear stress sensor technique which is based on hot-wire anemometry. The “surface hot wire” is a flush-mounted thermal resistive wire with a tiny slot underneath. The arrangement of this sensor guarantees an improved signal-to-noise ratio compared to a common surface hot film. The setup and the application of single sensors and of surface hot-wire arrays are shown. Some results are presented that were acquired in several experiments in the field of laminar-turbulent transition. Received: 26 May 2000/Accepted: 7 February 2001     

Simultaneous wall shear rate measurements by a three segmented electrodiffusion probe and laser-Doppler-anemometry

 In a fully developed turbulent channel flow, the instantaneous wall shear gradient was measured simultaneously by an electrodiffusion (ED) probe and by laser-Doppler anemometry (LDA). The LDA measurements were done in the viscous sublayer and the linear velocity profile was used to calculate the instantaneous shear gradient at the wall. The electrodiffusion probe requires presence of a suitable electrolyte (e.g. an iodine and potassium-iodine aqueous solution) for reduction of the species on the working electrode (cathode). This results in a concentration gradient of the ions due to the convective diffusion of the [I₃]^- to the cathode. Due to the limited molecular diffusion rate and due to the dimension of the electrodiffusion probe its frequency resolution is limited. It is shown in which limits the correction of the signal on the concentration boundary layer inertia is not necessary and the probe can be used for measuring instantaneous wall shear gradients in turbulent flows. The results are compared with those obtained by LDA and elucidate the capability of the ED probe for measuring instantaneous shear rates at low Reynolds numbers. Received: 25 October 1995/Accepted: 16 July 1996     

Turbulent Energy Scale-Budget Equations for nearly Homogeneous Sheared Turbulence

For moderate Reynolds numbers, the isotropic relation between second-order and third-order moments for velocity increments (Kolmogorov's equation) is not respected, reflecting a non-negligible correlation between the scales responsible for the injection, transfer and dissipation of the turbulent energy. For (shearless) grid turbulence, there is only one dominant large-scale phenomenon, which is the non-stationarity of statistical moments resulting from the decay of energy downstream of the grid. In this case, the extension of Kolmogorov's analysis, as carried out by Danaila, Anselmet, Zhou and Antonia, J. Fluid Mech. 391, 1999 359-369) is quite straightforward. For shear flows, several large-scale phenomena generally coexist with similar amplitudes. This is particularly the case for wall-bounded flows, where turbulent diffusion and shear effects can present comparable amplitudes. The objective of this work is to quantify, in a fully developed turbulent channel flow and far from the wall, the influence of these two effects on the scale-by-scale energy budget equation. A generalized Kolmogorov equation is derived. Relatively good agreement between the new equation and hot-wire measurements is obtained in the outer region (40 < x ⁺ ₃ < 150) of the channel flow, for which the turbulent Reynolds number is R _λ≈ 36.     

Lattice Boltzmann simulations of turbulent shear flow between parallel porous walls

The effects of two parallel porous walls are investigated, consisting of the Darcy number and the porosity of a porous medium, on the behavior of turbulent shear flows as well as skin-friction drag. The turbulent channel flow with a porous surface is directly simulated by the lattice Boltzmann method （LBM）. The Darcy-Brinkman- Forcheimer （DBF） acting force term is added in the lattice Boltzmann equation to simu- late the turbulent flow bounded by porous walls. It is found that there are two opposite trends （enhancement or reduction） for the porous medium to modify the intensities of the velocity fluctuations and the Reynolds stresses in the near wall region. The parametric study shows that flow modification depends on the Darcy number and the porosity of the porous medium. The results show that, with respect to the conventional impermeable wall, the degree of turbulence modification does not depend on any simple set of param- eters obviously. Moreover, the drag in porous wall-bounded turbulent flow decreases if the Darcy number is smaller than the order of O（10-4） and the porosity of porous walls is up to 0.4.     

Dynamic film thickness between bubbles and wall in a narrow channel

The present paper describes a novel technique to characterize the behavior of the liquid film between gas bubbles and the wall in a narrow channel. The method is based on the electrical conductance. Two liquid film sensors are installed on both opposite walls in a narrow rectangular channel. The liquid film thickness underneath the gas bubbles is recorded by the first sensor, while the void fraction information is obtained by measuring the conductance between the pair of opposite sensors. Both measurements are taken on a large two-dimensional domain and with a high speed. This makes it possible to obtain the two-dimensional distribution of the dynamic liquid film between the bubbles and the wall. In this study, this method was applied to an air–water flow ranging from bubbly to churn regimes in the narrow channel with a gap width of 1.5 mm.