Next Article in Journal
The Effect of Climate on the Solar Radiation Components on Building Skins and Building Integrated Photovoltaics (BIPV) Materials
Next Article in Special Issue
Experimental Study on Thermal Performance of a Loop Heat Pipe with Different Working Wick Materials
Previous Article in Journal
A Data-Mining Approach for Wind Turbine Fault Detection Based on SCADA Data Analysis Using Artificial Neural Networks
Previous Article in Special Issue
Experimentally Verified Flow Distribution Model for a Composite Modelling System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Curved Surface Minijet Impingement Phenomena Analysed with ζ-f Turbulence Model

1
Department of Fundamental Research in Energy Engineering, AGH University of Science and Technology, al. Mickiewicza 30, 30-059 Krakow, Poland
2
Institute of Energy, Faculty of Mechanical Engineering and Ship Technology, Gdansk University of Technology, ul. Narutowicza 11/12, 80-233 Gdansk, Poland
3
Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, van der Maasweg 9, 2629 HZ Delft, The Netherlands
*
Author to whom correspondence should be addressed.
Energies 2021, 14(7), 1846; https://doi.org/10.3390/en14071846
Submission received: 21 February 2021 / Revised: 19 March 2021 / Accepted: 20 March 2021 / Published: 26 March 2021
(This article belongs to the Special Issue Heat and Mass Transfer Issues in Mini Gaps)

Abstract

:
The jet impingement phenomenon plays an important role among the heat transfer intensification methods. Very often its application and analyses refer to simple flat surfaces, while there is a lack of information in the literature for cases addressing curved surfaces. In the present work, the single jet impingement on the non-flat (concave and convex) surface is studied for a wide range of geometries, which originate from the mini-jet heat-exchanger design. The numerical simulations were performed by an advanced ζ-f turbulence model implemented in the open-source OpenFOAM (ESI-OpenCFD Ltd, Bracknell, United Kingdom) code. Noticeable differences in the phenomena occurring on the convex and concave surfaces were identified in the stagnation zone. Besides, the existence and location of the secondary peak in the Nusselt number distribution differed between the cases. These distributions were influenced by the shape of geometry, which determined flow characteristics and resulting heat transfer performance. The origins of these differences were looked at in the turbulence characteristics close to the impinged surface of the stagnations zone and its vicinity, where turbulence kinetic energy and enstrophy were analysed. It was stated that the differences are already noticeable for the single jet impingement case, but they might sum up when multiple jets are considered. Therefore, here presented results would be important for the analysis of the overall unit of mentioned mini-jets heat-exchanger.

1. Introduction

1.1. Flat Surface Jet Impingement

Between numerous turbulence-related cases that are considered numerically, the jet impingement is an example of the phenomenon that despite the development of various computational methods and many scientists involved, still causes problems, even though its analyses started as early as just after the Second World War [1]. Uncertain accuracy of available computational methods, especially related to the thermal results, was commented by Zuckerman and Lior [2] or Dewan et al. [3].
Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS)-based methods, previously applied for analyses of particular, simple cases of low Reynolds number, have already become available also for complex flows, including for example additional heat transfer or crossflow phenomena. However, it does not mean that the Reynolds Averaged Navier-Stokes (RANS) methods should be forgotten. They still are important due to significantly lower computational time, sufficient accuracy and common availability [2,3]. Moreover their advanced examples are step by step introduced with an ultimate goal of reaching the highest accuracy and possible universality.
An overall summary of the DNS-based research of jet impingement was reported by Ries et al. [4]. Other DNS research examples were reported by Satake and Kunugi [5]: mechanism of eddies transport in the impingement zone, Reynolds number Re = 5300; by Hattori and Nagano [6]: turbulence statistics of velocity and heat fluxes, Re = 9120; by Dairay et al. [7]: unsteady processes impacting the spatial distribution of the heat transfer coefficient at the impinged wall, Re = 10,000; by Wilke and Sesterhenn [8]: properties of the flow at the inlet on the properties in the impinging area, Re = 8000; by Koido et al. [9]: 4 jets array impingement, Re = 1500.
While DNS cases of multiple jet impingement need to be limited by low Reynolds number or the number of jets, LES methods allow to overcome some of those limitations. The widely cited experimental work of Geers et al., [10,11,12,13] who investigated the array of 13 turbulent jets impinging the flat surface, was repeated numerically and described in the article of Draksler et al. [14], who performed the LES calculations of Reynolds number equal to 20,000. However, the utilised mesh consisted of more than 500 million elements. Nevertheless problems with the accuracy occur also when LES methods are used to analyse jet impingement—just to mention some findings from articles of Hadžiabdić and Hanjalić [15] and Uddin et al. [16].
Therefore, the more complex analysed problem is, the more useful are RANS based methods to examine it. Unfortunately, according to [2] only a few of them are capable to correctly represent the jet impingement, when it comes to thermal effects. The most recommended was the v 2 ¯ f turbulence model of Durbin [17]. Due to the mentioned issues with the thermal outcome of RANS based simulations of jet impingement, many researchers tried to propose new models, or improve the existing ones. General approach is to introduce additional equations to the system. One of the most recent works was published by Billard and Lawrence [18], who decided to modify the four equation v 2 ¯ f model of Durbin [17] and showed good accuracy of the jet impingement results with the reference data. Another example of the v 2 ¯ f evolution is the model called ζ-f proposed by Hanjalić et al. [19]. This model was designed to overcome the main drawback of the original Durbin’s model [17], the unstable behaviour due to the formulation of boundary condition for f variable, which determines the near-wall relaxation function. The same idea was behind the research of Lien and Kalitzin [20], who also tried to make the v 2 ¯ f model more robust. Finally Kenjeres et al. [21] proposed quasi-linear algebraic heat flux and buoyancy extended stress model, based on the [17]. The valuable benchmark for the single jet impingement analysis of the RANS method based calculations are analyses of Behnia et al. [22,23], based on Durbin’s model [17]. They tried to reproduce numerically the experiments by Cooper et al. [24] and Baughn and Shimizu [25]. Their results were also later verified experimentally by Yan [26]. All those publications were then included in the ERCOFTAC association database, as the reference data for scientists, used for example for validation purposes of the results from [14,18,19,20,21].

1.2. Curved Surface Jet Impingement

All aforementioned literature sources concerned flat surface impingement. However, situation changes when the non-flat impingement surfaces are considered. In their cases the available research data is rather limited. Recently the DNS analysis of concave surface jet impingement was published by Camerlengo et al. [27]. They conducted analysis of Reynolds number equal to Re = 3300 and the results exhibited an instantaneous peak values of the heat transfer approximately 30% lower than those found in the reference study of flat surface jet impingement. On the other hand negligible differences were found in the time-averaged integral heat transfer rate. More importantly, they underlined the general lack of the numerical simulations-based publications concerning the non-flat jet impingement. Some experimental works were published in [28,29]. One of the few articles trying to deliver universal correlation for the jet impingement on the concave surfaces was published by Ying et al. [30]. For that case, however, only single jet was taken into account. Few examples of jet arrays impinging on the non-flat surfaces were presented in the literature, however they were mostly related to the linear array of jets [31,32,33]. The comparison of the flat and non-flat jet impingement of the jets array was published by Kura et al. [34], with utilisation of RANS approach.

1.3. Motivation

