Receptivity of plane Poiseuille flow to local micro-vibration disturbance on wall

* Corresponding author.



Receptivity of plane Poiseuille flow to local micro-vibration disturbance on wall

Wei-dong Cao*,Qi-xiang Hu,Bing Liu

Research Center of Fluid Machinery Engineering and Technology,Jiangsu University,Zhenjiang 212013,PR China

Received 16 June 2014; accepted 23 February 2015 Available online 26 April 2015

* Corresponding author.

Abstract

The receptivity of plane Poiseuille flow to local single-period micro-vibration disturbances with different phases at the top and bottom walls was investigated through direct numerical simulation of three-dimensional incompressible Navier-Stokes equations.Results show that the disturbance presents a symmetrical distribution in the spanwise direction when the micro-vibration on the wall ends,and the initial disturbance velocities and spatial distribution of the disturbance structure are different at the top and bottom walls.The disturbance's velocity,amplitude,and high- and low-speed streaks increase with time,and the amplitude of streamwise disturbance velocity is larger than those of spanwise and vertical disturbance velocities.However,no significant Tollmien-Schlichting wave was found in the flow field.The number of disturbance vortex cores gradually increases with the disturbance area.High-speed disturbance fluid concentrates near the wall and its normal velocity largely points to the wall,while low-speed disturbance fluid largely deviates from the wall.Furthermore,the streamwise velocity profiles near the top and bottom walls both become plump because of the existence of the disturbances,and the streamwise velocity profiles show a trend of evolving into turbulent velocity profiles.The shear stress near the wall increases significantly.The local micro-vibration disturbance on the wall in plane Poiseuille flow can induce the development of a structure similar to turbulent spots.©2015 Hohai University.Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Keywords:Poiseuille flow; Disturbance; Velocity profile; Vortex; High- and low-speed streaks

This work was supported by the National Natural Science Foundation of China (Grant No.51179075)and the Natural Science Foundation of Jiangsu Province (Grant No.BK20131256).

E-mail address: cwd@ujs.edu.cn (Wei-dong Cao).Peer review under responsibility of Hohai University.

http://dx.doi.org/10.1016/j.wse.2015.04.006

1674-2370/©2015 Hohai University.Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

1.Introduction

Due to the significant influence of laminar-turbulent transition on heat transfer,mass transfer,momentum transfer,and wall friction near the wall,a large number of scholars have focused their research on the laminar-turbulent transition and control of the boundary layer in recent years.The three processes of transition include receptivity,linear disturbance evolution,and nonlinear evolution.Transition is a complex process relating to the Reynolds number,wall temperature,the wall shape,wall roughness,fluid compressibility,the pressure gradient,external noise,and external disturbance.Various flow structures are found under various conditions,especially in channel flow,pipes,and other simple shear flows.Transition at the initial stage is usually characterized by distortions in the local flow field caused by local micro-disturbance,and unstable high- and low-speed strip areas emerging near the wall,leading to complex vortices with gradually increasing vortex intensity downstream.The dynamic characteristics of disturbance are similar to those of coherent structures in fully developed turbulent boundary layers,whereas they are significantly different from the classical transition process,which is characterized by Tollmien-Schlichting waves showing linear growth,nonlinear instability,secondary instability,and,finally,a three-dimensional nonlinear effect.

