Evaluation AC6-14: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page SilverP:Evaluation AC6-14 to Evaluation AC6-14)
 
(74 intermediate revisions by the same user not shown)
Line 12: Line 12:
models to filter out all unsteadiness.
models to filter out all unsteadiness.
The survey axes, <math>{S*}</math>, at sections W0&nbsp;&ndash;&nbsp;W2
The survey axes, <math>{S*}</math>, at sections W0&nbsp;&ndash;&nbsp;W2
(see [[AC6-14#figure1|Fig. 1]]), are normalized
(see [[Description_AC6-14#figure1|Fig. 1]]), are normalized
by the throat radius, <math>{R_{throat}}</math>&nbsp;=&nbsp;0.05m, and the velocity is normalized by the bulk
by the throat radius, <math>{R_{throat}}</math>&nbsp;=&nbsp;0.05m, and the velocity is normalized by the bulk
velocity at the throat, <math>{W_{throat}}</math>.
velocity at the throat, <math>{W_{throat}}</math>.
Line 20: Line 20:


[[Evaluation_AC6-14#figure8|Figure 8]] compares the axial and tangential velocity distributions of the
[[Evaluation_AC6-14#figure8|Figure 8]] compares the axial and tangential velocity distributions of the
high-Reynolds number turbulence models, standard $k-\epsilon$, SST $k-\omega$,
high-Reynolds number turbulence models, standard <math>{k-\epsilon}</math>, SST <math>{k-\omega}</math>,
realizable $k-\epsilon$ and RNG $k-\epsilon$, with the experimental results at W1 and W2.
realizable <math>{k-\epsilon}</math> and RNG <math>{k-\epsilon}</math>, with the experimental results at W1 and W2.
It is shown that the turbulence models predict the mean velocity very similarly.
It is shown that the turbulence models predict the mean velocity very similarly.
It is known that a high level of swirl leads to an on-axis recirculation region due to the
It is known that a high level of swirl leads to an on-axis recirculation region due to the
high centrifugal force~\cite{Javadi2015}.
high centrifugal force&nbsp;[[Best_Practice_Advice_AC6-14#References|[1]]].
The models overestimate the size of that region, although the RNG $k-\epsilon$ and
The models overestimate the size of that region, although the RNG <math>{k-\epsilon}</math> and
standard $k-\epsilon$ models predict the recirculation region more realistically than
standard <math>{k-\epsilon}</math> models predict the recirculation region more realistically than
the other models.
the other models.
Javadi et al.~\cite{Javadietal2014} used the RNG $k-\epsilon$ model in the swirl generator
Javadi ''et&nbsp;al.''&nbsp;[[Best_Practice_Advice_AC6-14#References|[10]]] used the RNG <math>{k-\epsilon}</math> model in the swirl generator
over a wide range of runner rotational speeds, from part load to full load.
over a wide range of runner rotational speeds, from part load to full load.
They concluded that the applicability of the RNG $k-\epsilon$ model is restricted to a
They concluded that the applicability of the RNG <math>{k-\epsilon}</math> model is restricted to a
range close to the best efficiency point.
range close to the best efficiency point.
As mentioned before, the current flow field resembles the Francis turbine at 70\% load,
As mentioned before, the current flow field resembles the Francis turbine at 70% load,
which is far from the best efficiency point.
which is far from the best efficiency point.
The disintegration process of the vortex rope is highly unsteady and needs to be resolved
The disintegration process of the vortex rope is highly unsteady and needs to be resolved
Line 42: Line 42:
<div id="figure8"></div>
<div id="figure8"></div>
{|align="center"
{|align="center"
|[[Image:AC6-14_Fig8a.png|325px]]|[[Image:AC6-14_Fig8b.png|325px]]
|[[Image:AC6-14_Fig8a.png|325px]][[Image:AC6-14_Fig8b.png|325px]]
|-
|-
|align="center" width=650px|'''Figure 8:''' Axial and tangential velocity with high-Reynolds number turbulence models at W1. Black: RNG $k-\epsilon$. Blue: standard $k-\epsilon$. Green: realizable $k-\epsilon$. Red: SST $k-\omega$. Star markers (*): experimental axial velocity. Point markers ($\bullet$): experimental tangential velocity. Left: W1 and right: W2.
|align="center" width=650px|'''Figure 8:''' Axial and tangential velocity with high-Reynolds number turbulence models at W1. Black: RNG <math>{k-\epsilon}</math>. Blue: standard <math>{k-\epsilon}</math>. Green: realizable <math>{k-\epsilon}</math>. Red: SST <math>{k-\omega}</math>. Star markers (*): experimental axial velocity. Point markers (<math>{\bullet}</math>): experimental tangential velocity. Left: W1 and right: W2.
|}
|}




Figure \ref{fig7} compares the axial and tangential velocity distributions of the
[[Evaluation_AC6-14#figure9|Figure 9]] compares the axial and tangential velocity distributions of the
low-Reynolds number turbulence models, Launder-Sharma and Lien-cubic $k-\epsilon$, with
low-Reynolds number turbulence models, Launder-Sharma and Lien-cubic <math>{k-\epsilon}</math>, with
the experimental results at W1.
the experimental results at W1.
The results show that Lien-cubic $k-\epsilon$ underestimates the main character of the
The results show that Lien-cubic <math>{k-\epsilon}</math> underestimates the main character of the
flow, the on-axis recirculation region, considerably.
flow, the on-axis recirculation region, considerably.
It can be seen that Launder-Sharma $k-\epsilon$ model, on the other hand, overestimates
It can be seen that Launder-Sharma <math>{k-\epsilon}</math> model, on the other hand, overestimates
the recirculation region.
the recirculation region.
The LS model was initially proposed for predicting swirling flows.
The LS model was initially proposed for predicting swirling flows.
The model, as a modified version of the standard $k-\epsilon$ model, is one of the
The model, as a modified version of the standard <math>{k-\epsilon}</math> model, is one of the
earliest and most widely used models for resolving the near-wall flow behavior.
earliest and most widely used models for resolving the near-wall flow behavior.
Although some damping functions are added in the LS model to account for the viscous and
Although some damping functions are added in the LS model to account for the viscous and
Line 63: Line 63:




Figure \ref{fig8} compares the axial and tangential velocity distributions of the hybrid
<div id="figure9"></div>
models DDES-SA, IDDES-SA and DDES-SST, with the experimental results at sections W0-W2.
{|align="center"
Figure \ref{fig8}a shows that the velocity increases linearly from the hub to the shroud
|[[Image:AC6-14_Fig9.png|650px]]
downstream of the runner, at W0.
|-
|align="center" width=650px|'''Figure 9:''' Axial and tangential velocity with low-Reynolds number turbulence models compared with experimental results at W1. Solid line: LS with an inlet boundary condition extracted from high-Reynolds number simulation. Solid line with + markers: LS with constant inlet boundary condition. Dashed line: LC.  Star markers (*): experimental axial velocity. Point markers (<math>{\bullet}</math>): experimental tangential velocity.
|}
 
 
[[Evaluation_AC6-14#figure10|Figure 10]]
compares the axial and tangential velocity distributions of the hybrid
models DDES-SA, IDDES-SA and DDES-SST, with the experimental results at sections W0&nbsp;&ndash;&nbsp;W2.
[[Evaluation_AC6-14#figure10|Figure 10a]]
shows that the velocity increases linearly from the hub to the shroud downstream of the runner, at W0.
The blade wakes are responsible for the change in the slope of velocity components at
The blade wakes are responsible for the change in the slope of velocity components at
$0.25<S^*<0.4$.
<math>{0.25<S^*<0.4}</math>.
The wakes travel closer to the hub, i.e. the velocity components decrease in the inner
The wakes travel closer to the hub, i.e. the velocity components decrease in the inner
half of the annular cross-section.
half of the annular cross-section.
Line 83: Line 92:
As mentioned before, the vortex breakdown leads to a disintegration of the helical vortex.
As mentioned before, the vortex breakdown leads to a disintegration of the helical vortex.
The vortex rope forms slightly downstream the nozzle and immediately deflects outwards,
The vortex rope forms slightly downstream the nozzle and immediately deflects outwards,
see Figs. \ref{fig11} and \ref{fig13}.
see [[Evaluation_AC6-14#figure12|Figs 12]]&nbsp;and&nbsp;[[Evaluation_AC6-14#figure14|14]].
The on-axis axial velocity is almost zero and the tangential velocity increases linearly
The on-axis axial velocity is almost zero and the tangential velocity increases linearly
with the radius.
with the radius.
Line 90: Line 99:
The central recirculation region thus ends at the throat, yielding an on-axis stagnation
The central recirculation region thus ends at the throat, yielding an on-axis stagnation
region at W1.
region at W1.
The experimental results show $W/W_{throat}$=0.08 on-axis at W1 while numerical results
The experimental results show <math>{W/W_{throat}}</math>&nbsp;=&nbsp;0.08 on-axis at W1 while numerical results
predict $W/W_{throat}$=-0.010, -0.041 and -0.187 with DDES-SA, DDES-SST and IDDES-SA,
predict
<math>{W/W_{throat}}</math>&nbsp;=&nbsp;-0.010, -0.041 and -0.187 with DDES-SA, DDES-SST and IDDES-SA,
respectively.
respectively.
Figure \ref{fig8}c shows the axial and tangential velocity components at W2.
Figure [[Evaluation_AC6-14#figure10|10c]] shows the axial and tangential velocity components at W2.
The central cone-shaped stagnant region in the draft tube occupies 60\% of the
The central cone-shaped stagnant region in the draft tube occupies 60% of the
cross-section and is thus much wider than at W1.
cross-section and is thus much wider than at W1.
The velocity profiles at W2 follow the same patterns as at W1 but with a lower magnitude.
The velocity profiles at W2 follow the same patterns as at W1 but with a lower magnitude.


Figure \ref{fig9} shows the resolved root mean square of the axial and tangential velocity
 
fluctuations at W0 and W1, $w^{'}_{rms}$ and $u^{'}_{rms}$, respectively.
<div id="figure10"></div>
The results of DDES-SA and the RNG $k-\epsilon$ models are compared with the experimental
{|align="center"
|[[Image:AC6-14_Fig10a.png|325px]]||[[Image:AC6-14_Fig10b.png|325px]]
|-
|align="center"|(a) W0||align="center"|(b) W1
|-
|colspan=2 align="center"|[[Image:AC6-14_Fig10c.png|325px]]
|-
|align="center" colspan=2|(c) W2
|-
|align="center" colspan=2 width=650px|'''Figure 10:''' Axial and tangential velocity distributions of the hybrid URANS-LES models, compared with the experimental results. Black: DDES-SA. Red: IDDES-SA. Blue: DDES-SST. Star markers (*): experimental axial velocity. Point markers (<math>{\bullet}</math>): experimental tangential velocity.
|}
 
 
[[Evaluation_AC6-14#figure11|Figure 11]]
shows the resolved root mean square of the axial and tangential velocity
fluctuations at W0 and W1, <math>{w^{'}_{rms}}</math> and <math>{u^{'}_{rms}}</math>, respectively.
The results of DDES-SA and the RNG <math>{k-\epsilon}</math> models are compared with the experimental
results.
results.
The modeled part of the turbulence is different model to model,
The modeled part of the turbulence is different model to model,
$(u^{'}_{rms}, w^{'}_{rms})_{modeled}/W_{throat}=0.1$ at both W0 and W1 for DDES-SA,
<math>{(u^{'}_{rms}, w^{'}_{rms})_{modeled}/W_{throat}=0.1}</math>
and $(u^{'}_{rms}, w^{'}_{rms})_{modeled}/W_{throat}=0.02$ at W0 and $(u^{'}_{rms},
at both W0 and W1 for DDES-SA,
w^{'}_{rms})_{modeled}/W_{throat}=0.14$ at W1 for the RNG $k-\epsilon$ model.
and <math>{(u^{'}_{rms}, w^{'}_{rms})_{modeled}/W_{throat}=0.02}</math> at W0 and
<math>{(u^{'}_{rms}, w^{'}_{rms})_{modeled}/W_{throat}=0.14}</math> at W1 for the RNG <math>{k-\epsilon}</math> model.
Therefore, the total turbulence predicted by DDES-SA is slightly higher than the
Therefore, the total turbulence predicted by DDES-SA is slightly higher than the
experimental results for both components at W0 and is closer to the experimental results
experimental results for both components at W0 and is closer to the experimental results
at W1.
at W1.
The ratio of modeled to total part of the turbulence at the peaks are almost 50\% at W0
The ratio of modeled to total part of the turbulence at the peaks are almost 50% at W0
and 15\% at W1 for DDES-SA.
and 15% at W1 for DDES-SA.
The RNG $k-\epsilon$ model still underestimates the turbulence at W0 and W1.
The RNG <math>{k-\epsilon}</math> model still underestimates the turbulence at W0 and W1.
At W0, $w^{'}_{rms}$ has a single experimental peak at $S^*$=0.25, which is predicted at
At W0, <math>{w^{'}_{rms}}</math> has a single experimental peak at <math>{S^*}</math>&nbsp;=&nbsp;0.25, which is predicted at
$S^*$=0.3 by DDES-SA.
<math>{S^*}</math>&nbsp;=&nbsp;0.3 by DDES-SA.
Figure \ref{fig11} shows that the blade wakes are responsible for this peak.
 
 
<div id="figure11"></div>
{|align="center"
|[[Image:AC6-14_Fig11a.png|325px]]||[[Image:AC6-14_Fig11b.png|325px]]
|-
|align="center"|(a) W0||align="center"|(b) W1
|-
|align="center" width=650px colspan="2"|'''Figure 11:''' Axial , <math>{w^{'}_{rms}}</math> , and tangential, <math>{u^{'}_{rms}}</math> , root mean square velocity fluctuations with DDES-SA compared with the experimental results. Green: DDES-SA, <math>{w^{'}_{rms}}</math>. Blue: DDES-SA, <math>{u^{'}_{rms}}</math>. Red: RNG <math>{k-\epsilon}</math>, <math>{w^{'}_{rms}}</math>. Turquoise: RNG <math>{k-\epsilon}</math>, <math>{u^{'}_{rms}}</math>. Star markers (*): experimental <math>{w^{'}_{rms}}</math>. Point markers (<math>{\bullet}</math>):&nbsp;experimental&nbsp;<math>{u^{'}_{rms}}</math>.
|}
 
 
[[Evaluation_AC6-14#figure12|Figure 12]]
shows that the blade wakes are responsible for this peak.
The numerical result is quantitatively close to the experimental results.
The numerical result is quantitatively close to the experimental results.
The experimental tangential velocity fluctuation,  $u^{'}_{rms}$, increases moderately
The experimental tangential velocity fluctuation,  <math>{u^{'}_{rms}}</math>, increases moderately
from hub to shroud, while the numerical profile suggests a more pronounced effect of the
from hub to shroud, while the numerical profile suggests a more pronounced effect of the
wakes.
wakes.
At W1, the axial velocity fluctuation, $w^{'}_{rms}$, has two peaks, while the tangential
At W1, the axial velocity fluctuation, <math>{w^{'}_{rms}}</math>, has two peaks, while the tangential
velocity fluctuation, $u^{'}_{rms}$, has a single peak.
velocity fluctuation, <math>{u^{'}_{rms}}</math>, has a single peak.
The two peaks of the axial velocity fluctuation are predicted reasonably well, which
The two peaks of the axial velocity fluctuation are predicted reasonably well, which
confirms that DDES-SA is able to resolve the turbulent structures very well.
confirms that DDES-SA is able to resolve the turbulent structures very well.
These two peaks are related to the shear layer between the inner stagnant region and the
These two peaks are related to the shear layer between the inner stagnant region and the
outer flow, which is the source of most of the turbulence in the draft tube.
outer flow, which is the source of most of the turbulence in the draft tube.
The total $w^{'}_{rms}$ predicted by DDES-SA is slightly higher than the experimental
The total <math>{w^{'}_{rms}}</math> predicted by DDES-SA is slightly higher than the experimental
results, if the right peak of the profile at $S^*\sim1.4$ is considered as the reference.
results, if the right peak of the profile at <math>{S^*\sim1.4}</math> is considered as the reference.
The total $w^{'}_{rms}$ predicted by the RNG $k-\epsilon$ model agrees excellently with
The total <math>{w^{'}_{rms}}</math> predicted by the RNG <math>{k-\epsilon}</math> model agrees excellently with
the experimental results.
the experimental results.
The peak of the tangential velocity fluctuation can be related to the swirl in the on-axis
The peak of the tangential velocity fluctuation can be related to the swirl in the on-axis
recirculation region.
recirculation region.
DDES-SA is capable of predicting the high level of turbulence production of
DDES-SA is capable of predicting the high level of turbulence production of
$u^{'}_{rms}$, yielding an increase of $u^{'}_{rms}$ by a factor of seven between
<math>{u^{'}_{rms}}</math>, yielding an increase of <math>{u^{'}_{rms}}</math> by a factor of seven between
W0 and W1.
W0 and W1.
The RNG $k-\epsilon$ model still underpredicts the on-axis peak of $u^{'}_{rms}$.
The RNG <math>{k-\epsilon}</math> model still underpredicts the on-axis peak of <math>{u^{'}_{rms}}</math>.
 
 
 




Line 147: Line 184:
|}
|}


Figure \ref{fig11} shows the vorticity magnitude predicted by DDES-SA at three planes in
[[Evaluation_AC6-14#figure12|Figure 12]] shows the vorticity magnitude predicted by DDES-SA at three planes in
the draft tube.
the draft tube.
It is worth mentioning that the wakes are tenuous compared with the vortex rope in the
It is worth mentioning that the wakes are tenuous compared with the vortex rope in the
draft tube.
draft tube.
If DDES-SA is considered as the reference results, the RNG $k-\epsilon$ gives a smaller
If DDES-SA is considered as the reference results, the RNG <math>{k-\epsilon}</math> gives a smaller
inter-blade region with lower vorticity and a larger wake with higher vorticity compared
inter-blade region with lower vorticity and a larger wake with higher vorticity compared
with LS.
with LS.
Thus, the RNG $k-\epsilon$ model can be interpreted as an eddy-viscosity model which
Thus, the RNG <math>{k-\epsilon}</math> model can be interpreted as an eddy-viscosity model which
predicts better blade wakes.
predicts better blade wakes.


Figure \ref{fig12} shows the iso-surface of the pressure, representing the vortex rope
[[Evaluation_AC6-14#figure13|Figure 13]]
captured by the standard $k-\epsilon$, RNG $k-\epsilon$, LS and DDES-SA models.
shows the iso-surface of the pressure, representing the vortex rope
It can be seen that standard $k-\epsilon$ and RNG $k-\epsilon$ models predict very similar
captured by the standard <math>{k-\epsilon}</math>, RNG <math>{k-\epsilon}</math>, LS and DDES-SA models.
It can be seen that standard <math>{k-\epsilon}</math> and RNG <math>{k-\epsilon}</math> models predict very similar
vortex ropes.
vortex ropes.
The LS model predicts a vortex rope which is axially more compact than other
The LS model predicts a vortex rope which is axially more compact than other
Line 167: Line 205:




 
<div id="figure13"></div>
 
 
<div id="figure14"></div>
{|align="center"
{|align="center"
|[[Image:AC6-14_Fig13.png|x300px]]
|[[Image:AC6-14_Fig13a.png|x330px]]||[[Image:AC6-14_Fig13b.png|x330px]]||[[Image:AC6-14_Fig13c.png|x310px]]||[[Image:AC6-14_Fig13d.png|x310px]]
|-
|align="center"|a)||align="center"|b)||align="center"|c)||align="center"|d)
|-
|-
|align="center"|'''Figure 14:''' Iso-surface of ''q''-criterion, DDES-SA
|align="center" width=650px colspan="4"|'''Figure 13:''' Iso-surface of pressure, visualizing the vortex rope, a) standard <math>{k-\epsilon}</math>, b) RNG <math>{k-\epsilon}</math>, c) LS, d) DDES&#8209;SA.
|}
|}


Figure \ref{fig13} shows an iso-surface of the \textit{q}-criterion by DDES-SA.
 
[[Evaluation_AC6-14#figure14|Figure 14]]
shows an iso-surface of the <math>{q}</math>-criterion by DDES-SA.
The separations and the consecutive coherent structures form after the leading edge at
The separations and the consecutive coherent structures form after the leading edge at
the suction side close to the hub and move radially outward due to the pressure gradient
the suction side close to the hub and move radially outward due to the pressure gradient
which is caused by the swirl.
which is caused by the swirl.
As can be seen in Fig. \ref{fig11},
As can be seen in [[Evaluation_AC6-14#figure12|Fig. 12]], the structures are stronger at the outer part of the blade.
the structures are stronger at the outer part of the blade.
The blade wakes are deflected toward the hub due to the divergent part of the draft tube.
The blade wakes are deflected toward the hub due to the divergent part of the draft tube.
Figure \ref{fig13} also shows the disintegration of the vortex rope.
[[Evaluation_AC6-14#figure14|Figure 14]] also shows the disintegration of the vortex rope.
Further downstream, the wakes and the vortex rope come together, and the precessing
Further downstream, the wakes and the vortex rope come together, and the precessing
vortex rope ruptures the structures of the wakes.
vortex rope ruptures the structures of the wakes.
Line 190: Line 228:
nozzle in the draft tube.
nozzle in the draft tube.
This forms a small recirculation region nested downstream of the nozzle in the draft tube.
This forms a small recirculation region nested downstream of the nozzle in the draft tube.
<div id="figure14"></div>
{|align="center"
|[[Image:AC6-14_Fig13.png|x300px]]
|-
|align="center"|'''Figure 14:''' Iso-surface of <math>{q}</math>-criterion, DDES-SA
|}
<br/>
<br/>
----
----
{{ACContribs
{{ACContribs
| authors=A. Javadi, A. Bosioc, H Nilsson, S. Muntean, R. Susan-Resiga
| authors=A. Javadi<sup>a</sup>, A. Bosioc<sup>b</sup>, H Nilsson<sup>a</sup>, S. Muntean<sup>c</sup>, R. Susan-Resiga<sup>b</sup>
| organisation=Chalmers University of Technology
| organisation=<sup>a</sup>Chalmers University of Technology, G&ouml;teborg, Sweden; <sup>b</sup>University Polytehnica Timi&#351;oara, Timi&#351;oara, Romania; <sup>c</sup>Center for Advanced Research in Engineering Sciences, Romanian Academy, Timi&#351;oara Branch, Timi&#351;oara, Romania
}}
}}
{{ACHeader
{{ACHeader

Latest revision as of 18:49, 11 February 2017


Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

Swirling flow in a conical diffuser generated with rotor-stator interaction

Application Challenge AC6-14   © copyright ERCOFTAC 2024

Comparison of Test Data and CFD

The mean velocity field is determined by time averaging over five complete runner revolutions for the hybrid models and eight complete runner revolutions for the URANS models to filter out all unsteadiness. The survey axes, , at sections W0 – W2 (see Fig. 1), are normalized by the throat radius,  = 0.05m, and the velocity is normalized by the bulk velocity at the throat, . The axial axis is downward and the runner rotates in the positive direction according to the right-hand rule.


Figure 8 compares the axial and tangential velocity distributions of the high-Reynolds number turbulence models, standard , SST , realizable and RNG , with the experimental results at W1 and W2. It is shown that the turbulence models predict the mean velocity very similarly. It is known that a high level of swirl leads to an on-axis recirculation region due to the high centrifugal force [1]. The models overestimate the size of that region, although the RNG and standard models predict the recirculation region more realistically than the other models. Javadi et al. [10] used the RNG model in the swirl generator over a wide range of runner rotational speeds, from part load to full load. They concluded that the applicability of the RNG model is restricted to a range close to the best efficiency point. As mentioned before, the current flow field resembles the Francis turbine at 70% load, which is far from the best efficiency point. The disintegration process of the vortex rope is highly unsteady and needs to be resolved to capture its realistic representation. The turbulence models reproduce the width of the central stagnant region reasonably well at W2.


AC6-14 Fig8a.pngAC6-14 Fig8b.png
Figure 8: Axial and tangential velocity with high-Reynolds number turbulence models at W1. Black: RNG . Blue: standard . Green: realizable . Red: SST . Star markers (*): experimental axial velocity. Point markers (): experimental tangential velocity. Left: W1 and right: W2.


Figure 9 compares the axial and tangential velocity distributions of the low-Reynolds number turbulence models, Launder-Sharma and Lien-cubic , with the experimental results at W1. The results show that Lien-cubic underestimates the main character of the flow, the on-axis recirculation region, considerably. It can be seen that Launder-Sharma model, on the other hand, overestimates the recirculation region. The LS model was initially proposed for predicting swirling flows. The model, as a modified version of the standard model, is one of the earliest and most widely used models for resolving the near-wall flow behavior. Although some damping functions are added in the LS model to account for the viscous and wall effects, the improved near-wall behavior does not significantly increase the quality of the mean velocity profiles in the draft tube.


AC6-14 Fig9.png
Figure 9: Axial and tangential velocity with low-Reynolds number turbulence models compared with experimental results at W1. Solid line: LS with an inlet boundary condition extracted from high-Reynolds number simulation. Solid line with + markers: LS with constant inlet boundary condition. Dashed line: LC. Star markers (*): experimental axial velocity. Point markers (): experimental tangential velocity.


Figure 10 compares the axial and tangential velocity distributions of the hybrid models DDES-SA, IDDES-SA and DDES-SST, with the experimental results at sections W0 – W2. Figure 10a shows that the velocity increases linearly from the hub to the shroud downstream of the runner, at W0. The blade wakes are responsible for the change in the slope of velocity components at . The wakes travel closer to the hub, i.e. the velocity components decrease in the inner half of the annular cross-section. The hybrid models realistically predict the inner half of the flow, close to the hub. The axial velocity in the outer part, close to the shroud, is overestimated and the tangential velocity is underestimated, i.e. the predicted flow is too axial in the outer part, at W0. The effects of the wakes on the averaged velocity is overestimated by DDES-SA and IDDES-SA. Nevertheless, the wakes predicted by these two models are remarkably in spite of while a different wall treatment being used by IDDES-SA. DDES-SST presents weaker wakes compared with the other models while the results agree with the experimental results. Figure 8b shows the axial and tangential velocity components at W1, close to the throat. As mentioned before, the vortex breakdown leads to a disintegration of the helical vortex. The vortex rope forms slightly downstream the nozzle and immediately deflects outwards, see Figs 12 and 14. The on-axis axial velocity is almost zero and the tangential velocity increases linearly with the radius. In the outer part of the draft tube, both the axial and tangential velocity magnitudes increase at W1. The central recirculation region thus ends at the throat, yielding an on-axis stagnation region at W1. The experimental results show  = 0.08 on-axis at W1 while numerical results predict  = -0.010, -0.041 and -0.187 with DDES-SA, DDES-SST and IDDES-SA, respectively. Figure 10c shows the axial and tangential velocity components at W2. The central cone-shaped stagnant region in the draft tube occupies 60% of the cross-section and is thus much wider than at W1. The velocity profiles at W2 follow the same patterns as at W1 but with a lower magnitude.


AC6-14 Fig10a.png AC6-14 Fig10b.png
(a) W0 (b) W1
AC6-14 Fig10c.png
(c) W2
Figure 10: Axial and tangential velocity distributions of the hybrid URANS-LES models, compared with the experimental results. Black: DDES-SA. Red: IDDES-SA. Blue: DDES-SST. Star markers (*): experimental axial velocity. Point markers (): experimental tangential velocity.


Figure 11 shows the resolved root mean square of the axial and tangential velocity fluctuations at W0 and W1, and , respectively. The results of DDES-SA and the RNG models are compared with the experimental results. The modeled part of the turbulence is different model to model, at both W0 and W1 for DDES-SA, and at W0 and at W1 for the RNG model. Therefore, the total turbulence predicted by DDES-SA is slightly higher than the experimental results for both components at W0 and is closer to the experimental results at W1. The ratio of modeled to total part of the turbulence at the peaks are almost 50% at W0 and 15% at W1 for DDES-SA. The RNG model still underestimates the turbulence at W0 and W1. At W0, has a single experimental peak at  = 0.25, which is predicted at  = 0.3 by DDES-SA.


AC6-14 Fig11a.png AC6-14 Fig11b.png
(a) W0 (b) W1
Figure 11: Axial , , and tangential, , root mean square velocity fluctuations with DDES-SA compared with the experimental results. Green: DDES-SA, . Blue: DDES-SA, . Red: RNG , . Turquoise: RNG , . Star markers (*): experimental . Point markers (): experimental .


Figure 12 shows that the blade wakes are responsible for this peak. The numerical result is quantitatively close to the experimental results. The experimental tangential velocity fluctuation, , increases moderately from hub to shroud, while the numerical profile suggests a more pronounced effect of the wakes. At W1, the axial velocity fluctuation, , has two peaks, while the tangential velocity fluctuation, , has a single peak. The two peaks of the axial velocity fluctuation are predicted reasonably well, which confirms that DDES-SA is able to resolve the turbulent structures very well. These two peaks are related to the shear layer between the inner stagnant region and the outer flow, which is the source of most of the turbulence in the draft tube. The total predicted by DDES-SA is slightly higher than the experimental results, if the right peak of the profile at is considered as the reference. The total predicted by the RNG model agrees excellently with the experimental results. The peak of the tangential velocity fluctuation can be related to the swirl in the on-axis recirculation region. DDES-SA is capable of predicting the high level of turbulence production of , yielding an increase of by a factor of seven between W0 and W1. The RNG model still underpredicts the on-axis peak of .


AC6-14 Fig11.png
Figure 12: Vorticity magnitude with DDES-SA at three horizontal planes in the draft tube

Figure 12 shows the vorticity magnitude predicted by DDES-SA at three planes in the draft tube. It is worth mentioning that the wakes are tenuous compared with the vortex rope in the draft tube. If DDES-SA is considered as the reference results, the RNG gives a smaller inter-blade region with lower vorticity and a larger wake with higher vorticity compared with LS. Thus, the RNG model can be interpreted as an eddy-viscosity model which predicts better blade wakes.

Figure 13 shows the iso-surface of the pressure, representing the vortex rope captured by the standard , RNG , LS and DDES-SA models. It can be seen that standard and RNG models predict very similar vortex ropes. The LS model predicts a vortex rope which is axially more compact than other eddy-viscosity models. DDES-SA predicts the most detailed vortex rope. Thus, DDES-SA predicts a higher axial surge of the vortex rope.


AC6-14 Fig13a.png AC6-14 Fig13b.png AC6-14 Fig13c.png AC6-14 Fig13d.png
a) b) c) d)
Figure 13: Iso-surface of pressure, visualizing the vortex rope, a) standard , b) RNG , c) LS, d) DDES‑SA.


Figure 14 shows an iso-surface of the -criterion by DDES-SA. The separations and the consecutive coherent structures form after the leading edge at the suction side close to the hub and move radially outward due to the pressure gradient which is caused by the swirl. As can be seen in Fig. 12, the structures are stronger at the outer part of the blade. The blade wakes are deflected toward the hub due to the divergent part of the draft tube. Figure 14 also shows the disintegration of the vortex rope. Further downstream, the wakes and the vortex rope come together, and the precessing vortex rope ruptures the structures of the wakes. It is noteworthy that the inner part of the wake attaches and then detaches to/from the nozzle in the draft tube. This forms a small recirculation region nested downstream of the nozzle in the draft tube.


AC6-14 Fig13.png
Figure 14: Iso-surface of -criterion, DDES-SA



Contributed by: A. Javadia, A. Bosiocb, H Nilssona, S. Munteanc, R. Susan-Resigab — aChalmers University of Technology, Göteborg, Sweden; bUniversity Polytehnica Timişoara, Timişoara, Romania; cCenter for Advanced Research in Engineering Sciences, Romanian Academy, Timişoara Branch, Timişoara, Romania

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice


© copyright ERCOFTAC 2024

Retrieved from "http://kbwiki.ercoftac.org/w/index.php?title=Evaluation_AC6-14&oldid=32986"