Next Article in Journal
Efficiency Improvement of Miniaturized Heat Exchangers
Next Article in Special Issue
On the Use of a Domain Decomposition Strategy in Obtaining Response Statistics in Non-Gaussian Seas
Previous Article in Journal
Monolithic Solvers for Incompressible Two-Phase Flows at Large Density and Viscosity Ratios
Previous Article in Special Issue
A Coupled OpenFOAM-WRF Study on Atmosphere-Wake-Ocean Interaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Evaluation of Mesh Resolution for Large-Eddy Simulation of Internal Flows Using Openfoam †

by
Zahra Seifollahi Moghadam
*,
François Guibault
and
André Garon
Polytechnique Montréal, Mechanical Engineering Department, Montréal, QC H3T 1J4, Canada
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in The 15th OpenFOAM Workshop (OFW15), Washington, DC, USA, 22–25 June 2020.
Fluids 2021, 6(1), 24; https://doi.org/10.3390/fluids6010024
Submission received: 1 November 2020 / Revised: 23 November 2020 / Accepted: 31 December 2020 / Published: 5 January 2021
(This article belongs to the Special Issue Selected Papers from the 15th OpenFOAM Workshop)

Abstract

:
The central aim of this paper is to use OpenFOAM for the assessment of mesh resolution requirements for large-eddy simulation (LES) of flows similar to the ones which occur inside the draft-tube of hydraulic turbines at off-design operating conditions. The importance of this study is related to the fact that hydraulic turbines often need to be operated over an extended range of operating conditions, which makes the investigation of fluctuating stresses crucial. Scale-resolving simulation (SRS) approaches, such as LES and detached-eddy simulation (DES), have received more interests in the recent decade for understanding and mitigating unsteady operational behavior of hydro turbines. This interest is due to their ability to resolve a larger part of turbulent flows. However, verification studies in LES are very challenging, since errors in numerical discretization, but also subgrid-scale (SGS) models, are both influenced by grid resolution. A comprehensive examination of the literature shows that SRS for different operating conditions of hydraulic turbines is still quite limited and that there is no consensus on mesh resolution requirement for SRS studies. Therefore, the goal of this research is to develop a reliable framework for the validation and verification of SRS, especially LES, so that it can be applied for the investigation of flow phenomena inside hydraulic turbine draft-tube and runner at their off-design operating conditions. Two academic test cases are considered in this research, a turbulent channel flow and a case of sudden expansion. The sudden expansion test case resembles the flow inside the draft-tube of hydraulic turbines at part load. In this study, we concentrate on these academic test cases, but it is expected that hydraulic turbine flow simulations will eventually benefit from the results of the current research. The results show that two-point autocorrelation is more sensitive to mesh resolution than energy spectra. In addition, for the case of sudden expansion, the mesh resolution has a tremendous effect on the results, and, so far, we have not capture an asymptotic converging behavior in the results of Root Mean Square (RMS) of velocity fluctuations and two-point autocorrelation. This case, which represents complex flow behavior, needs further mesh resolution studies.

1. Introduction

Hydraulic-turbines are considered a highly reliable power source that can cover an extensive range of operating conditions in response to electricity demand. The importance of the large-eddy simulation (LES) for off-design operating conditions is due to the fact that RANS (Reynolds-averaged Navier–Stokes) studies are mainly capable of accurate flow simulations near the Best Efficiency Point, since optimal swirling flow occurs and the level of turbulence is low. However, at off-design operating conditions, swirling vortices are the main phenomena inside hydraulic turbine draft-tubes and the flow is highly turbulent. Therefore, it is necessary to use high fidelity methods, such as large-eddy simulation, to capture the turbulent fluctuations in the flow.
Analysis of flow phenomena for different components of hydraulic turbines has been performed by a large number of researchers in the last few decades. The URANS (Unsteady Reynolds-averaged Navier–Stokes) is the most common approach used by researchers for numerical analysis of the flow field in hydraulic turbines [1,2,3]. The limitation of URANS in predicting self-induced vortex rope was reported by Foroutan and Yavuzkurt (2012), and the importance of time-dependent boundary conditions with turbulent fluctuations at the inlet of the draft tube was shown [4]. SAS (Scale Adaptive Simulation), due to its capability to dynamically adjust the length scale in the turbulent flow and hence provide more detailed predictions of turbulent flow structures, is often employed for draft tube flow field behavior investigations [5,6]. For instance, in a study by Neto et al. (2012), SAS-SST (Shear Stress Transport) results of velocity components inside a draft tube were validated against experimental data for a Francis turbine at part load [7]. DES (Detached Eddy Simulation) as a Hybrid RANS/LES approach, applies the RANS approach for near wall modeling of turbulence and resolves the core flow region using LES. This approach was employed by Sentyabov et al. (2014) for the analysis of vortex core precession of a Kaplan turbine [8]. A comparison of URANS and DES potentials to predict turbulence statistics for the draft tube was presented by Paik et al. (2005). Although both methods agree in mean velocity field, the results of turbulence statistics show significant discrepancies for URANS [9]. DES was also used for the prediction of pressure pulsation in high-head Francis turbines by Minakov et al. (2015), and an accuracy of 10% was reported for the simulations [10].
A comparison of LES and SAS results for a sudden-expansion test case which resembles the vortex rope of a draft tube at part load was performed by Javadi and Nilsson (2014) [11]. Rotating stall mechanism of a pump-turbine was investigated by Pacot et al. (2014) using LES [12]. The results of efficiency prediction of a Kaplan turbine using Zonal LES (ZLES) is compared to SAS in a study by Morgut et al. (2015) [13]. LES with a cavitation model for Francis turbine is also performed at part load and high load in a study by Yang et al. (2016). In this study, however, the LES mesh requirements were not well-satisfied [14]. Comparison of LES results with DES and URANS results for the Francis-99 draft tube at part load, Best Efficiency Point (BEP), and high load was performed by Minakov et al. (2017) for mean velocity profiles and pressure fluctuations [15]. Compressible LES is performed by Trivedi and Dahlhaug (2018) for a whole configuration of the Francis-99 turbine to compare the high-amplitude stochastic fluctuations at speed no-load [16] and runaway [17].
In the papers featuring LES studies for hydraulic turbines, the effect of the mesh has barely been studied. Moreover, there is no consensus on the space and time resolution requirements for the LES or DES studies. The validation studies in most of the references are less than adequate, since most of the studies avoid to validate the turbulent characteristics of the flow. This points to a choice of academic test cases with proper Direct Numerical Simulation (DNS) or experimental validation data. In order to circumvent the complexities related to the geometry and other requirements for the boundary conditions of the flow simulation inside hydraulic turbine draft-tube, two academic test cases including internal non-swirling and swirling flows with available DNS and experimental data were chosen. First, a case of channel flow which has a simple geometry and includes complex near wall turbulent behavior was selected. Second, a sudden-expansion test case which resembles the swirling flow at part load operating condition of a hydraulic turbine draft-tube was investigated. In this paper, LES results are presented and analyzed for both test cases.

