DNS 1-3 Description: Difference between revisions

From KBwiki
Jump to navigation Jump to search
 
(39 intermediate revisions by 2 users not shown)
Line 1: Line 1:


{{DNSHeaderLib
=Flow in a 3D diffuser=
{{DNSHeader
|area=1
|area=1
|number=3
|number=3
Line 6: Line 7:
__TOC__
__TOC__
= Introduction =
= Introduction =
The 3D (Stanford) Diffuser is a well documented case with complex internal corner flow and 3D separation while having a relatively simple geometry. It has an inlet section, an expansion section and an outlet section (see Figure 1). The flow at the inlet is assumed to be a fully developed rectangular channel flow. At the outlet, standard Dirichlet condition for the pressure is prescribed. An inflow Reynolds number of 10000 is considered based on the duct height and the flow is considered to be incompressible. The following DNS data has been obtained using the in-house Finite Element Method (FEM) code Alya developed at BSC.
The 3D (Stanford) Diffuser is a well documented case with complex internal corner flow and 3D separation while having a relatively simple geometry. It has an inlet section, an expansion section and an outlet section (see Figure 1). The flow at the inlet is assumed to be a fully developed rectangular channel flow. At the outlet, standard Dirichlet condition for the pressure is prescribed. An inflow Reynolds number of 10000 is considered based on the duct height and the flow is considered to be incompressible. Such configuration is also of engineering relevance. It represents a diffuser situated between the compressor and the combustor chamber of a jet engine, whose task is to decelerate the flow from the compressor to adapt for the combustor section. The following DNS data has been obtained using the in-house Finite Element Method (FEM) code Alya developed at BSC.


= Review of previous studies =
= Review of previous studies =
This diffuser configuration has already been investigated in the framework of two ERCOFTAC-SIG15 Workshops and in the European ATAAC project. They studied both 3D diffuser configurations (denoted as SIG15 Case 13.2-1 and SIG15 Case 13.2-2, respectively) and were held in Austria (September, 2008) and Italy (September, 2009). Both in the workshops and in the ATAAC project a wide range of turbulence models in both LES and RANS frameworks as well as some novel Hybrid LES/RANS formulations have been employed.  
The physics of the Stanford diffuser were experimentally analyzed by [[lib:DNS_1-3_description#CFD_Methods|Cherry ''et al.''  (2008)]]. The sharp corners of the diffuser were smoothed with a fillet radius of 6.0. Experimental data was gathered in a recirculating water channel using the method of magnetic resonance velocimetry (MRV). In this turbulence causes a loss of net magnetization signal from the studied region. Such decrease in signal strength can be related to turbulent velocity statistics. This method was found to be accurate within 20% of the field of view and within 5% in  in regions of high turbulence. Three scans were performed using three different magnetic field gradient intensities. For each one, 30 scans were completed and averaged. Then, the three averaged data sets were averaged to obtain a final data set. Such dataset provided detailed experimental data comprising the mean velocity field on its three components, the streamwise Reynolds stresses and the pressure distribution along the bottom wall of the diffuser. This provided information on the mean flow configuration, however, little insight was given on the more complex time-motions of this flow.  


The workshop reports are published in the ERCOFTAC Bulletin Issues, see [[lib:DNS_1-3_description#3|Steiner ''et al.''  (2009)]], [[lib:DNS_1-3_description#3| Jakirlić ''et al.''  (2010b)]] and the [[UFR_4-16_Test_Case#CFD_Methods| CFD methods section of UFR4-16]]. The ATAAC reports can be found through the links through the links [https://kbwiki-images.s3.amazonaws.com/c/c8/ATAAC_D3-2-36_excerpt3DDiffuser.pdf ATAAC_D3-2-36_excerpt3DDiffuser.pdf] (excerpt from an ATAAC report) and
In the framework of two ERCOFTAC-SIG15 Workshops and in the European ATAAC project the flow configuration of the Stanford diffuser was thoroughly investigated. They studied both 3D diffuser configurations (denoted as SIG15 Case 13.2-1 and SIG15 Case 13.2-2, respectively) and were held in Austria (September, 2008) and Italy (September, 2009). Both in the workshops and in the ATAAC project a wide range of turbulence models in both LES and RANS frameworks as well as some novel Hybrid LES/RANS formulations have been employed. The workshop reports are published in the ERCOFTAC Bulletin Issues, see [[lib:DNS_1-3_description#3|Steiner ''et al.''  (2009)]], [[lib:DNS_1-3_description#3| Jakirlić ''et al.''  (2010)]] and the [[UFR_4-16_Test_Case#CFD_Methods| CFD methods section of UFR4-16]]. The ATAAC reports can be found through the links [https://kbwiki-images.s3.amazonaws.com/c/c8/ATAAC_D3-2-36_excerpt3DDiffuser.pdf ATAAC_D3-2-36_excerpt3DDiffuser.pdf] (excerpt from an ATAAC report) and
[https://kbwiki-images.s3.amazonaws.com/a/a1/ATAAC_finalWorkshop_ST04-Diffuser-ANSYS.pdf ATAAC_finalWorkshop_ST04-Diffuser-ANSYS.pdf] (PowerPoint presentation at ATAAC final workshop).
[https://kbwiki-images.s3.amazonaws.com/a/a1/ATAAC_finalWorkshop_ST04-Diffuser-ANSYS.pdf ATAAC_finalWorkshop_ST04-Diffuser-ANSYS.pdf] (PowerPoint presentation at ATAAC final workshop).


The only high-fidelity data available is the DNS performed by [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et al.''  (2010)]]. As with the present case, this DNS had a Re=10000 and was computed using a spectral element code with 11th order polynomials. Their calculations were performed on the Blue Gene/P at ALCF, using 32768 cores and 8 million core hours. Another computation was performed on the cluster “Ekman” at PDC, Stockholm, Sweden, using 2048 cores and a total of 4 million core hours. The flow was computed for 13 flow-through-times (based on the bulk inlet velocity and diffuser length) before gathering statistics, which were gathered over an additional 21 flowthrough-times.
The only high-fidelity data available is the DNS performed by [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et al.''  (2010)]]. The computational domain was designed to have a close agreement with the diffuser geometry, with the corners resulting from the diffuser expansion smoothly rounded with a radius of 6.0. As with the present case, this DNS had a Re=10000 and was computed using a spectral element code with 11th order polynomials. Their calculations were performed on the Blue Gene/P at ALCF, using 32768 cores and 8 million core hours. Another computation was performed on the cluster “Ekman” at PDC, Stockholm, Sweden, using 2048 cores and a total of 4 million core hours. The flow was computed for 13 flow-through-times (based on the bulk inlet velocity and diffuser length) before gathering statistics, which were gathered over an additional 21 flowthrough-times.


= Description of the test case =
= Description of the test case =
The diffuser studied is the [[UFR_4-16_Test_Case|UFR4-16 Test Case]], diffuser 1, provided in the ERCOFTAC database.  
The diffuser studied is the [[UFR_4-16_Test_Case|UFR4-16 Test Case]], diffuser 1, provided in the ERCOFTAC database.  
==Geometry and flow parameters==
==Geometry and flow parameters==
The diffuser shape, dimensions and the coordinate system are shown in [[UFR_4-16_Test_Case#figure3|Fig. 1]] (reproduced from [[UFR_4-16_Test_Case|UFR4-16 Test Case]]).  
The diffuser shape and dimensions as used in the experiment, and the coordinate system are shown in [[lib:DNS_1-3_description_#figure2|Fig. 2]] (reproduced from [[UFR_4-16_Test_Case|UFR4-16 Test Case]]).  


<div id="figure1"></div>
<div id="figure2"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:UFR4-16_figure3.png|740px]]
|[[Image:UFR4-16_figure3.png|740px]]
|-
|-
|'''Figure 1:''' Geometry of the 3-D diffuser 1 considered (not to  scale),  [[UFR_4-16_References#7|Cherry ''et&nbsp;al.'' (2008)]]; see also [[UFR_4-16_References#14|Jakirli&#x107; ''et&nbsp;al.'' (2010a)]].
|'''Figure 2:''' Geometry of the 3-D diffuser 1 considered (not to  scale),  [[UFR_4-16_References#7|Cherry ''et&nbsp;al.'' (2008)]]; see also [[UFR_4-16_References#14|Jakirli&#x107; ''et&nbsp;al.'' (2010a)]].
|}
|}


===Computational Domain===
===Computational Domain===
For the current diffuser, the upper-wall expansion angle is 11.3° and the side-wall expansion angle is 2.56°. The flow in the inlet duct (height h=1, width B=3.33) corresponds to fully-developed turbulent rectangular duct flow. The origin of coordinates is set at the entrance of the diffuser. The <math>L=15h</math> long diffuser section is followed by a straight outlet part (12.5h long). Downstream of this the flow goes through a 10h long contraction followed by a 5h straight duct in order to minimize the effect of the outlet to the diffuser.


A difference from the experiment of Cherry ''et&nbsp;al.'' (2008) and the DNS of Ohlsson ''et&nbsp;al.'' (2010) is that the present geometry does not have any smoothing (by means of a fillet or curvature radius) on the walls transitioning between diffuser and the straight duct parts. The geometry also includes a long inlet duct of 55h long in order to allow the flow in the inlet duct to fully develop. Before this, there is a section of 5h long with a small chevron at 3h from the inlet for 0.1h in order to trigger the turbulent transition on the square duct. An overview of the computational geometry with details on the square duct is shown on [[lib:DNS_1-3_description_#figure2|Fig. 2]].
For the current diffuser, the upper-wall expansion angle is 11.3° and the side-wall expansion angle is 2.56°. The flow in the inlet duct (height h=1, width B=3.33) corresponds to fully-developed turbulent rectangular duct flow. The origin of the x-coordinates (x=0) is set at the entrance of the diffuser and z=0 at the left sidewall. The <math>L=15h</math> long diffuser section is followed by a straight outlet part (12.5h long). Downstream of this the flow goes through a 10h long contraction followed by a 5h straight duct in order to minimize the effect of the outlet to the diffuser. A difference from the previous experimental and DNS works is that the geometry considered does keep the sharp angles on the walls transitioning between diffuser and the straight duct parts. The computational domain also includes a long inlet duct of 65h length in order to allow the flow in the inlet duct to fully develop. Before this, there is a section of 5h length with a small chevron placed 2h from the inlet in order to trigger the turbulent transition in the rectangular duct. An overview of the computational geometry with details of the rectangular duct is shown on [[lib:DNS_1-3_description_#figure3|Fig. 3]].
 
The bulk velocity in the inflow duct is <math>{U_\textrm{bulk}=U_\textrm{inflow}=1\ m/s}</math> in the x-direction resulting in the Reynolds number based on the inlet channel height (<math>{h}</math>) of <math>{Re=\frac{U_\textrm{bulk}h}{\nu} = 10000}</math>.


<div id="figure2"></div>
<div id="figure3"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:DNS_1_3_whole_domain.png|740px]]
|[[Image:DNS_1_3_whole_domain.png|740px]]
|[[Image:DNS_1_3_inlet_detail.png|440px]]
|-
|-
|'''Figure 2:''' Detail of the computational domain used for this test case (left) and detail of the inlet roughness element (right).
|[[Image:DNS_1_3_inlet_detail.png|740px]]
|-
|-
NEW FIGURES FROM BENET OR SAMUEL
|'''Figure 3:''' Computational domain used for this test case sideview top with detail of the inlet roughness element, topview bottom.
|}
|}
The bulk velocity in the inflow duct is <math>{U_\textrm{bulk}=U_\textrm{inflow}=1\ m/s}</math> in the x-direction resulting in the Reynolds number based on the inlet channel height of 10000.


==Boundary conditions==
==Boundary conditions==
The inflow is set to <math>U=1</math> at the inlet of the diffuser (x=-65h). As aforementioned, turbulence is triggered by creating a small discontinuity in the form of a small chevron in the duct. The entrance to the diffuser section at x=5h corresponds to that of a fully developed square duct.
The inflow is set to <math>U=1</math> at the inlet of the inlet duct (x=-70h), setting a Reynolds number of 10000, matching with the previous DNS and experimental data. Turbulence is triggered by creating a small discontinuity in the form of a small chevron in the duct. This method is preferred over using a precursor calculation of rectangular duct flow with streamwise periodicity conditions as [[lib:DNS_1-3_description#6| Nikitin (2008)]] argued that such conditions might not be suitable since they cause a spatial periodicity, which is not physical for turbulent flows. As a matter of fact, [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et&nbsp;al.''  (2010)]] use a 63h long inflow development duct, accounting even for the transition of the initially laminar inflow. The chevron method of this work is aimed to accelerate such transition. At the outlet (x=47.5h) <math>dU/dn = 0</math> and standard Dirichlet condition for the pressure is prescribed. The walls of the duct and the diffuser are set to no slip.


<div id="figure3"></div>
<div id="figure4"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:DNS_1_3_inlet_flow_profile.png|500px]]
|[[Image:DNS_1_3_inlet_flow_profile.png|500px]]
|-
|-
|'''Figure 3:''' Profile of the fully-developed turbulent inlet flow compared with the DNS of [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et&nbsp;al.''  (2010)]] and the experimental data of [[lib:DNS_1-3_description#CFD_Methods|Cherry ''et&nbsp;al.''  (2008)]].
|'''Figure 4:''' Profile of the fully-developed turbulent inlet flow compared with the DNS of [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et&nbsp;al.''  (2010)]] and the experimental data of [[lib:DNS_1-3_description#CFD_Methods|Cherry ''et&nbsp;al.''  (2008)]].
|}  
|}  


<div id="figure4"></div>
<div id="figure5"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:DNS_1_3_uplus_yplus_y.png|1000px]]
|[[Image:DNS_1_3_uplus_yplus_y.png|740px]]
|-
|-
|'''Figure 4:''' Axial velocity profile in semi-log coordinates corresponding to the fully-developed flow in the inflow duct compared with the DNS of [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et&nbsp;al.''  (2010)]].
|'''Figure 5:''' Axial velocity profile in semi-log coordinates corresponding to the fully-developed flow in the inflow duct compared with the DNS of [[lib:DNS_1-3_description#CFD_Methods|Ohlsson ''et&nbsp;al.''  (2010)]]. For the log law B is taken as 5.2.
|}
|}


<div id="figure5"></div>
The development of turbulence along the long inlet rectangular duct is shown in [[lib:DNS_1-3_description_#figure6|Fig. 6]] by plotting the evolution of <math>Re_{\tau}</math> defined as <math>{Re_{\tau} = \frac{u_\tau h/2}{\mu}}</math> for the top and bottom walls and <math>{Re_{\tau} = \frac{u_\tau B/2}{\mu}}</math> for the side walls. The straight line refers to <math>{Re_{\tau}}</math> of a fully developed rectangular duct flow. Eventually, the entrance to the diffuser section at x=0h corresponds to that of a fully developed rectangular duct flow.
 
<div id="figure6"></div>
{|align="center" width=750
{|align="center" width=750
|[[Image:DNS_1_3_retau.png|1000px]]
|[[Image:DNS_1_3_retau.png|740px]]
|-
|-
|'''Figure 5:''' Evolution of the Reynolds number based on the wall friction for the long inlet duct.
|'''Figure 6:''' Evolution of the Reynolds number based on the wall friction velocity for the long inlet duct.
|}  
|}  


As specified in [[UFR_4-16_Evaluation|UFR4-16 Evaluation]], the flow in a duct with rectangular cross-section is no longer unidirectional. It is characterized by a secondary motion with velocity components perpendicular to the axial direction, as shown in [[UFR_4-16_Evaluation#figure22|Fig. 22]] for the UFR4-16 and [[lib:DNS_1-3_description_#figure7|Fig. 7]] for the present case. Briefly, this secondary flow transporting momentum into the duct corners is characterized by jets directed towards the duct walls bisecting each corner with associated vortices at both sides of each jet. They have an influence on the development of the flow in the diffuser.


The outlet at x=47.5h is set to <math>dU/dn = 0</math> while the walls of the duct and the diffuser are set to no slip.
<div id="figure7"></div>
{|align="center" width=750
|[[Image:DNS_1_3_uquiver.png|740px]]
|-
|'''Figure 7:''' Velocity vectors in the y-z plane in the inflow duct at x/h = -2. Points have been sampled at each 10th grid node.
|}


==References==
==References==
Line 76: Line 83:
#<div id="2">'''Ohlsson, J., Schlatter, P., Fischer P.F. and  Henningson,  D.S.  (2010):''' DNS of separated flow in a three-dimensional diffuser  by  the  spectral-element method. ''J.&nbsp;Fluid&nbsp;Mech.'', Vol.&nbsp;650, pp.&nbsp;307&ndash;318</div>
#<div id="2">'''Ohlsson, J., Schlatter, P., Fischer P.F. and  Henningson,  D.S.  (2010):''' DNS of separated flow in a three-dimensional diffuser  by  the  spectral-element method. ''J.&nbsp;Fluid&nbsp;Mech.'', Vol.&nbsp;650, pp.&nbsp;307&ndash;318</div>
#<div id="3">'''Steiner, H., Jakirli&#x107;, S., Kadavelil, G., &#352;ari&#x107;,  S.,  Manceau,  R.  and Brenn. G. (2009):''' Report on 13<sup>th</sup> ERCOFTAC Workshop on Refined  Turbulence Modelling. September&nbsp;25&ndash;26, 2008, Graz University of Technology, ''ERCOFTAC Bulletin'', No.&nbsp;79, pp.&nbsp;24&ndash;29</div>
#<div id="3">'''Steiner, H., Jakirli&#x107;, S., Kadavelil, G., &#352;ari&#x107;,  S.,  Manceau,  R.  and Brenn. G. (2009):''' Report on 13<sup>th</sup> ERCOFTAC Workshop on Refined  Turbulence Modelling. September&nbsp;25&ndash;26, 2008, Graz University of Technology, ''ERCOFTAC Bulletin'', No.&nbsp;79, pp.&nbsp;24&ndash;29</div>
#<div id="4">'''Jakirli&#x107;, S., Kadavelil, G., Sirbubalo, E., von Terzi,  D.,  Breuer,  M. and Borello, D. (2010b):'''  14th  ERCOFTAC  SIG15  Workshop  on  Turbulence Modelling: Turbulent  Flow  Separation  in  a  3-D  Diffuser.  "Sapienza" University of Rome,  September  18,  2009,  ''ERCOFTAC  Bulletin'',  December Issue, No.&nbsp;85, pp.&nbsp;5&ndash;13</div>
#<div id="4">'''Jakirli&#x107;, S., Kadavelil, G., Sirbubalo, E., von Terzi,  D.,  Breuer,  M. and Borello, D. (2010):'''  14th  ERCOFTAC  SIG15  Workshop  on  Turbulence Modelling: Turbulent  Flow  Separation  in  a  3-D  Diffuser.  "Sapienza" University of Rome,  September  18,  2009,  ''ERCOFTAC  Bulletin'',  December Issue, No.&nbsp;85, pp.&nbsp;5&ndash;13</div>
#<div id="5">'''Cherry, E.M., Elkins, C.J. and Eaton, J.K. (2008):''' Geometric  sensitivity of three-dimensional separated flows. ''Int.&nbsp;J.&nbsp;of  Heat  and  Fluid  Flow'', Vol.&nbsp;29(3), pp.&nbsp;803&ndash;811</div>
#<div id="5">'''Cherry, E.M., Elkins, C.J. and Eaton, J.K. (2008):''' Geometric  sensitivity of three-dimensional separated flows. ''Int.&nbsp;J.&nbsp;of  Heat  and  Fluid  Flow'', Vol.&nbsp;29(3), pp.&nbsp;803&ndash;811</div>
#<div id="6">'''Nikitin, N. (2008):''' On the rate of spatial predictability in near-wall turbulence. ''J.&nbsp;of Fluid Mechanics'', Vol.&nbsp;614, pp.&nbsp;495&ndash;507</div>


<br/>
<br/>
Line 85: Line 93:
| organisation= Barcelona Supercomputing Center (BSC)
| organisation= Barcelona Supercomputing Center (BSC)
}}
}}
{{DNSHeaderLib
{{DNSHeader
|area=1
|area=1
|number=3
|number=3

Latest revision as of 09:21, 5 January 2023

Flow in a 3D diffuser

Front Page

Description

Computational Details

Quantification of Resolution

Statistical Data

Instantaneous Data

Storage Format

Introduction

The 3D (Stanford) Diffuser is a well documented case with complex internal corner flow and 3D separation while having a relatively simple geometry. It has an inlet section, an expansion section and an outlet section (see Figure 1). The flow at the inlet is assumed to be a fully developed rectangular channel flow. At the outlet, standard Dirichlet condition for the pressure is prescribed. An inflow Reynolds number of 10000 is considered based on the duct height and the flow is considered to be incompressible. Such configuration is also of engineering relevance. It represents a diffuser situated between the compressor and the combustor chamber of a jet engine, whose task is to decelerate the flow from the compressor to adapt for the combustor section. The following DNS data has been obtained using the in-house Finite Element Method (FEM) code Alya developed at BSC.

Review of previous studies

The physics of the Stanford diffuser were experimentally analyzed by Cherry et al. (2008). The sharp corners of the diffuser were smoothed with a fillet radius of 6.0. Experimental data was gathered in a recirculating water channel using the method of magnetic resonance velocimetry (MRV). In this turbulence causes a loss of net magnetization signal from the studied region. Such decrease in signal strength can be related to turbulent velocity statistics. This method was found to be accurate within 20% of the field of view and within 5% in in regions of high turbulence. Three scans were performed using three different magnetic field gradient intensities. For each one, 30 scans were completed and averaged. Then, the three averaged data sets were averaged to obtain a final data set. Such dataset provided detailed experimental data comprising the mean velocity field on its three components, the streamwise Reynolds stresses and the pressure distribution along the bottom wall of the diffuser. This provided information on the mean flow configuration, however, little insight was given on the more complex time-motions of this flow.

In the framework of two ERCOFTAC-SIG15 Workshops and in the European ATAAC project the flow configuration of the Stanford diffuser was thoroughly investigated. They studied both 3D diffuser configurations (denoted as SIG15 Case 13.2-1 and SIG15 Case 13.2-2, respectively) and were held in Austria (September, 2008) and Italy (September, 2009). Both in the workshops and in the ATAAC project a wide range of turbulence models in both LES and RANS frameworks as well as some novel Hybrid LES/RANS formulations have been employed. The workshop reports are published in the ERCOFTAC Bulletin Issues, see Steiner et al. (2009), Jakirlić et al. (2010) and the CFD methods section of UFR4-16. The ATAAC reports can be found through the links ATAAC_D3-2-36_excerpt3DDiffuser.pdf (excerpt from an ATAAC report) and ATAAC_finalWorkshop_ST04-Diffuser-ANSYS.pdf (PowerPoint presentation at ATAAC final workshop).

The only high-fidelity data available is the DNS performed by Ohlsson et al. (2010). The computational domain was designed to have a close agreement with the diffuser geometry, with the corners resulting from the diffuser expansion smoothly rounded with a radius of 6.0. As with the present case, this DNS had a Re=10000 and was computed using a spectral element code with 11th order polynomials. Their calculations were performed on the Blue Gene/P at ALCF, using 32768 cores and 8 million core hours. Another computation was performed on the cluster “Ekman” at PDC, Stockholm, Sweden, using 2048 cores and a total of 4 million core hours. The flow was computed for 13 flow-through-times (based on the bulk inlet velocity and diffuser length) before gathering statistics, which were gathered over an additional 21 flowthrough-times.

Description of the test case

The diffuser studied is the UFR4-16 Test Case, diffuser 1, provided in the ERCOFTAC database.

Geometry and flow parameters

The diffuser shape and dimensions as used in the experiment, and the coordinate system are shown in Fig. 2 (reproduced from UFR4-16 Test Case).

UFR4-16 figure3.png
Figure 2: Geometry of the 3-D diffuser 1 considered (not to scale), Cherry et al. (2008); see also Jakirlić et al. (2010a).

Computational Domain

For the current diffuser, the upper-wall expansion angle is 11.3° and the side-wall expansion angle is 2.56°. The flow in the inlet duct (height h=1, width B=3.33) corresponds to fully-developed turbulent rectangular duct flow. The origin of the x-coordinates (x=0) is set at the entrance of the diffuser and z=0 at the left sidewall. The long diffuser section is followed by a straight outlet part (12.5h long). Downstream of this the flow goes through a 10h long contraction followed by a 5h straight duct in order to minimize the effect of the outlet to the diffuser. A difference from the previous experimental and DNS works is that the geometry considered does keep the sharp angles on the walls transitioning between diffuser and the straight duct parts. The computational domain also includes a long inlet duct of 65h length in order to allow the flow in the inlet duct to fully develop. Before this, there is a section of 5h length with a small chevron placed 2h from the inlet in order to trigger the turbulent transition in the rectangular duct. An overview of the computational geometry with details of the rectangular duct is shown on Fig. 3.

The bulk velocity in the inflow duct is in the x-direction resulting in the Reynolds number based on the inlet channel height () of .

DNS 1 3 whole domain.png
DNS 1 3 inlet detail.png
Figure 3: Computational domain used for this test case sideview top with detail of the inlet roughness element, topview bottom.

Boundary conditions

The inflow is set to at the inlet of the inlet duct (x=-70h), setting a Reynolds number of 10000, matching with the previous DNS and experimental data. Turbulence is triggered by creating a small discontinuity in the form of a small chevron in the duct. This method is preferred over using a precursor calculation of rectangular duct flow with streamwise periodicity conditions as Nikitin (2008) argued that such conditions might not be suitable since they cause a spatial periodicity, which is not physical for turbulent flows. As a matter of fact, Ohlsson et al. (2010) use a 63h long inflow development duct, accounting even for the transition of the initially laminar inflow. The chevron method of this work is aimed to accelerate such transition. At the outlet (x=47.5h) and standard Dirichlet condition for the pressure is prescribed. The walls of the duct and the diffuser are set to no slip.

DNS 1 3 inlet flow profile.png
Figure 4: Profile of the fully-developed turbulent inlet flow compared with the DNS of Ohlsson et al. (2010) and the experimental data of Cherry et al. (2008).
DNS 1 3 uplus yplus y.png
Figure 5: Axial velocity profile in semi-log coordinates corresponding to the fully-developed flow in the inflow duct compared with the DNS of Ohlsson et al. (2010). For the log law B is taken as 5.2.

The development of turbulence along the long inlet rectangular duct is shown in Fig. 6 by plotting the evolution of defined as for the top and bottom walls and for the side walls. The straight line refers to of a fully developed rectangular duct flow. Eventually, the entrance to the diffuser section at x=0h corresponds to that of a fully developed rectangular duct flow.

DNS 1 3 retau.png
Figure 6: Evolution of the Reynolds number based on the wall friction velocity for the long inlet duct.

As specified in UFR4-16 Evaluation, the flow in a duct with rectangular cross-section is no longer unidirectional. It is characterized by a secondary motion with velocity components perpendicular to the axial direction, as shown in Fig. 22 for the UFR4-16 and Fig. 7 for the present case. Briefly, this secondary flow transporting momentum into the duct corners is characterized by jets directed towards the duct walls bisecting each corner with associated vortices at both sides of each jet. They have an influence on the development of the flow in the diffuser.

DNS 1 3 uquiver.png
Figure 7: Velocity vectors in the y-z plane in the inflow duct at x/h = -2. Points have been sampled at each 10th grid node.

References

  1. Ohlsson, J., Schlatter, P., Fischer P.F. and Henningson, D.S. (2009): DNS of three-dimensional separation in turbulent diffuser flows. In Advances in Turbulence XII, Proceedings of the 12th EUROMECH European Turbulence Conference, Marburg. Springer Proceedings in Physics, Vol. 132, ISBN 978-3-642-03084-0
  2. Ohlsson, J., Schlatter, P., Fischer P.F. and Henningson, D.S. (2010): DNS of separated flow in a three-dimensional diffuser by the spectral-element method. J. Fluid Mech., Vol. 650, pp. 307–318
  3. Steiner, H., Jakirlić, S., Kadavelil, G., Šarić, S., Manceau, R. and Brenn. G. (2009): Report on 13th ERCOFTAC Workshop on Refined Turbulence Modelling. September 25–26, 2008, Graz University of Technology, ERCOFTAC Bulletin, No. 79, pp. 24–29
  4. Jakirlić, S., Kadavelil, G., Sirbubalo, E., von Terzi, D., Breuer, M. and Borello, D. (2010): 14th ERCOFTAC SIG15 Workshop on Turbulence Modelling: Turbulent Flow Separation in a 3-D Diffuser. "Sapienza" University of Rome, September 18, 2009, ERCOFTAC Bulletin, December Issue, No. 85, pp. 5–13
  5. Cherry, E.M., Elkins, C.J. and Eaton, J.K. (2008): Geometric sensitivity of three-dimensional separated flows. Int. J. of Heat and Fluid Flow, Vol. 29(3), pp. 803–811
  6. Nikitin, N. (2008): On the rate of spatial predictability in near-wall turbulence. J. of Fluid Mechanics, Vol. 614, pp. 495–507




Contributed by: Oriol Lehmkuhl, Arnau Miro — Barcelona Supercomputing Center (BSC)

Front Page

Description

Computational Details

Quantification of Resolution

Statistical Data

Instantaneous Data

Storage Format


© copyright ERCOFTAC 2024