### Page Content

## B04: Numerical simulation of pulsation impinging jet arrays II

**Principal investigators: **

Prof. J. Sesterhenn (mail [1])

**Scientist:** Gabriele Camerlengo (mail [2])**phone ** (030) 314
24665

### Summary

A pulsating impinging jet has substantially higher heat transfer than a non pulsating jet. The extra heat transfer comes from a vortex system impinging on the plate. This vortex system was shown to be related to the eigenmodes of the jet but it is unclear how it behaves under cross flow and curvature like in the realistic case in the turbine blade. On top of that it is not the linear eigenmode rather than a saturated limit cycle which is of interest. We investigate this limit cycle and its receptivity in order to optimize saturation. The interaction of several jets under cross flow and curvature is investigated numerically for a line array.

## 2nd Funding period 2016 - 2020

- Fig. 1: Turbulent structures highlighted through the contours of the second invariant of the instantaneous velocity gradient tensor Q on a plane passing through the jet axis.
- [3]
- © TUB

**Analysis of an impinging jet with fully
turbulent inlet conditions**

By looking at Fig. 1, which compares the instantaneous values of
the second invariant of the velocity gradient tensor of the laminar
and turbulent inlet simulations, we note that the coherent structures
generated by Kelvin-Helmholtz instabilities are only observable in
the free jet region of the laminar inlet case. On the other hand, the
flow in the free jet shear layer appears strongly more chaotic in the
turbulent inlet case. Although the wall jet region appears fully
turbulent in both cases, we observe a recurrent aggregations of
vortices flowing downstream (vortex rings) only in the laminar inlet
case. As shown in Fig. 2, this reflects in the mean Nusselt number
distribution at the wall, where the characteristic shoulder is not
observable in the turbulent inlet case. We found that in the latter
case, the Nusselt number is approximately 20% higher in the region
*r/D*<1.2, whereas the laminar inlet case features a 5%
larger heat flux in *r/D>*1.8. The surprising result is that
both jets provide about the same heat flow rate (in the region
*r/D*<4), because the region where the turbulent inlet jet
has a much better cooling rate (approx. 20% higher) is much smaller,
too (approx. 2.5 times smaller).

- Fig. 2: Average Nusselt number distribution at the wall as a function of the dimensionless distance from the jet axis r/D (turbulent inlet case).
- [4]
- © TUB

Root mean squares (RMS) of the velocity in the shear layer are in good agreement when experimental and numerical data are compared. In the laminar inlet case, as a result of the laminar shear layer profile, natural modes of the impinging jet can develop undisturbed. In particular, modal structures appear in form of axis-symmetric ring vortices, which are accompanied with strong RMS of axial and radial velocity fluctuations in the free and wall jet regions. In the experiments, however, several different modes are induced by the turbulent inlet. They prevent the free development of natural modes so that the ring vortices in the experiment occur only in weakened and irregular form. As a result, the corresponding RMS in the shear layer areas are smaller.

- Fig. 3: Adjoint pressure component (modulus) of the most unstable mode.
- [5]
- © TUB

**Stability analysis of the impinging
jet**

In the 1st funding period we carried out a dynamic mode decomposition (DMD) of our DNS flow field. With this method, a dominant mode which has same non-dimensional frequency as the most effective inlet pulsation in the experiment could be detected. It is known that the application of the DMD in a linear system gives the eigenfunctions of the system. The linear stability analysis of the impinging jet system showed that global modes exist, whose the frequencies of correspond to the DMD frequencies. It is noteworthy that the sensitivity to external forcing is attained in the proximity of the jet inlet (Fig. 3). This feature can be gainfully exploited for the purpose of devising cooling efficiency enhancement methods.

- Fig. 4: Average Nusselt number Nu on the target plate as a function of the non-dimensional distance from jet axis r/D (curved plate case).
- [6]
- © TUB

**Effects of plate curvature**

