CFD Simulations AC2-07: Difference between revisions

From KBwiki
Jump to navigation Jump to search
m (Dave.Ellacott moved page Gold:CFD Simulations AC2-07 to CFD Simulations AC2-07 over redirect)
 
(18 intermediate revisions by 2 users not shown)
Line 1: Line 1:
{{AC|front=AC 2-07|description=Description_AC2-07|testdata=Test Case_AC2-07|cfdsimulations=CFD Simulations_AC2-07|evaluation=Evaluation_AC2-07|qualityreview=Quality Review_AC2-07|bestpractice=Best Practice Advice_AC2-07|relatedUFRs=Related UFRs_AC2-07}}
='''Confined double annular jet'''=
='''Confined double annular jet'''=


Line 14: Line 16:


   
   
<math>\kappa = \frac{1}{2}\left(u^{\prime\;2} +
<math>k = \frac{1}{2}\left(u^{\prime\;2} +
v^{\prime\;2} + w^{\prime\;2}\right)
v^{\prime\;2} + w^{\prime\;2}\right)
\qquad\qquad\qquad\qquad\qquad(2)</math>
\qquad\qquad\qquad\qquad\qquad(2)</math>
Line 23: Line 25:
w&prime;<sup>2</sup> is the tangential turbulence intensity.
w&prime;<sup>2</sup> is the tangential turbulence intensity.


Here we have assumed that v&prime;<sup>2</sup>&nbsp;»&nbsp;>>w&prime;<sup>2</sup>, so that equation (2) becomes:
Here we have assumed that v&prime;<sup>2</sup>&nbsp;=&nbsp;w&prime;<sup>2</sup>, so that equation (2) becomes:


   
   
<math>\kappa = \frac{1}{2}\left(u^{\prime\;2} + 2v^{\prime\;2}\right)
<math>k = \frac{1}{2}\left(u^{\prime\;2} + 2v^{\prime\;2}\right)
\qquad\qquad\qquad\qquad\qquad\qquad(3)</math>
\qquad\qquad\qquad\qquad\qquad\qquad(3)</math>
   
   


The assumption v&prime;<sup>2</sup>&nbsp;»&nbsp;w&prime;<sup>2</sup>
The assumption v&prime;<sup>2</sup>&nbsp;=&nbsp;w&prime;<sup>2</sup>
is more or less true for a fully developed pipe flow.
is more or less true for a fully developed pipe flow.


Line 62: Line 64:


   
   
<math>\[\varepsilon = C_\mu\cdot\frac{\mu}{\mu_{it}}\cdot\frac{\rho_{ref}^{K^2}}{\mu}\]</math>
<math>\varepsilon = C_\mu\cdot\frac{\mu}{\mu_{it}}\cdot\frac{{\rho_{ref}}k^2}{\mu}
\qquad\qquad\qquad\qquad\qquad(4)</math>
   
   


(4)
Where C<sub>&mu;</sub> is 0.09,


Where: - Cm is 0.09
&mu;<sub>it</sub>&nbsp;/&nbsp;&mu; is 1 (for external flows)


- mIt/m is 1. (For external flows)
&mu; is the kinematic viscosity (=1.57 .10-5 m2/s)


- m is the kinematic viscosity (=1.57 .10-5 m2/s)
&rho;<sub>ref</sub> is density of air (=1.23 kg/m3)
 
- rref is density of air (=1.23 kg/m3)


The real outlets of the combustion chamber are 4 gaps, positioned on the outer side of the combustion chamber. To have an axisymmetric grid, an axisymmetric equivalent of the experimental outlet is modeled. The total surface of the equivalent outlet of the numerical model is similar to the one of the empirical model (see Figure 5).
The real outlets of the combustion chamber are 4 gaps, positioned on the outer side of the combustion chamber. To have an axisymmetric grid, an axisymmetric equivalent of the experimental outlet is modeled. The total surface of the equivalent outlet of the numerical model is similar to the one of the empirical model (see Figure 5).
Line 154: Line 155:


