Next Article in Journal
Mixing of Particles in a Rotating Drum with Inclined Axis of Rotation
Next Article in Special Issue
Numerical Study on the Aerodynamic Characteristics of the NACA 0018 Airfoil at Low Reynolds Number for Darrieus Wind Turbines Using the Transition SST Model
Previous Article in Journal
Subspace Based Model Identification for an Industrial Bioreactor: Handling Infrequent Sampling Using Missing Data Algorithms
Previous Article in Special Issue
Numerical Analysis of the Flow around Two Square Cylinders in a Tandem Arrangement with Different Spacing Ratios Based on POD and DMD Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Study on CFD Turbulence Models for the Flow Field in Air Cooled Radiator

Changchun Institute of Optics, Fine Mechanics and Physics, Chinese Academy of Sciences, Changchun 130033, China
*
Author to whom correspondence should be addressed.
Processes 2020, 8(12), 1687; https://doi.org/10.3390/pr8121687
Submission received: 3 December 2020 / Revised: 17 December 2020 / Accepted: 18 December 2020 / Published: 21 December 2020
(This article belongs to the Special Issue Advancement in Computational Fluid Mechanics and Optimization Methods)

Abstract

:
This paper compares the performances of three Computational Fluid Dynamics (CFD) turbulence models, Reynolds Average Navier-Stokes (RANS), Detached Eddy Simulation (DES), and Large Eddy Simulation (LES), for simulating the flow field of a wheel loader engine compartment. The distributions of pressure fields, velocity fields, and vortex structures in a hybrid-grided engine compartment model are analyzed. The result reveals that the LES and DES can capture the detachment and breakage of the trailing edge more abundantly and meticulously than RANS. Additionally, by comparing the relevant calculation time, the feasibility of the DES model is proved to simulate the three-dimensional unsteady flow of engine compartment efficiently and accurately. This paper aims to provide a guiding idea for simulating the transient flow field in the engine compartment, which could serve as a theoretical basis for optimizing and improving the layout of the components of the engine compartment.

1. Introduction

The turbulent structure is very important to the flow field distribution of the power compartment of the loader. Its transportation velocity and distribution have important research values. The small-scale turbulent structure is transported at the speed of local flow field, while the large-scale turbulent structure maintains a fixed proportion between the transport speed and the average flow velocity [1].
Bull et al. [2] conducted experiments on subsonic boundary layers, and the results showed that the proportion between the transport velocity and the average flow velocity of the small-scale turbulent structure was 0.53, while that of large-scale turbulent structure was 0.825. Favre et al. [3,4,5] have extensively studied the spatio-temporal correlation characteristics of boundary layer turbulent fields. The results show that the high and intermediate frequency pulsations contain about 50% turbulent kinetic energy, which is mainly characterized by high viscous dissipation and strong space-time attenuation, while the low frequency pulsations contain another 50% turbulent kinetic energy, of which only 1% is lost due to viscosity. Demetriades [6] conducted an experiment on the compressible supersonic boundary layer, which shows that the spatial correlation between flow direction and radial direction is basically the same, and only after the distance of vortex structure migration exceeds its own scale does the correlation significantly decrease. Ganapathisubramani et al. [7] measured the boundary layer with particle image velocimetry (PIV), and the results showed that the flow direction and the spreading scale in the boundary layer increased with the increase of the distance to the wall, and the change rate was basically linear. In the external flow field, the trend of scale variation is the opposite. The lower speed zone contributes more to the attenuation of flow correlation Ruu. Kline et al. [8] observed that there were a large number of low-velocity strip structures in the boundary layer through the H2 bubble experiment. Falco et al. [9] studied the quasisturon motion in the turbulent boundary layer based on the smoke spectrum, determined the motion scale, and observed a large number of vortex structures. Moin and Kim [10] carried out a simulation study on turbulent channel flow based on large eddy simulation. Thereafter, numerical studies based on DNS were carried out in large quantities [11,12,13].
Different from traditional experimental method, numerical technology for the data of the turbulent flow field are measured and stored. The data processing means a great deal of diversity, and it can measure all direct physical quantities, such as pressure, velocity, and deals with predefined indirect quantities. For example, in the previous experiment, single point measurement was carried out based on thermocouple anemometer, so quadrant analysis could only be limited to a few measurement points, and it was concluded that up-jet and down-sweep were upstream and downstream relationships [14]. Numerical studies will not be subject to such restrictions, and not only can we obtain the three-dimensional structure of the quadrant component of specific instantaneous shear product, we can also determine the vortex structure [15,16,17]. Compared with the experiment, the biggest advantage of numerical study lies in the traceability of specific experimental phenomena, that is, it can trace forward and deduce the quasi-structure and quasi-motion at a certain time point. At the same time, it can also be studied by the means of defining phenomena and using conditional sampling [18]. Calautit J. K. et al. [19] used virtual wind tunnels to replace traditional wind tunnel tests to, respectively, conduct flow field and transmission of the power cabin of passenger cars and loaders. Lu P. et al. [20] adopted the method of one-dimensional and three-dimensional co-simulation to analyze the flow characteristics of the power cabin and the change rules of the internal thermal environment.
However, different turbulence models are used to predict uniformity, which results in different results. Therefore, hybrid meshing inside the engine compartment is used in this research, and a lot of energy and time is spent to improve the grid quality. By comparing it with the experimental results, the prediction performance of different turbulence models on the flow field and temperature field of the engine compartment is evaluated. The temperature field, pressure streamline structure, vorticity, and velocity distribution in the engine compartment are analyzed. Compared with the Reynolds-averaged Navier-Stokes equations (RANS) model, the Large Eddy Simulation (LES) and Detached Eddy Simulation (DES) models can capture transient vorticity characteristics such as the generation, development, and formation of vortices, as well as the trailing edge shedding and fracture. In addition, the DES model has been proven to be a feasible method to accurately and efficiently simulate three-dimensional unsteady flow in complex cabin space.

