Transcript Large Eddy Simulation of Impinging Jets with Heat Transfer
Accuracy
of present data
1
u rms
Grid dependency (I)
at r/D=0.5 u rms at r/D=1 U at r/D=1 No. cells on M-G level 1 (
): 194 x 98 x 194 (A) (
): 194 x 146 x 194 (B) (
): 194 x 226 x 194 (C)
N
N
(Re 1 / 2 ) (Re 3 / 4 ) (
U
3 / )
u rms
• • • •
Re=20000
LES: ratio≈141 ≈0.007D
( : average eddy size) DNS: ratio≈1700 • •
u rms U
The grid refinement study shows a monotonic behavior.
Doubling the number of cells in y does not affect the general flow character.
The character of the wall jet is only weakly influenced by four times higher wall normal resolution.
Additional studies have shown that the growth of instabilities at the nozzle outlet is very little influenced by the axial and radial grid resolution.
The present grid (B) provides a grid independent solution, within appropriate error limit, with respect to the wall-normal grid resolution. The important large-scale dynamics is not affected by the finite grid size. 2
•
Grid dependency (II)
As the grid spacing increases with the radius, the flow becomes successively less resolved.
This has a negative influence on the transition process in the wall jet. • • • Accurate results within this work is only obtained within a certain limit on r/D. This limit has not exactly been determined (significant cost to refine the grid) .
May contribute to the lack of a second maximum in Nu for H/D=2.
U at r/D=1 k at y/D=0.15
H/D=4 H/D=4 U No. cells on M-G level 1 (
): 194 x 194 x 194
(
): 170 x 194 x 170
k
3
Grid stretching in y (H/D=2)
No. cells on M-G level 1 (
): 194 x 98 x 194 (A) (
): 194 x 146 x 194 (B) (
): 194 x 226 x 194 (C) y/D
4
Physical
t
(Re 1 / 2 )
t
Numerical
h
/
U N LES
t
10
h
10
t
t
10 10
t h
/
U
10 (Re 1 / 2 )
Convergence error
Signal: U within the wall jet at r/D=1.5
(
): Mean velocity (
): RMS velocity
n
|| || ||
Q Q Q n n n
1
Q n
1
Q Q n n
1 || || || 1
Re=20000
LES: steps≈1400 DNS: steps≈17000 for 10 cycles • • • •
t
Convergence error below 10 -3 , fallen approx. three magnitudes accuracy of approx 0.1.
Doubling of the total sampling time Negligible change of 1 st order statistics.
Second order slightly higher change, not influencial.
5
Model errors and errors due to the boundary conditions
(
): LES (pipe) (
): LES (top-hat) (О): Cooper et al.
(
): Geers et al.
• •
U u rms
With correct boundary conditions the results correlate well with experimental data ("same" boundary conditions) .
The model/scheme provides a physically relevant solution.
6
PSD
Model error
Velocity spectrum at r/D=0.5, y/D=1 for H/D=4
• •
St
The flow is not fully developed but fairly close to that state.
The spectrum shows a behavior close to the -5/3 law within approx. one magnitude of frequencies.
7
Additional observations Additional to the figures shown:
• • • • • The frequency of the natural mode agrees with theory and experiments.
The frequency of the jet-column mode agrees with experimental data for H/D=4. Growth of instabilities exhibit physically correct behavior, i.e. the model/scheme provides appropriate amount of dissipation. From visual observations the dynamical behavior of the jet agrees with expermental observations.
The near-wall resolution is appropriate in terms of describing turbulent mechanisms (see e.g. near-wall peak of 1 st order statistics) .
8
SGS-models
9
Smagorinsky model
• • • Boussinesq hypothesis:
ij
1 3
ij
kk
2
S ij
The SGS-scales can be described by a length- and a time-scale
SGS
2 /
t
, The second important parameter at cut-off is the dissipation:
SGS
1 / 3 4 / 3 • • Assuming the equilibrium hypothesis Production = Dissipation: 2
SGS S ij S ij
SGS
(
C s
) 2 ( 2
S ij S ij
) 1 / 2
C
2 |
S
| 0 Considers the energetic action and neglects the structural information.
10
SSM
• Idea: The difference between the once and twice filtered field is related to the unresolved scales.
• •
ij
Simplest division: 1) the largest resolved scales, 2) the smallest resolved scales, 3) the unresolved scales.
AB
A B
if filter is a Reyn. Op. terms cancels out.
1 3
ij
ij
kk L ij
C b
(
u i u j
C ij
R ij
u i u j
) (
u i u j
1 3
u i
ij
u j
)
kk
C b
(
u i u
'
j
(
L ij
u
R
approx
.
ij
C ij j u i
' ) ),
C b
(
u i
'
u
'
j
),
u
1
G.I.
u
u
' • • •
in SSM : R ij
C ij
(
u i u j
u i u j
) L: large-scale interactions, C: cross-scale interactions, R: SGS interaction In Bardina model: L is computed exactly, C is modeled and R=0 non-dissipative combine linearly with the Smagorinsky:
R ij
2
C
2 |
S
|
S ij u i u j
(
u i
u i
' )(
u j
u
'
j
)
u i u j
u i u
'
j
u j u i
'
u i
'
u
'
j
11
Dynamic approach
• • • • • In order to adapt the model to the structure of the flow the model constant is dynamically calculated. Applicable to any model that explicitly uses an arbitrary constant.
Assumption: the SGS-stresses have the same asymptotic behavior at different filter widths.
Germano identity: the SGS and STS stress are related through the resolved stress
L
ij
T
ij
~
ij
,
where
ij
L
ij
C
ij
R
ij
u
i
u
j
u
i
u
j
T
ij
u
i
u
j
~ ~
u
i
u
j
,
(upper represents ~) L
ij
u
i
u
j
~
u
i
~
u
j
,
(upper represents ~)
Assume that the two SGS tensors can be modeled by the same constant
ij
1 3
ij
kk
C
ij T ij
1 3
ij T kk
C
ij
12
Dynamic approach
• • • • • • Using the same SGS model: assumption of scale-similarity. Smagorinsky:
ij
2
C
2 |
S
|
S ij
ij
2
C
~ 2 | ~
S
| ~
S ij
, ~ 2
L ij
1
L kk
ij
C
ij
C
ij
,
(overbar represents ~)
3
C
slowly varying function in space.
ij
~
ij
,
(overbar represents ~)
C is computed such as the introduced error becomes minimal.
C can take negative values account for backscatter.
Averaging or clipping of C to avoid numerical instability since the denominator may become zero.
13
LES
14
Filtering
u
(
x
)
G
(
x
x
' )
u
(
x
' )
d
3
x
'
G
*
u
Properties
• • • Linear:
A
B
A
B
Conservation of constant: (
x
Cummutation with derivatives:
x
' )
d
3
x
' 1 ,
C
x
,
G
*
u
C
x u
x u
0 15
LES-errors
16
Errors
• • • • A second order filter introduces commutation error of second order.
A fourth-order scheme using fourth-order commutative explicit filters with a filter width of at least twice the grid-spacing adds a numerically clean environment, suitable for SGS model evaluation.
Due to the effect from the modified wave-number, aliasing is less important at high wave-numbers for FD than for spectral methods.
However, the representation at high wave-numbers for FD methods is less accurate because of the large truncation errors.
17
Commutation errors
• • • Commutation errors are introduced 1.
Due to the non-linear term (also for uniform filters, meshes).
2.
Due to non-uniform filtering, meshes (also the linear terms).
The first commutation error is accounted for by the SGS model the second group can also be accounted for but this adds significant complexity to the model.
Error
x
,
G
*
u
x u
x u
The commutation errors can be made smaller by using high order schemes and by having a smooth varying grid.
18
Galilean invariance
19
Galilean invariance
• •
t
If the method is not GI additional terms (errors) are introduced under translation, that must be accounted for.
x
* *
t x
Vt
b
*
L ij
u i
*
u
*
j
u i
*
u
*
j
L ij
(
V i u j
V j u i
) (
V i u j
V j u i u
*
u
V
L ij
(
V i u
'
j
V i u
'
j
),
not GI
x k
*
t
*
x k
V k t
x k C ij
*
R ij
* *
L ij
C ij
R ij
C ij
* (
V i u
'
j
L ij
V i u
'
j
),
not GI
C ij p
*
p
ij
* *
L ij
C ij
*
R ij
*
L ij
C ij
R ij
ij
) The N-S are GI so also the filtered N-S with the tripple decomposition of the SGS stress.
20
Numerical Errors
21
• • •
Modified wavenumber
The formal order of accuracy of FD scheme describes the behavior for fairly well solved features. For scales that are described by only a few grid points the wave-# is replaced by a modified wave-#.
u
(
x
)
u
ˆ (
k
) exp(
ikx
)
u
ˆ (
k
) exp(
ikx
)
dk
With FD the derivative become:
du dx
For spectral methods k*=k
du dx
k
max
k
max
ik
ik u
ˆ (
k
*
u
ˆ (
k
) ) exp(
ikx
)
dk
exp(
ikx
)
dk
• •
k c =k * max
• • • •
k
max
k
* max /
x
/( 2
x
) Any UW-scheme is equiv. to the next higher CD-scheme.
Third order UW (Kawamura et al.):
k
* sin(
k
x
)( 4 3
x
cos(
k
x
))
i
4 sin 4 (
k x
x
/ 2 ) k*: the actual wave-# that is resolved on the smallest resolved scales. The energetic scales must be smaller than k c .
large influence Im(k*): dissipation (CDO6+FilterO4 = UWO3) The Im-part change the amplitude of u.
Due to the modified wave-# the effect from aliasing is of less importance for FD than for spectral. However, the high wave-# for FD is influenced by large amount of T.E. less accurate (Chow & Moin, JCP, No.184, 203).
22
• • • CDO2:
diff
(
e ikx
Modified wavenumber
)
e ikx
x
e ik
(
x
x
)
e
2
x ik
(
x
x
)
i
sin(
k
x
)
x e ikx
ie ikx
k
* sin(
k
x
x
) For low wave-# (smooth functions) (Ferziger & Peric) :
k
* sin(
k
x
)
x
k
k
3 (
x
) 2
smal l k,
x
6
k *
k
For UW-schemes the modified wave-# is complex, type of approximation (Ferziger & Peric) .
k
* 1
e
ik
x
x
23
• • • • • • • •
Numerical dissipation
Sagaut: The numerical and the subgrid dissipation is correlated in space.
A SGS viscosity model: second order dissipation, associated with a spectrum: (
k
/
k c
) 2
E
(
k
) A nth-order numerical dissipation is associated with a spectrum: (
k
/
k
Third order dissipative scheme: the numerical dissipation will be largest for high wave-#. SGS-dissipation be largest for small wave-#.
c
)
n
The numerical dissipation increases exp. as k k c , i.e. as k
k
max /
E
(
k
x
) Numerical dissipative schemes may give to steep spectrum at large k: se pipe simulation.
The dissipative role of the SGS terms can be replaced by a numerical one if it does not affect the energy level of the larger scales.
For fine enough resolution the effects of the SGS scales (backscatter etc) does not affect the large scales.
With implicit modeling the effects from the SGS scales cannot be accounted for in a physical manner. ‘Non-physical’ related backscatter may occur due to the dispersive behavior of the scheme, i.e. not monotone.
24
• • • • • •
Truncation error
For linear problems the truncation error acts as a source of the discretization error, which is convected and diffused by the operator, in this case the NS-operator. The T.E. cannot be computed exactly since the exact solution is not known.
PDE
FDE
d h
Exact analysis is not possible for non-linear problems. Small enough errors linearization. An approximation to the T.E. can be found from successive refinement of the grid. Guide to where grid refinement is needed.
For sufficiently fine grids the T.E. is prop. to the leading term in the Taylor series.
The error can be estimated as,
d h
h
F
p
1
Fh p
log
Fh h
log
F
2
Fh Fh
where F is the grid coarsening ratio and p the order of the scheme.
When solutions on several grids are available, an the convergence is monotonic, one can use Richardson extrapolation to find a better approximation to the solution on the finest grid h:
h
d h
25
• •
Truncation error
For the present UWO3 scheme: (
uu
x
)
d
uu
x
1 12 |
u
|
h
3
u
xxxx
(
h
4 ) The fourth order derivatives are dissipative in nature.
26