{|border="1" cell padding="20" cell spacing="5"
{|border="1" align="center" width="550"
!Segment!!Number pf points!!colspan="3"| Clustering
|+ align="bottom" | <b>Table 9. Clustering of numerical grid.</b>
!align="center" rowspan="2"| Segment
!align="center" rowspan="2"| Number of points
!align="center" colspan="3"| Clustering
|-
|-
! !!  !!Init spacing at start!!Init spacing at end!! Number of constant cells
!align="center" | Init spacing at start
!align="center" | Init spacing at end
!align="center" | Number of constant cells
|-
|-
|1 || 150 || 0.0005 || 0.0 || 120
|align="center"|1 ||align="center"| 150 ||align="center"| 0.0005 ||align="center"| 0.0 ||align="center"| 120
|-
|-
|2 || 240 || 0.0095 || 0.05 || 50
|align="center"|2 ||align="center"| 240 ||align="center"| 0.0095 ||align="center"| 0.05 ||align="center"| 50
|-
|-
|3 || 209 || 0.0005 || 0.05 || 0
|align="center"|3 ||align="center"| 209 ||align="center"| 0.0005 ||align="center"| 0.05 ||align="center"| 0
|-
|-
|4 || 33 || 0.0005 || 0.0005 || 17
|align="center"|4 ||align="center"| 33 ||align="center"| 0.0005 ||align="center"| 0.0005 ||align="center"| 17
|-
|-
|5 || 49 || 0.0 || 0.005 || 20
|align="center"|5 ||align="center"| 49 ||align="center"| 0.0 ||align="center"| 0.005 ||align="center"| 20
|-
|-
|6 || 204 || 0.0095 || 0.05 || 50
|align="center"|6 ||align="center"| 204 ||align="center"| 0.0095 ||align="center"| 0.05 ||align="center"| 50
|-
|-
|7 || 150 || 0.0005 || 0.0 || 120
|align="center"|7 ||align="center"| 150 ||align="center"| 0.0005 ||align="center"| 0.0 ||align="center"| 120
|-
|-
|8 || 81 || 0.0 || 0.0005 || 20
|align="center"|8 ||align="center"| 81 ||align="center"| 0.0 ||align="center"| 0.0005 ||align="center"| 20
|-
|-
|9 || 33 || 0.0005 || 0.0005 || 15
|align="center"|9 ||align="center"| 33 ||align="center"| 0.0005 ||align="center"| 0.0005 ||align="center"| 15
|-
|-
|10 || 65 || 0.0005 || 0.0005 || 40
|align="center"|10 ||align="center"| 65 ||align="center"| 0.0005 ||align="center"| 0.0005 ||align="center"| 40
|-
|-
|11 || 33 || 0.0005 || 0.0005 || 15
|align="center"|11 ||align="center"| 33 ||align="center"| 0.0005 ||align="center"| 0.0005 ||align="center"| 15
|-
|-
|12 || 81 || 0.0005 || 0.0005 || 40
|align="center"|12 ||align="center"| 81 ||align="center"| 0.0005 ||align="center"| 0.0005 ||align="center"| 40
|}
|}
Table 9: Clustering of numerical grid


=='''Physical Modeling'''==
=='''Physical Modeling'''==
Line 216: Line 218:
The fluid modeling is not a critical element in this test case. The fluid model is an incompressible perfect gas.
The fluid modeling is not a critical element in this test case. The fluid model is an incompressible perfect gas.


The Reynolds number based on the diameter of the outer jet, 156 mm and on the equivalent velocity of 2 m/s (equivalent velocity corresponding to a circular jet of diameter 156 mm) is 6.104. These Reynolds number is constant for all turbulent models.
The Reynolds number based on the diameter of the outer jet, 156 mm and on the equivalent velocity of 2 m/s (equivalent velocity corresponding to a circular jet of diameter 156 mm) is 6.104. This Reynolds number is constant for all turbulent models.
 