Ellingsen and Palm (1975)first put forward a probable growth mechanism in the sense of a non-natural transition.Asmicro-disturbance of spanwise vortices occurs in the shear layer,the streamwise disturbance velocity grows linearly with time,and high- and low-speed streaks appear in non-viscous ideal fluid.Landahl (1980)also obtained a similar conclusion that constantly lifting fluid particles with horizontal momentum in shear flow cause the disturbance of stream velocity,which means that there is a transient growth mechanism where non-viscous growthcoexists withviscous damping.Blackwelder (1983)found that vortex dynamics of reverse streamwise vortices in the turbulent boundary layer and in the laminar-turbulent transition zone of the boundary layer are similar and their scales are also similarly represented by the viscous length.Acarlar and Smith (1987)found that artificial disturbance near the wall of the laminar boundary layer led to lateral vortex line deformation and the disturbance developed into a hairpin vortex.Henningson et al.(1987)investigated the characteristics of vortex structures in Poiseuille flow and the boundary layer through numerical methods by designing an initial local disturbance near the wall.Testing results of the hot-wire anemometer proved that the vortex structures exist as one kind of multi-eddy structure.Haidari and Smith (1994)examined the generation and growth of single hairpin vortices created by controlled surface fluid injection within a laminar boundary layer over a range of the Reynolds number.Rosenfeld et al.(1999)proposed a general model to describe the evolution of local three-dimensional disturbance,which showed that the size of disturbance was much smaller than the characteristic length of the external shear flow.Using this model,the spanwise vortex was pulled to the outer region due to the effect of the disturbance jet,and the disturbance extended in the normal direction,while it was also stretched by the shear flow of the outer region and re-rotated into the inner region; these processes promoted the growth and enlargement of streamwise vortices,which directly caused the increase of the normal velocity as well as the intensity of upper jet-flow (Rosenfeld et al.,1999).Andersson et al.(1999)found that,owing to the effect of viscous dissipation,there was considerable linear growth in three-dimensional disturbance before attenuation.Streamwise vortices can induce instantaneous maximum growth in space in a non-parallel flat plate boundary layer.If the amplitude of streaks reaches a sufficiently large value,secondary instability will occur and induce the primeval breakdown and transition (Andersson et al.,1999).Li (2001)experimentally investigated the generation mechanism of streamwise vortices in the transition region and examined the physical process,caused by the axial vortex instability induced by the interaction of the L-vortex and secondary vortex rings.Svizher and Cohen (2001)used a continuous injection to generate hairpin vortices in subcritical plane Poiseuille flow.Zhang and Tang (2006)simulated the generation and development of turbulent spots with the local streamwise velocity pulse of fluid as the initial disturbance near the wall in channel flow,and analyzed the characteristics of nonlinear evolution of turbulent spots.Lu et al.(2008)researched the evolutionary mechanisms and characteristics of the vortex structure stimulated by the local constant wall velocity pulse during a period of time in the boundary layer with the pressure gradient.Lee and Wu (2008)presented direct comparisons of experimental results of transition in wallbounded flows obtained by flow visualizations,hot-film measurement,and particle-image velocimetry,along with a brief mention of relevant theoretical progresses,based on a critical review of about 120 selected publications.Despite somewhat different initial disturbance conditions used in experiments,the flow structures were found to be practically the same.

Although there has been some progress,through experiments and numerical simulations,improving our understanding of the temporal and spatial evolution characteristics of disturbance in laminar shear flow,further research is needed to investigate the effect and mechanisms of basic flow deformation on the disturbance growth,as well as the receptivity problems induced by initial disturbance differences at the walls in simple shear flow.Currently,wall disturbance is a common method of investigating receptivity problems such as the evolution of disturbance and mean flow profile changes in Poiseuille flow in experiments and numerical calculations.In this study,local single-period micro-vibrations with opposite phases were applied to the top and bottom walls in plane Poiseuille flow.

2.Numerical methods

2.1.Governing equations and numerical method

The governing equations adopted were the incompressible,non-dimensional Navier-Stokes equations,and continuity equation: where Re is the Reynolds number; U is the numerical solution of Poiseuille basic flow,andwhere u0,v0,and w0are the streamwise,vertical,and spanwise velocities,respectively; u is the three-dimensional disturbance velocity vector,andwhere u0,v0,and w0are the streamwise disturbance velocity,vertical disturbance velocity,and spanwise disturbance velocity,respectively; and p0is the three-dimensional disturbance pressure.In this paperhere U∞is twice the maximum velocity of Poiseuille basic flow,h is 0.5 times the width of the two-dimensional channel,and y is the kinematic viscosity.The procedures of direct numerical simulation of Eqs.(1)and (2)were as follows: the third-order mixed explicit-implicit scheme was used for time discretization,the Fourier spectral expansion was used in the spanwise direction,the fifth-order upwind compact finite difference scheme was used for the nonlinear terms,the five-point central non-equidistant difference scheme was used for the Helmholtz equations,the fifthorder symmetrical compact finite difference scheme was used for the viscous terms,and the time step was 0.01.The numerical method is described in detail in Lu et al.(2006).

2.2.Computation domain and boundary conditions

Owing to the limitation of the computing capacity,the computational domain wasandin the streamwise,vertical,and spanwise directions,respectively.The numbers of mesh points in the x and z directions were 400 and 32,respectively,with uniform distribution.The number of mesh points in the y direction was 80,with non-uniformly distributed mesh points near the wall deliberately refined.The node coordinate y(j)in the y direction can be expressed as