2. Simulation and Modeling of Turbulent Flow

In LES, the velocity field ( U i ) is decomposed into a filtered (or resolved) component ( U i ¯ ) and a residual (or sub-grid scale, SGS) component ( u ). By formally applying the filtering operation to the continuity equation and Navier–Stokes equations, it is possible to derive conservation laws for the filtered flow variables. Due to the linearity of the continuity equation, applying the filtering is straightforward. The form of the equation remains unchanged.
U i ¯ x i = 0 .
The resulting Navier–Stokes equations are of the same form as the original Navier–Stokes equations and also a residual stress tensor or SGS stress tensor ( τ i j r ) appears:
U j ¯ t + U i ¯ U j ¯ x i = ν 2 U j ¯ x i x i 1 ρ p ¯ x j τ i j r x i .
The equation of motions must be closed by modeling the SGS stress tensor. This is usually performed using an eddy viscosity model. A one-equation model and the WALE (Wall-Adapting Local Eddy-viscosity) model are used in the present research to model eddy-viscosity. In the WALE model, the eddy-viscosity reads as:
ν r = ( C w Δ ) 2 ( S i j d S i j d ) 3 / 2 ( S ¯ i j S ¯ i j ) 5 / 2 + ( S i j d S i j d ) 5 / 4 ,
where
S i j d = 1 2 ( g ¯ i j 2 + g ¯ j i 2 ) 1 3 δ i j g ¯ k k 2 g ¯ i j = U i ¯ x j g ¯ i j 2 = g ¯ i k g ¯ k j ,
and the constant is in the range 0.55 C w 0.60 . A one-equation model can be used to model SGS turbulent kinetic energy:
k s g s t + ( U ¯ j k s g s ) x j = x j ( ν + ν s g s ) k s g s x j + P k s g s ϵ ,
ν s g s = c k Δ k s g s 1 / 2 ,
ϵ = C ϵ k s g s 3 / 2 Δ ,
P k s g s = 2 ν s g s S ¯ i j S ¯ i j .
Very close to the wall, the van Driest damping function is used to correct the behavior for the ν s g s . This function has the following form:
D = 1 e y + A + ,
where A + = 26 . The final length scale is given by:
Δ = m i n ( κ y C S D , Δ g ) ,
where Δ g is a geometric-based filter length, such as the element cube-root volume delta.
Several approaches are presented in the literature to evaluate the resolution of LES. Following Gant (2010) [18] and Celik (2006) [19], some of these methods can be used by a prior RANS simulation or one LES calculation. For example, a good practice is to obtain the integral length scale L i n t from a prior RANS simulation as L i n t = k 3 / 2 ϵ and, therefore, to estimate the initial grid resolution for LES [20]. Therefore, they can also be called single-grid estimators, which means that, by running one LES calculation, they can be estimated. The multi-grid estimator methods require a number of LES calculations and some sort of Richardson Extrapolation. In this research, two-point autocorrelation is used for the evaluation of mesh resolutions in LES.
Correlation means the tendency of two values to change together, either in the same or opposite way. If we think of u ( x ) as a component of the velocity along a line in statistically homogeneous turbulence, the autocovariance or two-point correlation in space is:
R ( r ) u ( x + r ) u ( x ) ,
where r is the separation distance between two points. As the two points move closer to each other, R increases. If the points move further away, R will go to zero. R ( r ) is often normalized by the root-mean-square of velocity. Integral length scale L i n t can be computed from the autocovariance which is the integral of R ( r ) over the separation distance r, i.e.,
L i n t = 1 u 2 0 R ( r ) d r .
The two point correlations are calculated by Davidson (2009, 2011) [20,21]. In this study, the separation distance is varied with the grid resolution. The number of cells needed for the autocovariance to approach zero is verified, thereby checking the grid resolution. It is recommended that in a good LES, the integral length scale should cover 8–16 cells. In this study, it is reported that the energy spectra and the ratio of the resolved turbulent kinetic energy to the total one is not a good measure of LES resolution. Two-point correlations are suggested as the best measures for estimating LES resolution.
The ways to create turbulent inlet boundary conditions vary and include the two main categories of synthesis inlets and precursor simulation methods which are briefly described in the following sections.
Divergence-free synthetic eddy method (DFSEM) was proposed by Poletto et al. (2013) [22] to reduce near-inlet pressure fluctuations and the development length in the original SEM method by Reference [23]. The steps in DFSEM are the following:
  • User selection of inlet surface Ω .
  • User definition of average velocity u ( x ) , Reynolds stresses and turbulence length-scales σ ( x ) , for x Ω .
  • Eddy bounding box taken as: max { x + σ } , min { x σ } for x Ω .
  • Definition of the number of eddies.
  • Assigning random positions x k and intensities α k to all the eddies.
  • Eddies being convected through the eddy box, x k = x k + U b * Δ t , where U b = Ω u d s / Ω d s is the bulk velocity.
  • u ( x ) calculated and superimposed to u to generate the inlet condition.
  • Repeat steps 6–7 for all the subsequent time steps.