The idea behind the following article comes from the publications [35,36,37]. They presented the innovative, cylindrical in shape, heat exchanger and its experimental analysis. Novelty came from the fact, that both sides of heat transfer wall (cylindrical) were surrounded by perforated pipes. The drilled orifices (of diameter equal to 0.001 m) generated jets, which hit the heat transfer surfaces, intensifying the process. Schematic view of heat exchanger is shown in Figure 1. The resulted macroscale heat transfer coefficients turned out to be significantly higher than in the case without jets, and, in addition, higher than in comparable heat exchangers of competing design. On the other hand, device-scale pressure losses, which were caused by those additional parts, increased just slightly.
The presented results [35,36,37], significant from the engineering point of view, refer to the overall output. However, the description of the small-scale turbulent phenomena occurring near the heat transfer curved wall due to the jets’ impingement, being the source of the heat transfer intensification, are not yet available. Such information will be of great importance, for the purposes of further optimisation of the similar devices. It requires distinction of numerical model, which will provide the best possible accuracy of the results. Previous studies, which have been already published [34], were based on the Intermittency Transition Model described by Menter et al. [38].
The results exhibited some issues referring to an accuracy and further research was necessary to make a proper choice. It was already performed and published in [39,40]. In both articles, which concerned single jet impingement on both flat and non-flat surfaces, the superior performance of the Hanjalic’s ζ-f model over others investigated was proven. Additionally this model allows one to predict the turbulent-laminar transition effects that might occur when the jets array is considered [41]. The main conclusion coming from [39,40] was, that the hydrodynamic phenomena occurring in both cases (flat and non-flat geometries impingement) are very similar, so the turbulence model predicting the heat and mass transfer well for the flat case would also be suitable for the non-flat ones.
Presented paper concentrates on analysis of non-flat, namely concave and convex, impinging jet, from the hydrodynamic and heat transfer point of view. Special attention is devoted towards their relation in the impinged surface proximity. Following scheme, presented in Figure 2, simply defines the problem, which authors tried to solve and the main outcomes of this paper.

2. Geometrical Model and Boundary Conditions

Two non-flat surface types were considered, concave and convex. Both of particular curvature radius R. The schematic view of important geometrical parameters is presented in Figure 3. Other necessary variables are as follows: the diameter D of the jet generating nozzle and distance H between the nozzle exit and stagnation point (the central point of the impinged surface). The analyses regarded the axisymmetric cases, so geometry is shown in such manner.
Applied boundary conditions were as follows:
  • at the inlet: the constant bulk velocity, determined from the Reynolds number values (namely Re = 5000; Re = 10,000; Re = 23,000) and the constant bulk temperature equal to 293 K,
  • at impinged wall: the constant heat flux equal to 1000 W/m2,
  • at the upper wall: the constant heat flux equal to 0 W/m2,
  • at the outlet: the gauge pressure equal to 0 Pa.
In Table 1, the summary of all investigated cases is presented. Variables are shown in a non-dimensional way, an approach typical in the literature regarding the jet-impingement. The inlet Reynolds number based on inlet bulk velocity ub:
Re = u b D ν ,
where D is the jet issuing part diameter, m and ν is the kinematic viscosity, m2/s.
As the inlet diameter size might strongly impact the results of the heat transfer at the impinged surface, which was proven in [42,43,44], it was kept constant throughout the research and equal to 0.01 m. This information is very often missed in published sources, which can rise misunderstandings. Moreover without it, particular results should not be used as the data validation.
Combining all the above cases, the 54 simulations were performed. The Reynolds number equal to 23,000 was selected as a reference case, widely used as the jet impingement benchmark [14,18,19,20,21]. Other values of the Reynolds number were adapted in the basis of [35,36,37] and conditions occurring in the minijet heat exchanger. The values of H/D and R/D were chosen arbitrarily. Nevertheless they reflect the possible geometries. For example curvatures of R/D ≤ 4 were not possible to be applied for the jet impingement cases, as the geometry would be too narrow. Also the H/D ratios higher than 4 would significantly increase the overall size, which is unwelcomed in the case of the heat exchangers.
In Figure 3, the placement and definition of all boundary conditions are also shown. At the inlet, the fully developed profile was applied, using the mapping method [45], which required prescription of the bulk velocity ub, and the bulk temperature Tb. Turbulent inlet conditions were then concluded with every consecutive iteration by mapping the downstream flow back to the inlet until the flow was fully developed. The final turbulence intensity was afterwards equal to 5.8%, while the ratio between the centreline velocity and bulk velocity uc/ub was equal to 1.2. It is important to note, that such intensity was relatively high, and strongly impacted the results in the stagnation zone. As shown in [22,23,44], assumed inlet conditions (both regarding velocity profile and turbulence parameters) affect significantly the heat transfer rates at the impinged surface. Such information very often is not presented in the literature, which precludes exact representation of boundary conditions, because behind the statement fully-developed profile many possible configurations can be hidden. For the purposes of the following article, aforementioned parameters were therefore kept unchanged. In the outlet, the zero gauge pressure boundary condition was applied in all cases. At the surface, the constant heat flux q was applied, equal to 1000 W/m2, being more physical boundary condition than fixed temperature [23].
The working fluid was the air, in an incompressible form. All its properties were kept constant, with values listed in Table 2. The reason was mostly due to the chosen heat flux applied at the impinged surface, equal to 1000 W/m2. The resulted temperature gradient was relatively small, and temperature dependent parameters of the working fluid would not have impacted the accuracy of the results.

3. Numerical Model and Solver Configuration

OpenFOAM was chosen to perform the numerical analyses, mainly due to its open-source character and possibility to extend its functionalities [45]. To investigate the jet impingement phenomenon, it was necessary to include the conservation laws of mass, momentum and energy in the model:
u j ¯ x j = 0 ,
ρ x j u i ¯ u j ¯ = p ¯ x j + x j 2 μ S i j ρ u i u j ¯ ,
x j u j ¯ Θ = x j a Θ x j u j θ ¯ .
where ui,j are the components of velocity vector, m/s; ρ is the fluid density, kg/m3; p is the pressure, Pa; μ is the fluid dynamic viscosity, Pa·s; Sij is the tensor of strain rate, 1/s; u i u j ¯ is the term for the Reynolds stress tensor, m2/s2; Θ is the mean fluid temperature, K; θ is the fluctuation of fluid mean temperature, K; a is the thermal diffusivity, m2/s and u j θ ¯ is the turbulent heat flux vector, K∙m/s. Turbulent component of the heat flux were obtained with the use of eddy diffusivity concept [45]. As the RANS based turbulence modelling approach was applied, averaged quantities were employed as the variables.
The single jet impingement was analysed, and axisymmetric cases were considered, with symmetry axis presented in Figure 3. Single jet impingement can be thought as the phenomenon that with time-constant boundary conditions reaches steady state. Therefore, the steady-state SIMPLE algorithm was applied. Second order numerical schemes were used [45].

Selection of Turbulence Model