3.Results and discussion

For Poiseuille basic flow,we have u0¼u0(y)¼0.5(1y)2,v0¼0,and w0¼0.Thus,the governing equations can be simplified as follows:

Inflow boundary conditions were x¼0,u¼0,and vp0/vx¼0.Outflow boundary conditions were x¼80,the non-reflecting boundary condition,and vp0/vx¼0.Boundary conditions at the top wall were y¼1,vu/vy¼0,and p0¼0.Boundary condition at the bottom wall were y¼1,vu/vy¼0,and p0¼0.

To perform approximate simulations of the wall micro-jet or local forced vibration,initial disturbance was set as a local single-period micro-vibration on the walls at y¼1 and y¼1.The mesh deformation was not considered in this study.The vertical disturbance velocity at mesh points in the circle region withis supposed to bewhere r is the distance from the mesh point to the point (10,1,0)at the top wall or the point (10,1,0)at the bottom wall,T is the vibration period,t is the time of vibration,T¼20,and t < 20.The amplitude of the small disturbance velocity is commonly set to be about one percent of the maximum mean basic flow velocity,and a relatively small amplitude of the disturbance velocity can lead to a relatively long period of evolution of the disturbance.Thus,the amplitude of the disturbance velocity was tested and set at 0.007 5.As the interior between the top and bottom walls is the flow region,the actual disturbance phase at the top wall was opposite to that at the bottom wall in this study.The computational domain is shown in Fig.1,where points (10,1,0)and (10,1,0)are the centre points of initially circular disturbance.

The disturbance amplitude is defined as

qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

The disturbance amplitude and the maximum absolute values of three-dimensional disturbance velocities are shown in Fig.2.

It can be seen from Fig.2 that when t < 10,the disturbance amplitude A gradually increases to 0.02 and approaches the value ofWhen 10 < t < 20,the disturbance amplitude A remains at 0.02.As the disturbance in Poiseuille flow evolves and develops constantly,the disturbance amplitude A shows a linear growth trend whendemonstrates a same growth trend,whilestill remain at a lower level.However,when the disturbance amplitude is greater than 0.06 and 40 < t < 90,the disturbance amplitude A shows a nonlinear growth trend and even reaches the maximum value of Poiseuille basic flowalso increase rapidly,to over 0.1,particularly when t > 80,presents an accelerating growth trend.It can be prehat such a large disturbance velocity will cause strong shear layers and complex vortices,and the velocity profiles of original Poiseuille basic flow will be greatly modified.

Fig.1.Computation domain.

Simulation results show that the absolute values of streamwise and vertical disturbance velocities on the plane z¼0 are at their maximum.Fig.3 shows the distribution of the streamwise and vertical disturbance velocities on the plane z¼0.When t¼20,there is an obvious clockwise vortex core in the disturbance flow field at the point (x¼14.2,y¼0),which is defined as the original vortex core hereafter,on the left and right sides of the original vortex core,there is arelatively obvious vortex core,and the main range of the disturbance is 10 < x < 18.When t¼40,the original vortex core moves to the point (x¼18.4,y¼0),two vortex cores occur on the right side of the original one,while there is still only one vortex core on its left side,and the main range of the disturbance is 14 < x < 26.When t¼60,the original vortex core moves to the point (x¼22.6,y¼0),three vortex cores are generated on the right side of the original one,while there is still only one vortex core on its left side,and the main range of disturbance is 17 < x < 33.When t¼90,the original vortex core moves to the point (x¼29.4,y¼0),five vortex cores are generated on the right side of the original one,and two are generated on its left side,and the main range of disturbance is 22 < x < 45.The disturbance vortex core's moving velocity is about 0.45 times the maximum velocity of Poiseuille basic flow.Singer (1996)concluded that the streamwise velocity of turbulent spots at the initial stage in the boundary layer was 0.5 times the maximum velocity of basic flow.When 20 < t < 60,the expansion speed of the main disturbance area is about 0.4 times the maximum velocity of Poiseuille basic flow,almost equal to the streamwise expansion speed of turbulent spots in the plane Couette flow simulated by Lundbladh and Johansson (1991).When 60 < t < 90,the disturbance is in the rapidly growing stage,and the mean expansion speed of the main disturbance area seems to show an accelerating trend,reaching 0.46 times the maximum velocity of Poiseuille basic flow.