The fluctuating velocity field is:
u ( x ) = 1 N k = 1 N q σ ( | r k | ) | r k | × α k ,
where:
  • N: the number of eddies introduced into the DFSEM;
  • r k = x x k σ k with σ k being the eddy length scale for the k t h eddy;
  • q σ ( | r | k ) is a suitable shape function; and
  • α i k are random numbers with zero averages which represent eddy intensity.
The fluctuating velocity field, taking into account the turbulence anisotropy, is:
u β ( x ) = 1 N k = 1 N σ β k 1 ( d k ) 2 ϵ β j l r j k α l k ,
where:
  • d k = ( r j k ) 2 ; and
  • ϵ β j l is the Levi-Civita symbol.
The expression for the Reynolds stresses is also expressed as follows:
u β u γ = 1 N k = 1 N σ β k σ γ k ϵ β j l ϵ γ m n 1 ( d k ) 2 2 r j k r m k ( α l k ) ( α n k ) .
Therefore, σ β k and σ γ k are the eddy length scales for the k t h eddy in β and γ directions.
This method in comparison with original SEM methods is able to provide a divergence-free velocity field and also to reproduce all possible state of Reynolds stress anisotropy as a function of the characteristic ellipsoid eddy shapes described by the aspect ratio Γ = σ x σ y = σ x σ z .

3. Verification and Validation Results

3.1. Channel Flow Test Case