In order to resemble the internal configuration of a gas turbine blade, a DNS of a jet impinging on flat plate has been performed. The analysis focused on assessing how good is the common hypothesis by which curved plates can be approximated as flat when estimating the heat transfer efficiency through them. It has been found that integral heat flux (integrated on the plate up to a distance from the jet axis of 5 jet diameters) differs from the flat plate case (reference case) by just 0.03% (Fig. 4). On the other hand, peak frequencies of the instantaneous Nusselt number appear reduced by 40% and 50% when compared to the flat plate case. A dynamic mode decomposition (DMD) of the system showed, as observed in the reference case, that the modes oscillating at the aforementioned peak frequencies are those responsible for the characteristic heat flux profile at the wall. It can be consequently affirmed that if only averaged properties are of interest, it is possible to model curved impingement plates as flat surfaces. On the contrary, when the dynamic response of the system is addressed in order to implement, for instance, dynamic heat transfer enhancement techniques, the presence of the curved plate can not be disregarded.

### [Translate to English:] Analysis of a linear array of impinging jets in cross flow

- [7]
- © Copyright??

[Translate to English:] DNS of an infinite linear
array of impinging jets is carried out with the addition of cross
flow. This configuration mimics a line of jets in the spanwise
direction of the blade which follows the first line and is therefore
subject to a cross flow induced by the previous lines of jets.
Reynolds and Mach numbers of the simulation are respectively 8000 and
0.8, whereas the nozzle-to-plate distance is, as in the previous case,
*5D*. Similarly to the experiments carried out in
**B03**, the blowing ratio is set to 5. Also in this
case, the DMD of the flow indicates that two dominant modes exist. By
looking at Fig. 5 we observe that the first and second mode are not
perfectly toroidal because of the effect of the cross flow.
Nevertheless, the cross flow does not prevent the formation of
Kelvin-Helmholtz instabilities, which are particularly correlated with
the first dominant mode. As in the absence of cross flow, we expect
that pulsating the jet with a frequency close to that of the dominant
mode is optimal. Importantly, the latter reduces by 15% when compared
to the case without cross flow. The frequency associated to the second
mode is affected as well, reducing in the same proportion. Analogously
to the curved plate case, we conclude that the cross flow plays a role
that is not negligible for the correct design of dynamic internal
cooling techniques for gas turbine applications.

### Acknowledgement

Simulations were performed at the High Performance Computing Center Stuttgart (HLRS) under the grant number JetCool/44127.

### Publications

Camerlengo, G. und J. Sesterhenn: DNS study of the
turbulent inflow effects on the fluid dynamics and heat transfer of a
compressible impinging jet flow. In: Nagel, W. E., D. H. Kröner und
M. M. Resch (ed.): *High Performance Computing in Science and
Engineering '19*, accepted 2019. Springer International
Publishing, Cham.

Camerlengo, G., J. Sesterhenn, F. Giannetti, V. Citro und P.
Luchini: On the stability of subsonic impinging jets. In: Paolone, A.
et al. (ed.): *Proceedings of the XXIV Conference – The Italian
Association of Theoretical and Applied Mechanics (AIMETA 2019)*,
accepted 2019. Springer International Publishing, Cham.

Camerlengo, G., Borello, D., Salvagni, A., & Sesterhenn, J.
(2019). Effects of Wall Curvature on the Dynamics of an Impinging Jet
and Resulting Heat Transfer. In *Active Flow and Combustion Control
2018* (pp. 355-366). Springer International Publishing, Cham.

Wilke, R., & Sesterhenn, J. (2017). Statistics of fully
turbulent impinging jets. *Journal of Fluid Mechanics*,
*825*, 795-824.

### Summary

An effective cooling of turbine components subject to high thermal stresses is vital for the success of new engine and combustion concepts. Therefore efficient cooling mechanisms have to be developed and optimized. A promising approach is the use of pulsating impinging jets. The impinging jets generates vortices which encapsulate the hot fluid at the wall and effectively transport the pockets of hot fluid away from the wall. This project investigates the underlying mechanisms of pulsating impinging jets by direct numerical simulation and optimizes them with the aim of improved cooling efficiency.

- Fig. 1: Computational domain with iso-surfaces at Ma=0.2 coloured with pressure and at Q=10^5 m^2s^-4 coloured with radial velocity. Re=8000.
- [8]
- © TUB

In order to investigate the heat transfer of a confined round impinging jet (not pulsating) as a reference, two direct numerical simulations are performed at Reynolds numbers Re=3300 and Re=8000 using a grid of 512 × 512 × 512 respectively 1024 × 1024 × 1024 points.

