Turbulent Vortex Shedding From Triangle Cylinder Using · PDF fileturbulent vortex shedding...
Transcript of Turbulent Vortex Shedding From Triangle Cylinder Using · PDF fileturbulent vortex shedding...
Proceedings of ASME FEDSM’00ASME 2000 Fluids Engineering Division Summer Meeting
June 11-15, 2000 Boston, Massachusetts
FEDSM2000-11172
TURBULENT VORTEX SHEDDING FROM TRIANGLE CYLINDER USING THETURBULENT BODY FORCE POTENTIAL MODEL
Xing Zhang Blair Perot
Department of Mechanical and Industrial Engineering, University of Massachusetts, AmherstAmherst, MA 01003, U.S.A.
ABSTRACTNumerical simulation of the turbulent flow around a
triangular cylinder at a Reynolds number of 45,000 ispresented in this paper. Both steady and unsteady vortex-shedding results are presented. A body force potentialmodel is used to model the turbulent motion. Thisapproach is able to model non-equilibrium turbulenceaccurately at a cost and complexity comparable to k-εmodels. The numerical method used in this calculation isan unstructured staggered mesh scheme. The property thatthis method conserves kinetic energy both locally withincells and globally makes it a good choice for performingturbulence modeling.
For the unsteady solution, the Strouhal number andtime-averaged velocity profile agree well withexperiments. However, the steady solution that wasobtained by using a symmetric boundary condition at thecenterline leads to poor predictions of the time-averagedmean velocity profile.
INTRODUCTIONThe flow around a triangle provides an example of
bluff body flow with fixed separation points. If theReynolds number is not too small the flow is inherentlyunsteady and a Von Karman vortex street appears with awell-defined frequency. If the Reynolds number issufficiently high the flow will be turbulent and aturbulence model must be included to model the turbulentfluctuations. In the case of turbulent vortex shedding, wehave the option of including the large scale vortexshedding in the turbulence model and calculating a steadymean flow, or of solving for the large scale vortex
shedding by numerical scheme while only including thesmall scale turbulence in the model. The former approachis less expensive, but we show here less likely to giveaccurate predictions. This is hypothesized to be due to thefact that the large-scale vortex structures do not behavelike equilibrium turbulence.
Sjunnesson (1991) measured the flow of a triangularcylinder in a duct. Their experimental study wasmotivated by the application to flame holders. Johanssonet al. (1993) carried out numerical simulation of this flowusing a k-ε model. Durbin (1994) also carried out the
same simulation using a k-ε- 2v model. Franke et al.(1991) compared the ability of different models to predictturbulent vortex shedding from a rectangular cylinder.Franke’s conclusion is that some k-ε models do notpredict the right shedding frequency and Reynolds stresstransport models can produce results in good agreementwith the experiments. The turbulent potential model is asimplified Reynolds stress transport model, which has theability of modeling non-equilibrium turbulence with thecomputing cost and complexity comparable to k-ε model.
TURBULENCE MODEL AND NUMERICAL SCHEMEThe primary difficulty of modeling unsteady
turbulent vortex shedding is thought to be that theturbulence is not in equilibrium with the mean flow,because the large vortices move and decay at the sametime-scale as the turbulence. The most commonconstitutive relation, the eddy viscosity hypothesis (orlinear Boussinesq hypothesis) is probably incorrect inthis case.
1 Copyright © 2000 by ASME
In the past, avoiding an algebraic constitutive relationfor the Reynolds stresses required solving coupledtransport equations for the Reynolds stress themselves.Recently, a new modeling approach, the turbulencepotential model has been developed. It is capable ofmodeling the complex turbulent physics associated withseparation and unsteady flow. This turbulent potentialmodel is well suited to vortex shedding problem because itdoes not require a constitutive relation relating theReynolds stress tensor to the mean flow. The modelhypothesizes evolution equations for the scalar and vectorpotentials of the turbulent body force (the divergence ofthe Reynolds stress tensor). It has the accuracy of aReynolds stress model, at a cost comparable to moderntwo equation models. The governing equations of theturbulence potential model will not be presented here.Their initial development is described in Perot (1999).
Our numerical method uses an unstructured staggeredmesh scheme which can conserve mass, momentum, andkinetic energy to machine precision. The turbulencequantities are advected using an unwinding scheme toguarantee positivity constraints. The model integrates upto the wall, so wall functions are not used, but the firstgrid point should be in the laminar sub-layer to obtainaccurate predictions. The details for this numericalmethod, including accuracy analysis and conservationproperty are discussed in Perot & Zhang (1999).
FLOW OVER A TRIANGLEIn order to compare with the experimental data, we
select a computational domain that is the same as theconfiguration of Sjunnesson’s experiment. The mesh isgenerated by TRIANGLE – an automatic 2D-Delaunaymesh maker. There are approximately 25,000 triangles inour calculation (see Figure 1).
In the present calculation, the inlet mean stream-wisevelocity is a constant value, the vertical velocity is zero.For turbulent kinetic energy and dissipation rate, we usethe same conditions described in Johnasson’s paper.
Figure 1. Computational domain and mesh.
0.17=inU m/s2)05.0( inin Uk =
l2.0
16.0 23in
in
k=ε
The total mass flow was 6.01 =m& kgs-1 in theirexperiment, and the inlet velocity is evaluated based onthat value. These values are also used as the initial valuefor the whole domain. l is the height of the duct. Azero gradient boundary condition is used for all thevariables at the outlet. Slip-wall boundary conditions areused for the duct wall.
In the steady calculation, we use half of the domainmentioned above and imposed a symmetric boundarycondition along the centerline.
RESULTS Calculation of 2D unsteady turbulent flow around a
triangle cylinder with Reynolds number
000,45Re ==νHU in is presented, where H is the height
of the triangle. Unsteady behavior is due to vorticesalternately shedding from the upper and lower edges of thecylinder, forming a Von Karmann vortex street behind thetriangle.
No special triggering measure is taken to start thevortex shedding, the unsteadiness in the computationalresult evolved naturally. It was triggered by the machineerror and asymmetry of the mesh.
To illustrate the periodicity of the flow, the streamfunction of a point about one triangle height behind thetriangle near the centerline is shown in Figure 2. It can beseen that an almost perfect periodicity exists. Theshedding frequency is 109.3 (s-1). The correspondingStrouhal number defined by,
inU
fHSr =
is 0.257, which should be compared with experimentaldata of 0.25 and the computed value of 0.27 inJohnnasson (1991). Figure 3 shows an instantaneousvelocity vector plot, we can see that the center of a vortexis rolled up at the lower edge and a new vortex isbeginning to roll up at the upper edge. The vortex streetcan also be seen in the instantaneous vorticity contoursplot shown in Figure 4.
2 Copyright © 2000 by ASME
t (sec)
S
0.018 0.02 0.022 0.024 0.026 0.028 0.03 0.032 0.034 0.036-0.45
-0.3
-0.15
0
0.15
0.3
0.45
Figure 2. The stream-function of one point about one cylinderheight behind the triangle near the centerline.
Figure 3 Instantaneous velocity vector plot.
Figure 4. Instantaneous vorticity contours plot.
3 Copyright © 2000 by ASME
Although the instantaneous flow is asymmetric, thetime-averaged fields are always symmetric or anti-symmetric. Figure 5 shows the stream-wise velocity atdifferent locations behind the triangle. The calculatedvelocity profiles are in reasonable agreement with theexperiment. However, it is hypothesized that theboundary layer on the triangle is not fully resolved due tomesh size restrictions. The computed boundary layer ismuch thicker than the real one, thus close to the back ofthe triangle, the fluid is slowed down and drivenbackwards more than it should be. This would explain themean velocity profile close to the centerline at x=15mmwhere the velocity is under-predicted. Figure 6 shows themean stream-wise velocity at the centerline. The lengthof recirculation zone is accurately predicted, while the
location of the maximum negative velocity is slightlyupstream compare with the experiments. The magnitudeof the maximum negative velocity is also a little lowerthan the experiment data.
Figure 7 shows comparison of the mean stream-wisevelocity contour plots between the time-averaged and thesteady solution. The contour levels in each plot are thesame. The predicted “steady-state” recirculation zone ismuch longer than the time-averaged unsteady solution.The reason for this is that the unsteady flow increases themomentum exchange between the wake and itssurrounding, thus reducing the recirculation zone. Theturbulence model does not adequately represent themomentum exchange due to these very large eddies
U [m/s]
y [
m]
-10 10 30-0.06
-0.04
-0.02
0
0.02
0.04
0.06
U [m/s]
y [
m]
-10 10 30-0.06
-0.04
-0.02
0
0.02
0.04
0.06
U [m/s]
y [
m]
-10 10 30-0.06
-0.04
-0.02
0
0.02
0.04
0.06
(a) (b) (c)
U [m/s]
y [
m]
-10 10 30-0.06
-0.04
-0.02
0
0.02
0.04
0.06
U [m/s]
y [
m]
-10 10 30-0.06
-0.04
-0.02
0
0.02
0.04
0.06
(d) (e)Figure 5. Mean stream-wise velocity behind the triangle: , calculations; *, experiments.
(a) 15mm, (b) 38mm, (c) 61mm, (d) 150mm, (e) 376mm
4 Copyright © 2000 by ASME
Figure 7. Mean stream-wise velocity contours of the time-averaged and steady solution
X [m]
U
[m/s
]
0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 0.36 0.4-24
-16
-8
0
8
16
24
32
calculationsexperiments
Figure 6. Mean stream-wise velocity at centerline.
CONCLUSIONSIn this paper, numerical simulation of flow past
triangular cylinder at high Reynolds number (45,000) ispresented. The instantaneous flow situation is verycomplex due to the presence of vortex shedding andturbulence.
The calculation was performed using an unstructuredstaggered mesh scheme. A turbulent potential model isused to model the small-scale fluctuation motion.
The capability of the turbulent potential model topredict turbulent vortex shedding has been demonstratedin this calculation. Computed Strouhal number and meanvelocity profiles down stream of the triangle cylinder arein agreement with experiment data. In addition, it hasbeen shown that statistical unsteadiness produced byvortex shedding must be resolved in order to simulate theflow correctly. The steady state computation of this flowwill lead to poor predictions.
5 Copyright © 2000 by ASME
ACKNOWLEDGMENTS The financial support of the Office of Naval Research(grant number: N00014-99-1-0194) is gratefullyacknowledged.
REFERENCESJ. B. Perot, Turbulence Modeling Using Body Force
Potentials. Physics of Fluids, 11 (9), 1999
J. B. Perot and X. Zhang, Reformulation of theunstructured staggered mesh method as a classic finitevolume method. Second International Symposium onFinite Volumes for Complex Applications: Problems andPerspectives, July 19-22,1999, Duisberg, Germany
Stefan H. Johansson, Lars Davidson and Erik Olsson,Numerical Simulation of Vortex Shedding PastTriangular Cylinders at High Reynolds Number Using ak-ε Turbulence Model. International Journal ForNumerical Methods In Fluids, VOL. 16, 859-878 (1993)
P. A. Durbin, Turbulence modeling for separated flow.Annual Research Briefs-1994, Center for TurbulenceResearch, Stanford University
A. Sjunnesson, C. Nelsons and E. Max, LDAmeasurement of velocities and turbulence in a bluff bodystabilized flame. Volvo Flygmotor AB, Trollhattan, 1991
R. Franke and W. Rodi, Calculation of vortexshedding past a square cylinder with various turbulencemodels. Proc. 8th Symp. On Turbulent Shear Flows,Munich, 1991, pp. 20.1.1- 20.1.6
.
6 Copyright © 2000 by ASME