The turbulent flow in the channel is simulated at a moderate Reynolds number. The geometry of the channel is shown in Figure 1. The flow in the channel is considered statistically homogeneous in the spanwise direction (z) and statistically developing in the streamwise (x) and channel-height (y) directions. Domain size for the channel is (20 π , 2, π ) [m]. The physical modeling is done based on friction Reynolds number which is R e u τ = u τ δ ν = 395 , where δ = 1 [m] is a characteristic length (Channel half-height) and friction velocity, u τ = τ w ρ is assumed to be equal to 1 [ m s ]. This research took advantage of OpenFOAM built-in mesh generators. Capabilities and meshing strategies of OpenFOAM are presented by Concli et al. (2020) [24]. The initial number of nodes for the mesh was chosen as (500, 46, 82) with a total of 1,866,000 cells. Therefore, the wall units which are the normalized distance of the first mesh vertex next to the wall are ( x + , y + , z + ) = ( 49.36 , 1.1 , 15.1 ) . This orthogonal mesh results in a maximum aspect ratio of 23. Two sets of boundary conditions are considered for the generation of the inflow turbulence, periodic and Divergence Free Synthetic Eddy Method (DFSEM). The boundaries in the spanwise direction are considered periodic. k-equation is used as the SGS model.
In the initial verification study, we want to investigate the influence of the inlet eddy sizes on the development of the flow inside the channel, and we want to verify that we can obtain the mean velocity and Reynolds stress field independent of the size of eddies generated at the inlet and transported into the domain. Therefore, the effect of the size of the inlet eddies will be evaluated based on the development length of the mean streamwise velocity profile inside the channel. Our hypothesis is that further increment in the size of the inlet eddies will not further exhibit changes in the flow development length. Starting the simulation by considering 1 cell per eddy simulation, and then increasing to 3 and 5, we expect that simulating each eddy with an adequate number of 3 cells would allow capturing non-decaying turbulence fluctuations in the domain.
The results of this verification study are provided in Figure 2 and Figure 3. In these figures, the downstream development of the profiles of the mean streamwise velocity and u u component of the Reynolds stress is presented, respectively, for several cross sections and for different values of the parameter nCellsPerEddy representing the eddy sizes at the inlet. Moreover, the results for the periodic boundary condition case is also provided in the same figure, as a reference for the comparison.
The results of these simulations show enormous changes in the profiles of the mean streamwise velocity and u u component of the Reynolds stress, which are shown in Figure 2 and Figure 3. In these figures, we can see how the flow develops across the channel using different conditions for the inlet eddies. Changing the value of the number of mesh cells per eddy to a value of 3, which means that each eddy is represented by at least 3 mesh cells in the domain, although requiring smaller time steps, improves the results of the development of the flow in the channel. The influence of further increasing this value to 5 does not exhibit a further improvement in the flow characteristics and only imposes more expensive computations in terms of time steps. In the aforementioned figures, the results at each streamwise section is averaged in the spanwise direction and time; therefore, the notation and { Y Z } index in the axis titles imply averaging in time and in the { Y Z } plane, respectively.
In Figure 3, the cross sections with the values of x / δ = 0.06 and x / δ = 62.7 represent the center of the first cell immediately after and before the inlet and outlet boundaries. It seems that a minimum number of cells in the streamwise direction and away from the boundaries is required for the attenuation of the effect of boundary conditions. Moreover, it seems that, although both boundary conditions with 3 and 5 mesh cells are representing shorter development length in comparison with the one with one mesh cell per eddy, which is evident by the results at downstream locations of x / δ of 12.5, 50.2, and 56.5 being almost on top of each other, the results of the 5 cells per eddy is worse than with 3 cells per eddy. A possible explanation could be that there should be an optimum number of cells per eddies for a given mesh, and increasing this value does not necessarily provide more accurate results. In other words, the optimum value for this parameter could be mesh-dependent. This hypothesis needs further investigation and analysis which is beyond the scope of this research.
For mesh resolution analysis, two other meshes are generated. The details of these meshes are presented in Table 1.
The results of the verification analysis for the three mesh resolutions including components of mean velocity and RMS of velocity fluctuations are presented in Figure 4. In this figure, the results are averaged in both { Y Z } and { X Y } planes, meaning both span-wise and stream-wise directions. The results show different flow behavior when moving farther from the wall in the outer layer ( y + > 50 ). The finest mesh gives results that are in very good agreement with DNS across the whole region. We can conclude that grid refinement has a positive effect which is reflected in the gradual increase in the accuracy of the obtained results. The results of the RMS of velocity fluctuations show that all three components converge towards the DNS data as the mesh gets refined. Closer to the core region of the channel ( y / δ > 0.5 ), mesh 2 shows generally better prediction for the three components of RMS of velocities.
The second verification study for the mesh resolution is done for the two-point autocorrelation of velocity fluctuations. The results of this analysis are presented in Figure 5 for y + = 5 . It should be noted that the results of autocorrelation in the x direction are averaged in a cross section plane (z-plane) and the results of autocorrelation in the z direction are averaged in the stream-wise plane (x-plane). In general, the profiles indicate that the chosen size of the domain is large enough to accommodate all the relevant turbulent structures, since the profiles reach near zero value or an asymptotic value in each computed direction. These profiles show higher sensitivity of this indicator to the mesh resolution, especially in the x-direction. The results computed with DFSEM boundary conditions, overall, show better performance than the periodic case, almost for all the cases, even for the coarsest mesh. In the integral length scale results in the Table 2, we do not observe a consistent trend in the results. These behaviors in the results show that further mesh refinement is still required to get a converged value for the integral length scales, especially in the streamwise direction, since the difference in the results for Mesh 3 and 2 is still big.
The third verification study for the mesh resolution is done based on the energy spectra. The objective of this verification analysis is to evaluate the influence of the mesh resolution on the energy spectra and to compare the spectra in the inertial subrange with the Kolmogorov spectrum. The results of this verification analysis are presented in Figure 6. These results show a much less sensitive behavior of the energy spectra to the mesh resolution; therefore, this parameter does not prove to be useful as an indicator of the LES mesh resolution.

3.2. Sudden-Expansion Test Case