Figure 1 shows the computational domain including a snapshot of the turbulent structures of the impinging jet at Reynolds number Re=8000. The two walls are isothermal. All other boundaries are non-reflecting. The computation is done on a grid with more than 1 billion points. Much more than a quarter of a million time steps are required to reach convergence in turbulent statistics.

- Fig. 2: Mean radial distributions of the local Nusselt number Nu (top) and turbulent heat flux (bottom). Re=3300.
- [9]
- © TUB

Each configuration features two annular regions with local maxima of heat transfer at the impinging plate (Figure 2). The periodical appearance and disap-pearance of the vortex pairs (Figure 3) lead to a high averaged turbulent heat flux exactly where the averaged Nusselt number reaches its local maxima.

- Fig. 3: Snapshots at Re=3300. Slice through the jet axis: temperature. Impingement plate: Nusselt number.
- [10]
- © TUB

Given these results, the aim for future work is to enhance the vortices by applying an pulsating inlet condition and thereby the wall normal turbulent heat flux concluding in a more efficient cooling of the impinging plate.

### Acknowledgement

The simulations were performed on the national supercomputer Cray XE6 at the High Performance Computing Centre Stuttgart (HLRS) under the grant number GCS-NOIJ/12993.

### Publications

Wilke, R. und J. L. Sesterhenn: Direct Numerical Simulation of Heat Transfer of a Round Subsonic Impinging Jet. In: Notes on Numerical Fluid Mechanics and Multidisciplinary Design, Bd. 127, S. 147–159. Springer, 2015.

Wilke, R. und J. L. Sesterhenn: Numerical Simulation of Impinging
Jets. In: High Performance Computing in Science and Engineering ‘14,
S. 275–287. Springer International Publishing, 2015.

Wilke, R.
und J. L. Sesterhenn: Numerical Simulation of Subsonic and Supersonic
Impinging Jets. Akzeptiert f ¨ ur High Performance Computing in
Science and Engineering ‘15), 2016.

Haucke, F., W.
Nitsche, R. Wilke und J. L. Sesterhenn: Experimental and Numerical
Investigations Regarding Pulsed Impingement Cooling. In: Deutscher
Luft- und Raumfahrt Kongress, Rostock, Germany, 2015.

Hossbach, S., R. Wilke und J. L. Sesterhenn: Identification of Material Vortices. In: Turbulence, Heat and Mass Transfer 8, isbn 978-1-56700-428-8, 2015.

## Speaker

**Prof. Dr.-Ing. Dieter Peitsch**

e-mail query [11]

## Office

Steffi Stehrsec. F 2

Room 107

Marchstr. 12

10623 Berlin

+49 (30) 314 22954

e-mail query [12]

parameter/en/id/122151/?no_cache=1&ask_mail=Xy5ImwA

AwtXWQNlwgnWfT4VGZSndTaTuDZhllMJToEwoowf06e8pwg%3D%3D&a

mp;ask_name=JOERN%20SESTERHENN

parameter/en/id/122151/?no_cache=1&ask_mail=Xy5ImwA

BBoR2b8LjTBMB0GZ6JpszmImOeQse5z0ME%2FuYAWEVWLsICw%3D%3D

&ask_name=ROBERT%20WILKE

rojekte/B04_Website_Fig1.png

rojekte/B04_Website_Fig2.png

rojekte/B04_Website_Fig3.png

rojekte/B04_Website_Fig4.png

rojekte/B04/B04_Fig5.png

dungen_allg/b04_1_Domain.png

dungen_allg/b04_2_Nu.png

ldungen_allg/b04_3_Snapshots.png

/parameter/en/id/122151/?no_cache=1&ask_mail=Xy5Imw

AC07LWigCUliHSMRQrfesink26cx9MJJPapGJj4%2B%2FsCGO6sA%3D

%3D&ask_name=Prof.%20Dr.-Ing.%20Dieter%20Peitsch

/parameter/en/id/122151/?no_cache=1&ask_mail=Xy5Imw

AC6fidtLk%2BDzn5QkbA6pk%2BA8botpQepEWOrdio9D%2BI8jR%2BT

w%3D%3D&ask_name=Ansprechpartner