Based on the conclusions from [39,40,44] the ζ-f model was chosen to be used in the following paper, after being successfully implemented in the OpenFOAM software. The governing equation for the model including the turbulence kinetic energy is presented in Equation (6) and the governing equation for its dissipation is shown in Equation (7):
ζ = v 2 ¯ k ,
x j u j ¯ k = G ε + x j ν + ν t σ k k x j ,
x j u j ¯ ε = C ε 1 G C ε 2 ε τ + x j ν + ν t σ ε ε x j ,
x j u j ¯ ζ = f ζ k G + x j ν + ν t σ ζ ζ x j ,
L 2 2 f f = 1 τ C 1 + C 2 G ε ζ 2 3 ,
where G is the turbulence kinetic energy production term, m2/s3, Cε1, Cε2, C1, C’2, σk, σε, σζ are the variables specific for this particular model, τ is the model time scale and L is its length scale, which are used to limit the model and enhance the quality of the results obtained with it.
Being a low-Reynolds type model, it avoids any wall-functions that might perturb the performance. Its main feature, in comparison to the precursor, the v 2 ¯ f model [17], has better robustness, achieved due to reformulation of the elliptic relaxation function f, defined by Equation (9), especially by modifying its wall-boundary value and assuring the better reproduction of the turbulence kinetic energy production in the near wall region. This region is especially important when investigating the jet impingement, because there the heat transfer strongly depends on the near-wall turbulence. Finally, instead of using the v 2 ¯ as the variable, which is the wall-normal velocity scale, ζ is introduced, defined by Equations (5) and (8), and obtained by normalization of the velocity scale with turbulence kinetic energy k.
As it was mentioned in the Motivation section, the numerical model was validated [39,40] and exemplary results are shown in Figure 4, presenting the comparison of the local heat transfer results for various RANS turbulence models and experimental data obtained for the flat surface impinging jet. The green lines in Figure 4 indicate the results obtained with ζ-f model [19], their difference lies in the formulation of the turbulence intensity at the inlet [44] and indicates the borderlines of the Nusselt number results that can be obtained with this model. The high accuracy of the ζ-f model, in comparison with the reference data, can be seen.
In Figure 5 the exemplary results of the mesh independence tests are presented, for two of the cases mentioned in Table 1, represented by Nusselt number values. The mesh was generated in the same manner as in [39,40,44]. The most important was to maintain dense enough mesh near the heat transfer wall, as well as near the nozzle wall. While typically it is recommended for the low-Reynolds RANS turbulence models to keep the non-dimensional wall normal distance y+ at a value of approximately 1, for the purposes of this research that value was below 0.5, with decent number of mesh elements in the viscous sublayer. As can be seen, despite the mesh size, obtained results remained almost the same, as long as the indicated requirements were maintained. It is worth mentioning, that the designed mesh which fulfilled the criteria for the highest Reynolds number cases was also used in less demanding cases, where Reynolds number had lower values.
For all analysed cases, the mesh consisted only of hexahedral elements to minimise the impact of non-orthogonality and skewness on the results. Also the boundary layer mesh was optimised to limit the aspect ratio of its elements to avoid the results distortion.

4. Results

Prior to the results description, an information regarding the way of their presentation is necessary to be provided. The location of some parameters will be shown in a non-dimensional way, by using the distance x/D, as indicated in Figure 6, where x is the length of the arc between the stagnation point and any selected point on the surface and D is the nozzle diameter. It may be argued, that instead of the arc length, the relative length of the straight segment connecting the points, also shown in Figure 6, should be used. However, for the range of the distances analysed in the following article, the difference between those two, was negligible.
Due to the shape of the surfaces, difference in the distance between them and the nozzle exit occurs (particular H/D value is maintained exactly constant only at the symmetry axis). However, the discrepancies between these distances in the whole stagnation zone are very small. In addition, it makes the comparison of results between flat and non-flat surface jet impingement possible.
Lastly, the h variable is also shown in Figure 6, which would be used in Section 4.3 to provide results in particular, various non-dimensioned distance h/H from the impinged surface.
The Results section is divided into three main parts:
  • overall thermal results of the jet impingement for all 54 cases,
  • influence of the impinged surface type on the flow behaviour near it,
  • impact of the impinged surface type on the impinging jet.

4.1. Thermal Performance of Convex and Concave Cases

In Figure 7, Figure 8 and Figure 9 the local values of the Nusselt number at the impinged surface for all 54 cases are plotted. Results related with the convex geometry are indicated by the dashed line and these related to the concave with the solid line. First conclusion that can be stated after analysis of the charts from Figure 7, Figure 8 and Figure 9 is related to the difference in the values of Nusselt number in the stagnation point and close to it. For all convex cases, those values were noticeably higher than for concave ones. As seen in Figure 7, Figure 8 and Figure 9, the reason lies in the stagnation zone, where the visible differences of Nusselt number values occur. Therefore, it is the area that has to be analysed in detail to explain the nature of such phenomenon.
With increasing H/D ratio, the characteristic feature of the local Nusselt number curve for some jet impingement cases is vanishing—the secondary peak. Its origin was described in detail in [46]. They related the skewness of the pressure temporal distribution with the rebound of primary vortices observed for a single vortex ring operating near the location of the secondary peak. They have stated, that the rebound is the event causing the mean heat transfer enhancement to stop. While the RANS based turbulence models do not allow the extraction of so detailed data, still their ability to predict the secondary peak on the Nusselt number curve indicates their suitability for trustworthy jet impingement analyses—which was another reason to choose the ζ-f model for the presented research [23]. This ability in reference to particular RANS models was also reviewed in [2].
In Figure 10, selected Nusselt number curves are plotted also in two-dimensional plot. Chosen geometry configurations reflects two outermost cases. As indicated in Figure 7, Figure 8 and Figure 9, differences between Nusselt number occurs due to the both hydrodynamic and geometric variations.
In general, the shape of all the curves, with mentioned secondary peak less visible with increasing H/D, is similar to findings from experimental research of Lee and Lee [47], and O’Donovan and Murray [48]. While they focused on the flat surface jet impingement, similarities between their results and the ones presented here can be observed. When it comes to shape of particular Nusselt number curves in relation to specific H/D ratios and chosen Reynolds number values—it is another suggestion, that the phenomena behind the jet impingement on flat and curved surfaces share the same origin and the turbulence model adequate for flat cases would also perform well for non-flat types analyses. An influence of the inlet Reynolds number on the Nusselt number is also visible in Figure 7, Figure 8 and Figure 9. For lower value of H/D, the Reynolds number required for the secondary peak to occur was also lower. This statement is true for both, convex and concave geometries. Nevertheless, the secondary peak, if occurs, is more visible for convex geometries. This secondary maximum of Nusselt number tends to move away from the stagnation point with increasing Reynolds number.
In Figure 11 the stagnation point Nusselt number value are presented, while in Figure 12 the Nusselt number values averaged over the surface, up to x/D = 2, are presented. This distance was chosen, as further away from the stagnation zone the stream start to change the behaviour to the typical wall jet, independent from the jet impingement phenomenon [18,19,22,25,26]. Generally, for the convex geometry with the increased curvature ratio R/D the stagnation point Nusselt number decreases, while for the concave geometry the behaviour is opposite. Such results are in agreement with the preliminary research presented in [41]. The source of such behaviour is related to the size and characteristics of the stagnation zone, varying between convex and concave impinged geometries. This is more visible for local values from Figure 11. Averaged values, shown in Figure 12, tend to differ less. It can be stated, that the averaged values remain more or less similar. Such tendency, however, occurs not for all cases and less clearly—the reason is due to averaging, which might diminish the differences. Similar situation happened, when analysing the area averaged Nusselt number values of jets array impingement in [41]. Despite the impinged surface type, discrepancies were small. Nevertheless, from Figure 11 and Figure 12 it can be concluded, that jet impingement in the convex geometry environment assures slightly higher heat transfer rates.
In Table 3, the location x/D of the secondary Nusselt number maximum is presented, for the cases in which its occurrence was the most distinguishable. The differences of about 5% occurred between convex and concave geometries. For other cases the maxima were less visible therefore they were omitted in Table 3.