Fig.2.Disturbance amplitude and maximum absolute values of threedimensional disturbance velocities.

As can be seen from Fig.3(a),the induced disturbance flow field is generally anti-symmetrical in the y direction at t¼20.However,there are exceptions,e.g.,the amplitude of disturbance velocity at the point (x¼12.8,y¼0.64)is different from that at the point (x¼12.8,y¼0.64),as shown in the block diagram in Fig.3(a).

Fig.3.Streamwise and vertical disturbance velocities for plane z¼0 at different times.

High-speed fluid (u0> 0)and low-speed fluid (u0< 0)are largely concentrated in the region near y¼1 or y¼1 as the disturbance moves downstream with time.Most high-speed fluid sweeps downwards along the wall,while low-speed fluid deviates from the wall upwards.These characteristics are similar to the bursting characteristics of coherent structures in classical turbulence examined by Kline et al.(1967).The disturbance flow field becomes complex at t¼60 and 90: it isno longer anti-symmetrical in the y direction,high-intensity vortices are distributed densely,and high- and low-speed streaks are distributed in a staggered manner.At t¼90,a large amount of high-speed fluid (u0> 0)is concentrated near the walls at y¼±1.This situation rapidly changes the shear stress.

Disturbance vortices play a vital role in energy generation and transportation as well as the formation of the laminarturbulent transition.Fig.4 shows the iso-surfaces with the disturbance vortex amplitude greater than 1.There are different scales of disturbance vortices in the flow field,mainly distributed in the region near the plane y¼1 or y¼1,and symmetrically distributed with respect to the plane z¼0.However,the area size,shape,and amplitude of disturbance vortices near the planes y¼1 and y¼1 are different,due to different levels of receptivity of the flow field to different initial disturbances.The distributions of vortices at t¼60 are focused on the top and bottom walls and the vortex area at the bottom wall is greater than that at the top wall.When t¼90,large vortices have reached the central area.The vertical velocity is relatively large at y¼0,as shown in Fig.3,but the vortices' amplitude is not large (Fig.4).Micro-complex flow structure is mainly concentrated near the wall,and the nearwall region is the most active area,where the flow instability,transition,and turbulence occur.

Fig.5 shows the distribution of streamwise disturbance velocity on the planes y¼0.1 and y¼0.1.Due to the periodic boundary conditions in the spanwise direction,the streamwise disturbance velocity distribution is symmetrical with respect to z¼0.The streamwise disturbance velocity generally shows a certain fluctuation in the flow direction,and the area with positive streamwise disturbance velocity is larger than the area with negative streamwise disturbance velocity.The disturbance velocity amplitude changes greatly and there are many local deformations.The streamwise disturbance velocity at x¼29 is a trough on the plane y¼0.1,and the distances between peaks are not equal.These facts indicate that Tollmien-Schlichting wave characteristics are nonexistent here,which is consistent with the direct numerical simulation results from Singer (1996).

Fig.4.Iso-surfaces with disturbance vortex amplitude greater than 1.0.

Fig.5.Streamwise disturbance velocity on different planes.

The solid lines in Fig.6 represent the velocity profile of the original Poiseuille basic flow.The dashed lines and dotted lines represent the velocity profiles of the Poiseuille basic flow added by the local mean value of the streamwise disturbance velocity.As can be seen from Fig.6,when t¼60,the three velocity profiles are similar to one another,and in fact there is almost no difference between them; when t¼90,the difference between the three velocity profiles can be easily discerned.On the one hand,because of the presence of disturbance velocity,the shear stress closest to the wall significantly increases.On the other hand,under the conditions of different initial disturbances at the top and bottom walls in Poiseuille basic flow,the wall mean velocity profiles all become plump after a period of evolution.Although the disturbance is only in its initial stage,the streamwise velocity profile has a tendency to evolve into the turbulent mean velocity profile.Furthermore,the velocity profile for 25 < x < 40 is plumper than the other two velocity profiles and the average shear stress near the top and bottom walls for 25 < x < 40 is greater than that in Poiseuille basic flow.

4.Conclusions

The receptivity of plane Poiseuille flow to local singleperiod micro-vibrations with different phases at the top and bottom walls was investigated with direct numerical simulation.The main conclusions are as follows:

Fig.6.Profiles of streamwise velocity at different times.