2. Numerical Simulation

2.1. Turbulence Models

The average time domain of the Navier-Stokes equation to obtain the continuity and momentum equations for incompressible flow is as follows:
ρ t + ( ρ u i ) x i = 0 ( ρ u i ) t + ( ρ u i u j ) x j = p x i + x j ( μ u i x j ) + τ i j x j
where ρ denotes the density and u i denotes velocity. p is the pressure, and μ is the coefficient of viscosity.
Where the Reynolds stress τ i j is:
τ i j = ρ u i u j ¯ .
The most widely used turbulence model is the standard k-ε model. The model is based on the kinetic energy k equation, and then an equation of turbulent dissipation rate ε is introduced. The ε is defined as:
k = 1 2 m v 2 ,
ε = μ ρ u l x k u l x k ¯ .
The turbulent viscosity μ t can be expressed as a function of k and ε:
μ t = ρ C μ k 2 ε .
Therefore, the transport equation of the standard k-ε model is [21]:
ρ k t + ρ k u i x i = x j μ + μ t σ k k x j + G k + G b ρ ε Y M + S k ρ ε t + ρ k u i x i = x j μ + μ t σ ε ε x j + C 1 ε ε k G k + C 3 ε G b C 2 ε ρ ε 2 k + S ε
where the kinetic energy caused by the average velocity gradient is:
G k = μ t u i x j + u i x i u i x j .
The turbulent energy generated by the buoyancy effect G b is [22]:
G b = β g i u t Pr t T x i .
The following values are used for standard k-ε model: C1ε = 1.44, C2ε = 1.92, C3ε = 0.09, σk = 1.0, and σε = 1.3.
The basic idea of RANS is to perform time averaging on the N-S equations to transform the unsteady turbulence problem into a steady problem study, at the cost of additional unknown numbers, which are also in the same form as stress, called Reynolds stress. Reynolds stress also needs to be characterized by a model, which is the so-called turbulence model. However, due to the time averaging of the problem, the information contained in the equation itself has been partially lost. It is actually very difficult to give a Reynolds stress model and it is also difficult to apply it to all flows.
LES is a commonly used turbulence simulation method. LES is based on the self-similarity theory and uses a subgrid-scale to simulate small vortices, while the large vortices are based on the geometric calculation. The disadvantage of this model is that it is difficult to calculate the region near the wall. In the LES simulation method, the transport of fluid momentum, mass, energy, and other physical quantities are mainly affected by large-scale vortices. Small-scale vortices with isotropic motion are not affected by geometric and boundary conditions. Therefore, LES uses the N-S equation to calculate the turbulent motion larger than the grid-scale. The LES model obtains a filtered momentum equation by filtering out vortices smaller than the filtered grid in the Fourier equation or spatial domain N-S equation. The filtered N-S equation and the continuity equation are transformed into [22]:
ρ t + ( ρ u i ¯ ) x i = 0 t ( ρ u i ¯ ) + x j ( ρ u i u j ¯ ) = p ¯ x i + x j μ u i ¯ x j τ i j x j .
The subgrid-scale stress tensor that results from filtering the Navier-Stokes equations consists of three terms:
τ i j = L i j + R i j + C i j , L i j = u i ¯ u ¯ j ¯ u ¯ i u ¯ j , C i j = u i ¯ u j + u i u j ¯ ¯ , R i j = u i u j ¯
where u ¯ i denotes the filtered velocity component and u i = u i u i ¯ denotes the SGS component of u i . L i j is the Leonard term, C i j is the cross term, and R i j is the SGS Reynolds stress.
The sub-grid stress term τ i j is obtained from the vortex viscous model equation:
τ i j 1 3 τ k k δ i j = 2 μ t S i j ¯
where μ t is the turbulent dynamic viscosity, and S i j ¯ is the rate-of-strain tensor, which is defined as:
S i j ¯ = 1 2 ( u i ¯ x j + u j ¯ x i ) .
Because the turbulence scale decreases rapidly with the increase of Reynolds number, LES simulation cannot be used in engineering calculation at a high Reynolds number. According to the characteristics of the RANS/LES model, a kind of large-scale separated flow in engineering is reasonably arranged by “model allocation.” Since the LES model is not an ideal model for solving near-wall problems, the DES turbulence model was developed. The DES turbulence model combines RANS with LES methods. The basic idea of DES is to use the RANS model within the boundary layer and the LES side in the separation area, which saves computational resources and ensures computational accuracy. The DES model selected SST (Shear Stress Transfer) as the RANS model in this study, and its flow energy diffusion term Yk is modified to [23]:
Y k = ρ β k w F D E S F D E S = max ( L t C D E S Δ max , 1 )
where CDES is a calibration constant in the DES model with a value of 0.61, and Δ max is the largest grid gap in the three directions ( Δ x ,   Δ y ,   Δ z ) in the grid.