4.2. Influence of Geometry Type on Stagnation Zone Flow Behaviour

In Figure 13, the illustrative close-up of the stagnation zones for selected geometry case is presented, for both convex and concave surfaces, for H/D = 2 and R/D = 4. Reynolds number was equal to 23,000. The zoom includes also small part of the nozzle. Analysis of this area, where relatively low velocity occurred, confirms already stated conclusions, that the stagnation zone is slightly larger in concave geometry. The reason is related strictly to the geometry and the flow distribution. In convex geometry it is easier for the stream to leave the impinging area and move further. In concave geometry, due to its shape, the outflow is being slightly limited by the shape of curvature. Stagnation zone is also characterised by the high-pressure gradient, a result of the velocity decrease.
In Figure 14, local values of turbulence kinetic energy, k, are shown, for both concave and convex geometries, with an example case, where H/D = 2 and R/D = 4. The values were obtained at the first mesh nodes near the heat transferring wall. The reason for choosing these particular values was shown in the authors’ publication [49]. As Nusselt number distribution is influenced mostly by the phenomena in the near-wall region, it should be referred to corresponding values of other variables from this zone, despite the fact, that for example global maxima of turbulence kinetic energy appear significantly above that region.
As can be seen in Figure 14, with increasing Reynolds number, the second peak in the turbulence kinetic energy distribution becomes more evident, and eventually higher, than the first one. Values of turbulence kinetic energy for convex cases are higher until reaching the first peak, afterwards they get lower than for concave cases. After reaching the first minima, they again exceed the values of concave cases. Eventually, after the second peak, they reach lower values. The distributions for convex and concave geometries intersperse.
The location of both inflection points is strongly related to the location of peaks in the Nusselt number values of the same cases, as indicated in Figure 14c by the green lines. For all concave cases, those peak points are moved slightly away from the stagnation points. Away from the stagnation zone, since the distance of about x/D ≥ 2.5, the wall jet creation can be noticed from the plots of Figure 14. The relaminarization processes can be assumed there, as the turbulence kinetic energy is decreasing and reaching relatively small values.
In Figure 15, the distribution of turbulence kinetic energy dissipation, ε, is compared with the distribution of local enstrophy e. It is described by Equation (10), on the basis of vorticity (velocity curl), with its unit being 1/s2 [45]:
e = 0.5 × u 2 ,
Both variables, e and ε are shown for both concave and convex geometries, using an example case, where H/D = 2, R/D = 4 and Reynolds number was equal to 23,000. One noticeable difference is related with the values in the stagnation zone, where the turbulence dissipation ε reached relatively high values, which is related to the rapid flow suppression. Local enstrophy, on the other hand, related to the change in velocity field, according to its definition from Equation (10), does not reach high values in the stagnation region.
The visible connection between those two variables can be concluded from the Figure 15. It confirms the findings of Zhu and Antonia [50]. Away from the stagnation point regions of high local enstrophy are associated with those of high turbulence kinetic energy dissipation. It is a valuable information, especially, if the analyses of the jets array impingement would be performed, since enstrophy can possibly be used for the optimisation process. Moreover, both variables reach high values in near-wall region, indicating the importance of analyses of that region and its impact on thermal performance of various jet impingement cases.
In Figure 16, the plots showing distribution of the turbulence kinetic energy dissipation, ε, and distribution of the local enstrophy, e, are presented. Both values are again shown for concave as well as convex geometries, with an example cases, where H/D = 2, R/D = 4 and Reynolds number was equal to 23,000. The values were obtained at the first mesh nodes near the heat transfer wall. The curves related to the dissipation, ε, share the same tendency of maxima location and shape, as the results of turbulence kinetic energy, shown in Figure 14c. However, analysis of the enstrophy curve leads to the conclusion, that slight difference between it and the dissipation can be seen. The first peaks of both concave and convex cases, noticeably higher, than the second, can be related to the end of the stagnation zone and start of the wall jet formulation. Even more distinguished change of location occurs for the second maxima. This time, however, they occur before the peak of turbulence kinetic energy, and appear close to the location of the secondary peak of Nusselt number. Therefore, while qualitative comparison of images from Figure 15 suggest close relation between the distribution of ε and e, plots shown in Figure 16 reveal more detailed comparison and some discrepancies. Finally, in Figure 17 the values of the v 2 ¯ , calculated with Equation (5), are shown. They were obtained along the centreline of the impinging jet and h/D = 0 corresponds then with the stagnation point. The “wall-normal” velocity scale v 2 ¯ is a good indicator of the turbulence transported with the impinging jet. With higher values of the Reynolds number, the jet core is preserved longer, which is indicated by more evident curvatures close to the wall. The location of the v 2 ¯ maxima near the impinged surface is different with different surface types—for convex geometry they are closer to the wall than for concave one, which is in agreement with the conclusions coming from Figure 13.

4.3. Impact of Surface Shape on Impinging Jet Hydrodynamics

The source of the heat transfer intensification away from the stagnation point is attributed to the generation of the boundary layer due to the wall jet formulation and its interference with the big eddies, which are generated at a distance x/D corresponding with the secondary Nusselt number maxima. It is related with the inlet conditions, as shown in Figure 7, Figure 8 and Figure 9. But even stronger influence of those occurs with the heat transfer intensification at the stagnation point and near it. As shown previously, the highest values of Nusselt number are obtained exactly there. However, as indicated in Figure 13, the stagnation zone is an area with negligible flow and high-pressure gradients. Thus the energy transfer occurs there, caused directly by the dissipation of eddies, which are transported with the inflow. It is assumed [5,50], that they are the main reason of such high heat transfer. In Figure 18, the velocity profiles in particular distance values from the heated and impinged surface, are shown, for the case, in which Re = 5000, H/D = 1 and R/D = 4. As it can be seen, while close to the nozzle exit conditions remain the same, the jet core (in which the velocity is stable along the jet normal direction) downwards starts to differ between the cases. The black arrow placed near the plots (Figure 18, Figure 19, Figure 20, Figure 21, Figure 22 and Figure 23) shows the flow direction, as the plots are oriented according to it. The ratio h/H indicates the particular distance above the impinged surface.
In Figure 19, the values of turbulence kinetic energy, corresponding with the locations and cases from Figure 18, are shown. The maxima, visible in the first 3 plots, indicate the regions, in which the shear occurs due to the jet hitting initially still fluid. It can be noticed, that the values of turbulence kinetic energy at the stagnation zone and close to it increase on consecutive figures—which proves the statement about the energy transfer due to the dissipation of eddies, which are transported with the inflow. In Figure 20, similarly to the situation from Figure 19, the values of enstrophy flux are presented. As for the case of turbulence kinetic energy, the shear regions are evident. However, when the jet approach the heated surface, vorticity values are low and increase while moving away from the stagnation point. It can be attributed with the already mentioned fact that almost no direct flow occurs in the stagnation zone.
In Figure 21, Figure 22 and Figure 23 the mentioned parameters are shown, for the case of H/D = 2 and R/D = 8. Reynolds number was equal to Re = 23,000.
As can be seen, the differences between the cases are less visible now. The reason is due to relatively high R/D ratio. As proven by authors in [34] such convex and concave geometries tend to differ less from the flat case impingement, and from themselves. Analysis of the velocity profiles indicates, that the jet core is preserved for a longer period of the jet—it leads then to the heat transfer intensification in the stagnation zone, compared to cases where Reynolds number was equal to 5000. The black arrows also show the flow direction.