(1)Although small-amplitude disturbance waves in Poiseuille flow at Re¼5 000 are attenuated based on linear stability theory,the disturbance growth mechanism described in this paper is similar to the bursting characteristics of coherent structures in classical turbulence.

(2)High- and low-speed streaks,vortices,propagation velocity and other characteristics of the disturbance structures are similar to the characteristics of turbulent spots and coherent structures in the turbulent boundary layer.

(3)The flow field of plane Poiseuille flow induced by local single-period micro-vibrations with different phases at the top and bottom walls makes the streamwise velocity profiles plump at the top and bottom walls.The streamwise velocity profiles even show a trend of evolving into a turbulent mean velocity profile.

References

Acarlar,M.S.,Smith,C.R.,1987.A study of hairpin vortices in a laminar boundary layer,Part 1: Hairpin vortices generated by a hemisphere protuberance.J.Fluid Mech.175(1),1e41.http://dx.doi.org/10.1017/S0022112087000272.

Andersson,P.,Berggren,M.,Henningson,D.S.,1999.Optimal disturbances and bypass transition in boundary layers.Phys.Fluids 11(1),134e151.

Blackwelder,R.F.,1983.Analogies between transitional and turbulent boundary layers.Phys.Fluids 26(8),2807e2815.http://dx.doi.org/10.1063/1.864047.

Ellingsen,T.,Palm,E.,1975.Stability of linear flow.Phys.Fluids 18(1),487e502.http://dx.doi.org/10.1063/1.861156.

Haidari,A.H.,Smith,C.R.,1994.The generation and regeneration of single hairpin vortices.J.Fluid Mech.277(1),135e162.http://dx.doi.org/10.1017/S0022112094002715.

Henningson,D.,Spalart,P.,Kim,J.,1987.Numerical simulations of turbulent spots in plane Poiseuille and boundary-layer flow.Phys.Fluids 30(10),2914e2918.http://dx.doi.org/10.1063/1.866067.

Kline,S.L.,Reynolds,W.C.,Schraub,F.A.,Runstadler,P.W.,1967.The structure of turbulent boundary layers.J.Fluid Mech.30(4),741e773.http://dx.doi.org/10.1017/S0022112067001740.

Landahl,M.T.,1980.A note on an algebraic instability of inviscid parallel shear flows.J.Fluid Mech.98(2),243e251.http://dx.doi.org/10.1017/S0022112080000122.

Lee,C.B.,Wu,J.Z.,2008.Transition in wall-bounded flows.Appl.Mech.Rev.61(3),030802.http://dx.doi.org/10.1115/1.2909605.

Li,C.B.,2001.On the formation of the streamwise vortex in a transitional boundary layer.Acta Phys.Sin.50(1),182e190 (in Chinese).

Lu,C.G.,Cao,W.D.,Qian,J.H.,2006.A study on numerical method of Navier-stokes equations and nonlinear evolution of the coherent structure in a laminar boundary layer.J.Hydrodynamics (Ser.B)18(3),372e377.http://dx.doi.org/10.1016/S1001-6058(06)60019-X.

Lu,C.G.,Cao,W.D.,Zhang,Y.M.,Peng,J.T.,2008.Large eddies induced by local impulse at wall of boundary layer with pressure gradients.Prog.Nat.Sci.18(7),873e878.http://dx.doi.org/10.1016/j.pnsc.2008.02.007.

Lundbladh,A.,Johansson,A.V.,1991.Direct simulation of turbulent spots in plane Couette flow.J.Fluid Mech.229(1),499e516.http://dx.doi.org/10.1017/S0022112091003130.

Rosenfeld,M.,Cohen,J.,Levinski,V.,1999.The effect of rotation on the growth of hairpin vortices in shear flows.In: the First International Symposium on Turbulence and Shear Flow Phenomena.Begell House Inc.,New York.

Singer,B.A.,1996.Characteristics of a young turbulent spot.Phys.Fluids 37(2),509e521.http://dx.doi.org/10.1063/1.868804.

Svizher,A.,Cohen,J.,2001.The evolution of hairpin vortices in subcritical air channel flow.In: 54th annual meeting of the division of fluid dynamics.American Physical Society,San Diego.

Zhang,L.,Tang,D.B.,2006.Nonlinear evolution of turbulent spots in the near-wall shear flow.Sci.China (Ser.G)36(1),103e112.http://dx.doi.org/10.1007/s11433-006-0158-4.