2.2. Geometric Model Construction and Mesh Generation

In this paper, different computational fluid dynamics (CFD) turbulence models are used to simulate the flow field of a loader’s engine compartment. In order to facilitate the calculation, the model is appropriately simplified, and the components such as the bucket, front axle, and wheels that have less influence on the flow field of the engine compartment are omitted. A 3D simulation model of the engine compartment is established to analyze its internal flow field. The simplified geometric model consists of an outer wind tunnel, an inner wind tunnel, and a loader model as shown in Figure 1 and Figure 2.
The GAMBIT software is used to perform hybrid meshing on the radiator group. Due to the irregular structure of the fan rotation domain and engine compartment, the mesh was built by tetrahedron and hexahedral. The grid of the air intake, inlet fan, and air duct area of the engine compartment are refinement. In order to reduce the total mesh size of the model, a two-layer wind tunnel nested structure is used. The size-Function and Hex-core functions are used to reduce the total mesh number in the outer and inner wind tunnels. The grid distribution in the engine compartment is dense. The outer wind tunnel and the inner wind tunnel mesh are sparse, and the virtual wind tunnel model mesh is shown in Figure 3. The minimum grid size of the simulation model is 1 mm. The total number of grids in the entire calculation model is 800,000.

2.3. Boundary Condition Setting

The incompressible N-S equation is solved by ANSYS Fluent software. During the calculation, the fan rotates along the mesh section at a certain time step (0.01 s). When the change of mass flow rate in two consecutive circulatory iterations and all normalized residuals in mass and momentum conservation equations were less than 10−4 for the differences, the solution was assumed to be converged [24,25,26,27]. The fan speed is the same as the experimental value. The mesh interface is set between the various modules of the engine compartment. The boundary conditions and settings of the CFD solver are shown in Table 1. The radiators in the engine compartment are simulated by a porous medium model [28].