5. Conclusions

In the article, the thermal and hydrodynamic results of non-flat surface single air jet impingement were presented. They were obtained with the use of tRANS-based, turbulence model of Hanjalic, ζ-f, [23], which in the authors’ opinion is capable of predicting properly the phenomena occurring in the analysed cases. The ζ-f turbulence model performs well in the stagnation region, where many RANS-based turbulence models tend to overpredict the generation of the turbulence kinetic energy [2]. According to the author’s knowledge, such cases were analysed for the first time with that model. It was implemented and validated in the open-source software—OpenFOAM. The model’s main features, which are relative robustness and convergence and more importantly the ability to indicate the secondary peak in the local Nusselt number distribution at the heated and impinged surface made it possible to calculate 54 various cases, with variable non-dimensional geometry constraints: R/D and H/D, as well as the inlet Reynolds number.
Noticeable differences between convex and concave surfaces impingement were identified, regarding not only stagnation point or, more generally, stagnation zone thermal performance, but also the existence of the secondary peak in relation to the mentioned variables. The shape of particular geometry impacts the size of stagnation zone, which results in the flow characteristics when the wall jet starts to formulate after impingement.
Distribution of local Nusselt number was also related to the turbulence kinetic energy, k, values near the heated wall, where differences, especially connected with the Reynolds number at the inlet, were shown.
In addition, the trial to interconnect the rate of turbulence dissipation, ε, and enstrophy, e, in the analysed models were done and presented. Enstrophy, related to the vorticity and velocity distribution can be taken as a measure of flow energy dissipation. While the big-scale comparison of ε and e showed similarities, the near-wall comparison indicated differences, related mostly with the impact of wall jet formulation on the enstrophy field.
It is worth to mention, that while enstrophy can provide a meaningful added value, it is not a common practice to use it when simulating jet impingement with the use of RANS models. Authors did not find any, related to the results above, examples of such approach. It is assumed, that this variable might be even more useful, when applied to the analyses of jets array impingement and used to identify the zones of particular thermal behaviour.
To sum up, only selected results were presented. The number of investigated cases made it impossible to include all possible variations. Nevertheless, the presented cases indicate the most important outcome from the research—altogether illustrating the most important differences between the convex and concave surfaces impingement.
The presented results will be of great importance after the next step of the research project—investigation of the jets arrays impingement—which is already started, to be able to optimise the occurring flow patterns.

Author Contributions

Conceptualization, E.F.-W., T.K. and J.W.; methodology, T.K., E.F.-W., S.K.; software, T.K.; validation, T.K., J.W., E.F.-W., S.K.; formal analysis, T.K.; investigation, T.K.; resources, T.K., J.W., E.F.-W.; data curation, T.K.; writing—original draft preparation, E.F.-W., T.K., J.W. and S.K.; writing—review and editing, T.K., E.F.-W., J.W., S.K.; visualization, T.K., E.F.-W.; supervision, E.F.-W. and J.W.; project administration, T.K., E.F.-W., J.W.; funding acquisition, T.K., J.W., E.F.-W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Higher Education (for T. Kura Grant AGH, while for J. Wajs Grant PG) and Project AGH Initiative for Excellence—Research University for E. Fornalik-Wajs.

Acknowledgments

This work was supported by the PLGrid Infrastructure.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Latin symbols
athermal diffusivity, m2/s
C1ζ-f model constant, –
C’2ζ-f model constant, –
Cε1ζ-f model constant, –
Cε2ζ-f model constant, –
cpspecific heat, J/(kg·K)
Ddiameter, m
eenstrophy, 1/s2
felliptic relaxation of the ζ-f model, 1/s
Gturbulence production rate, m2/s3
Hdistance between the pipe exit and impinged surface, m
hdistance between the stagnation point and the potential jet core termination, m
iiteration
kturbulence kinetic energy, m2/s2
Lturbulent length scale, 1/m
ppressure, Pa
pgaugegauge outlet pressure, Pa
Sij strain rate tensor, 1/s
Ttemperature, K
Tbbulk temperature, K
uvelocity, m/s
ubbulk velocity, m/s
uvelocity fluctuation, m/s
v 2 ¯ root mean square of the fluctuation of the velocity normal to streamlines, m2/s2
xradial distance from the stagnation point, m
ypoint to wall distance, m
y+dimensionless wall distance, –
Greek symbols
αef effective heat transfer coefficient, W/(m2·K)
εdissipation of the turbulence kinetic energy rate, m2/s3
ζ v 2 ¯ normalised by k, –
Θmean temperature, K
θtemperature fluctuation, K
λthermal conductivity of the fluid, W/(m·K)
μdynamic viscosity, Pa·s
νkinematic viscosity, m2/s
νteddy viscosity, m2/s
ρdensity, kg/m3
τturbulent time scale, 1/s
σε ζ-f model constant, –
σkζ-f model constant, –
σζζ-f model constant, –