=='''Initial Solution'''==
=='''Initial Solution'''==
Line 232: Line 233:
Site Design and Implementation: [[Atkins]] and [[UniS]]
Site Design and Implementation: [[Atkins]] and [[UniS]]
        Top              Next
{{AC|front=AC 2-07|description=Description_AC2-07|testdata=Test Case_AC2-07|cfdsimulations=CFD Simulations_AC2-07|evaluation=Evaluation_AC2-07|qualityreview=Quality Review_AC2-07|bestpractice=Best Practice Advice_AC2-07|relatedUFRs=Related UFRs_AC2-07}}

Latest revision as of 15:34, 11 February 2017

Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice

Confined double annular jet

Application Challenge 2-07 © copyright ERCOFTAC 2004


Computational Domain and Boundary Conditions

Experimental measurements show that the overall axisymmetry of the mean flow in the combustion chamber is very good. The whole combustion chamber is modeled as an axisymmetric model, because the boundary of the measured grid is not known.

The measured velocity at the outlet of the primary and the secondary ducts of the burner are imposed as the inflow condition. Figure 9 until 16 shows the LDA measurements of the mean axial and radial velocity profiles at the exit of the primary and second ducts. The turbulent kinetic energy and the turbulent dissipation rate cannot be measured directly and are imposed.

The turbulent kinetic energy are derived from the axial and radial turbulence intensity:



Where u′2 is the axial turbulence intensity, v′2 is the radial turbulence intensity and w′2 is the tangential turbulence intensity.

Here we have assumed that v′2 = w′2, so that equation (2) becomes:



The assumption v′2 = w′2 is more or less true for a fully developed pipe flow.


Image189.jpg Image191.jpg


Figure 9: Axial velocity at primary duct. Figure 10: Axial velocity at secondary duct.


Image193.jpg Image195.jpg


Figure 11: Radial velocity at primary duct. Figure 12: Radial velocity at secondary duct


Image197.jpg Image199.jpg


Figure 13: Turbulent kinetic energy at primary duct. Figure 14: Turbulent kinetic energy at secondary duct.


Image201.jpg Image203.jpg


Figure 15: Turbulent dissipation rate at primary duct. Figure 16: Turbulent dissipation rate at secondary duct.


The value of the turbulent dissipation is specified by the ratio of the turbulent viscosity to the laminar viscosity (equation 4) and is on the basis of the k profile:



Where Cμ is 0.09,

μit / μ is 1 (for external flows)

μ is the kinematic viscosity (=1.57 .10-5 m2/s)

ρref is density of air (=1.23 kg/m3)

The real outlets of the combustion chamber are 4 gaps, positioned on the outer side of the combustion chamber. To have an axisymmetric grid, an axisymmetric equivalent of the experimental outlet is modeled. The total surface of the equivalent outlet of the numerical model is similar to the one of the empirical model (see Figure 5).


Image156.gif



Figure 17: Axisymmetric equivalent of the outlets of the combustion chamber.

The numerical model is very sensitive for negative mass flow at the outlet. The geometry to the outlet has to be carefully dimensioned in order to avoid backflows at the downstream boundary. Increasing the surface of the outlet solves this problem.

Figure 18 shows the stream-function contours calculated from the numerical velocity calculated with the Spalart Allmaras turbulence model (later more about the turbulence models).


Image208.jpg


Figure 18: The results of the stream function contours.

The solid boundaries are modeled as adiabatic walls with no-slip condition for the velocity field, except the symmetric axis.

The static pressure is prescribed to be atmospheric at the outlet.

Discretization and Grid Resolution

The empirical study of the structure of the measured flow identifies several zones, depending on the axial distance from the nozzle. The nozzle region can be separated in 4 zones, which we consider below. The different zones are illustrated in the following Figure 19.


Image157.gif


Figure 19 : Velocity vectors on some grid lines. Their amplitude gives the color.


The different zones are:

- For : there are two recirculation regions, one around the axis, and the other annular, between the two streams.

- For : there is one recirculation region, around the axis, with around it a complex annular jet, with a local minimum of its velocity.

- For : there is one recirculation region, around the axis, with around it an annular jet. This axisymmetric jet behaves as if there was only one stream.