3. Results and Analysis

3.1. Experiment Analysis

The credibility of the computational model is verified by experiments. There are 9 monitoring points set at the exit of the engine compartment. When the fan speed reaches 1000 rpm, they collect the wind speed value. The monitoring point position is shown in Figure 4. Figure 5 shows the comparison of the three simulation results at 9 monitoring points with the experimental data collected at 1000 rpm. Compared with the experimental results, the error of the calculation results of each simulation model is less than 15%, and the simulation model can solve the actual fluid flow of the loader’s engine compartment relatively correctly.

3.2. Pressure Distribution

The arrangement of various heat dissipating components in the loader plays an important role in the pressure distribution of the engine compartment. Figure 6 shows the evolution of pressure in the engine compartment with fan speed. As the fan speed increases, the velocity of the fluid in the cabin gradually increases. A low-pressure region appears at the fan shaft center and a high-pressure region appears at the outer edge of the fan. At low speeds, the surface unevenness of the radiator is consistent with the dynamic pressure of the fan outlet. Due to the large space in the upper part of the engine, the fluid inflow is large at high speeds. Under the commonality of the viscous force and pressure gradient, the radiator gradually approaches the equilibrium pressure.
There is a slight difference in the pressure gradient between the RANS, LES, and DES models. Compared with the other two simulation models, the LES model has a significantly low-pressure region at the fan shaft center when fan speeds are high. The pressure trends described in the three turbulence models are the same, indicating that all simulations can effectively capture the mainstream structure.

3.3. Velocity Distribution

The evolution of the fluid flow field in the engine compartment is shown in Figure 7. All simulation models are good at predicting the magnitude and level of flow. However, after a closer look at the YZ section, it can be seen that there is a recirculation zone generated by the intake air in the LES and DES models, while it remains substantially stationary in the RANS model. Figure 8 shows the flow field distribution of the XY section. The RANS model shows multiple small turbulences between the engine and fan, while the LES and DES models are relatively smooth in this area.
According to Figure 7 and Figure 8, as the rotational speed increases, the high-speed zone on the outer edge of the fan matches the pressure distribution described above. A high flow rate region occurs at the boundary of the radiator, and the seal between the radiators of each layer allows the cold air to pass sufficiently. There is a retardation zone between adjacent radiators, and the efficient radiator has a relatively small pressure drop retardation zone. In the simulation of the engine compartment, it can be seen that the cold air mainly flows from above and below the power compartment.

3.4. Vortex Structure

The vortex was shown by a special isosurface (called vortex core). Based on the Q criterion, the Q criterion index is used to describe it. The recipe is defined as follows [29,30]:
Q = 1 2 ( Ω 2 + S 2 )
where Ω and S represent the vorticity tensor and strain rate tensor, respectively.
Figure 9 is a vortex core diagram of the rotating region. As we can see in the standard transient structure for the Q standard, the vortex structure among RANS, LES, and DES models are completely different. LES and DES can reflect more detailed vortex cores.
It can be seen from the figure that the vortex core of the fan is mainly composed of two parts; the vortex core is generated from the blade root vortex and tip vortex formed by the leading edge of the blade. The vortex sheets remain attached on the pressure side while separated at the tip. A high-velocity vortex core appears on the windward surface of the fan. As the rotational speed increases, the fan blade surface gradually forms a vortex. Due to the faster transition between DES and LES, the vortex is also present at the root of the blade, which is gradually developed more abundantly and orderly.

3.5. Prediction Performance of the Three Turbulent Models