References

  1. Freidman, S.J.; Mueller, A.C. Heat Transfer to Flat Surfaces. General Discussion on Heat Transfer. Proc. Inst. Mech. Eng. 1951, c164, 138–142. [Google Scholar]
  2. Zuckerman, N.; Lior, N. Jet Impingement Heat Transfer: Physics, Correlations, and Numerical Modelling. Adv. Heat Transf. 2006, 39, 565–631. [Google Scholar]
  3. Dewan, A.; Dutta, R.; Srinivasan, B. Recent Trends in Computation of Turbulent Jet Impingement Heat Transfer. Heat Transf. Eng. 2012, 33, 447–460. [Google Scholar] [CrossRef]
  4. Ries, F.; Li, Y.; Rißmann, M.; Klingenberg, D.; Nishad, K.; Böhm, B.; Dreizler, A.; Janicka, J.; Sadiki, A. Database of Near-Wall Turbulent Flow Properties of a Jet Impinging on a Solid Surface under Different Inclination Angles. Fluids 2018, 3, 5. [Google Scholar] [CrossRef] [Green Version]
  5. Satake, S.; Kunugi, T. Direct numerical simulation of an impinging jet into parallel disks. Int. J. Numer. Methods Heat Fluid Flow 1998, 8, 768–780. [Google Scholar] [CrossRef]
  6. Hattori, H.; Nagano, Y. Direct numerical simulation of turbulent heat transfer in plane impinging jet. Int. J. Heat Fluid Flow 2004, 25, 749–758. [Google Scholar] [CrossRef]
  7. Dairay, T.; Fortuné, V.; Lamballais, E.; Brizzi, L.-E. Direct numerical simulation of a turbulent jet impinging on a heated wall. J. Fluid Mech. 2015, 764, 362–394. [Google Scholar] [CrossRef]
  8. Wilke, R.; Sesterhenn, J. Statistics of fully turbulent impinging jets. J. Fluid Mech. 2017, 825, 795–824. [Google Scholar] [CrossRef] [Green Version]
  9. Koide, T.; Tsujimoto, K.; Shakouchi, T.; Ando, T. DNS analysis of multiple impinging jets. J. Fluid Sci. Technol. 2014, 9, 27. [Google Scholar] [CrossRef] [Green Version]
  10. Geers, L.F.G.; Hanjalic, K.; Tummers, M.J. Experimental investigation of impinging jet arrays. Exp. Fluids 2004, 36, 946–958. [Google Scholar] [CrossRef]
  11. Geers, L.F.G.; Hanjalic, K.; Tummers, M.J. Particle imaging velocimetry-based identification of coherent structures in normally impinging multiple jets. Phys. Fluids 2005, 17, 1–13. [Google Scholar] [CrossRef] [Green Version]
  12. Geers, L.F.G.; Hanjalic, K.; Tummers, M.J. Wall imprint of turbulent structures and heat transfer in multiple impinging jet arrays. J. Fluid Mech. 2006, 546, 255–284. [Google Scholar] [CrossRef] [Green Version]
  13. Geers, L.F.G.; Tummers, M.J.; Bueninck, T.J.; Hanjalic, K. Heat transfer correlation for hexagonal and in-line arrays of impinging jets. Int. J. Heat Mass Transf. 2008, 51, 5389–5399. [Google Scholar] [CrossRef]
  14. Draksler, M.; Končar, B.; Cizelj, L.; Ničeno, B. Large Eddy Simulation of multiple impinging jets in hexagonal configuration—Flow dynamics and heat transfer characteristics. Int. J. Heat Mass Transf. 2017, 109, 16–27. [Google Scholar] [CrossRef]
  15. Hadžiabdić, M.; Hanjalić, K. Vortical structures and heat transfer in a round impinging jet. J. Fluid Mech. 2008, 596, 221–260. [Google Scholar] [CrossRef]
  16. Uddin, N.; Neumann, S.; Weigand, B. LES simulations of an impinging jet: On the origin of the second peak in the Nusselt number distribution. Int. J. Heat Mass Transf. 2013, 57, 356–368. [Google Scholar] [CrossRef]
  17. Durbin, P.A. Separated flow computations with the kεv2 model. AIAA J. 1995, 33, 659–664. [Google Scholar] [CrossRef]
  18. Billard, F.; Laurence, D. A robust kεv2/k elliptic blending turbulence model applied to near-wall, separated and buoyant flows. Int. J. Heat Fluid Flow 2012, 33, 45–58. [Google Scholar] [CrossRef]
  19. Hanjalic, K.; Popovac, M.; Hadziabdic, M. A robust near-wall elliptic-relaxation eddy-viscosity turbulence model for CFD. Int. J. Heat Fluid Flow 2004, 25, 1047–1051. [Google Scholar] [CrossRef]
  20. Lien, F.; Kalitzin, G. Computations of transonic flow with the v2f turbulence model. Int. J. Heat Fluid Flow 2001, 22, 53–61. [Google Scholar] [CrossRef]
  21. Kenjeres, S.; Gunarjo, S.B.; Hanjalic, K. Contribution to elliptic relaxation modelling of turbulent natural and mixed convection. Int. J. Heat Fluid Flow 2005, 26, 569–586. [Google Scholar] [CrossRef]
  22. Behnia, M.; Parneix, S.; Durbin, P.A. Prediction of heat transfer in an axisymmetric turbulent jet impinging on a flat plate. Int. J. Heat Mass Transf. 1998, 41, 1845–1855. [Google Scholar] [CrossRef]
  23. Behnia, M.; Parneix, S.; Shabany, Y.; Durbin, P.A. Numerical study of turbulent heat transfer in confined and unconfined impinging jets. Int. J. Heat Fluid Flow 1999, 20, 1–9. [Google Scholar] [CrossRef]
  24. Cooper, D.; Jackson, D.C.; Launder, B.E.; Liao, G.X. Impinging jet studies for turbulence model assessment: Part I, flow field experiments. Int. J. Heat Mass Transf. 1993, 36, 2675–2684. [Google Scholar] [CrossRef]
  25. Baughn, J.W.; Shimizu, S. Heat transfer measurements from a surface with uniform heat flux and an impinging jet. J. Heat Transf. 1989, 111, 1096–1098. [Google Scholar] [CrossRef]
  26. Yan, X. A Preheated-Wall Transient Method Using Liquid Crystals for the Measurement of Heat Transfer on External Surfaces and in Ducts. Ph.D. Thesis, University of California, Davis, CA, USA, 1993. [Google Scholar]
  27. Camerlengo, G.; Borello, D.; Salvagni, A.; Sesterhenn, J. Effects of Wall Curvature on the Dynamics of an Impinging Jet and Resulting Heat Transfer. In Active Flow and Combustion Control. Notes on Numerical Fluid Mechanics and Multidisciplinary Design; Springer: Cham, Switzerland, 2018; pp. 355–366. [Google Scholar]
  28. Hashiehbaf, A.; Baramade, A.; Agrawal, A.; Romano, G. Experimental investigation on an axisymmetric turbulent jet impinging on a concave surface. Int. J. Heat Fluid Flow 2015, 53, 167–182. [Google Scholar] [CrossRef]
  29. Yang, G.; Choi, M.; Lee, J.S. An experimental study of slot jet impingement cooling on concave surface: Effects of nozzle configuration and curvature. Int. J. Heat Mass Transf. 1999, 42, 2199–2209. [Google Scholar] [CrossRef]
  30. Zhou, Y.; Lin, G.; Bu, X.; Bai, L.; Wen, D. Experimental study of curvature effects on jet impingement heat transfer on concave surfaces. Chin. J. Aeronaut. 2017, 30, 586–594. [Google Scholar] [CrossRef]
  31. Craft, T.J.; Iacovides, H.; Mostafa, N.A. Modelling of three-dimensional jet array impingement and heat transfer on a concave surface. Int. J. Heat Fluid Flow 2008, 29, 687–702. [Google Scholar] [CrossRef]
  32. Fenot, M.; Dorignac, E.; Vullierme, J.-J. An experimental study on hot round jets impinging a concave surface. Int. J. Heat Fluid Flow 2008, 29, 945–956. [Google Scholar] [CrossRef]
  33. Lee, D.H.; Kim, S.J.; Kim, Y.H.; Park, H.J. Heat transfer with fully developed slot jets impinging on confined concave and convex surfaces. Int. J. Heat Mass Transf. 2015, 88, 218–223. [Google Scholar] [CrossRef]
  34. Kura, T.; Fornalik-Wajs, E.; Wajs, J. Thermal and hydraulic phenomena in boundary layer of minijets impingement on curved surfaces. Arch. Thermodyn. 2018, 39, 147–166. [Google Scholar]
  35. Wajs, J.; Mikielewicz, D.; Fornalik-Wajs, E. Cylindrical Jet Heat Exchanger Dedicated to Heat Recovery, Especially From Low Temperature Waste Sources. Patent PL224494, 8 July 2013. [Google Scholar]
  36. Wajs, J.; Mikielewicz, D.; Fornalik-Wajs, E.; Bajor, M. Recuperator with microjet technology as a proposal for heat recovery from low-temperature sources. Arch. Thermodyn. 2015, 36, 49–64. [Google Scholar] [CrossRef] [Green Version]
  37. Wajs, J.; Mikielewicz, D.; Fornalik-Wajs, E.; Bajor, M. High Performance Tubular Heat Exchanger with Minijet Heat Transfer Enhancement. Heat Transf. Eng. 2019, 40, 772–783. [Google Scholar] [CrossRef]
  38. Menter, F.R.; Smirnov, P.E.; Liu, T.; Avancha, R. A One-Equation Local Correlation-Based Transition Model. FlowTurbul. Combust. 2015, 95, 583–619. [Google Scholar] [CrossRef]
  39. Kura, T.; Fornalik-Wajs, E.; Wajs, J.; Kenjeres, S. Turbulence models impact on the flow and thermal analyses of jet impingement. Matec Web Conf. 2018, 240, 01016. [Google Scholar] [CrossRef]
  40. Kura, T.; Fornalik-Wajs, E.; Wajs, J.; Kenjeres, S. Local Nusselt number evaluation in the case of jet impingement. J. Phys. Conf. Ser. 2018, 1101, 012018. [Google Scholar] [CrossRef] [Green Version]
  41. Ramadani, K.; Kelterer, M.E.; Pecnik, R.; Sanz, W. Application of Zeta-f-Turbulence Model to Steady Transitional Flow. Ercoftac. Bull. 2009, 80, 35–40. [Google Scholar]
  42. Garimella, S.V.; Nenaydykh, B. Nozzle-Geometry Effects in Liquid Jet Impingement Heat Transfer. Int. J. Heat Mass Transf. 1996, 39, 2915–2923. [Google Scholar] [CrossRef]
  43. Lee, D.H.; Song, J.; Jo, M.C. The effects of nozzle diameter on impinging jet heat transfer and fluid flow. J. Heat Transf. 2004, 126, 554–557. [Google Scholar] [CrossRef]
  44. Kura, T.; Wajs, J.; Fornalik-Wajs, E.; Kenjeres, S.; Gurgul, S. Thermal and Hydrodynamic Phenomena in the Stagnation Zone—Impact of the Inlet Turbulence Characteristics on the Numerical Analyses. Energies 2021, 14, 105. [Google Scholar] [CrossRef]
  45. OpenFOAM Extended Code Guide. Available online: https://www.openfoam.com/documentation/guides/latest/doc/> (accessed on 18 February 2021).
  46. Aillaud, P.; Duchaine, F.; Gicquel, L.Y.M.; Didorally, S. Secondary peak in the Nusselt number distribution of impinging jet flows: A phenomenological analysis. Phys. Fluids 2016, 28, 095110. [Google Scholar] [CrossRef]
  47. Lee, J.; Lee, S.J. Stagnation Region Heat Transfer of a Turbulent Axisymmetric Jet Impingement. Exp. Heat Transf. 1999, 12, 137–156. [Google Scholar] [CrossRef]
  48. O’Donovan, T.S.; Murray, D.B. Jet impingement heat transfer–Part I: Mean and root-mean-square heat transfer and velocity distributions. Int. J. Heat Mass Transf. 2007, 50, 3291–3301. [Google Scholar] [CrossRef]
  49. Kura, T.; Fornalik-Wajs, E.; Wajs, J.; Kenjeres, S. Heat transfer intensification by jet impingement—Numerical analysis using RANS approach. E3S Web Conf. 2019, 108, 01025. [Google Scholar] [CrossRef]
  50. Zhu, Y.; Antonia, R.A. Correlation between the Enstrophy and the Energy Dissipation Rate in a Turbulent Wake. In Proceedings of the Advances in Turbulence VI. Proceedings of the Sixth European Turbulence Conference, Lausanne, Switzerland, 2–5 July 1996; pp. 507–510.