- For : there is one axisymmetric jet, with a complex shape possessing a local minimum at the center.

- For , we are outside the nozzle region, and the burner behaves as if there was only a simple axisymmetric jet.

The numerical model is a structured mesh of one block, divided in 2 zones: An area with a high spatial clustering of grid points in the recirculation and the vortices regions, and a mild clustering at the end of the combustion chamber. Stretching is applied towards the boundary of the outlet. A study of the grid sensibility is done and discussed in paragraph 4.2. The details of the final grid, which is used for the numerical simulations, are resumed in Table 9. The details of the grid, and the number of points of the segment are shown in Figure 20 and 21. The way of clustering of the segments is given in Table 9. The total number of grid points in the longitudinal direction is 353 and in the transversal direction 289. The number of grid levels is (5/5/0).

Image158.gif



Image229.jpg



Figure 20: Grid of numerical model (grid level 2/2/0)


Image159.gif



Image233.jpg


Figure 21: Detail grid at inlet of numerical model (grid level 1/1/0)


Table 9. Clustering of numerical grid.
Segment Number of points Clustering
Init spacing at start Init spacing at end Number of constant cells
1 150 0.0005 0.0 120
2 240 0.0095 0.05 50
3 209 0.0005 0.05 0
4 33 0.0005 0.0005 17
5 49 0.0 0.005 20
6 204 0.0095 0.05 50
7 150 0.0005 0.0 120
8 81 0.0 0.0005 20
9 33 0.0005 0.0005 15
10 65 0.0005 0.0005 40
11 33 0.0005 0.0005 15
12 81 0.0005 0.0005 40

Physical Modeling

The confined double annular burner is accurately modeled by a time-marching finite volume Navier-Stokes solver. The three dimensional compressible Reynolds averaged continuity, turbulent Navier-Stokes with different turbulence models, and energy equations are solved using the structured meshes.

The system of equations is integrated in time using an explicit four-stage Runga-Kutta scheme and convergence towards steady state is accelerated using multigrid, combined with residual smoothing and local time stepping.

The convective fluxes are treated through a second-order central Jameson scheme with second- and fourth-order scalar dissipation, and a pseudo-compressibility method and a second-order dual-time stepping procedure are used.

The Full multigrid strategy is based on the V cycle. The system of equations is integrated in time using an explicit four–stage Runga-Kutta scheme and convergence towards steady-state is accelerated using multigrid combined with residual smoothing and local time stepping with a Courant-Friedrichs-Lewy number of three.

The preconditioning parameters β and the characteristic velocity are suggested as three and the maximum velocity in the flow yield (It is the maximum flow velocity at the inlet of the primary duct).

Turbulence are modeled by the next models:

- Baldwin-Lomax algebraic turbulence model

- Spalart-Allmaras turbulence model

- k-ε turbulence model of Chien (Low Reynolds)

- k-ε turbulence model of Launder-Sharma (Low Reynolds)

- k-ε turbulence model of Yang-Shih (Low Reynolds)

All the quoted turbulence models were used for the CFD study, and compared with the experimental results. It will be discussed in paragraph 4: “Evaluation-comparison of test data and CFD”.

The fluid modeling is not a critical element in this test case. The fluid model is an incompressible perfect gas.

The Reynolds number based on the diameter of the outer jet, 156 mm and on the equivalent velocity of 2 m/s (equivalent velocity corresponding to a circular jet of diameter 156 mm) is 6.104. This Reynolds number is constant for all turbulent models.

Initial Solution

The initial solution is generated from a uniform initial field. A constant pressure field higher than the imposed pressure at the outlet is chosen, and a constant non-zero axial velocity is imposed to prevent a negative mass flows at the outlet of the combustion chamber.

© copyright ERCOFTAC 2004



Contributors: Charles Hirsch; Francois G. Schmitt - Vrije Universiteit Brussel

Site Design and Implementation: Atkins and UniS


Front Page

Description

Test Data

CFD Simulations

Evaluation

Best Practice Advice