As shown in Figure 10, the simulation data are in good agreement with the experimental data, and the wind speed at each monitoring point increases gradually. A low-speed zone occurs at the detection point 5, which is consistent with the unevenness of the fan flow rate. The prediction speed of the Standard k-ε model and LES model is relatively larger than the DES model. However, the DES model is more accurate than the Standard k-ε model and LES model, whose absolute error is under 20%.
The calculation time of the three simulation models is shown in Table 2. The RANS simulation is about 2.8 times faster than the LES simulation. In industrial design, the RANS model is used for CFD simulation. Although the calculation cost of the RANS model is low, its accuracy is lower than that of the LES model. The LES model has high precision, but the calculation time is long and it is difficult to apply to actual engineering. Besides, the DES model, which combines the advantages of LES and RANS models, is an ideal choice for CFD simulation with short computation time and high accuracy.

4. Conclusions

In this paper, the hybrid grid model is used to simulate the unsteady flow of the engine power compartment. The ability of the RANS, DES, and LES models to predict the turbulence distribution is compared and analyzed. The engine compartment propulsion simulation using RANS, DES, and LES is proposed. Scholars can use the RANS model to estimate the flow velocity and structure reasonably and accurately. However, LES or DES models should be used when it is necessary to obtain more accurate turbulence distributions with complex geometric structures. Under normal working conditions, the natural wind energy in the environment improves the heat dissipation performance of the engine compartment of the loader. In this study, considering the influence of natural wind on heat dissipation can improve the accuracy of performance prediction. The DES model combines the advantages of the RANS model and LES model, which can accurately simulate a complex eddy current field. At the same time, the LES model can simulate large-scale eddy currents in the boundary layer. The DES model has proven to be an accurate and efficient model for simulating three-dimensional unsteady turbulent flow in complex channels. Meanwhile, theoretical research on this aspect will be the focus of CFD simulation in the future.

Author Contributions

Conceptualization, C.Y.; methodology, C.Y.; software, X.X.; validation, K.S.; formal analysis, K.S.; investigation, X.X.; resources, M.S.; data curation, M.S.; writing—original draft preparation, C.Y.; writing—review and editing, C.Y.; visualization, Y.L.; supervision, X.X.; project administration, X.X.; funding acquisition, X.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

kTurbulent kinetic energy [m2/s2]εTurbulent dissipation rate [m2/s3]
NRotating speed [rpm]ωRelative velocity [m/s]
uiTime-averaged velocity [m/s]νTurbulent viscosity [m2/s]
XiCoordinate [m]ρdensity [kg/m3]
tTime [s]μCircumferential velocity [m/s]
PFluid pressure [Pa]μtTurbulence viscosity coefficient
τijReynolds stress

