UFR 4-13 Evaluation: Difference between revisions

From KBwiki
Jump to navigation Jump to search
(New page: {{UFR|front=UFR 4-13|description=UFR 4-13 Description|references=UFR 4-13 References|testcase=UFR 4-13 Test Case|evaluation=UFR 4-13 Evaluation|qualityreview=UFR 4-13 Quality Review|bestp...)
 
m (Dave.Ellacott moved page Gold:UFR 4-13 Evaluation to UFR 4-13 Evaluation over redirect)
 
(6 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{UFR|front=UFR 4-13|description=UFR 4-13 Description|references=UFR 4-13 References|testcase=UFR 4-13 Test Case|evaluation=UFR 4-13 Evaluation|qualityreview=UFR 4-13 Quality Review|bestpractice=UFR 4-13 Best Practice Advice|relatedACs=UFR 4-13 Related ACs}}
{{UFR|front=UFR 4-13|description=UFR 4-13 Description|references=UFR 4-13 References|testcase=UFR 4-13 Test Case|evaluation=UFR 4-13 Evaluation|qualityreview=UFR 4-13 Quality Review|bestpractice=UFR 4-13 Best Practice Advice|relatedACs=UFR 4-13 Related ACs}}


Line 19: Line 18:
The calculations presented here were performed using STAR-CD, a transient, multi-phase and multidimensional CFD code, now widely used for internal combustion engine simulations. They were performed by Imperial College within a European Research program. Validation was based on the IMFT square piston test case.
The calculations presented here were performed using STAR-CD, a transient, multi-phase and multidimensional CFD code, now widely used for internal combustion engine simulations. They were performed by Imperial College within a European Research program. Validation was based on the IMFT square piston test case.


 
The first simulations presented here where performed with an operational guillotine, producing a compression ratio of about 10. The fine 2D grid of 8000 cells was applied with RNG k-ε and the Launder Gibson-Reynolds Stress Transport Model (LG-RSTM). Previous results obtained with coarse mesh with RNG turbulence model are also presented.
 
The first simulations presented here where performed with an operational guillotine, producing a compression ratio of about 10. The fine 2D grid of 8000 cells was applied with RNG k-<span><font face="Symbol">e </font></span>and theLaunder Gibson-Reynolds Stress Transport Model (LG-RSTM). Previous results obtained with coarse mesh with RNG turbulence model are also presented.
 
 
 
During the intake at 90 Crank Angle Degrees (CA) After Top Dead Centre (ATDC), the axial velocity using the RNG k-<span><font face="Symbol">e </font></span>model, is qualitatively close to the experimental data. During the compression and especially at the end of the compression (at 270 CA. ATDC) when of the tumble vortex breaks down, The RNG k-<span><font face="Symbol">e </font></span>shows a poor comparison with experimental data. Finer meshing produces better prediction.


 
During the intake at 90 Crank Angle Degrees (CA) After Top Dead Centre (ATDC), the axial velocity using the RNG k-&epsilon; model, is qualitatively close to the experimental data. During the compression and especially at the end of the compression (at 270 CA. ATDC) when of the tumble vortex breaks down, The RNG k-&epsilon; shows a poor comparison with experimental data. Finer meshing produces better prediction.


The flow field predicted with the Launder Gibson Reynolds Stress Transport Model (LG-RSTM) are shown in the Figures 8, 9 and 10. The predicted axial velocity during the intake (at 90 deg ATDC) is less accurate than k-epsilon results. On the contrary, during compression, the LG-RSTM model shows better results than the k-epsilon, cf. Figure 7. A fairly large discrepancy is observed near the wall in the case of the Reynolds stress model; this problem is shared by the other numerical results but to a lesser extent. This defect occurs also during the intake and the compression strokes.
The flow field predicted with the Launder Gibson Reynolds Stress Transport Model (LG-RSTM) are shown in the Figures 8, 9 and 10. The predicted axial velocity during the intake (at 90 deg ATDC) is less accurate than k-epsilon results. On the contrary, during compression, the LG-RSTM model shows better results than the k-epsilon, cf. Figure 7. A fairly large discrepancy is observed near the wall in the case of the Reynolds stress model; this problem is shared by the other numerical results but to a lesser extent. This defect occurs also during the intake and the compression strokes.
Line 33: Line 26:
LES calculations were carried out with an open guillotine at a compression ratio with a longer intake duct and random stream-wise fluctuations imposed in the inlet. A coarse 3D grid of about 90K cells was used for a one-cycle simulation. Undeniably, LES calculations give better agreement with experiments during both the intake and the compression even with coarse mesh and a single cycle, see Figures 6,7 and 11.
LES calculations were carried out with an open guillotine at a compression ratio with a longer intake duct and random stream-wise fluctuations imposed in the inlet. A coarse 3D grid of about 90K cells was used for a one-cycle simulation. Undeniably, LES calculations give better agreement with experiments during both the intake and the compression even with coarse mesh and a single cycle, see Figures 6,7 and 11.


 


 
 
 
 
 
 
 
 
 
 
 
<span><font face="&quot;Times New Roman&quot;"><font size="11.0pt"><span class="msoChangeProp"><br style="page-break-before: always" clear="all" /></span></font></font></span>
<center> </center>
<center></center>
{| align="left"
| width="74" height="11" |
|-
|
|
[[Image:U4-13d32_files_image064.gif]]
[[Image:U4-13d32_files_image064.gif]]
|}
 
<center> </center>
<center> </center>
<center> </center>
<center> </center>
<center> </center>


 


 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
{| align="left"
| width="127" height="11" |
|-
|
|
[[Image:U4-13d32_files_image066.gif]]
[[Image:U4-13d32_files_image066.gif]]
|}


 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 


<br style="page-break-before: always" clear="all" /><br clear="ALL" />
<br style="page-break-before: always" clear="all" /><br clear="ALL" />
Line 195: Line 55:
Figure (10)- LGRST Predicted flow field
Figure (10)- LGRST Predicted flow field


with compression tumble CR=10 and at 310 ATDC<br /><br style="page-break-before: always" clear="all" />
with compression tumble CR=10 and at 310 ATDC<br />
 
<br style="page-break-before: always" clear="all" />


<span style="position: absolute; z-index: 6; left: 0px; margin-left: 309px; margin-top: 104px; width: 148px; height: 33px"> </span>
<span style="position: absolute; z-index: 6; left: 0px; margin-left: 309px; margin-top: 104px; width: 148px; height: 33px"> </span>
Line 206: Line 68:
<div class="shape" style="padding: 3.6pt 7.2pt 3.6pt 7.2pt">
<div class="shape" style="padding: 3.6pt 7.2pt 3.6pt 7.2pt">


<span><font size="11.0pt">Fig (11-a)</font></span>
Fig (11-a)


</div>
</div>
Line 215: Line 77:


[[Image:U4-13d32_files_image071.gif]]
[[Image:U4-13d32_files_image071.gif]]
 


<span style="position: absolute; z-index: 7; left: 0px; margin-left: 315px; margin-top: 40px; width: 148px; height: 33px"> </span>
<span style="position: absolute; z-index: 7; left: 0px; margin-left: 315px; margin-top: 40px; width: 148px; height: 33px"> </span>
Line 227: Line 87:
<div class="shape" style="padding: 3.6pt 7.2pt 3.6pt 7.2pt">
<div class="shape" style="padding: 3.6pt 7.2pt 3.6pt 7.2pt">


<span><font size="11.0pt">Fig (11-b)</font></span>
Fig (11-b)


</div>
</div>
Line 236: Line 96:


[[Image:U4-13d32_files_image072.gif]]
[[Image:U4-13d32_files_image072.gif]]
 


Fig (11) Predicted flow field in longitudinal mid-plane during induction and exhaust for non compressing tumble case CR = 4, obtained using low Reynolds number subgrid model
Fig (11) Predicted flow field in longitudinal mid-plane during induction and exhaust for non compressing tumble case CR = 4, obtained using low Reynolds number subgrid model
 


<br style="page-break-before: always" clear="all" />
<br style="page-break-before: always" clear="all" />
Line 253: Line 109:
The simulations also showed that using the 3rd order TTGC (a finite element two-step Taylor-Galerkin convective scheme) , allowed for a given mesh resolution, to resolve more flow details near the cut-off scale than a 2nd order centred Lax-Wendroff type scheme, a definitive advantage for LES simulations. The mean velocity profiles obtained by phase averaging six consecutive 3D engine cycles with a standard Smagorinsky LES turbulence model combined to wall laws and using a LW scheme showed encouraging reproduction of mean experimental profiles ( cf. Figure 13).
The simulations also showed that using the 3rd order TTGC (a finite element two-step Taylor-Galerkin convective scheme) , allowed for a given mesh resolution, to resolve more flow details near the cut-off scale than a 2nd order centred Lax-Wendroff type scheme, a definitive advantage for LES simulations. The mean velocity profiles obtained by phase averaging six consecutive 3D engine cycles with a standard Smagorinsky LES turbulence model combined to wall laws and using a LW scheme showed encouraging reproduction of mean experimental profiles ( cf. Figure 13).


 
 
 
 
</div>


<br style="page-break-before: auto" clear="all" />
<br style="page-break-before: auto" clear="all" />
Line 267: Line 114:
<div class="Section4">
<div class="Section4">


[[Image:U4-13d32_files_image074.jpg]]
<center>[[Image:U4-13d32_files_image074.jpg]]</center>


<center>At 90°CA BTDC</center>
<center>At 90&deg;CA BTDC</center>


<br clear="all" />
<br clear="all" />


<center></center>


[[Image:U4-13d32_files_image076.jpg]]
<center>[[Image:U4-13d32_files_image076.jpg]]</center>


 


<center>at TDC </center>
<center>at TDC </center>
 


</div>
</div>
Line 291: Line 134:
<center>Figure 12: mean velocity profiles at the compression phase, Moureau ''et al''. 200</center>
<center>Figure 12: mean velocity profiles at the compression phase, Moureau ''et al''. 200</center>


 


<center>[[Image:U4-13d32_files_image078.jpg]]</center>
<center>[[Image:U4-13d32_files_image078.jpg]]</center>


<center>Figure 13: instantaneous velocity fields of the compression, 45° Crank Angles (CA) Before TDC (BTDC), obtained with Lax-Wendroff (left) and TTGC (right) schemes on 2D meshes.</center>
<center>Figure 13: instantaneous velocity fields of the compression, 45&deg; Crank Angles (CA) Before TDC (BTDC), obtained with Lax-Wendroff (left) and TTGC (right) schemes on 2D meshes.</center>


<center>The velocity normalization is the same. (Moureau ''et al''. 2004)</center>
<center>The velocity normalization is the same. (Moureau ''et al''. 2004)</center>
Line 307: Line 149:


{{UFR|front=UFR 4-13|description=UFR 4-13 Description|references=UFR 4-13 References|testcase=UFR 4-13 Test Case|evaluation=UFR 4-13 Evaluation|qualityreview=UFR 4-13 Quality Review|bestpractice=UFR 4-13 Best Practice Advice|relatedACs=UFR 4-13 Related ACs}}
{{UFR|front=UFR 4-13|description=UFR 4-13 Description|references=UFR 4-13 References|testcase=UFR 4-13 Test Case|evaluation=UFR 4-13 Evaluation|qualityreview=UFR 4-13 Quality Review|bestpractice=UFR 4-13 Best Practice Advice|relatedACs=UFR 4-13 Related ACs}}
[[Category:Underlying Flow Regime]]

Latest revision as of 14:32, 12 February 2017

Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References




Compression of vortex in cavity

Underlying Flow Regime 4-13               © copyright ERCOFTAC 2004


Evaluation

Comparison of CFD calculations with Experiments

The calculations presented here were performed using STAR-CD, a transient, multi-phase and multidimensional CFD code, now widely used for internal combustion engine simulations. They were performed by Imperial College within a European Research program. Validation was based on the IMFT square piston test case.

The first simulations presented here where performed with an operational guillotine, producing a compression ratio of about 10. The fine 2D grid of 8000 cells was applied with RNG k-ε and the Launder Gibson-Reynolds Stress Transport Model (LG-RSTM). Previous results obtained with coarse mesh with RNG turbulence model are also presented.

During the intake at 90 Crank Angle Degrees (CA) After Top Dead Centre (ATDC), the axial velocity using the RNG k-ε model, is qualitatively close to the experimental data. During the compression and especially at the end of the compression (at 270 CA. ATDC) when of the tumble vortex breaks down, The RNG k-ε shows a poor comparison with experimental data. Finer meshing produces better prediction.

The flow field predicted with the Launder Gibson Reynolds Stress Transport Model (LG-RSTM) are shown in the Figures 8, 9 and 10. The predicted axial velocity during the intake (at 90 deg ATDC) is less accurate than k-epsilon results. On the contrary, during compression, the LG-RSTM model shows better results than the k-epsilon, cf. Figure 7. A fairly large discrepancy is observed near the wall in the case of the Reynolds stress model; this problem is shared by the other numerical results but to a lesser extent. This defect occurs also during the intake and the compression strokes.

LES calculations were carried out with an open guillotine at a compression ratio with a longer intake duct and random stream-wise fluctuations imposed in the inlet. A coarse 3D grid of about 90K cells was used for a one-cycle simulation. Undeniably, LES calculations give better agreement with experiments during both the intake and the compression even with coarse mesh and a single cycle, see Figures 6,7 and 11.


U4-13d32 files image064.gif


U4-13d32 files image066.gif




U4-13d32 files image068.jpg

Figure (8) - LGRST Predicted flow field and axial velocity with

compression tumble CR=10 and at 90 ATDC

 

U4-13d32 files image069.gif

Figure (9)- LGRST Predicted flow field and axial velocity

with compression tumble CR=10 and at 270 ATDC

 

U4-13d32 files image070.gif

Figure (10)- LGRST Predicted flow field

with compression tumble CR=10 and at 310 ATDC


Fig (11-a)

 

U4-13d32 files image071.gif

Fig (11-b)

 

U4-13d32 files image072.gif

Fig (11) Predicted flow field in longitudinal mid-plane during induction and exhaust for non compressing tumble case CR = 4, obtained using low Reynolds number subgrid model


More recently, LES developments were reported in Moureau et al. 2004 and Devesa et al. 2004 with the square piston as a validation test case. The unstructured parallel AVBP CFD code was used to predict cyclic variability occurring in Internal Combustion (IC) engines.

An Arbitrary Lagrangian Eulerian (ALE) method was implemented with 2nd order finite volume and 3rd order finite element convective schemes with moving solid boundaries. LES modelling validation was performed by comparing simulations results obtained on 2D and 3D meshes and using two convective schemes with experimental data on the square piston experiment.

The simulations clearly showed that even if the mean flow is essentially 2D during the intake phase, only a full 3D simulation is able to capture the breakdown of the tumble during piston compression as shown in Figure 12.

The simulations also showed that using the 3rd order TTGC (a finite element two-step Taylor-Galerkin convective scheme) , allowed for a given mesh resolution, to resolve more flow details near the cut-off scale than a 2nd order centred Lax-Wendroff type scheme, a definitive advantage for LES simulations. The mean velocity profiles obtained by phase averaging six consecutive 3D engine cycles with a standard Smagorinsky LES turbulence model combined to wall laws and using a LW scheme showed encouraging reproduction of mean experimental profiles ( cf. Figure 13).



U4-13d32 files image074.jpg
At 90°CA BTDC



U4-13d32 files image076.jpg


at TDC


Figure 12: mean velocity profiles at the compression phase, Moureau et al. 200


U4-13d32 files image078.jpg
Figure 13: instantaneous velocity fields of the compression, 45° Crank Angles (CA) Before TDC (BTDC), obtained with Lax-Wendroff (left) and TTGC (right) schemes on 2D meshes.
The velocity normalization is the same. (Moureau et al. 2004)

© copyright ERCOFTAC 2004



Contributors: Afif Ahmed - RENAULT


Front Page

Description

Test Case Studies

Evaluation

Best Practice Advice

References