The Dellenback Abrupt Expansion test case resembles the swirling flow inside the draft-tube of hydraulic turbines at part load condition. Measurements of mean and fluctuating velocities were performed in a water flow with a laser Doppler anemometer [25]. The measurements were taken at the cross-sections shown in Figure 7. The simulation are performed at Reynolds number 30,000 and swirl number 0.6. The swirl number may be physically interpreted as the ratio of axial fluxes of swirl to linear momentum divided by a characteristic radius.
S = 0 R i n V θ V z r 2 d r R i n 0 R i n V z 2 r d r | z / D i n = 2.00 .
The largest uncertainties in the experiments were reported to be about 2% in Reynolds number, 8% in swirl number, and 1% in probe volume positioning. Uncertainties in mean and RMS velocities stemming from the many possible biases and broadening errors were estimated to be about ±3% and ±10%, respectively.
The mesh is shown in Figure 8. This mesh has 1,567,944 cells. This mesh corresponds to a maximum aspect ratio of 21, maximum non-orthogonality of 28, average non-orthogonality of 5, and and maximum skewness of 0.7, which are within the acceptable range. For inlet boundary conditions, velocity profiles from the measurements are specified as axisymmetric profiles using the profile1DfixedValue boundary condition in OpenFOAM. The k ω SST model is used as the turbulence model for the RANS and WALE model is used as the SGS model for the LES. For the case of simulations at swirl number 0.6, no inlet turbulence generation method is used, as it is stated in Reference [26], and the inlet turbulence is only important when the swirl in the flow is low.
The validation study is performed to compare the results of the mean velocity field and RMS of velocity fluctuations against experimental data. These results for the axial and tangential components of the mean velocity at the aforementioned cross sections are presented and compared with experimental data in Figure 9 to Figure 10. The same comparison is also drawn for the axial and tangential components of the RMS of the velocity fluctuations in Figure 11 to Figure 12. The results of the LES WALE model show great improvement over the results of the RANS k- ω SST model. These results show a huge improvement in the mean velocity results compared to the RANS model, especially immediately after the expansion where z / D = 0.25 to z / D = 1.5 . LES also captured the RMS of velocity components with a good agreement with experimental data.
The verification analysis is performed to investigate the impact of the mesh resolution on the LES results. To this aim, two other meshes one coarser and one finer than the one for which results were presented in the previous paragraph were generated. The details about these meshes are given in Table 3. The results of these simulations are provided in Figure 13 and Figure 14 for the mean velocities for axial and tangential components and in Figure 15 and Figure 16 for RMS of velocity fluctuations, again for axial and tangential components. These results show higher sensitivity of the velocity profiles to mesh resolution for both mean and RMS values. A higher mesh resolution does not necessarily give less deviation from experimental data, due to the complex nature of the flow, especially for RMS of velocity fluctuations. However, in general, results for Mesh 3 appear to be slightly more consistent with experimental data.
The results of the two-point correlations of the components of the velocity fluctuations at the centerline of the cylinder downstream of the expansion are calculated for the three meshes. These results are presented in Figure 17. This figure shows that the complex vortical structures also represent a complex interactions which does not linearly go to zero as it changes randomly further downstream of the expansion.
The third verification analysis of the mesh resolution is performed for the results of the energy spectra of pressure fluctuations. The objective of this analysis, as for the case for the channel flow, is to investigate the effect of the mesh resolution on the convergence behavior of the energy spectra and also to compare the spectra in the inertial subrange with the Kolmogorov spectrum. Time history of the pressure fluctuations at several points at the wall is recorded and monitored during the simulations. Power spectral density of the pressure fluctuations at one point at z / D = 2 is calculated for the three meshes, and results are presented in Figure 18. This case also shows less sensitivity of the results to the mesh resolution, and this criteria can again not be used as an indicator of the mesh resolution adequacy.

4. Conclusions

This work is dedicated to developing an expertise on how to use, validate, and verify the OpenFOAM CFD code for the large-eddy simulation of the flow types similar to the ones which occur at off-design operating condition, such as part load inside hydraulic turbines. Two main test cases, a case of turbulent channel flow and a case of sudden-expansion, were considered in this study. The first one brings the classical problem of near-wall turbulence complexities, and the second one resembles the swirling flow inside a hydraulic turbine draft-tube at part load off-design operating condition. The calculated profiles of two-point autocorrelation and hence the integral length scales showed that, even for a simple case of channel flow, obtaining mesh-independent results is still challenging, and further mesh-dependency analysis is required. The verified framework for LES on the channel flow was then applied to study the flow in a sudden expansion pipe flow at a Reynolds number of 30,000 and swirl number of 0.6. As with this case, the validation analysis of the LES WALE results versus RANS k- ω SST results show a large improvement of the LES results over the RANS ones for the mean and RMS of axial and tangential velocities. These improvements were also more evident at the cross sections which were immediately after the expansion. The influence of mesh resolution was also studied and showed a higher sensitivity of the results of this case to the mesh resolution and characteristics. The energy spectrum of the pressure fluctuations also did not exhibit much sensitivity to the mesh resolution, as it was also evident in the channel flow test case.

Author Contributions

Conceptualization, Z.S.M., F.G. and A.G.; methodology, Z.S.M., F.G. and A.G.; validation, Z.S.M., F.G. and A.G.; formal analysis, Z.S.M., F.G. and A.G.; investigation, Z.S.M., F.G. and A.G.; writing—original draft preparation, Z.S.M.; writing—review and editing, Z.S.M., F.G. and A.G.; visualization, Z.S.M.; supervision, F.G. and A.G.; project administration, F.G. and A.G.; funding acquisition, F.G. and A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Natural Sciences and Engineering research Council of Canada through Discovery Grants (RGPIN-04935-2015) and a Grant from the Collaborative research and Training Experience (CREATE-481695-2016) program in Simulation-based Engineering Science (Génie Par la Simulation).

