Next Article in Journal
FSNB-YOLOV8: Improvement of Object Detection Model for Surface Defects Inspection in Online Industrial Systems
Next Article in Special Issue
A Validated Thermal Computational Fluid Dynamics Model of Wine Warming in a Glass
Previous Article in Journal
A Method to Estimate Dynamic Pore Water Pressure Growth of Saturated Sand-Gravel Materials
Previous Article in Special Issue
Asymmetry Propagation in a Pipe Flow Downstream of a 90° Sharp Elbow Bend
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Evaluation of Turbulent Supersonic Jet Impinging on Inclined Plate

by
Antonio Mezzacapo
and
Giuliano De Stefano
*
Engineering Department, University of Campania Luigi Vanvitelli, Via Roma 29, 81031 Aversa, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(17), 7910; https://doi.org/10.3390/app14177910
Submission received: 7 August 2024 / Revised: 26 August 2024 / Accepted: 2 September 2024 / Published: 5 September 2024
(This article belongs to the Special Issue Applied Computational Fluid Dynamics and Thermodynamics)

Abstract

:
A computational fluid dynamics investigation of a turbulent supersonic jet impacting a solid flat plate is conducted utilizing the OpenFOAM software. The research focuses on simulating the three-dimensional mean compressible flow for jet impingement on an inclined plate by analyzing the complex flow field and the surface distribution of pressure. Various simulations are carried out at a jet Mach number of 2.2 maintaining a constant nozzle-to-plate distance while varying the angle of inclination of the plate. In contrast to earlier numerical studies, this work employs a modern turbulence modeling technique known as detached eddy simulation (DES), along with a traditional unsteady Reynolds-averaged Navier–Stokes model. Making a comparison with experimental findings, the current analysis reveals that both turbulence modeling techniques effectively predict the mean pressure distribution on the plate. However, the DES approach offers deeper insights into the turbulent flow field, showing notable consistency with the experiments. The complex compressible flow patterns are simulated with higher accuracy compared to the traditional approach. Enhanced turbulence resolution is attained by utilizing the same computational grid with a limited increase in computational complexity.

1. Introduction

Despite decades of research and progress, accurate and efficient computational fluid dynamics (CFD) simulation of wall-bounded compressible turbulent flows still represents a challenging task, e.g., [1]. Within this framework, impinging supersonic jet flows continue to receive great attention, due to their technological importance as well as high complexity, in both theoretical and applied engineering research. In particular, the impingement of a turbulent supersonic jet on an inclined flat plate appears in various aerospace engineering applications, including rocket launchers, vertical/short-takeoff-and-landing aircraft operation, and jet-engine thrust control [2,3]. Moreover, this particular jet flow configuration has been also employed in gas-liquid separators, heat exchanges, and mixing devices. In all these cases, the jet flow characteristics play a crucial role in enhancing mixing, as well as heat and mass transfer between the fluid phases, e.g., [4].
Generally, using applied CFD techniques makes it less expensive to examine a wide range of flow parameters and configurations, compared to experimental works. In this particular case, notwithstanding the geometric simplicity, the interaction between the jet flow and the inclined solid surface results in highly intricate fluid dynamic phenomena, such as flow separation, reattachment, vortex dynamics, and complex interactions between shocks and boundary layers. The significant computational complexity has directed research efforts primarily towards the examination of perpendicular jet impingement, e.g., [5]. However, recent advancements in computational power accessible to industrial research, coupled with the evolution of sophisticated numerical methods for applied CFD, have led to an increased focus on supersonic jet flow dynamics, including more complex configurations, e.g., [6]. Specifically, engineering research has taken into account the variations in nozzle-to-wall distance, nozzle exit diameter, and impingement angle [7]. In fact, these geometric parameters, in conjunction with the nozzle pressure ratio, play a crucial role in determining the characteristics of impinging jet flow. For instance, an increase in the angle of inclination of the plate results in heightened surface pressure oscillations, due to more pronounced flow separation and reattachment phenomena, e.g., [8].
The present study specifically examines the scenario of near-wall impingement, characterized by a ratio of L/ D j < 5 , where L denotes the distance from the nozzle to the plate (measured along the centerline of the nozzle) and D j refers to the diameter of the nozzle exit section. In this context, the pressure exerted on the plate is significantly augmented, owing to the intricate flow structures that emerge from the intensified interactions between shock waves and expansion waves. Indeed, experimental research indicates that the pressure load resulting from the supersonic jet striking an inclined solid surface can exceed that observed at normal incidence, primarily due to the presence of multiple shock waves between the nozzle exit and the plate [9]. Furthermore, various studies have emphasized that the compressible flow structures formed at the zone of impingement are also affected by the fluid turbulence, e.g., [10].
In the initial experimental studies focused on near-wall impingement, pressure taps were employed to assess the pressure distribution on the plate. Typically, these measurements yield pressure fields characterized by low spatial resolution, which restricts further analysis, particularly in small-scale experiments. Additionally, in certain instances, the use of physical probes for measurement is not feasible, as they may significantly modify the flow field. In response to these limitations, the pressure-sensitive paint (PSP) technique was introduced to ascertain the pressure distribution with relatively high spatial resolution, thus serving as an effective tool for the examination of jet impingement [11]. For example, Nakai and collaborators [12] explored the behavior of supersonic jets impacting inclined flat plates at various nozzle-to-plate distances ( L / D j ), impingement angles ( θ ), and nozzle pressure ratios ( P R ). The latter parameter is defined as P R = P j / P , where P j and P represent nozzle exit and ambient pressures, respectively. By integrating the high-resolution pressure maps obtained from PSP technology with Schlieren imaging, a methodology was developed to anticipate the overall impinging jet configuration. These experimental findings revealed the existence of various flow patterns, which are dependent on the characteristics of the impinging jet [13]. Additionally, Ito et al. [14] conducted experiments on jet impingement at elevated plate angles, systematically varying the position and inclination of the plate within the ranges of 1 L / D j 4.5 and 60° ≤ θ ≤ 90°, respectively.
It is crucial to note that CFD analysis effectively replicates the surface pressure distribution over the plate with the desired spatial resolution. The precise assessment of the pressure field greatly enhances flow diagnostics and, thus, industrial design, particularly in predicting and mitigating pressure loading and improving system reliability, which is especially vital for aerospace and defense applications. For example, considering the particular application to the takeoff of an aircraft (or a missile), the local increase in velocity results in a decrease in local pressure near the impingement area, which can lead to diminished lift. This reduction in lift may result in either a failed takeoff or increased fuel consumption. In recent years, several experimental studies have been conducted to explore the lift loss induced by this phenomenon, utilizing a simplified test configuration that includes a baffle plate surrounding a single air jet. For instance, a series of experiments examining various nozzle-to-ground heights were presented in [15], indicating that a closer proximity to the ground surface correlates with a greater loss of lift.
In this demanding context, the primary objective of this work is to conduct a computational investigation of the impinging jet flow field utilizing a modern turbulence model. Unlike earlier numerical investigations, such as those presented in [16], which relied on conventional unsteady Reynolds-averaged Navier–Stokes (RANS) models, this research also employs a hybrid turbulence modeling procedure known as detached-eddy simulation (DES). Indeed, traditional RANS simulations are limited in their ability to capture spectral content in the flow field, even when sufficient spatio-temporal resolution is achieved, as they compute the (statistical) mean flow, thereby eliminating all turbulence information from the resolved flow field through the averaging process. Although multi-scale RANS models that can accurately predict flow separation details and anisotropic turbulence are still in development, e.g., [17], the growing accessibility of high-performance computing (HPC) resources for industrial research enables more precise evaluation of turbulent supersonic jet impingement. Specifically, the current DES formulation, which is based on the hybridization of the two-equation eddy-viscosity k ω shear-stress transport (SST) model [18], akin to large-eddy simulation (LES), facilitates the emergence of a resolved turbulent spectrum, allowing for partial resolution of turbulent fluctuations. To examine the capability of the two different techniques to reproduce the complex flow physics arising from the jet impingement, various computations are conducted for varying inclination angles.
Furthermore, contrasting with similar investigations where commercial proprietary software is usually employed, the use of an open-source CFD solver allows us to overcome a couple of significant limitations. In fact, these solvers do not provide access to the underlying code, so users cannot modify the algorithms or scrutinize the underlying physical-mathematical model. Moreover, even a basic-level solver license can be prohibitively expensive, restricting the possibility of performing high-fidelity simulations. The ability of the selected CFD solver to correctly simulate wall-confined supersonic flows is tested by comparing the present results against both numerical data obtained with a commercial code, and experimental findings. Particular attention is devoted to the very complex flow structures created by the jet impact on the inclined plate, and the consequent pressure distribution on the solid surface, while examining the effect of varying the plate inclination.
The remainder of the paper is organized as follows. Section 2 introduces the overall computational model, including the meshing procedure and the main numerical settings, while the two different turbulence models employed are given in Section 3. The results of the various calculations are presented and discussed in Section 4. Finally, Section 5 provides some conclusions and perspectives of this study.