References

  1. Del Álamo, J.C.; Jiménez, J. Estimation of turbulent convection velocities and corrections to Taylor’s approximation. J. Fluid Mech. 2009, 640, 5–26. [Google Scholar]
  2. Bull, M. Wall-pressure fluctuations associated with subsonic turbulent boundary layer flow. J. Fluid Mech. 1967, 28, 719–754. [Google Scholar]
  3. Favre, A. Review on space-time correlations in turbulent fluids. J. Appl. Mech. 1965, 32, 241–257. [Google Scholar]
  4. Favre, A.; Gaviglio, J.; Dumas, R. Further space-time correlations of velocity in a turbulent boundary layer. J. Fluid Mech. 1958, 3, 344–356. [Google Scholar]
  5. Favre, A.; Gaviglio, J.; Dumas, R. Structure of Velocity Space-Time Correlations in a Boundary Layer. Phys. Fluids 1967, 10, S138–S145. [Google Scholar]
  6. Demetriades, A. Turbulence correlations in a compressible wake. J. Fluid Mech. 1976, 74, 251–267. [Google Scholar]
  7. Ganapathisubramani, B.; Hutchins, N.; Hambleton, W.T.; Longmire, E.K.; Marusic, I. Investigation of large-scale coherence in a turbulent boundary layer using two-point correlations. J. Fluid Mech. 2005, 524, 57–80. [Google Scholar]
  8. Kline, S.J.; Reynolds, W.C.; Schraub, F.A.; Runstadler, P.W. The structure of turbulent boundary layers. J. Fluid Mech. 1967, 30, 741–773. [Google Scholar]
  9. Falco, R.E. Coherent motions in the outer region of turbulent boundary layers. Phys. Fluids 1977, 20, S124. [Google Scholar]
  10. Moin, P.; Kim, J. Numerical investigation of turbulent channel flow. J. Fluid Mech. 1982, 118, 341–377. [Google Scholar]
  11. Kim, J.; Moin, P.; Moser, R. Turbulence statistics in fully developed channel flow at low Reynolds number. J. Fluid Mech. 1987, 177, 133–166. [Google Scholar] [CrossRef] [Green Version]
  12. Moser, R.D.; Kim, J.; Mansour, N.N. Direct numerical simulation of turbulent channel flow up to Reτ = 590. Phys. Fluids 1999, 11, 943–945. [Google Scholar] [CrossRef]
  13. Abe, H.; Kawamura, H.; Matsuo, Y. Direct numerical simulation of a fully developed turbulent channel flow with respect to the Reynolds number dependence. J. Fluids Eng. Trans. ASME 2001, 123, 382–393. [Google Scholar] [CrossRef]
  14. Wallace, J.M.; Eckelmann, H.; Brodkey, R.S. The wall region in turbulent shear flow. J. Fluid Mech. 1972, 54, 39–48. [Google Scholar] [CrossRef]
  15. Lozano-Durán, A.; Flores, O.; Jiménez, J. The three-dimensional structure of momentum transfer in turbulent channels. J. Fluid Mech 2012, 694, 100–130. [Google Scholar] [CrossRef]
  16. Robinson, S.K. Coherent motions in the turbulent boundary layer. Annu. Rev. Fluid Mech. 1991, 23, 601–639. [Google Scholar] [CrossRef]
  17. Adrian, R.J. Hairpin vortex organization in wall turbulence. Phys. Fluids 2007, 19, 041301. [Google Scholar] [CrossRef] [Green Version]
  18. Jiménez, J. Coherent structures in wall-bounded turbulence. J. Fluid Mech. 2018, 842, 1. [Google Scholar] [CrossRef] [Green Version]
  19. Calautit, J.K.; Chaudhry, H.N.; Hughes, B.R.; Sim, L.F. A validated design methodology for a closed-loop subsonic wind tunnel. J. Wind Eng. Ind. Aerodyn. 2014, 125, 180–194. [Google Scholar] [CrossRef]
  20. Lu, P.; Gao, Q.; Wang, Y. The simulation methods based on 1D/3D collaborative computing for the vehicle integrated thermal management. Appl. Therm. Eng. 2016, 104, 42–53. [Google Scholar] [CrossRef]
  21. Galvan, S.; Reggio, M.; Guibault, F. Assessment study of k-ε turbulence models and near-wall modeling for steady state swirling flow analysis in draft tube using fluent. Eng. Appl. Comput. Fluid Mech. 2011, 5, 459–478. [Google Scholar] [CrossRef] [Green Version]
  22. Yang, X.; Gupta, S.; Kuo, T.W.; Gopalakrishnan, V. RANS and large Eddy simulation of internal combustion engine flows-a comparative study. J. Eng. Gas Turbines Power 2014, 136, 051507. [Google Scholar] [CrossRef]
  23. Xu, C.; Chen, L.; Lu, X. Large-eddy and detached-eddy simulations of the separated flow around a circu-lar cylinder. J. Hydrodyn. 2007, 19, 559–563. [Google Scholar] [CrossRef]
  24. Wood, H.G.; Yan, Q.; Wei, W. Parametric analysis and optimization of inlet deflection angle in Torque converters. J. Fluid Eng. 2015, 137, 031101. [Google Scholar]
  25. Chen, X.; Liu, Y. Finite Element Modeling and Simulation with ANSYS Workbench; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  26. Yu, C.; Qin, S.; Liu, Y.; Chai, B. Heat exchange performance optimization of a wheel loader cooling system based on computational fluid dynamic simulation. Adv. Mech. Eng. 2018, 10, 1687814018803984. [Google Scholar] [CrossRef] [Green Version]
  27. Yu, C.; Qin, S.; Chai, B.; Huang, S.; Liu, Y. The Effect of Compressible Flow on Heat Transfer Performance of Heat Exchanger by Computational Fluid Dynamics (CFD) Simulation. Entropy 2019, 21, 829. [Google Scholar] [CrossRef] [Green Version]
  28. Yang, L.; Sicheng, Q. The optimization design of Off-Highway machinery radiator based on genetic algorithm and ε-NTU. Acta Tech. CSAV 2017, 62, 465–476. [Google Scholar]
  29. Ji, B.; Wang, J.; Luo, X.; Miyagawa, K.; Xiao, L.Z.; Long, X.; Tsujimoto, Y. Numerical simulation of cavitation surge and vortical flows in a diffuser with swirling flow. J. Mech. Sci. Technol. 2016, 30, 2507–2514. [Google Scholar] [CrossRef]
  30. Hunt, J.C.R.; Wray, A.A.; Moin, P. Eddies, streams, and convergence zones in turbulent flows. Stud. Turbul. Using Numer. Simul. Databases 1988, 2, 193–208. [Google Scholar]