Acknowledgments

This research was supported by the Natural Sciences and Engineering research Council of Canada through Discovery Grants (RGPIN-04935-2015) and a Grant from the Collaborative research and Training Experience (CREATE-481695-2016) program in Simulation-based Engineering Science (Génie Par la Simulation).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gentner, C.; Sallaberger, M.; Widmer, C.; Bobach, B.J.; Jaberg, H.; Schiffer, J.; Senn, F.; Guggenberger, M. Comprehensive experimental and numerical analysis of instability phenomena in pump turbines. IOP Conf. Ser. Earth Environ. Sci. 2014, 22, 032046. [Google Scholar] [CrossRef] [Green Version]
  2. Nicolle, J.; Giroux, A.M.; Morissette, J.F. CFD configurations for hydraulic turbine startup. IOP Conf. Ser. Earth Environ. Sci. 2014, 22, 032021. [Google Scholar] [CrossRef]
  3. Hosseinimanesh, H.; Devals, C.; Nennemann, B.; Reggio, M.; Guibault, F. A Numerical Study of Francis Turbine Operation at No-Load Condition. J. Fluids Eng. 2016, 139, 011104. [Google Scholar] [CrossRef]
  4. Foroutan, H.; Yavuzkurt, S. Simulation of flow in a simplified draft tube: Turbulence closure considerations. IOP Conf. Ser. Earth Environ. Sci. 2012, 15, 022020. [Google Scholar] [CrossRef]
  5. Krappel, T.; Riedelbauch, S.; Jester-Zuerker, R.; Jung, A.; Flurl, B.; Unger, F.; Galpin, P. Turbulence Resolving Flow Simulations of a Francis Turbine in Part Load using Highly Parallel CFD Simulations. IOP Conf. Ser. Earth Environ. Sci. 2016, 49, 062014. [Google Scholar] [CrossRef] [Green Version]
  6. Pasche, S.; Avellan, F.; Gallaire, F. Part Load Vortex Rope as a Global Unstable Mode. J. Fluids Eng. 2017, 139, 051102. [Google Scholar] [CrossRef]
  7. Neto, A.D.A.; Jester-Zuerker, R.; Jung, A.; Maiwald, M. Evaluation of a Francis turbine draft tube flow at part load using hybrid RANS-LES turbulence modelling. IOP Conf. Ser. Earth Environ. Sci. 2012, 15, 062010. [Google Scholar] [CrossRef]
  8. Sentyabov, A.V.; Gavrilov, A.A.; Dekterev, A.A.; Minakov, A.V. Numerical investigation of the vortex core precession in a model hydro turbine with the aid of hybrid methods for computation of turbulent flows. Thermophys. Aeromech. 2014, 21, 707–718. [Google Scholar] [CrossRef]
  9. Paik, J.; Sotiropoulos, F.; Sale, M.J. Numerical Simulation of Swirling Flow in Complex Hydroturbine Draft Tube Using Unsteady Statistical Turbulence Models. J. Hydraul. Eng. 2005, 131, 441–456. [Google Scholar] [CrossRef]
  10. Minakov, A.; Platonov, D.; Dekterev, A.; Sentyabov, A.; Zakharov, A. The analysis of unsteady flow structure and low frequency pressure pulsations in the high-head Francis turbines. Int. J. Heat Fluid Flow 2015, 53, 183–194. [Google Scholar] [CrossRef]
  11. Javadi, A.; Nilsson, H. A comparative study of scale-adaptive and large-eddy simulations of highly swirling turbulent flow through an abrupt expansion. IOP Conf. Ser. Earth Environ. Sci. 2014, 22, 022017. [Google Scholar] [CrossRef] [Green Version]
  12. Pacot, O.; Kato, C.; Avellan, F. High-resolution LES of the rotating stall in a reduced scale model pump-turbine. IOP Conf. Ser. Earth Environ. Sci. 2014, 22, 022018. [Google Scholar] [CrossRef] [Green Version]
  13. Morgut, M.; Jošt, D.; Nobile, E.; Škerlavaj, A. Numerical investigation of the flow in axial water turbines and marine propellers with scale-resolving simulations. J. Phys. Conf. Ser. 2015, 655, 012052. [Google Scholar] [CrossRef]
  14. Yang, J.; Zhou, L.; Wang, Z. The numerical simulation of draft tube cavitation in Francis turbine at off-design conditions. Eng. Comput. 2016, 33, 139–155. [Google Scholar] [CrossRef]
  15. Minakov, A.; Platonov, D.; Sentyabov, A.; Gavrilov, A. Francis-99 turbine numerical flow simulation of steady state operation using RANS and RANS/LES turbulence model. J. Phys. Conf. Ser. 2017, 782, 012005. [Google Scholar] [CrossRef] [Green Version]
  16. Trivedi, C. Compressible Large Eddy Simulation of a Francis Turbine During Speed-No-Load: Rotor Stator Interaction and Inception of a Vortical Flow. J. Eng. Gas Turbines Power 2018, 140, 112601. [Google Scholar] [CrossRef]
  17. Trivedi, C.; Dahlhaug, O.G. Interaction between trailing edge wake and vortex rings in a Francis turbine at runaway condition: Compressible large eddy simulation. Phys. Fluids 2018, 30, 075101. [Google Scholar] [CrossRef]
  18. Gant, S.E. Reliability Issues of LES-Related Approaches in an Industrial Context. Flow Turbul. Combust. 2010, 84, 325–335. [Google Scholar] [CrossRef]
  19. Celik, I.; Klein, M.; Freitag, M.; Janicka, J. Assessment measures for URANS/DES/LES: An overview with applications. J. Turbul. 2006, 7, N48. [Google Scholar] [CrossRef]
  20. Davidson, L. Large Eddy Simulations: How to evaluate resolution. Int. J. Heat Fluid Flow 2009, 30, 1016–1025. [Google Scholar] [CrossRef] [Green Version]
  21. Salvetti, M.V. (Ed.) Quality and Reliability of Large-Eddy Simulations II; Number 16 in ERCOFTAC Series; OCLC: 837896204; Springer: Dordrecht, The Netherlands, 2011. [Google Scholar]
  22. Poletto, R.; Craft, T.; Revell, A. A New Divergence Free Synthetic Eddy Method for the Reproduction of Inlet Flow Conditions for LES. Flow Turbul. Combust. 2013, 91, 519–539. [Google Scholar] [CrossRef]
  23. Jarrin, N.; Prosser, R.; Uribe, J.C.; Benhamadouche, S.; Laurence, D. Reconstruction of turbulent fluctuations for hybrid RANS/LES simulations using a Synthetic-Eddy Method. Int. J. Heat Fluid Flow 2009, 30, 435–442. [Google Scholar] [CrossRef] [Green Version]
  24. Concli, F.; Schaefer, C.T.; Bohnert, C. Innovative Meshing Strategies for Bearing Lubrication Simulations. Lubricants 2020, 8, 46. [Google Scholar] [CrossRef] [Green Version]
  25. Dellenback, P.A.; Metzger, D.E.; Neitzel, G.P. Measurements in turbulent swirling flow through an abrupt axisymmetric expansion. AIAA J. 1988, 26, 669–681. [Google Scholar] [CrossRef]
  26. Schlüter, J.U.; Pitsch, H.; Moin, P. Large-Eddy Simulation Inflow Conditions for Coupling with Reynolds-Averaged Flow Solvers. AIAA J. 2004, 42, 478–484. [Google Scholar] [CrossRef]