2. Computational Model

The flow configuration that was investigated comprises a supersonic air jet originating from an under-expanded diverging nozzle, and impacting on an inclined flat plate. The complex compressible flow field was simulated for varying plate inclination angles, once both the nozzle pressure ratio and the nozzle-to-plate distance were fixed. The evolving fluid was considered as an ideal gas, while the Sutherland law was used to model the dependence of dynamic viscosity μ on the temperature. A number of different computations were conducted by employing the open-source CFD software OpenFOAM 9, which is based on the finite volume (FV) discretization scheme.

2.1. Flow Domain and Boundary Conditions

The flow geometry consists of a diverging axisymmetric nozzle, generating the turbulent supersonic jet, together with a large receiver chamber containing the inclined flat plate. The large computational domain was prescribed in order to accurately capture most of the complex flow field features. It extends up to 100 D j , 64 D j , 64 D j , and 71 D j in the longitudinal, upper, lower, and lateral directions, respectively, where the nozzle exit diameter D j is assumed as the reference length. The nozzle geometry is designed to ensure the desired Mach number at the exit section that is M j = 2.2 . The rectangular flat plate has the following extents: 10 D j (length), 10 D j (width), and 0.28 D j (thickness). Some important data are summarized in Table 1, while a schematic view of the computational domain is presented in Figure 1, where the X-axis direction is aligned with the nozzle centerline. The flow region of particular interest between the nozzle exit and the plate surface is depicted in the zoomed-in picture. Note, that the size of the present computational domain is actually reduced compared to reference numerical simulations in [16].
At the inlet section, corresponding to the nozzle throat, since the nozzle geometry is given and known, pressure, temperature, and velocity are directly imposed using Dirichlet boundary conditions, with the thermodynamic state being determined by means of isentropic flow relations. Both the impingement plate and nozzle surfaces are treated as adiabatic walls so that corresponding homogeneous Neumann boundary conditions are adopted for the thermal field. In addition, no-slip boundary conditions are imposed on the walls. Importantly, the wave transmissive boundary condition is imposed for both pressure and velocity, in order to avoid the propagation of transmissive and reflective acoustic waves through the computational domain. Following [19], such a particular boundary condition enables the CFD solution to capture the propagation of acoustic and pressure waves with minimal artificial reflections, ensuring that unbounded flow domains are realistically represented. Note, that in this study in contrast with the pure RANS analysis performed in [16], where the symmetry boundary condition at the plane Z = 0 was imposed for the mean flow solution, the (partial) resolution of the fluctuating turbulent flow component makes this condition inadequate. Indeed, owing to the adoption of a hybrid RANS/LES formulation, the jet flow was simulated in the entire computational domain, without resorting to any simplification. The different types of boundary conditions are summarized in Table 2, while the associated values for the field variables are given in Table 3. Note, that the far-field values for the turbulent kinetic energy (k) and the specific dissipation rate ( ω ) were determined by setting the initial level of turbulence intensity to 1%.

2.2. Spatial Discretization

The discretization of the computational domain was performed by using hexahedral elements in the proximity of the plate, aligned with the solid surface, together with tetrahedral elements in the rest of the domain. The sketch of the overall computational grid is reported in Figure 2, where a couple of side views with different zoom levels illustrate the flow region of interest, including the diverging nozzle. Owing to the FV approach, the compressible flow governing equations were numerically solved while inherently satisfying the conservation principles over each computational cell, e.g., [20], which is particularly important for non-uniform grid arrangement [21]. The particular meshing strategy allowed the smooth transition from structured to unstructured meshes, with the computational complexity of the overall CFD model being alleviated. The mesh resolution at the wall was prescribed by controlling the size of the first computational cell away from the solid surface, while ensuring that the non-dimensional height was corresponding to y + 35 , in inner units. In fact, to save computational resources, the turbulence modeling procedure was supplied with wall functions, which allowed us to avoid using more refined grids in the wall region.
The robustness of the present computational model was first assessed by performing a numerical grid convergence analysis for a given inclination angle, namely, θ = 30 ° . Three different grid resolutions were prescribed on the plate surface, say N y × N z , with the resulting overall mesh sizes listed in Table 4. Note, that owing to the prescribed Courant–Friedrichs–Lewy (CFL) condition, the time integration step was consistently reduced by refining the spatial grid. Based on the CFD results that were obtained, the fine resolution could be considered sufficient to accurately predict the surface pressure load on the plate, as demonstrated in Section 4. Also, looking at the relative errors for the pressure peak value measured against the experimental findings, that are given in Table 5, this resolution was needed for a correct estimation of the complex compressible flow field.

2.3. Numerical Settings

In the framework of the OpenFOAM package, the rhoPimpleFoam solver was employed for the current numerical investigations. This transient, density-based, compressible flow solver has been validated over a wide range of low to medium supersonic Mach numbers, showing a good ability to capture shock structures in several different cases. For instance, it was successfully used for freely exhausting jets [22], and transonic flow over airfoils [23].
Due to the hybrid nature of the present DES formulation, suitable numerical schemes should have low dissipation and dispersion errors over the resolved range of scales, while being dissipative for non-resolvable small scales [24]. Therefore, the numerical method used for the gradient discretization was the second-order cellLimited Gauss linear scheme, with a correction factor of 0.5. In fact, due to the non-orthogonality of the mesh away from the plate, such a parameter needed to be used for an accurate estimation of gradients and diffusive fluxes. The divergence terms, including the turbulence quantities, were discretized using the second-order Gauss limitedLinear scheme. Moreover, particular attention was paid to the discretization of the velocity field. Due to the high compressibility nature of the flow, the total variation diminishing (TVD) scheme Gauss limitedLinearV was chosen for the calculation of the velocity divergence. Such a particular scheme, which prevents the creation of unphysical local extrema, allows a corrected shock wave capturing effect in the proximity of the wall, resulting in a more accurate estimation of the pressure distribution on the plate [25]. Also, a linear scheme was used as the default interpolation scheme. Finally, the second-order implicit Adams–Bashforth scheme was selected for the time integration, as it has been shown to provide a good combination of accuracy and stability, compared to other methods such as the Crank–Nicolson scheme, which shows an oscillating behaviour [26]. The time step size was allowed to be adjusted such that the maximum CFL number did not exceed the value of 0.8, which was empirically found to be sufficient for numerical stability.

3. Turbulence Modeling

Two different turbulence modeling procedures were employed to assess the predictive capabilities of the present CFD model to simulate the jet flow impingement on the inclined flat plate, namely, unsteady RANS and DES. The compressible flow governing equations, based on the k- ω formulation, can be formally written the same way for both approaches [27]. Specifically, the balance equations for mass, momentum, and energy read:
ρ t + x i ρ u i = 0 ,
t ρ u i + x j ρ u j u i = p x i + x j t j i + ρ τ j i ,
t ρ e + 1 2 u i u i + k + x j ρ u j h + 1 2 u i u i + k = x j u i t i j + ρ τ i j + μ Pr L + μ T Pr T h x j + μ + σ * ρ k ω k x j ,
where ρ , ρ u i , and ρ E represent the balanced variables per unit volume. The energy Equation (3) ensures conservation of the total energy that is E = e + 1 2 u i u i + k , including internal energy e, kinetic-energy, and turbulent kinetic-energy. Note, that in the present formulation, the energy transport equation is based on the enthalpy variable h = e + p / ρ and, once this equation is solved, the temperature field can be directly estimated.
In the above equations, t i j and τ i j are the viscous and turbulent stresses, μ T the turbulent eddy-viscosity, P r L and P r T the molecular and turbulent Prandtl numbers. Here, according to the eddy-viscosity assumption, the unknown turbulent stresses are approximated as:
τ i j = 2 μ T S i j 1 3 S l l δ i j 2 3 ρ k δ i j ,
where S i j = 1 2 u i x j + u j x i is the resolved strain-rate tensor, and δ i j stands for the Kronecker delta, with μ T being determined in terms of the resolved turbulence variables.

3.1. Unsteady RANS Simulation

In this work, the high Reynolds number version of the k- ω turbulence model proposed by Wilcox [27] was utilized. The turbulence modeling procedure involves the solution of two additional evolution equations for the turbulence variables k and ω , specifically:
t ( ρ k ) + x j ρ u j k = P k D k + x j μ + σ * ρ k ω k x j ,
t ( ρ ω ) + x j ρ u j ω = P ω β ρ ω 2 + σ d ρ ω k x j ω x j + x j μ + σ ρ k ω ω x j .
Compared to the original version [18], such a model variant uses the additional cross-diffusion term at the right-hand side of the ω equation, as well as a stress limiter. These extra terms are applied using switches so that the model formulation maintains its simplicity, without the need for blending functions. In the above equations, the production of turbulent energy is represented by P k = ρ τ i j u i x j , whereas the corresponding dissipation rate is given by D k = β * ρ k ω , while P ω = α ω k P k . Moreover, α , β , β * , σ , σ * , and σ d are constant coefficients for which standard values were selected. For a comprehensive presentation of this two-equation eddy-viscosity model and relevant discussion, the interested reader is referred to [28].

3.2. Detached-Eddy Simulation

The DES method represents a more sophisticated hybrid turbulence modeling formulation that can be adopted in conjunction with various RANS models. In this work, the k ω SST-based DES model is actually employed. The k ω SST model was among the different pure RANS models tested in the previous CFD study used as reference [16]. Therein, a correct surface pressure distribution for the impinging jet flow at low incidence (θ = 30°) was attained, due to the proven ability of this model to provide accurate mean flow results, especially for wall-bounded separated flows [18].
In the current approach, the evolution equations for the two turbulence variables are slightly modified as follows:
    t ( ρ k ) + x j ρ u j k = P k F D E S D k + x j μ + σ k μ t k x j ,
t ( ρ ω ) + x j ρ u j ω = P ω β ρ ω 2 + 2 ( 1 F 1 ) σ ω 2 ρ ω k x j ω x j + x j μ + σ ω μ t ω x j ,
with F D E S given by:
F DES = max 1 F 2 L RANS C DES Δ , 1 ,
where the switching functions F 1 and F 2 are the same as in the standard k ω SST model [18]. In the above definition, L R A N S = k / ( β * ω ) represents the resolved turbulence length scale, Δ is the largest (linear) dimension of the FV cell, with C DES standing for a calibration coefficient.
In practice, RANS and LES solutions are coupled by introducing a hybrid turbulent length scale, which is modified with respect to the pure RANS approach, where it holds D k = ρ k 3 / 2 / L R A N S . This way, a single turbulence closure acts as either a sub-grid scale model, where the local grid resolution is fine enough for LES, or a RANS model, where that is not true. The two different regions are actually separated by a transition zone where the eddy-viscosity level is sufficiently low to allow the rapid regeneration of turbulent eddies. Basically, DES aims at the prediction of separated flows in fluids engineering research at a manageable cost. More details about the DES formulation can be found in the original work by Strelets [29].

4. Results

As a preliminary step, the grid convergence analysis presented in Section 2.2 is demonstrated in Figure 3 by comparing the pressure distributions on the plate mid-section for the different mesh resolutions with each other, as well as against the experimental data in [12]. The three different FV grids are made of 0.66 (coarse), 2.7 (medium), and 6.8 (fine) million cells for the whole computational domain. In this plot, the spatial coordinate along the plate is non-dimensionalized by the nozzle exit radius ( R j D j / 2 ), whereas the pressure variable is non-dimensionalized by the nozzle chamber (stagnation) pressure that is P c = 100 kPa. By inspection of these data, the numerical resolution corresponding to the finest grid can be regarded, a posteriori, as adequate for capturing most of the complex flow features. Therefore, the various calculations for the two different turbulence modeling approaches were conducted using this particular grid.
The pressure contours over the plate surface predicted by the two different methods are presented in Figure 4. These results have been obtained once the statistically steady state is reached, for different inclination angles that are θ = 30°, 40°, and 50°, respectively. Note, that practically since the near-wall region is modeled using the RANS approach for both turbulence models the shape of the pressure pattern of impingement is pretty similar for the two solutions, particularly at low plate incidences. As the inclination angle increases, slight differences start appearing. In fact, the higher the inclination angle, the stronger the jet deviation imposed by the plate (where θ = 90° would correspond to perpendicular jet impingement), and the more complex is the separated compressible flow field. This way, the LES capability of the hybrid model leads to a better description of the resolved pressure field with increasing plate incidence. Apparently, at higher inclinations, the DES solution shows a more elongated impingement zone over the plate, with the stagnation flow region being slightly larger, compared to RANS.
The present qualitative results show a good agreement with reference numerical data in [16], which were obtained using a commercial code. Expectedly, since the plate center is perfectly aligned to the nozzle axis, the maximum pressure level is achieved at the mid-section of the plate, as illustrated in Figure 4. This way, when performing a quantitative analysis and comparison of pressure load using different methods, the plate section at Z = 0 is selected. Figure 5 shows the pressure distribution on the mid-section of the plate for the two different solutions, and varying plate incidence. Making a comparison with reference experimental data [12], both models result in quite correctly reproducing the pressure peak and its location on the plate. At the lowest inclination angle, the numerical solutions are able to predict the pressure oscillations away from the high-pressure region, even though the values and locations of local extrema are not well captured.
The two different CFD solutions are further analyzed by looking at the resolved velocity and temperature fields, for varying plate inclination. Figure 6 shows the normalized speed contours in the mid-plane Z = 0 for the space region between the nozzle exit and the plate surface, while the corresponding thermal contours are depicted in Figure 7. In any case, a stagnation zone is visible in the impingement region, together with the consequent boundary layers development in both directions along the solid surface. It is worth noting that both solutions are able to capture the jet boundary and the jet shock away from the wall with acceptably high fidelity. However, a slower and more consistent velocity field along the plate is obtained with DES, as the RANS computation seems to overpredict the fluid velocity in the wall jet region. Moreover, the DES approach is able to more clearly reveal the presence of the expansion waves that are produced at the exit of the under-expanded nozzle and their reflection from the jet boundary. Such waves turn into compression waves that coalesce to form a barrel-shaped shock. Increasing the plate angle, in accordance with the experimental findings, the flow structure changes and the flow complexity grows. Apparently, as the jet spreads along the plate, the fluid velocity tends to decrease, while the boundary layer thickness increases. Moreover, the stagnation zone tends to shift slightly upwards. This flow scenario is in agreement with previous numerical investigations, where more sophisticated, and more costly, LES techniques were adopted [30]. However, the drift of the stagnation bubble should be further investigated, due to its crucial impact on the evaluation of heat transfer and wall shear stress over the plate. By inspection of the temperature contours in Figure 7, it seems that unsteady RANS tends to over-predict the temperature in the jet core region, as well as in the wall proximity, while the thermal field exhibits a more uniform behavior for DES. Moreover, the temperature in the far-field region appears slightly higher for the RANS solution.
In accordance with previous investigations, the pure RANS solution is shown to provide too slow decay of the mean velocity profile along the jet axis, including the core region, which can be attributed to the underprediction of the turbulence mixing in the shear layer. The nearby temperature field is overestimated, leading to the misrepresentation of the heat transfer rate in the impingement zone. On the contrary, the hybrid RANS/LES approach produces faster mixing, due to the capability to take into account the vortex breakup process in the shear layer of the jet and in the near-wall region, which leads to a more correct estimation of the thermal field in the proximity of the plate [31].
A more detailed description of the complex flow structures arising from the jet impingement is presented in the following, using the numerical Schlieren technique. In Figure 8, the various flow scenarios associated with different plate inclinations are visualized by reporting the contour maps of the density gradient modulus, making a direct comparison against experimental Schlieren images provided by Nakai et al. [12]. For the numerical data, the contour level varies between zero (white color) and 0.053 (yellow color). Indeed, based on the experimental findings, it was found that the complex flow patterns can be classified into different types, namely: Type I (with/without bubble), Type II (with/without bubble), and Type III, as discussed in [13,14]. The actual flow type can be predicted in advance, once the plate angle, the nozzle-to-plate distance, and the shock cell length are given and known. Apparently, at θ = 30°, both solutions correctly resemble the experimental flow field, with no bubble being present. In this case, the flow can be classified as Type III. A good agreement with experiments is also found at θ = 40°, where the plate shock forming at the upper side of the plate is predicted by both solutions, confirming what was found in [16]. According to these results, the flow can be classified as Type II (with bubble). Note, however, that the plate shock size seems thicker for DES than for unsteady RANS. Lastly, at θ = 50°, the flow type is still of Type II (with bubble). In fact, the intermediate tail shock and the stagnation bubble are clearly recognized in the present Schlieren pictures, consistent with the pressure peak reported in Figure 5. The presence of the bubble on the upper part of the impingement zone is more evident here than for previous numerical investigations [16]. Note, that the current results are fully consistent with the jet flow classifications reported in [14].
Overall, both approaches seem to correctly capture the essential flow structures, compared to experiments. However, DES is shown to predict some fluid flow structures more accurately, both near and away from the plate. For example, jet boundary and intermediate tail shock are captured with high detail by the hybrid model, whereas the traditional RANS model shows a lower level of resolution. Importantly, enhanced turbulence resolution is attained with DES by utilizing the same computational grid and, thus, with a limited increase in computational complexity.
The DES approach, which partly resolves the turbulent fluctuations, is able to provide an estimation for the turbulent Reynolds stresses. A qualitative analysis of the normal stresses u 2 and v 2 is presented in Figure 9, where the stress level is non-dimensionalized by the inlet velocity squared U * 2 . The axial component is higher in magnitude, in agreement with previous numerical and experimental investigations [32,33]. This fact can be attributed to the compressible effects at the exit of the under-expanded supersonic nozzle. Moreover, the amplitude of the velocity fluctuations in the axial direction tends to increase along the plate, as the inclination gets higher. Practically, a rise in the plate incidence leads to an increase in the level of disturbance for the impinging jet flow, resulting in a higher magnitude for the transversal component v 2 along the plate.

5. Conclusions and Perspectives

The open-source software OpenFOAM was utilized for the computational evaluation of a turbulent supersonic jet impinging on an inclined flat plate at various incidences. The proposed applied CFD model was based on the unsteady RANS approach using wall functions, thus avoiding excessive resolution requirements in the wall region. In addition, the modern DES technique, which integrates RANS and LES methods, was tested for this particular application in fluids engineering, with the present results being validated against reference numerical data and experimental findings. According to the analyses carried out, both qualitatively and quantitatively, the current models are shown to correctly reproduce the surface pressure distribution on the inclined plate, depending on the plate inclination. Although more sophisticated turbulence models can be used, this study demonstrates that simpler RANS-based models provide a way to create computationally less expensive simulations, where the mean flow can be predicted with engineering accuracy and moderate computational complexity. However, the DES method was able to better reproduce the very complex structure of the impinging jet flow, combining RANS-like turbulence resolution in the wall proximity with partial capture of the turbulent fluctuations away from the plate.
As far as perspectives are concerned, the deeper analysis of the temperature field will be the next step of this continuing research, in order to thoroughly investigate the heat transfer at the plate. The importance of such an analysis is actually twofold. On one hand, in comparison to other conventional systems, jet impingement cooling has proven to attain higher heat transfer rates. In this context, the computational study contributes to determining the jet flow conditions (including the geometry) ensuring the desired thermal energy drain. On the other hand, in a number of industrial applications, the increasing temperature of the plate has to be controlled, to avoid the undesired deterioration of the thermophysical properties of the solid surface coating. Further investigations will be performed by testing innovative unsteady RANS modeling approaches based on adaptive numerical methods [34], along with enhanced DES formulations, like delayed detached-eddy simulation (DDES) [1,35]. Future research may also focus on applied aeroacoustics, specifically examining the acoustic waves generated by fluid-structure interaction, which can be estimated using the pressure fluctuations predicted by the LES component of the hybrid solution.

Author Contributions

Conceptualization, G.D.S.; methodology, G.D.S.; validation, A.M. and G.D.S.; investigation, A.M. and G.D.S.; resources, A.M. and G.D.S.; data curation, A.M.; writing—original draft preparation, A.M. and G.D.S.; writing—review and editing, A.M. and G.D.S.; visualization, A.M.; supervision, G.D.S. All authors read and agreed to the published version of the manuscript.

Funding

The authors acknowledge financial support under the National Recovery and Resilience Plan (NRRP), Mission 4, Component 2, Investment 1.1, Call for tender No.104 published on 2.2.2022 by the Italian Ministry of University and Research (MUR), funded by the European Union–NextGenerationEU– Project 2022B2X937 Title “Next Generation Space Propulsion Design Techniques”–CUP E53D23003080006- Grant Assignment Decree No.961 adopted on 30.6.2023 by MUR.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors thank Kozo Fujii, Tokyo University of Science, Japan, for kindly providing the experimental images reported in Figure 8. The authors acknowledge the CINECA award under the ISCRA initiative (Project HP10BB0905) for the availability of high performance computing resources.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CFDcomputational fluid dynamics
CFLCourant–Friedrichs–Lewy (condition)
DESdetached-eddy simulation
FVfinite volume (method)
HPChigh-performance computing
LESlarge-eddy simulation
PSPpressure-sensitive paint (method)
RANSReynolds-averaged Navier–Stokes
SSTshear-stress transport
TVDtotal variation diminishing

References

  1. Ge, X.; De Stefano, G.; Hussaini, M.Y.; Vasilyev, O.V. Wavelet-based adaptive eddy-resolving methods for modeling and simulation of complex wall-bounded compressible turbulent flows. Fluids 2021, 6, 331. [Google Scholar] [CrossRef]
  2. Lubert, C.P.; Gee, K.L.; Tsutsumi, S. Supersonic jet noise from launch vehicles: 50 years since NASA SP-8072. J. Acoust. Soc. Am. 2022, 151, 752–791. [Google Scholar] [CrossRef] [PubMed]
  3. Zhu, H.; Guo, H.; Sun, J.; Tian, H.; Cai, G. Research progress on active secondary jet technology in supersonic flow field of aerospace propulsion systems. Fluids 2023, 8, 313. [Google Scholar] [CrossRef]
  4. Li, M.; Karami, S.; Sandberg, R.; Soria, J.; Ooi, A. Heat transfer investigation on under-expanded supersonic impinging jets. Int. J. Heat Fluid Flow 2023, 101, 109132. [Google Scholar] [CrossRef]
  5. Wilke, R.; Sesterhenn, J. Statistics of fully turbulent impinging jets. J. Fluid Mech. 2017, 825, 795–824. [Google Scholar] [CrossRef]
  6. Chin, C.; Li, M.; Harkin, C.; Rochwerger, T.; Chan, L.; Ooi, A.; Risborg, A.; Soria, J. Investigation of the flow structures in supersonic free and impinging jet flows. J. Fluids Eng. 2013, 135, 031202. [Google Scholar] [CrossRef]
  7. Nguyen, T.; Maher, B.; Hassan, Y. Flowfield characteristics of a supersonic jet impinging on an inclined surface. AIAA J. 2020, 58, 1240–1254. [Google Scholar] [CrossRef]
  8. Worden, T.; Gustavsson, J.; Shih, C.; Alvi, F.S. Acoustic measurements of high-temperature supersonic impinging jets in multiple configurations. In Proceedings of the 19th AIAA/CEAS Aeroacoustics Conference, Berlin, Germany, 27–29 May 2013. [Google Scholar]
  9. Lamont, P.J.; Hunt, B.L. The impingement of underexpanded, axisymmetric jets on perpendicular and inclined flat plates. J. Fluid Mech. 1980, 100, 471–511. [Google Scholar] [CrossRef]
  10. 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]
  11. Peng, D.; Liu, Y. Fast pressure-sensitive paint for understanding complex flows: From regular to harsh environments. Exp. Fluids 2020, 61, 8. [Google Scholar] [CrossRef]
  12. Nakai, Y.; Fujimatsu, N.; Fujii, K. Experimental study of underexpanded supersonic jet impingement on an inclined flat plate. AIAA J. 2006, 44, 2691–2699. [Google Scholar] [CrossRef]
  13. Nakai, Y.; Fujimatsu, N.; Fujii, K. Flow classification of the under-expanded super sonic jet impinging on a flat plate. In Proceedings of the 33rd AIAA Fluid Dynamics Conference, Orlando, FL, USA, 23–26 June 2003. [Google Scholar]
  14. Ito, M.; Oyama, A.; Fujii, K.; Hayashi, A.K. Flow field analysis of jet impinging on an inclined flat plate at high plate angles. In Proceedings of the 44th AIAA Aerospace Sciences Meeting, Reno, NV, USA, 9–12 January 2006. [Google Scholar]
  15. Bhargav, V.N.; Mehta, Y.; Kumar, R.; Alvi, F.S. Effect of relative jet momentum in supersonic dual impinging jets. Exp. Fluids 2022, 63, 163. [Google Scholar] [CrossRef]
  16. Dharavath, M.; Chakraborty, D. Numerical simulation of supersonic jet impingement on inclined plate. Def. Sci. J. 2013, 4, 355–362. [Google Scholar] [CrossRef]
  17. Ge, X.; Vasilyev, O.V.; De Stefano, G.; Hussaini, M.Y. Wavelet-based adaptive unsteady Reynolds-averaged Navier–Stokes computations of wall-bounded internal and external compressible turbulent flows. In Proceedings of the 56th AIAA Aerospace Sciences Meeting, Kissimme, FL, USA, 8–12 January 2018. [Google Scholar]
  18. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef]
  19. You, R.G.Y.; New, T.H.; Chan, W.L. Large-eddy simulations of a supersonic impinging jet using OpenFOAM. Computation 2024, 12, 124. [Google Scholar] [CrossRef]
  20. Denaro, F.M.; De Stefano, G. A new development of the dynamic procedure in large-eddy simulation based on a Finite Volume integral approach. Application to stratified turbulence. Theoret. Comput. Fluid Dyn. 2011, 25, 315–355. [Google Scholar] [CrossRef]
  21. Iannelli, P.; Denaro, F.M.; De Stefano, G. A deconvolution-based fourth-order finite volume method for incompressible flows on non-uniform grids. Int. J. Numer. Meth. Fluids 2003, 43, 431–462. [Google Scholar] [CrossRef]
  22. Gutiérrez Marcantoni, L.F.; Tamagno, J.P.; Elaskar, S. rhoCentralRfFoam: An OpenFOAM solver for high speed chemically active flows–Simulation of planar detonations. Comput. Phys. Commun. 2017, 219, 209–222. [Google Scholar] [CrossRef]
  23. Vitulano, M.C.; De Tavernier, D.A.M.; De Stefano, G.; von Terzi, D.A. Numerical analysis of transonic flow over the FFA-W3-211 wind turbine tip airfoil. Wind Energ. Sci. Discuss. 2024, 2024, 1–19. [Google Scholar] [CrossRef]
  24. Montecchia, M.; Brethouwer, G.; Wallin, S.; Johansson, A.V.; Knacke, T. Improving LES with OpenFOAM by minimising numerical dissipation and use of explicit algebraic SGS stress model. J. Turbul. 2019, 20, 697–722. [Google Scholar] [CrossRef]
  25. Sweby, P.K. High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM J. Numer. Anal. 1984, 21, 995–1011. [Google Scholar] [CrossRef]
  26. Seng, S.; Monroy, C.; Malenica, S. On the use of Euler and Crank–Nicolson time-stepping schemes for seakeeping simulations in OpenFOAM. In Proceedings of the VII International Conference on Computational Methods in Marine Engineering, Nantes, France, 15–17 May 2017. [Google Scholar]
  27. Wilcox, D.C. Turbulence Modeling for CFD, 3rd ed.; DCW Industries: La Canada, CA, USA, 2006. [Google Scholar]
  28. Wilcox, D.C. Formulation of the k-ω turbulence model revisited. AIAA J. 2008, 46, 2823–2838. [Google Scholar] [CrossRef]
  29. Strelets, M. Detached eddy simulation of massively separated flows. In Proceedings of the 39th AIAA Aerospace Sciences Meeting, Reno, NV, USA, 8–11 January 2001. [Google Scholar]
  30. Gojon, R.; Bogey, C.; Marsden, O. Large-eddy simulation of supersonic planar jets impinging on a flat plate at an angle of 60 to 90 degrees. In Proceedings of the 21st AIAA/CEAS Aeroacoustics Conference, Dallas, TX, USA, 22–26 June 2015. [Google Scholar]
  31. DeBonis, J.R. Prediction of turbulent temperature fluctuations in hot jets. AIAA J. 2018, 56, 3097–3111. [Google Scholar] [CrossRef]
  32. de Cacqueray, N.; Bogey, C.; Bailly, C. Investigation of a high-Mach-number overexpanded jet using large-eddy simulation. AIAA J. 2011, 49, 2171–2182. [Google Scholar] [CrossRef]
  33. Nguyen, D.T.; Maher, B.; Hassan, Y. Effects of nozzle pressure ratio and nozzle-to-plate distance to flowfield characteristics of an under-expanded jet impinging on a flat surface. Aerospace 2019, 6, 4. [Google Scholar] [CrossRef]
  34. De Stefano, G.; Vasilyev, O.V. Hierarchical adaptive eddy-capturing approach for modeling and simulation of turbulent flows. Fluids 2021, 6, 83. [Google Scholar] [CrossRef]
  35. Salomone, T.; Piomelli, U.; De Stefano, G. Wall-modeled and hybrid large-eddy simulations of the flow over roughness strips. Fluids 2023, 8, 10. [Google Scholar] [CrossRef]
Figure 1. Sketch of computational domain, with zoomed-in picture of the space region of primary interest.
Figure 1. Sketch of computational domain, with zoomed-in picture of the space region of primary interest.
Applsci 14 07910 g001
Figure 2. Sketch of the computational grid, with zoomed-in pictures of the space region of primary interest.
Figure 2. Sketch of the computational grid, with zoomed-in pictures of the space region of primary interest.
Applsci 14 07910 g002
Figure 3. Surface pressure distribution on the plate for varying numerical resolution, compared to experimental data [12].
Figure 3. Surface pressure distribution on the plate for varying numerical resolution, compared to experimental data [12].
Applsci 14 07910 g003
Figure 4. Normalized surface pressure contours for unsteady RANS (top) and DES (bottom), for different plate angles that are θ = 30°, 40°, and 50° (from left to right). The pressure level is non-dimensionalized by the stagnation pressure P c .
Figure 4. Normalized surface pressure contours for unsteady RANS (top) and DES (bottom), for different plate angles that are θ = 30°, 40°, and 50° (from left to right). The pressure level is non-dimensionalized by the stagnation pressure P c .
Applsci 14 07910 g004
Figure 5. Surface pressure distribution on the plate for different inclinations that are θ = 30°, 40°, and 50° (from top to bottom), compared to experimental data [12].
Figure 5. Surface pressure distribution on the plate for different inclinations that are θ = 30°, 40°, and 50° (from top to bottom), compared to experimental data [12].
Applsci 14 07910 g005
Figure 6. Normalized speed contours in the plane Z = 0 for unsteady RANS (left) and DES (right), for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom). The speed level is non-dimensionalized by the inlet velocity U*.
Figure 6. Normalized speed contours in the plane Z = 0 for unsteady RANS (left) and DES (right), for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom). The speed level is non-dimensionalized by the inlet velocity U*.
Applsci 14 07910 g006
Figure 7. Normalized temperature contours in the plane Z = 0 for unsteady RANS (left) and DES (right), for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom). The temperature level is non-dimensionalized by the inlet temperature T*.
Figure 7. Normalized temperature contours in the plane Z = 0 for unsteady RANS (left) and DES (right), for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom). The temperature level is non-dimensionalized by the inlet temperature T*.
Applsci 14 07910 g007
Figure 8. Schlieren pictures for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom): unsteady RANS (left) and DES (middle) solutions compared to experimental images (right) [12].
Figure 8. Schlieren pictures for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom): unsteady RANS (left) and DES (middle) solutions compared to experimental images (right) [12].
Applsci 14 07910 g008
Figure 9. DES solution: turbulent stresses u 2 (left) and v 2 (right) for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom). The stress level is non-dimensionalized by the inlet velocity squared U * 2 .
Figure 9. DES solution: turbulent stresses u 2 (left) and v 2 (right) for different plate angles that are θ = 30°, 40°, and 50° (from top to bottom). The stress level is non-dimensionalized by the inlet velocity squared U * 2 .
Applsci 14 07910 g009
Table 1. Jet flow configuration parameters.
Table 1. Jet flow configuration parameters.
ParameterSymbolValue
Nozzle exit Mach number M j 2.2
Nozzle exit diameter D j 7.08 mm
Nozzle throat diameterD*5.0 mm
Nozzle pressure ratioPR = P j / P 7
Nozzle-to-plate distance L / D j 2
Plate inclination angle θ 30°, 40°, 50°
Table 2. Boundary conditions: selected types in OpenFOAM.
Table 2. Boundary conditions: selected types in OpenFOAM.
PatchVelocityPressureTemperature
InletFixedValueFixedValueFixedValue
OutflowwaveTransmissivewaveTransmissiveadvective
Nozzle surfacenoSlipzeroGradientzeroGradient
Plate surfacenoSlipzeroGradientzeroGradient
Table 3. Boundary conditions: physical parameters.
Table 3. Boundary conditions: physical parameters.
ParameterSymbolValue
Inlet velocityU*315 m/s
Inlet pressureP*52.80 kPa
Inlet temperatureT* 248 K
Farfield (ambient) pressure P 1.34 kPa
Farfield temperature T 298 K
Farfield turbulent energy k 14.88 m2/s2
Farfield dissipation rate ω 1.40 s−1
Table 4. Grid convergence test: spatial resolutions.
Table 4. Grid convergence test: spatial resolutions.
GridOverall SizeSurface Resolution  N y × N z
Coarse0.66 M42 × 19
Medium2.7 M79 × 35
Fine6.8 M100 × 50
Table 5. Grid convergence test: relative error against experiment.
Table 5. Grid convergence test: relative error against experiment.
Grid P / P c Relative Error
Coarse0.12417.2%
Medium0.1369.9%
Fine0.1421.4%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mezzacapo, A.; De Stefano, G. Computational Evaluation of Turbulent Supersonic Jet Impinging on Inclined Plate. Appl. Sci. 2024, 14, 7910. https://doi.org/10.3390/app14177910

AMA Style

Mezzacapo A, De Stefano G. Computational Evaluation of Turbulent Supersonic Jet Impinging on Inclined Plate. Applied Sciences. 2024; 14(17):7910. https://doi.org/10.3390/app14177910

Chicago/Turabian Style

Mezzacapo, Antonio, and Giuliano De Stefano. 2024. "Computational Evaluation of Turbulent Supersonic Jet Impinging on Inclined Plate" Applied Sciences 14, no. 17: 7910. https://doi.org/10.3390/app14177910

APA Style

Mezzacapo, A., & De Stefano, G. (2024). Computational Evaluation of Turbulent Supersonic Jet Impinging on Inclined Plate. Applied Sciences, 14(17), 7910. https://doi.org/10.3390/app14177910

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