Figure 1. Virtual wind tunnel model. (a) The whole model; (b) the wheel loader.
Figure 1. Virtual wind tunnel model. (a) The whole model; (b) the wheel loader.
Processes 08 01687 g001
Figure 2. Simplified engine compartment.
Figure 2. Simplified engine compartment.
Processes 08 01687 g002
Figure 3. Engine compartment virtual wind tunnel model meshing (x = 0 section).
Figure 3. Engine compartment virtual wind tunnel model meshing (x = 0 section).
Processes 08 01687 g003
Figure 4. Flow rate collection points of experiment.
Figure 4. Flow rate collection points of experiment.
Processes 08 01687 g004
Figure 5. Simulation calculation and experimental data comparison (at 2000 rpm).
Figure 5. Simulation calculation and experimental data comparison (at 2000 rpm).
Processes 08 01687 g005
Figure 6. 2D contours of pressure magnitude at YZ plane.
Figure 6. 2D contours of pressure magnitude at YZ plane.
Processes 08 01687 g006
Figure 7. 2D contours of velocity magnitude at YZ plane.
Figure 7. 2D contours of velocity magnitude at YZ plane.
Processes 08 01687 g007
Figure 8. 2D contours of velocity magnitude at XY plane.
Figure 8. 2D contours of velocity magnitude at XY plane.
Processes 08 01687 g008
Figure 9. Vortex structures of the fan.
Figure 9. Vortex structures of the fan.
Processes 08 01687 g009
Figure 10. Prediction performance of three models.
Figure 10. Prediction performance of three models.
Processes 08 01687 g010
Table 1. CFD (computational fluid dynamics) model description.
Table 1. CFD (computational fluid dynamics) model description.
Analysis TypeTransient State
Turbulence ModelLES, DES, and Standard k ε model
Pressure-Velocity CouplingSIMPLEC
Transient FormulationSecond-order implicit
Fan StatusVaried from 0 rpm to 2000 rpm
Vent StatusOutflow
Environment Temperature25 °C
Inlet velocity1.5 m/s
Table 2. Comparison of computation time with different models (3.4 GHz, 2 CPU).
Table 2. Comparison of computation time with different models (3.4 GHz, 2 CPU).
ModelCPU Time (h)
RANS12.3
DES22.1
LES34.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, C.; Xue, X.; Shi, K.; Shao, M.; Liu, Y. Comparative Study on CFD Turbulence Models for the Flow Field in Air Cooled Radiator. Processes 2020, 8, 1687. https://doi.org/10.3390/pr8121687

AMA Style

Yu C, Xue X, Shi K, Shao M, Liu Y. Comparative Study on CFD Turbulence Models for the Flow Field in Air Cooled Radiator. Processes. 2020; 8(12):1687. https://doi.org/10.3390/pr8121687

Chicago/Turabian Style

Yu, Chao, Xiangyao Xue, Kui Shi, Mingzhen Shao, and Yang Liu. 2020. "Comparative Study on CFD Turbulence Models for the Flow Field in Air Cooled Radiator" Processes 8, no. 12: 1687. https://doi.org/10.3390/pr8121687

APA Style

Yu, C., Xue, X., Shi, K., Shao, M., & Liu, Y. (2020). Comparative Study on CFD Turbulence Models for the Flow Field in Air Cooled Radiator. Processes, 8(12), 1687. https://doi.org/10.3390/pr8121687

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