Figure 1. The geometry of the channel.
Figure 1. The geometry of the channel.
Fluids 06 00024 g001
Figure 2. Mean velocity profiles at selected streamwise locations using various inlet conditionsas (a) 1 cell Per Eddy (b) 3 cells Per Eddy (c) 5 cells Per Eddy (d) Periodic.
Figure 2. Mean velocity profiles at selected streamwise locations using various inlet conditionsas (a) 1 cell Per Eddy (b) 3 cells Per Eddy (c) 5 cells Per Eddy (d) Periodic.
Fluids 06 00024 g002
Figure 3. u u profiles at selected streamwise locations using various inlet conditions as (a) 1 cell Per Eddy (b) 3 cells Per Eddy (c) 5 cells Per Eddy (d) Periodic.
Figure 3. u u profiles at selected streamwise locations using various inlet conditions as (a) 1 cell Per Eddy (b) 3 cells Per Eddy (c) 5 cells Per Eddy (d) Periodic.
Fluids 06 00024 g003
Figure 4. Comparison of RMS of velocity fluctuations for different mesh resolutions for (a) stream-wise, (b) wall-normal, (c) span-wise and (d) shear stress.
Figure 4. Comparison of RMS of velocity fluctuations for different mesh resolutions for (a) stream-wise, (b) wall-normal, (c) span-wise and (d) shear stress.
Fluids 06 00024 g004
Figure 5. Spatial autocorrelation of velocity components at y + = 5 : stream-wise velocity along (a) x and (b) z directions, wall-normal velocity along (c) x and (d) z directions, span-wise velocity along (e) x and (f) z directions.
Figure 5. Spatial autocorrelation of velocity components at y + = 5 : stream-wise velocity along (a) x and (b) z directions, wall-normal velocity along (c) x and (d) z directions, span-wise velocity along (e) x and (f) z directions.
Fluids 06 00024 g005
Figure 6. Energy spectra at two y + values of (a) 15 and (b) 150.
Figure 6. Energy spectra at two y + values of (a) 15 and (b) 150.
Fluids 06 00024 g006
Figure 7. Measurement cross-sections. Numbers refer to Z/D, where D is the inlet diameter, and Z = 0 at the abrupt expansion.
Figure 7. Measurement cross-sections. Numbers refer to Z/D, where D is the inlet diameter, and Z = 0 at the abrupt expansion.
Fluids 06 00024 g007
Figure 8. Numerical mesh for abrupt expansion test case.
Figure 8. Numerical mesh for abrupt expansion test case.
Fluids 06 00024 g008
Figure 9. Mean axial velocity at several cross sections (dashed line: k- ω Shear Stress Transport (SST) Reynolds-averaged Navier–Stokes (RANS), solid line: Wall-Adapting Local Eddy-viscosity (WALE) large-eddy simulation (LES), dots: experiments).
Figure 9. Mean axial velocity at several cross sections (dashed line: k- ω Shear Stress Transport (SST) Reynolds-averaged Navier–Stokes (RANS), solid line: Wall-Adapting Local Eddy-viscosity (WALE) large-eddy simulation (LES), dots: experiments).
Fluids 06 00024 g009
Figure 10. Mean tangential velocity at several cross sections (dashed line: k- ω SST RANS, solid line: WALE LES, dots: experiments).
Figure 10. Mean tangential velocity at several cross sections (dashed line: k- ω SST RANS, solid line: WALE LES, dots: experiments).
Fluids 06 00024 g010
Figure 11. RMS of axial velocity at several cross sections (solid line: WALE LES, dots: experiments).
Figure 11. RMS of axial velocity at several cross sections (solid line: WALE LES, dots: experiments).
Fluids 06 00024 g011
Figure 12. RMS of tangential velocity at several cross sections (solid line: WALE LES, dots: experiments).
Figure 12. RMS of tangential velocity at several cross sections (solid line: WALE LES, dots: experiments).
Fluids 06 00024 g012
Figure 13. Mean axial velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Figure 13. Mean axial velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Fluids 06 00024 g013
Figure 14. Mean tangential velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Figure 14. Mean tangential velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Fluids 06 00024 g014
Figure 15. RMS of axial velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Figure 15. RMS of axial velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Fluids 06 00024 g015
Figure 16. RMS of tangential velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Figure 16. RMS of tangential velocity at several cross sections for different mesh resolutions (red: Mesh 1, green: Mesh 2, blue: Mesh 3, dots: experiments).
Fluids 06 00024 g016
Figure 17. Two-point auto-correlation of the (a) stream-wise, (b) wall-normal, (c) span-wise components of the velocity fluctuations along the centerline.
Figure 17. Two-point auto-correlation of the (a) stream-wise, (b) wall-normal, (c) span-wise components of the velocity fluctuations along the centerline.
Fluids 06 00024 g017
Figure 18. Power spectral density of the pressure fluctuations at one point on the wall at z / D = 2 .
Figure 18. Power spectral density of the pressure fluctuations at one point on the wall at z / D = 2 .
Fluids 06 00024 g018
Table 1. Details of the numerical grids.
Table 1. Details of the numerical grids.
NameNumber of CellsTotal Size Δ x + Δ z + y +
Mesh 1 397 × 38 × 65 980,59062.5119.091.34
Mesh 2 500 × 46 × 82 1,886,00049.6315.131.11
Mesh 3 630 × 58 × 103 3,763,62039.3912.050.91
Table 2. Integral length scales.
Table 2. Integral length scales.
y + L x L z
Mesh: 123PeriodicDNSMesh: 123PeriodicDNS
u u 50.740.750.510.900.530.020.010.010.020.02
u u 1500.560.380.300.500.800.040.040.030.050.08
v v 50.740.750.510.900.130.020.010.010.020.01
v v 1500.560.380.300.500.150.040.040.030.050.07
w w 50.740.750.510.900.200.020.010.010.020.02
w w 1500.560.380.300.500.080.040.040.030.050.16
Table 3. Details of the numerical grids for sudden expansion test case.
Table 3. Details of the numerical grids for sudden expansion test case.
NameNumber of Cells Δ t y min + y max + y average +
Mesh 1787,5120.00050.19535.425.29
Mesh 21,567,9440.00050.36028.704.21
Mesh 33,108,3970.00010.41254.176.38
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Seifollahi Moghadam, Z.; Guibault, F.; Garon, A. On the Evaluation of Mesh Resolution for Large-Eddy Simulation of Internal Flows Using Openfoam. Fluids 2021, 6, 24. https://doi.org/10.3390/fluids6010024

AMA Style

Seifollahi Moghadam Z, Guibault F, Garon A. On the Evaluation of Mesh Resolution for Large-Eddy Simulation of Internal Flows Using Openfoam. Fluids. 2021; 6(1):24. https://doi.org/10.3390/fluids6010024

Chicago/Turabian Style

Seifollahi Moghadam, Zahra, François Guibault, and André Garon. 2021. "On the Evaluation of Mesh Resolution for Large-Eddy Simulation of Internal Flows Using Openfoam" Fluids 6, no. 1: 24. https://doi.org/10.3390/fluids6010024

APA Style

Seifollahi Moghadam, Z., Guibault, F., & Garon, A. (2021). On the Evaluation of Mesh Resolution for Large-Eddy Simulation of Internal Flows Using Openfoam. Fluids, 6(1), 24. https://doi.org/10.3390/fluids6010024

Article Metrics

Back to TopTop