Figure 1. The design idea of the heat exchanger of Wajs et al. [35,36,37]. Red arrows indicate the hot fluid flow, while blue arrows indicate a flow of cold fluid. Orifices, in which the jets are generated, are visible on the right side of the image. The solid wall, which constricts the cold fluid, is not shown, for clarity.
Figure 1. The design idea of the heat exchanger of Wajs et al. [35,36,37]. Red arrows indicate the hot fluid flow, while blue arrows indicate a flow of cold fluid. Orifices, in which the jets are generated, are visible on the right side of the image. The solid wall, which constricts the cold fluid, is not shown, for clarity.
Energies 14 01846 g001
Figure 2. Definition of the problem.
Figure 2. Definition of the problem.
Energies 14 01846 g002
Figure 3. Geometric constraints used in the simulation and the placement of boundary conditions for (a) convex and (b) concave geometries.
Figure 3. Geometric constraints used in the simulation and the placement of boundary conditions for (a) convex and (b) concave geometries.
Energies 14 01846 g003
Figure 4. Comparison of the flat surface jet impingement, local heat transfer results, represented by Nusselt number, in relation to the RANS turbulence model. Green lines indicate the results obtained with the ζ-f, depending on the turbulence boundary conditions at the inlet. Points indicate the reference experimental data from the ERCOFTAC database. Reynolds number Re = 23,000, H/D = 2.
Figure 4. Comparison of the flat surface jet impingement, local heat transfer results, represented by Nusselt number, in relation to the RANS turbulence model. Green lines indicate the results obtained with the ζ-f, depending on the turbulence boundary conditions at the inlet. Points indicate the reference experimental data from the ERCOFTAC database. Reynolds number Re = 23,000, H/D = 2.
Energies 14 01846 g004
Figure 5. Mesh independence test example. (a) H/D = 1, R/D = 4, Re = 5000, convex surface. (b) H/D = 1, R/D = 8, Re = 23,000, concave surface.
Figure 5. Mesh independence test example. (a) H/D = 1, R/D = 4, Re = 5000, convex surface. (b) H/D = 1, R/D = 8, Re = 23,000, concave surface.
Energies 14 01846 g005
Figure 6. Methodology for determination of the reference point x/D, (a) convex and (b) concave geometries.
Figure 6. Methodology for determination of the reference point x/D, (a) convex and (b) concave geometries.
Energies 14 01846 g006
Figure 7. Local Nusselt number values, obtained at the heated surface, for Reynolds number Re = 5000. (a) R/D = 4, (b) R/D = 6, (c) R/D = 8.
Figure 7. Local Nusselt number values, obtained at the heated surface, for Reynolds number Re = 5000. (a) R/D = 4, (b) R/D = 6, (c) R/D = 8.
Energies 14 01846 g007
Figure 8. Local Nusselt number values, obtained at the heated surface, for Reynolds number Re = 10,000. (a) R/D = 4, (b) R/D = 6, (c) R/D = 8.
Figure 8. Local Nusselt number values, obtained at the heated surface, for Reynolds number Re = 10,000. (a) R/D = 4, (b) R/D = 6, (c) R/D = 8.
Energies 14 01846 g008
Figure 9. Local Nusselt number values, obtained at the heated surface, for Reynolds number Re = 23,000. (a) R/D = 4, (b) R/D = 6, (c) R/D = 8.
Figure 9. Local Nusselt number values, obtained at the heated surface, for Reynolds number Re = 23,000. (a) R/D = 4, (b) R/D = 6, (c) R/D = 8.
Energies 14 01846 g009
Figure 10. Local Nusselt number values: (a) convex case, H/D = 1, R/D = 4, (b) concave case, H/D = 4, R/D = 8.
Figure 10. Local Nusselt number values: (a) convex case, H/D = 1, R/D = 4, (b) concave case, H/D = 4, R/D = 8.
Energies 14 01846 g010
Figure 11. Stagnation point Nusselt number values, x/D = 0. Top: Reynolds number Re = 23,000, middle: Re = 10,000, bottom: Re = 5000.
Figure 11. Stagnation point Nusselt number values, x/D = 0. Top: Reynolds number Re = 23,000, middle: Re = 10,000, bottom: Re = 5000.
Energies 14 01846 g011
Figure 12. Nusselt number values averaged over the surface, up to x/D = 2. Top: Reynolds number Re = 23,000, middle: Re = 10,000, bottom: Re = 5000.
Figure 12. Nusselt number values averaged over the surface, up to x/D = 2. Top: Reynolds number Re = 23,000, middle: Re = 10,000, bottom: Re = 5000.
Energies 14 01846 g012
Figure 13. Stagnation zone velocity field: (a) convex and (b) concave geometry. H/D = 2; R/D = 4. Re = 23,000.
Figure 13. Stagnation zone velocity field: (a) convex and (b) concave geometry. H/D = 2; R/D = 4. Re = 23,000.
Energies 14 01846 g013
Figure 14. Local values of turbulence kinetic energy, extracted from the first mesh nodes along the impinged wall. H/D = 2 and R/D = 4. (a) Re = 5000, (b) Re = 10,000 and (c) Re = 23,000. Green lines indicate the second maximum of Nusselt number for the same case.
Figure 14. Local values of turbulence kinetic energy, extracted from the first mesh nodes along the impinged wall. H/D = 2 and R/D = 4. (a) Re = 5000, (b) Re = 10,000 and (c) Re = 23,000. Green lines indicate the second maximum of Nusselt number for the same case.
Energies 14 01846 g014
Figure 15. Relation between the (a,c) turbulence kinetic energy (TKE) dissipation and (b,d) local enstrophy. H/D = 2; R/D = 4, Re = 23,000.
Figure 15. Relation between the (a,c) turbulence kinetic energy (TKE) dissipation and (b,d) local enstrophy. H/D = 2; R/D = 4, Re = 23,000.
Energies 14 01846 g015
Figure 16. Local values of (a) dissipation of turbulence kinetic energy, ε, and (b) enstrophy, e, extracted from the first mesh nodes along the impinged wall. H/D = 2 and R/D = 4, Re = 23,000.
Figure 16. Local values of (a) dissipation of turbulence kinetic energy, ε, and (b) enstrophy, e, extracted from the first mesh nodes along the impinged wall. H/D = 2 and R/D = 4, Re = 23,000.
Energies 14 01846 g016
Figure 17. Values of the v 2 ¯ obtained along the centreline of the impinging jet, R/D = 4.
Figure 17. Values of the v 2 ¯ obtained along the centreline of the impinging jet, R/D = 4.
Energies 14 01846 g017
Figure 18. Velocity profiles at particular values of h/H from the heated surface. H/D = 1 and R/D = 4, Re = 5000.
Figure 18. Velocity profiles at particular values of h/H from the heated surface. H/D = 1 and R/D = 4, Re = 5000.
Energies 14 01846 g018
Figure 19. Turbulence kinetic energy profiles at particular values of h/H from the heated surface. H/D = 1 and R/D = 4, Re = 5000.
Figure 19. Turbulence kinetic energy profiles at particular values of h/H from the heated surface. H/D = 1 and R/D = 4, Re = 5000.
Energies 14 01846 g019
Figure 20. Enstrophy flux profiles at particular values of h/H from the heated surface. H/D = 1 and R/D = 4, Re = 5000. Each plot has different scale of Y axis.
Figure 20. Enstrophy flux profiles at particular values of h/H from the heated surface. H/D = 1 and R/D = 4, Re = 5000. Each plot has different scale of Y axis.
Energies 14 01846 g020
Figure 21. Velocity profiles at particular values of h/H from the heated surface H/D = 2 and R/D = 8, Re = 23,000.
Figure 21. Velocity profiles at particular values of h/H from the heated surface H/D = 2 and R/D = 8, Re = 23,000.
Energies 14 01846 g021
Figure 22. Turbulence kinetic energy profiles at particular values of h/H from the heated surface. H/D = 2 and R/D = 8, Re = 23,000.
Figure 22. Turbulence kinetic energy profiles at particular values of h/H from the heated surface. H/D = 2 and R/D = 8, Re = 23,000.
Energies 14 01846 g022
Figure 23. Enstrophy flux profiles at particular values of h/H from the heated surface. H/D = 2 and R/D = 8, Re = 23,000. Each plot has different scale of Y axis.
Figure 23. Enstrophy flux profiles at particular values of h/H from the heated surface. H/D = 2 and R/D = 8, Re = 23,000. Each plot has different scale of Y axis.
Energies 14 01846 g023
Table 1. Geometrical variables used in the paper, as well as inlet Reynolds number (based on the bulk inlet velocity ub).
Table 1. Geometrical variables used in the paper, as well as inlet Reynolds number (based on the bulk inlet velocity ub).
Surface TypeH/DR/DReynolds Number
Concave
convex
1
2
4
4
6
8
5000
10,000
23,000
Table 2. Properties of the working fluid, air.
Table 2. Properties of the working fluid, air.
PropertyValue
density, kg/m31.225
kinematic viscosity, m2/s1.5·10−5
specific heat capacity, J/(kg∙K)1.005·103
Pr (Prandtl number)0.713
Prt (turbulent Prandtl number)1
Table 3. Location of secondary maximum of local Nusselt number.
Table 3. Location of secondary maximum of local Nusselt number.
CaseConvexConcave
x/D
Re = 23,000; H/D = 1; R/D = 41.892.02
Re = 23,000; H/D = 1; R/D = 61.871.97
Re = 23,000; H/D = 1; R/D = 81.891.99
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kura, T.; Fornalik-Wajs, E.; Wajs, J.; Kenjeres, S. Curved Surface Minijet Impingement Phenomena Analysed with ζ-f Turbulence Model. Energies 2021, 14, 1846. https://doi.org/10.3390/en14071846

AMA Style

Kura T, Fornalik-Wajs E, Wajs J, Kenjeres S. Curved Surface Minijet Impingement Phenomena Analysed with ζ-f Turbulence Model. Energies. 2021; 14(7):1846. https://doi.org/10.3390/en14071846

Chicago/Turabian Style

Kura, Tomasz, Elzbieta Fornalik-Wajs, Jan Wajs, and Sasa Kenjeres. 2021. "Curved Surface Minijet Impingement Phenomena Analysed with ζ-f Turbulence Model" Energies 14, no. 7: 1846. https://doi.org/10.3390/en14071846

APA Style

Kura, T., Fornalik-Wajs, E., Wajs, J., & Kenjeres, S. (2021). Curved Surface Minijet Impingement Phenomena Analysed with ζ-f Turbulence Model. Energies, 14(7), 1846. https://doi.org/10.3390/en14071846

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop