Next Article in Journal
Investigation of Hydrodynamically Dominated Membrane Rupture, Using Smoothed Particle Hydrodynamics–Finite Element Method
Next Article in Special Issue
Impact of Longitudinal Acceleration and Deceleration on Bluff Body Wakes
Previous Article in Journal
Instability and Convection in Rotating Porous Media: A Review
Previous Article in Special Issue
Stability Analysis on Nonequilibrium Supersonic Boundary Layer Flow with Velocity-Slip Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Turbulence Modeling Effects on the CFD Predictions of Flow over a Detailed Full-Scale Sedan Vehicle

1
North Carolina Motorsports and Automotive Research Center, UNC Charlotte, 9201 University City Blvd, Charlotte, NC 28223, USA
2
Hyundai America Technical Center, Inc., 6800 Geddes Road, Superior Township, MI 48198, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Current address: AXISCADES Engineering Technologies Ltd., Peoria, IL 61614, USA.
Fluids 2019, 4(3), 148; https://doi.org/10.3390/fluids4030148
Submission received: 3 July 2019 / Revised: 27 July 2019 / Accepted: 30 July 2019 / Published: 1 August 2019
(This article belongs to the Special Issue Turbulence and Transitional Modeling of Aerodynamic Flows)

Abstract

:
In today’s road vehicle design processes, Computational Fluid Dynamics (CFD) has emerged as one of the major investigative tools for aerodynamics analyses. The age-old CFD methodology based on the Reynolds Averaged Navier–Stokes (RANS) approach is still considered as the most popular turbulence modeling approach in automotive industries due to its acceptable accuracy and affordable computational cost for predicting flows involving complex geometries. This popular use of RANS still persists in spite of the well-known fact that, for automotive flows, RANS turbulence models often fail to characterize the associated flow-field properly. It is even true that more often, the RANS approach fails to predict correct integral aerodynamic quantities like lift, drag, or moment coefficients, and as such, they are used to assess the relative magnitude and direction of a trend. Moreover, even for such purposes, notable disagreements generally exist between results predicted by different RANS models. Thanks to fast advances in computer technology, increasing popularity has been seen in the use of the hybrid Detached Eddy Simulation (DES), which blends the RANS approach with Large Eddy Simulation (LES). The DES methodology demonstrated a high potential of being more accurate and informative than the RANS approaches. Whilst evaluations of RANS and DES models on various applications are abundant in the literature, such evaluations on full-car models are relatively fewer. In this study, four RANS models that are widely used in engineering applications, i.e., the realizable k ε two-layer, Abe–Kondoh–Nagano (AKN) k ε low-Reynolds, SST k ω , and V2F are evaluated on a full-scale passenger vehicle with two different front-end configurations. In addition, both cases are run with two DES models to assess the differences between the flow predictions obtained using RANS and DES.

1. Introduction

Since the last twenty years, Computational Fluid Dynamics (CFD) has played an important role in the early design stages of modern-day vehicles. In practice, Reynolds-Averaged Navier–Stokes (RANS) modeling still remains a widely-used approach in automotive applications of CFD, due to its acceptable accuracy, affordable cost, and fast turnaround time. However, due to the well-known drawbacks of RANS models in predicting complex separated flows, variants of the hybrid Detached Eddy Simulation (DES) approach are seen to gain popularity currently in many automotive applications. As demonstrated by many researchers [1,2,3], DES models generally offer an advantage over RANS models in better predicting the flow field and integral aerodynamic quantities like drag, lift, and moment coefficients. However, in spite of showing a superior flow prediction capability than the RANS approach, DES still cannot fully correctly capture the flow details. Additionally, the DES approach is still cost prohibitive for many applications; for example, for the recently developed DrivAer vehicle, the DES simulation is 17-times more expensive than the RANS computations [2]. In present day high-fidelity CFD analyses, simulations with 120 million cells are very common; even a RANS-based well-converged simulation of this size takes more 3000 core-hours for completion (see [4]), implying that a DES would have taken more than 50,000 core-hours for a single case. This makes the use DES analysis impractical for motorsports applications for sure and may be for the passenger vehicle industries as well. Subsequently, RANS approaches still remain a popular method in situations like the early design cycle of a new vehicle model where hundreds of design options need to be assessed.
In RANS approaches, the entire turbulence spectrum is approximated through a set of transport equations to solve the time-averaged mean flow field. These equations were established during the past several decades based on empiricism, mathematical analysis, and physical reasoning of the key mechanisms observed in canonical flows. RANS models have enjoyed remarkable success in various industrial applications. However, the inability of RANS models to predict complex separated flows correctly has motivated researchers to develop and apply higher order models in industrial applications. One promising model is the DES approach. DES is a hybrid modeling approach that combines the best of both RANS and LES: an RANS simulation in the viscosity-dominated boundary layers and an LES in the unsteady separated regions.
However, one major problem of the DES approach is how to deal with the spurious buffer layer between the RANS and Large Eddy Simulation (LES) regions [5,6]. Spalart et al. [7] proposed the Delayed Detached Eddy Simulation (DDES) to prevent DES from a too early switch to LES mode. Another modification is called the Improved Delayed Detached Eddy Simulation (IDDES), proposed by Shur et al. [8]. IDDES allows RANS to be used in a much thinner near-wall region and could provide the so-called Wall-Modeled LES (WMLES) capabilities to the DES formulation.
The existing literature contains a vast body of works on the evaluation of RANS and DES models for various flow configurations. Due to the scope of this study, only a few notable validation studies related to automotive flows will be discussed in this paper.
Due to the complexity and proprietary nature of the works associated with production passenger vehicles, the evaluations of turbulence models on full-scale vehicles are rare. Instead, researchers have mainly focused on simplified, generic car geometries, such as the Ahmed body [9] and the more recently introduced DrivAer models [10]. The Ahmed body was introduced in 1984 to analyze the key flow field features of ground vehicles [9]. The Ahmed body features a flat front with rounded corners and a sharp slanted rear upper surface with various angles. Since ground vehicles can be viewed as bluff bodies moving in close vicinity to the road, broad aerodynamic features of the Ahmed body resemble a ground vehicle, e.g., the flow separation on the rear window or the wake structure behind the body. The Ahmed body has been extensively studied by using both experimental and numerical methods. Until recently, it has been the most widely-used object for evaluation works on turbulence models related to automotive flows.
The main model variable for the Ahmed body geometry is the rear slant angle. The flow patterns on the rear slant and behind the body vary as the slant angle changes. Existing literature contains a vast amount of works that investigated the influence of the slant angle on the flow patterns. It was first shown by Ahmed et al. [9] in their introductory work on the Ahmed body that 30 degrees is the critical slant angle at which the topology of the flow is suddenly altered and the drag is reduced noticeably. Above this angle, the flow fully separates over the slant starting from its sharp edge due to the strong adverse pressure gradient between the slant and the roof. For such configurations, RANS models succeed in predicting the correct separation region, which agreed well with experimental data [2]. Below this critical angle, especially at 25 degrees, the flow is most complex. The flow still separates from the slant; however, the pressure gradient between the slant area and the side walls is still strong enough, which generates streamwise counter-rotating vortices at the lateral slant edges. These vortices induce a downward motion over the slant, mainly in the downstream part of the slant. Consequently, the flow that was partially detached from the upstream part of the slant can reattach further downstream, forming a confined recirculation zone over the rear slant. Moreover, the wake structures, which are dominated by the interactions between the separation zone and strong counter-rotating vortices, become very complex in this configuration. These vortices are also responsible for the high induced drag. Subsequently, the 25 degree slant angle model has become the most extensively studied variant of the Ahmed body geometry as it poses a significant challenge for CFD analyses. As demonstrated by [11,12,13,14], the majority of RANS models fail to predict the flow field past this 25 degree model correctly. These simulations either predict no separation or they cannot predict the correct size of the separation zone. Not only for RANS, the CFD simulation of the 25 degree Ahmed model is also challenging for the higher fidelity models, such as LES and hybrid LES/RANS [15]. For very small slant angles, i.e., smaller than 12 degrees, the flow is fully attached on the rear slant and only separates over the rear end of the body. For such simple flow patterns, most RANS models again perform very well in terms of both integrated quantities and flow structures [16]. To sum up, for the Ahmed car body, the performance of various turbulence models, especially RANS models, is highly influenced by the geometry and flow topology present. This raises the potential uncertainty in predicting automotive flows when using numerical approaches, especially RANS approaches.
A new realistic generic car model, called the DrivAer model, was proposed recently by the Institute of Aerodynamics and Fluid Mechanics of the Technical University of Munich (TUM), in cooperation with BMW and Audi AG. The DrivAer model is available in three typical car configurations, i.e., estate, fastback, and notchback. All three of these configurations have been experimentally investigated to provide data for numerical comparisons [10,17]. TUM researchers also investigated the fastback configuration numerically using the SST k ω turbulence model [18]. These studies found that although there are good agreements for drag coefficients between the experiment and simulation, observable discrepancies in pressure distributions on the top and bottom of the car exist between the CFD simulation and experiment results. A recent work by Peters et al. [19], which presented higher order structured finite difference solutions using overset grids and NASA’s overset grid solver Overflow, also observed similar pressure distribution differences when using the SST k ω model. Moreover, several widely-used RANS and hybrid LES/RANS models were also applied to the DrivAer model [1,2,20,21]. It was found that RANS models were unable to predict the force coefficients correctly, especially the lift coefficient. In addition, these studies indicated that, while inaccuracies still exist when compared to the experimental data, hybrid approaches predict both force coefficients and flow fields with better experimental correlations, but at the expense of much higher computational cost.
However, the evaluations of turbulence models on full-car geometries are quite rare. A recent work presented a study of RANS turbulence models on the prediction of aerodynamic characteristics of a NASCAR Gen-6 Cup race car [22]; this work showed that discernible differences in the predicted flow features, especially in the wake region, exist between the results from different RANS models. In this backdrop, the goal of the current study is to provide an evaluation of several RANS and DES approaches on a full-scale passenger vehicle. The RANS models investigated are the ones that are commonly used in various industrial external flow applications, i.e., the realizable k ε two-layer, Abe–Kondoh–Nagano (AKN) k ε low-Reynolds, SST k ω , and V2F models; further details about these models will be provided in the proceeding sections. The DES variants studied are the DDES and IDDES approaches. The prediction differences between these turbulence models were analyzed for two vehicle configurations in terms of integral aerodynamic forces, pressure, and velocity fields.

2. Simulation Details

2.1. Governing Equations

For incompressible flows, a two-equation RANS SST approach is based on the continuity, three momentum, and two additional transport equations. With x i denoting the coordinate axis in the i-direction and using the Einstein notation, the time-averaged continuity and momentum equations (Navier–Stokes equations) can be written as:
U i x i = 0
ρ U j U i x j = ρ f i ¯ + x j P δ i j + 2 μ S i j ρ u i u j ¯
where the symbols U i , P, μ , and ρ stand for the i-component of velocity, pressure, dynamic viscosity, and density, respectively; δ i j and f i ¯ are the Kronecker delta and resultant of body forces, and S i j ( 1 / 2 ) ( U i / x j + U j / x i ) is the rate of the mean strain tensor. The additional terms ρ u i u j ¯ in Equation (2) are known as the Reynolds stress terms, which are due to the fluctuating velocity field. Note that, using the Reynolds decomposition, a time-dependent quantity a ˜ can be written as a ˜ ( t ) = A + a ( t ) , where A is the time-averaged value of a ˜ ( t ) and a ( t ) is the fluctuating component. The Reynolds stress terms have a total of nine components, of which six of them are independent. Resolving the mean flow-field around any object requires solving the system of equations represented by Equations (1) and (2), which is now an unclosed system as it contains ten unknowns, yet only four equations. This is termed as the so called turbulence closure issue. The nonlinear Reynolds stress terms require additional modeling for closure, which can be achieved using a number of turbulence modeling approaches; for details, see [23]. These turbulence models differ from one another with respect to the choice of transported turbulence variables, turbulence viscosity formulation, or some additional functions embedded in the equations. The most popular turbulence viscosity relation is the Eddy Viscosity Model (EVM), which is based on an analogy between the molecular gradient-diffusion process and turbulent motion in which the Reynolds stress tensor is modeled as a function of mean flow quantities by the concept of turbulent eddy viscosity, ν t , the so-called Boussinesq approximation (to be defined later in Equation (5)).

2.2. Turbulence Models

Four commonly-used RANS models and two DES variants were investigated in this study. The RANS models include: the Realizable k ε (RKE) model of Shih et al. [24], the Abe–Kondoh–Nagano low-Re k ε or AKN model [25], the υ 2 ¯ f k ε or V2F model of Durbin [26], and the SST k ω model (SST) of Menter [27,28]. The DES variants include Delayed Detached Eddy Simulation (DDES) as proposed by Spalart and coworkers [7] and the Improved Delayed Detached Eddy Simulation (IDDES) model proposed by Shur et al. [8] and Gritskevich et al. [29].
As the names suggest, the RKE, AKN, and V2F turbulence models are formulated based on the Standard k ε model (SKE) of Launder and co-workers [30,31]. The k ε models use two additional transport equations for solving a fluid flow problem, one for the turbulent kinetic energy or TKE, denoted by k, and the other for the rate of TKE dissipation, denoted by ε . One of the main flaws of the standard k ε model is that it may predict an anomalously large growth of the Turbulent Kinetic Energy (TKE) in stagnation point flows. The usual explanation for the stagnation point anomaly is that the linear eddy viscosity formulation gives an erroneous difference in normal stresses. One way of circumventing the stagnation anomaly is by imposing a “realizability” constraint on the normal stresses [32,33]. All three k ε models investigated in this study applied the realizability; the AKN and V2F models achieved the realizability through a turbulent time scale, while RKE achieved it by using a damping function. The other difference between these three k ε variants is the way these models are applied to the viscosity-affected layer, which includes the viscous sub-layer, buffer layer, and the overlap region. AKN uses a low-Reynolds number approach by applying damping functions to some of the coefficients in the model. RKE adopts a two-layer approach [34], which solves for turbulent kinetic energy k, but prescribes turbulent dissipation rate ε algebraically with distance from the wall in the viscosity-dominated near-wall flow regions. On the other hand, the V2F model is valid up to the solid walls, as it is designed to handle wall effects in turbulent boundary layers and to accommodate non-local effects [26,35].

2.2.1. Realizable k ε Model

The Realizable k ε (RKE) model of Shih et al. [24] is reported in the literature to perform better, compared to the SKE, in flows involving complex geometries or over bluff bodies where flow separation and recirculation occur. The RKE improves the standard k ε model of Launder and coworkers [30,31] by introducing a new realizable formulation of the eddy viscosity and a new modeled dissipation rate equation, as given in Equation (4). This dissipation transport equation is based on the dynamic equation of the mean-squared vorticity fluctuation at large turbulence Reynolds number. The transport equations and algebraic relations relevant to the RKE model are given below.
U j k x j = τ i j U i x j ε + x j ν + ν t σ k k x j
U j ε x j = C 1 ε k τ i j U i x j C 2 ε 2 k + ν ε + x j ν + ν t σ ε ε x j
where:
τ i j = u i u j ¯ = 2 ν t S i j 2 3 k δ i j
In Equation (5), ν t is the turbulent eddy viscosity, and S i j is the mean rate of the strain tensor defined as:
S i j 1 2 U i x j + U j x i .
The eddy viscosity formulation, calculated the same as the SKE model and given in Equation (7), is based on the realizability constraints, which ensures the positivity of normal Reynolds stresses and the Schwarz inequality for turbulent shear stresses.
ν t = C μ k 2 ε
The other difference between the RKE and SKE turbulence eddy viscosity formulation is that the closure coefficient C μ in Equation (7) is a variable given by Equation (8) as opposed to a constant as in the SKE model.
C μ = 1 A 0 + A s κ U * ε
In Equation (8), A 0 is a constant, while the other variable coefficient A s is defined as:
A s = 6 cos ϕ
where ϕ = 1 3 cos 1 6 W , W = ( S i j S j i S k i ) / S ˜ and S ˜ = S i j S j i .

2.2.2. Abe–Kondoh–Nagano k ε Model

The AKN model is shown to work well with low-Reynolds number complex flows [25,36]. The main improvement of the AKN is the usage of the Kolmogorov velocity scale u η , u η ( v ε ) 1 / 4 , instead of the friction velocity, u τ , to account for the near-wall and low-Reynolds number effects in the attached and detached flows [25]; u τ τ w / ρ , where τ w is the shear stress at the wall. The k-transport equation for the AKN model is the same as that of the RKE model as given in Equation (3), while the ε equation is different as given below:
U j ε x j = C ε 1 ε k τ i j U i x j C ε 2 f ε ε 2 k + x j ν + ν t σ ε ε x j
The eddy viscosity, ν t , in the AKN model is defined as:
ν t = C μ f μ k 2 ε .
The damping functions f μ and f ε , defined below, ensure two things. Firstly, the effect of the presence of multiple length scales in shear flows is taken care of, and secondly, the solution meets the requirement of the near-wall limiting turbulence behavior.
f μ = 1 e y + 14 2 1 + 5 R t 3 / 4 e R t 200 2 ,
f ε = 1 e y + 3.4 2 1 0.3 e R t 6.5 2 ,
where R t = k 2 / ν ε , and the normalized wall distance y + is defined as y + u τ y / ν .

2.2.3. υ 2 f k ε Model

The V2F model, developed by Durbin [26], is similar to the SKE, but incorporates the near-wall turbulence anisotropy and non-local pressure-strain effects. The model solves four transport equations, those for k, ε , velocity scale υ 2 , and elliptic relaxation factor f. This is a general low-Reynolds number model, which does not require any wall function and is valid all the way up to the wall. The model was originally developed for attached or mildly separated boundary layers; however, subsequent studies showed that it can also accurately simulate separation-dominated flows (see [37]).
With T as the turbulence length scale, the turbulent viscosity, ν t , is defined in the V2F model as:
ν t = C μ υ 2 ¯ T
Note that the value of the coefficient C μ in Equation (14) is different from the value of the same symbol as used in the SKE model. The additional transport equations that the V2F model solves for the velocity scale υ 2 and elliptic relaxation factor f are:
υ 2 ¯ t + U j υ 2 ¯ x j = k f υ 2 ¯ ε k + x j ν + ν t σ υ 2 ¯ x j
L 2 f 2 x j f = C 1 1 T υ 2 ¯ k 2 3 C 2 P k ε
where L is the turbulence length scales and P k represents the production of turbulence kinetic energy due to the mean flow velocity gradients. The turbulence length and velocity scales, L and T, respectively, are defined as follows:
L = C L max k 3 / 2 ε , C η ν 3 ε 1 / 4
T = max k ε , C T ν ε 1 / 2

2.2.4. SST k ω Model

Another widely-used two-equation model is the standard k ω model (SKW) of Wilcox [38,39]. However, a well-known flaw of the standard k ω model is its sensitivity to free-stream/inlet conditions. Menter [27] addressed this problem in the SST model, which effectively blends a k ε model in the far-field with a k ω model near the wall. Thus, the SST model does not have the problem of sensitivity to free-stream/inlet conditions and is valid up to the wall. A modification is also made to the definition of the eddy-viscosity, which accounts for the effect of the transport of principal turbulent shear stress [27]. The relevant equations for the SST model are given below.
ν t = a 1 k max ( a 1 ω , S F 2 )
U j k x j = P k β * k ω + x j ν + σ k ν t k x j
U j ω x j = γ S 2 β ω 2 + x j ν + σ ω ν t ω x j + 2 1 F 1 σ ω 2 1 ω k x i ω x i
F 1 = tanh min max k 0.09 ω d , 500 ν d 2 ω , 2 k C D k ω d 2 4
F 2 = tanh max 2 k β * ω d , 500 ν d 2 ω 2
C D k ω = max 1 ω k x i ω x i , 10 20
γ = F 1 γ 1 + ( 1 F 1 ) γ 2
γ 1 = β 1 β * σ ω 1 κ 2 β *
γ 2 = β 2 β * σ ω 2 κ 2 β * .
In the above equations, α , β , β * , σ k , σ ω , and σ ω 2 are turbulence model closure coefficients, and F 1 and F 2 are blending functions. F 1 blends the near-wall value or the SKW model value of a model closure coefficient with its SKE or far away from the wall value. That means:
σ k = F 1 × σ k 1 + ( 1 F 1 ) × σ k 2 .
In Equation (28) the second subscripts 1 and 2 imply the value of σ k as obtained from the SKE and SKW models, respectively.

2.2.5. DES Models

Besides the RANS approaches that model all scales, another approach to model turbulence is called the Scale-Resolving Simulation (SRS), which resolves the large scales of turbulence and models the small-scale motions only. A common scale-resolving simulation is called the Large Eddy Simulation (LES). The problem with LES is that it is still too computationally costly for automotive flows. The Detached Eddy Simulation (DES), proposed by Spalart [6,40], is a hybrid modeling approach that combines the best of both RANS and LES: an RANS simulation in the boundary layers and an LES in the unsteady separated regions. Two variants of the DES were investigated in this study: the Delayed Detached Eddy Simulation (DDES) and Improved Delayed Detached Eddy Simulation (IDDES).
A number of studies revealed a major flaw of the DES methodology that these simulations are prone to producing an artificial separation when refining the maximum cell edge length inside the boundary layer below a critical value (see Menteret al. [41]). This phenomenon is known as the Grid-Induced Separation (GIS) as this fictitious separation depends on the grid spacing and not the flow physics. GIS can in principle be avoided by shielding the RANS model from the DES formulation for wall boundary layers. Later, Spalart et al. [7] proposed a generic formulation of the shielding function, which depends only on the eddy-viscosity and the wall distance, and termed this implementation as the Delayed Detached Eddy Simulation (DDES). Essentially, the DDES model incorporates a delay factor that enhances the ability of the model to distinguish between the LES and RANS regions on meshes where spatial refinements could give rise to ambiguous behavior. The Improved Delayed Detached Eddy Simulation (IDDES) model consists of a combination of various new and existing techniques aimed at providing a more flexible and convenient DDES for high-Reynolds number flows; see [8,29]. In addition to providing a shielding against GIS, IDDES allows the model to run in Wall-Modeled LES (WMLES) mode. However, the model formulation is relatively complex for both DDES and IDDES, and, as such, for brevity, the details of the relevant equations are omitted; the interested reader is referred to the original publications of Shur et al. [5] and and Gritskevich et al. [29].
Please note that all RANS models investigated in this paper used the all- y + wall treatment, which is a hybrid treatment that emulates the low- y + wall treatment for fine meshes (when the wall-node y + is less than five), and the high- y + wall treatment for coarse meshes (when the wall-node y + is greater than 30). Additionally, this all- y + approach makes use of a blending function to produce reasonable answers for meshes of intermediate resolution, that is when the wall-cell centroid falls within the buffer region of the boundary layer ( 5 < y + < 30 ).

2.3. Vehicle Model

The vehicle model used in this study was a 2014 Hyundai Veloster with detailed underhood and underbody components as shown in Figure 1. The length, width, and height dimensions of this vehicle were 4220, 1800, and 1400 mm, respectively. The front ground clearance was kept at 143 mm for all simulations. The geometry preparation process was performed using a commercial pre-processor code, ANSA Version 15.02. This involved the removal of unnecessary vehicle components that did not interact with the external airflow, e.g., the vehicle interior parts. Modifications of some lesser important details were also made to further simplify the model, such as removal of bolts and nuts, closure of small holes and gaps, etc. All parts that had any impact on the external flow field were retained. The cleaned geometry was surface meshed in ANSA, and the mesh in NASTRAN format was imported into a commercial CFD code, Star-CCM+ Version 12.02 by CD-adapco, for volume-meshing, simulation, and post-processing.

2.4. Cases Investigated

Two cases with different front-end configurations were investigated in this study. Case 1 was the baseline case with a vertical Condenser, Radiator, and Fan Module (CRFM) and no air-duct. Case 2 differed from Case 1 in two aspects: (1) it had an inclined CRFM and included an air-duct, as shown in Figure 2; and (2) the side grille openings on either side of the primary grille opening were sealed off in this case; see Figure 3. The combination of air-duct, inclined CRFM, and side grille opening sealing were investigated to achieve a simultaneous improvement of the radiator performance and cooling drag reduction. The benefits of such design changes have been well demonstrated by the authors in a previous study [42].

2.5. Boundary and Operating Conditions

The vehicle model was simulated in a sufficiently large virtual wind tunnel. The diagram of the virtual tunnel is shown in Figure 4, where L represents the overall length of the vehicle. Velocity inlet and pressure outlet boundary conditions were applied to the inlet and outlet, respectively. A moving wall boundary condition was applied on the ground with a tangential speed the same as the inlet air speed. The top and side walls were modeled as symmetry planes. The inlet air speed was set at a typical highway speed at 70 mph or 31.293 m/s.
The condenser and radiator were modeled as porous mediums. The cooling fan rotation was simulated using the Moving Reference Frame (MRF) method, which can provide reasonable pressure jump across the fan and swirl components downstream without experimental data as input. The wheel rotation was simulated by setting the tangential velocity the same as the inlet air speed. In addition, the MRF method was applied to the regions between the spokes to simulate the pumping effect.

2.6. Mesh and Solver

The volume mesh primarily consisted of hexahedral cells generated by the trimmed mesher. In addition, to capture flow features with a high velocity gradient close to the solid boundaries properly, 12 prism layers with the first node height of 0.005 mm and a total thickness of 6 mm were applied on the external surfaces of the vehicle and the ground. To achieve a good balance between the accuracy and computational cost, mesh refinements were chosen to use volume sources that followed the external shape of the vehicle. This technique can generate a high quality volume mesh with a number of cells much fewer than that of traditional cubic block volume sources. Furthermore, volume mesh refinements were also applied on regions with complex flow patterns or separation, e.g., under-hood, under-body, wheels, and wake region. The cell size of the core mesh around the vehicle was 6 mm. Using the above mesh setting, a high-quality volume mesh with 75 million cells was generated, of which only 0.003% of the cells showed volume change less than 1e-4 and 0.003% of the faces had face validity less than 0.95. All simulations presented in this study used the same mesh setting. The mesh on the central plane of the vehicle is shown in Figure 5. Using the aforementioned prism layer setting, the y + values on the external surfaces of the vehicle in all the cases were mostly less than 0.5, as can be seen in Figure 6.
The finite volume method was used for the spatial and temporal discretization of the governing equations. The integral form of the conservation equations of mass and momentum were solved by the segregated flow solver in a sequential manner, i.e., the non-linear governing equations were solved iteratively one after the other for the solution variables such as the three velocity components and the pressure. The pressure-velocity coupling was achieved with a predictor-corrector approach called SIMPLE (Semi-Implicit Method for Pressure-Linked Equations). For RANS models, the discretization scheme for computing the convection flux on a cell face was Hybrid MUSCL (Monotonic Upwind Scheme for Conservation Laws) third-order/central differencing; for DES variants, the discretization scheme was a hybrid second-order-upwind bounded-Central. A second-order upwind scheme was used for computing the convection flux on a cell face for each turbulence model solver. The interested reader is directed to the Star-CCM+ user manual for further details.

3. Results and Discussions

All the simulations were run on a cluster with 96 computing cores. These simulations were run until the flow reached an acceptable convergence judged by two criteria. Firstly, the normalized residuals of all transported quantities were steady, and secondly, the drag force coefficient, averaged over last 500 iterations and 0.4 s (3 flow units) for RANS and DES, respectively, did not show fluctuations over three counts. Based on these two stopping criteria, the RANS and DES models took around 28 and 250 h, respectively, to be completed. The DES simulations were run with a Courant-Friedrichs-Lewy or CFL number of one for a period, T, of 15 convective flow units ( T = 15 × L / U = 2 s , where L is the car length and U is the freestream velocity). The CFL number, C, is defined as, C | U x | Δ t / Δ x , where U x is the local cell velocity, and Δ x and Δ t are local grid and time resolutions (time step), respectively. In order to have a comprehensive understanding of the predictive differences between the six models, the proceeding discussions will cover analyses of force coefficients, velocity, pressure field, wake structures, and some other flow fields that are crucial to the flow features.
A mesh independence study was conducted first for the baseline case with the RKE and IDDES models to further eliminate the effects of the mesh on simulation results. For the RKE model, the fine mesh contained 100 million cells, which were generated by reducing the base size from 24 mm to 20 mm. A change of drag coefficient within 1% between 75 and 100 million cells was achieved in this process. However, for the IDDES model, simply reducing the base size will increase the computational cost dramatically, not only due to the increased cell number, but also due to the smaller time step resulting from reducing the core mesh size close to the vehicle. Thus, a different strategy, where only the cells in the wake region and underbody area were halved, was used to refine the mesh for the IDDES model. As suggested by [43], the recommended mesh size in the detached flow regions for a DES-type simulation should be close to the Taylor microscale given by λ = 15 1 / 2 A 1 / 2 R l 1 / 2 l , where A is an undetermined constant set to be 0.5, l is the bulk integral scale (taken as the vehicle length, L), and R l the integral scale Reynolds number (here taken as the Reynolds number based on the free-stream velocity and vehicle overall length L). For the Reynolds number used in this study, the Taylor microscale λ was estimated to be 7.9 mm. Therefore, the mesh strategy used to refine the IDDES case with a mesh size of 6 mm in the wake region and underbody was reasonably fine enough. The fine mesh of the IDDES model (shortened to IDDES fine in this study) contained 87 million cells by using this strategy. Further to eliminate the numerical uncertainty, the IDDES fine case was run with a CFL number of 0.5. It was found that the change of drag coefficient was within 1%; thus, the coarser mesh with a CFL number of one is sufficient enough for performing the intended investigations.

3.1. Drag and Lift Forces

Due to the proprietary nature of the data, the drag coefficient values presented in this paper were normalized using an arbitrary reference area, such that the baseline car C D obtained from the RKE model matched with the one available in the public domain. Data downloaded from the website Edmunds.com,“features & specs 2014 Hyundai Veloster”, https://www.edmunds.com/hyundai/veloster/2014/features-specs/, accessed 19 June 2017, indicates that the drag coefficient of the 2014 model Hyundai Veloster is 0.320. Thus, the C D for the baseline case obtained from the RKE model was set to 0.320, and the predictions from all other approaches were nominalized accordingly.
Table 1 shows the percent deviation of the predicted C D values from the wind-tunnel test data for all investigated models; the data presented in this table correspond to the baseline vehicle model. Firstly, it can be seen that C D prediction from the RKE agreed very well with the experimental data. However, all other models predicted significantly higher drag coefficient, with the highest C D value coming from the AKN model. Unfortunately, while both the DDES and IDDES models, which apply the RANS SST model in the boundary layers, predicted a slightly better drag coefficient compared with that of the SST model, they failed to predict well-correlated C D values. Furthermore, for the DES models, while the C D prediction of the IDDES was slightly better than that of the DDES, the IDDES model with a fine mesh produced negligible improvement. The fact that even the IDDES fine case failed to predict well correlated C D values indicates the monumental challenge that is still encountered by the CFD simulation of automotive flows. Although the RKE model predicted a significantly better C D value than all other investigated models, it should be cautioned that this does not suggest the superior performance of the RKE model. For the RANS and DES models, the drag prediction may contain error cancellation of flow fields in different regions around the vehicle.
Table 2 summarizes drag and lift coefficients corresponding to both vehicle configurations as obtained from the studied turbulence models. These include the coefficients of drag, total lift, front and rear lift, and the front to total lift ratio, known as the lift balance or percent front. Data presented in Table 2 are graphically summarized in Figure 7, Figure 8 and Figure 9 to illustrate the predictive differences between all models in terms of drag, lift, and percent front. It can be seen that the RKE model predicted significantly smaller C D values for both cases compared to the other models. For instance, the RKE predicted a 11.1% less C D value than the AKN for Case 1 and a 17.1% smaller C D value than the DDES model for Case 2. Moreover, while all models predicted a higher C D value for the baseline case than that for the air-duct case, compared to the RANS predictions, the two DES variants predicted a much smaller difference between the two cases; see Table 3. For instance, while the IDDES predicted a Δ C D of three counts between the two car configurations, the RKE model predicted a 30 count difference. While the variation of Δ C D predictions between the RANS was within ± 5 counts, the spread of Δ C L predictions was much larger; note that, for any force coefficient C F , Δ C F is defined as the C F value from Case 1 minus the corresponding value from Case 2. Clearly, the predictive disagreements between the RANS models were most visible with the front lift ratio. Further, comparing the two IDDES simulations with different mesh strategies for the baseline case, while there was no difference in C D prediction, a noticeable discrepancy can be seen in C L prediction. In short, for all the turbulence models and mesh strategies studied, the prediction of drag was more consistent than that of lift. Moreover, even for industrial automotive flows that feature complex geometry and exhibit massive separations, a certain model could predict the C D value that correlates very well with the experimental data.
The accumulated drag coefficients along the vehicle centerline for the baseline case as obtained from all turbulence models are shown in Figure 10. It can be seen, firstly, that before the airflow passed through the CRFM ( 0 < X / L < 0.06 ), all the models predicted almost identical results; then, in the underhood region from the CRFM to the cowl area ( 0.06 < X / L < 0.25 ), the RKE predicted a significantly lower drag value compared to all other models. Nevertheless, the results from the RKE in this region still followed the same trend as all other models. Further, all the models predicted almost identical trends, though with different magnitudes, in the region between the cowl area and the rear spoiler where the flow was mostly attached ( 0.25 < X / L < 0.85 ). However, at the rear end of the vehicle, the predictions from all the models started to diverge from each other significantly, especially that of the DDES and the RKE. In short, for the vehicle model used in this study, all the turbulence models investigated could predict consistent results even in the underhood compartment where the flow was featured with separations and complex interactions between the cooling airflow and underbody airflow. However, no consistent results between all the models were observed for the prediction of the wake region where massive separation occurred.
The delta accumulated drag coefficients (between Case 1 and Case 2 vehicle configurations) obtained using all turbulence models are shown in Figure 11. Firstly, it is not surprising to observe significant prediction fluctuations in the front part of the vehicle since the only difference between the two vehicle configurations was the front end geometry. For this part of the vehicle, although all turbulence models predicted similar trends, noticeable differences in magnitude indicated a high sensitivity of the C D prediction to the choice of turbulence models. Further, the noticeable fluctuations in the rear end of the vehicle, especially for the RKE, DDES, and IDDES models, highlighted the strong influence of the predictions of upstream flow on the downstream flow. Moreover, the prediction trend of the RKE model differed significantly from all other models in two areas: the windshield and the rear end. The causes of prediction differences will be qualitatively and quantitatively analyzed in the following sections.

3.2. Pressure and Velocity Fields

3.2.1. Baseline Case

As stated in the previous section, the front and rear ends of the vehicle contributed most to the C D difference for all the turbulence models studied. To better analyze the C D prediction differences for all the models, pressure and velocity fields at the front and rear ends of the vehicle were investigated.
For the baseline case, since the RKE predicted a significantly lower C D value, the pressure coefficient on the back of the vehicle from this model is shown first in Figure 12. It can be seen that the most noticeable higher pressure area was right above the exhaust pipe. This area was the saddle of attachment of the wake flow as supported by the streamline structure in Figure 13 for the RKE model (Top left). To better illustrate the differences between the RKE and other models for the baseline case flow prediction, contours of Δ C P , where Δ C P represents C P from a model relative to C P obtained from the RKE, are shown in Figure 14. Firstly, all other models predicted overall lower pressure when compared with that of the RKE, especially the three RANS models. This is consistent with the results of accumulated drag as shown in Figure 10. More specifically, the Δ C P predicted by the AKN model might be explained by the key differences between the RKE and the AKN model on how these two models solve the viscosity-dominated near-wall flow regions. The RKE model used a two-layer approach (the Star-CCM+ implementation of this approach was after [34]), and the AKN model adopted a low-Reynolds number approach with the introduction of the Kolmogorov velocity scale, instead of the friction velocity, to account for the near-wall and low-Reynolds number effects in both attached and detached flow. On the other hand, both of V2F and SST models can resolve down to the viscosity-dominated near-wall regions. Thus, more evident prediction differences in certain areas, such as the area around the right taillight and both sides of the vehicle, are expected. Further, the delta pressure distributions predicted by the two DES variants showed a completely different pattern compared to those from the other three RANS models. Interestingly, they both predicted less discrepancies against that of the RKE model, with slightly higher pressure on the upper part, but lower pressure on the lower part.
The vector scenes with constrained streamlines presented insights into the wake structure, as demonstrated in Figure 13. Firstly, the RKE model predicted two strong counter-rotating vortices with the saddle-like attachment area right above the exhaust pipes (see Figure 12). Secondly, compared to the RKE model, the other three RANS models predicted significantly different wake structures: (1) no counter-rotating vortices existed; (2) the lower vortex was much weaker and the saddle of attachment area at a lower location; (3) the wake area was much smaller (see the areas on the right of the red spline line in the corresponding images). On the other hand, again interestingly, two DES variants predicted a mostly similar wake structure as that of the RKE model, e.g., counter-rotating vortices, saddle of attachment location, and similar wake area. However, some differences can still be observed: (1) the top vortex core was at a further downstream location; (2) wake structures were more complex, especially as the DDES predicted a tertiary vortex near the exhaust pipes, and the IDDES predicted a tertiary vortex under the spoiler.
To further investigate the prediction differences between all the models studied, streamwise velocity profiles for the baseline case are shown at four selected locations. These lines are marked with red lines in Figure 15 and are located (1) on the hood, (2) behind the antenna, (3) on the rear end of the spoiler, and (4) on the rear glass. The velocity profiles are shown in Figure 16. It can be seen that, for Locations 1 and 3, where the flow was attached, similar trends can be observed in the predictions from all models except for the AKN model at Location 1 on the hood. At Location 3, right on the rear end of the spoiler, where the flow was mostly attached and fully turbulent, all the models showed similar predictions in the viscosity-dominated near-wall region. This implies that the wall treatments used by different models worked rather well in the attached flow region. However, above the viscosity-dominated region, the models studied start to show differences in terms of boundary layer growth. Clearly, the RKE showed the smallest boundary layer thickness, and the IDDES predicted the thickest boundary layer at Location 3. The unexpected streamwise velocity profile from the AKN at Location 1 looked quite suspicious, and indicated that the low-Reynolds number approach implemented in the AKN model may not work well even in the attached flow region. On the other hand, for Locations 2 and 4, which were in wake regions, all models studied differed from each other significantly. In short, for the models studied, while a consistent streamwise velocity profile could be observed in attached flow regions except for the AKN model, the significant different predictions of the streamwise velocity in the wake regions further showed a strong dependence of the turbulence models on the prediction of complex automotive flows.
For the sake of a more comprehensive analysis of the flow fields, three plane sections were selected for flow visualizations. As shown in Figure 17, these planes are: (1) through the underhood area; (2) at the rear of the vehicle; (3) inside the wake region. Streamwise velocity scenes on these three selected plane sections from all models for the baseline case are shown in Figure 18. It can be seen that, on Plane 1 at X / L = 0.15 , all models predicted a similar streamwise velocity distribution. However, compared to the other three RANS models, the RKE model predicted a larger negative streamwise velocity area in the upper part of the wake region on planes X / L = 1.0 1.1 (see the area bound by the black lines). Moreover, the RKE model predicted a smaller negative streamwise velocity area in the lower part of the wake region on plane X / L = 1.1 (see the yellow dashed lines). Surprisingly, the predictions from the two DES variants were quite similar to that of the RKE model on both selected planes. The similarity of the streamwise velocity profile on these two planes between the predictions from the DES variants and the RKE model supports that the RKE model performed better than the other three RANS models for the baseline case. This can be further confirmed by the delta streamwise velocity against the IDDES fine model on these planes, as shown in Figure 19. It can be seen that the RKE predicted closer velocity profiles in the wake region to that of the IDDES fine model than the other three RANS models. Moreover, the IDDES case with the coarse mesh predicted very similar results as that from the fine mesh, which indicates that the coarse mesh was sufficient enough for this study.

3.2.2. Comparison between the Baseline Case and the Air-Duct Case

This section presents analyses of how the selected turbulence models predicted the flow differences between two car configurations. This was done by presenting delta values of pressure and velocity; delta pressure or velocity means the quantity for Case 2 (the air-duct case) subtracted from that of Case 1 (the baseline case).
Scenes of delta pressure coefficient ( C p of Case 2 relative to the baseline Case 1), as obtained from the selected models, are shown in Figure 20. Firstly, as expected, all the models predicted a higher pressure coefficient in the front end part of the vehicle for the baseline case, since the larger front grille opening in the baseline configuration will result in more airflow entering the underhood compartment. However, noticeable differences can be observed for the pressure distribution predictions on the rest of the car. In particular, the prediction from the RKE model differed significantly from that of all other models; the RKE model predicted noticeable lower pressure on the rest of the car. This is consistent with the delta accumulated drag plot discussed earlier in Figure 11. Further, the two DES variants predicted lower pressure on the rear bottom part of the vehicle. These observations indicate that the change of upstream geometry may result in noticeably different downstream predictions of the flow field depending on the choice of turbulence models.
Scalar scenes of delta mean streamwise velocity (the mean streamwise velocity of Case 2 relative to the baseline Case 1), as obtained from the selected models, are shown in Figure 21; see Figure 17 for the location of these planes. The most significant differences in delta mean streamwise velocity at X / L = 0.15 corresponded to RKE predictions near the far left and right sides. These two area denote the region of the flow most influenced by the tire rotation effects. This indicates that the prediction of flows subject to rotation was different for the RKE model compared to the rest RANS model. Surprisingly, the RKE deltas were much more aligned with the predictions from the scale-resolved DDES or IDDES models, which again strengthens the earlier argument that the RKE-based RANS approach produced a more realistic flow field compared to the other steady state RANS models for this particular vehicle geometry. Another interesting point is that the RKE model predicted larger delta pressure distributions on most part of the vehicle, as shown in Figure 20. Furthermore, the delta streamwise velocity on the two planes in the wake regions predicted by the RKE model was much smaller than that of the other three RANS models. On the other hand, the two DES variants predicted a slightly larger deviation than that of the RKE model, but much smaller than that of the other three RANS models.

4. Conclusions

This study presented a numerical evaluation of various commonly-used RANS models and hybrid RANS/LES models on a full-scale passenger vehicle with two different front-end configurations. The effect of turbulence modeling on the prediction of drag and lift forces, pressure, and velocity field were investigated. It was found that, firstly, the RKE model predicted a well-correlated C D value for the baseline case; however, all the models investigated still needed significant improvement in order to predict the lift values consistently. Secondly, despite the DES variants not showing superior performance over RANS models for predicting the absolute C D value, they were able to predict more complex flow structures with a higher fidelity. For example, in Figure 13, the DDES model predicted a tertiary vortex near the exhaust pipes, and the IDDES model predicted a tertiary vortex under the spoiler. The similarity of the streamwise velocity profile and turbulent kinetic energy distribution on the selected planes between the predictions from the RKE model and the DES variants supports that the RKE model performed better than the other three RANS models for the baseline case.
In summary, for automotive flows that feature complex geometry and exhibit massive separation, the success of CFD simulation still heavily relies on the choice of mesh strategy and turbulent models, which in turn depend on the experience of the engineer. RANS models might still be a good choice for the prediction of drag values due to their acceptable accuracy, low computational cost, and fast turnaround time. However, for applications that require predictions of the detailed flow features, such as aeroacoustic, soiling, and water management, the transient DES type approaches would be a better choice.

Author Contributions

Conceptualization, C.Z. and M.U.; methodology, C.Z. and M.U.; software, C.Z.; validation, C.Z. and M.U.; formal analysis, C.Z., C.P.B. and M.U.; investigation, C.Z. and M.U.; resources, L.F. and M.U.; data curation, C.Z., C.P.B., and M.U.; writing, original draft preparation, C.Z., C.P.B., and M.U.; writing, review and editing, C.Z., C.P.B., L.F., and M.U.; visualization, C.Z.; supervision, M.U.; project administration, L.F. and M.U.

Funding

This research received no external funding.

Acknowledgments

Lead author Chuhui Zhang was supported by the Department of Mechanical Engineering and Engineering Science at UNC Charlotte as a Teaching Assistant during the execution of the project. The authors also acknowledge with thanks computational supports from UNC Charlotte College of Engineering MOSAIC computing and UNC Charlotte University Research Computing (URC).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ashton, N.; Revell, A. Comparison of RANS and DES Methods for the DrivAer Automotive Body; Technical Report, SAE Technical Paper; SAE: Warrendale, PA, USA, 2015. [Google Scholar]
  2. Ashton, N.; West, A.; Lardeau, S.; Revell, A. Assessment of RANS and DES methods for realistic automotive models. Comput. Fluids 2016, 128, 1–15. [Google Scholar] [CrossRef] [Green Version]
  3. Guilmineau, E.; Deng, G.; Leroyer, A.; Queutey, P.; Visonneau, M.; Wackers, J. Assessment of hybrid RANS-LES formulations for flow simulation around the Ahmed body. Comput. Fluids 2017, 176, 302–319. [Google Scholar] [CrossRef]
  4. Fu, C.; Bounds, C.P.; Selent, C.; Uddin, M. Turbulence modeling effects on the aerodynamic characterizations of a NASCAR Generation 6 race car subject to yaw and pitch changes. Proc. Inst. of Mech. Eng. Part D J. Automob. Eng. 2019. [Google Scholar] [CrossRef]
  5. Shur, M.; Spalart, P.; Strelets, M.; Travin, A. Detached-eddy simulation of an airfoil at high angle of attack. In Engineering Turbulence Modelling and Experiments 4; Elsevier: Amsterdam, The Netherlands, 1999; pp. 669–678. [Google Scholar]
  6. Spalart, P.R. Young-Person’s Guide to Detached-Eddy simulation Grids; Technical Report; NASA Langley Research Center: Hampton, VA, USA, 2001. [Google Scholar]
  7. Spalart, P.R.; Deck, S.; Shur, M.L.; Squires, K.D.; Strelets, M.K.; Travin, A. A new version of detached-eddy simulation, resistant to ambiguous grid densities. Theor. Comput. Fluid Dynam. 2006, 20, 181. [Google Scholar] [CrossRef]
  8. Shur, M.L.; Spalart, P.R.; Strelets, M.K.; Travin, A.K. A hybrid RANS-LES approach with delayed-DES and wall-modeled LES capabilities. Int. J. Heat Fluid Flow 2008, 29, 1638–1649. [Google Scholar] [CrossRef]
  9. Ahmed, S.; Ramm, G.; Faltin, G. Some Salient Features of the Time-Averaged Ground Vehicle Wake; Technical Report, SAE Technical Paper; SAE: Warrendale, PA, USA, 1984. [Google Scholar]
  10. Heft, A.; Indinger, T.; Adams, N. Investigation of unsteady flow structures in the wake of a realistic generic car model. In Proceedings of the 29th AIAA Applied Aerodynamics Conference, Honolulu, HI, USA, 27–30 June 2011; p. 3669. [Google Scholar]
  11. Han, T. Computational analysis of three-dimensional turbulent flow around a bluff body in ground proximity. AIAA J. 1989, 27, 1213–1219. [Google Scholar] [CrossRef]
  12. Gilliéron, P.; Chometon, F. Modelling of stationary three-dimensional separated air flows around an Ahmed reference model. In ESAIM: Proceedings; EDP Sciences: Les Ulis, France, 1999; Volume 7, pp. 173–182. [Google Scholar]
  13. Guilmineau, E. Computational study of flow around a simplified car body. J. Wind Eng. Ind. Aerodyn. 2008, 96, 1207–1217. [Google Scholar] [CrossRef]
  14. Ashton, N.; Revell, A. Key factors in the use of DDES for the flow around a simplified car. Int. J. Heat Fluid Flow 2015, 54, 236–249. [Google Scholar] [CrossRef]
  15. Krajnović, S.; Davidson, L. Large-Eddy Simulation of the Flow around Simplified Car Model; Technical Report, SAE Technical Paper; SAE: Warrendale, PA, USA, 2004. [Google Scholar]
  16. Minguez, M.; Pasquetti, R.; Serre, E. High-order large-eddy simulation of flow over the “Ahmed body” car model. Phys. Fluids 2008, 20, 095101. [Google Scholar] [CrossRef]
  17. Heft, A.I.; Indinger, T.; Adams, N.A. Introduction of a New Realistic Generic Car Model for Aerodynamic Investigations; Technical Report, SAE Technical Paper; SAE: Warrendale, PA, USA, 2012. [Google Scholar]
  18. Heft, A.I.; Indinger, T.; Adams, N.A. Experimental and numerical investigation of the DrivAer model. In ASME 2012 Fluids Engineering Division Summer Meeting Collocated with the ASME 2012 Heat Transfer Summer Conference and the ASME 2012 10th International Conference on Nanochannels, Microchannels, and Minichannels; American Society of Mechanical Engineers: New York, NY, USA, 2012; pp. 41–51. [Google Scholar]
  19. Peters, B.C.; Uddin, M.; Bain, J.; Curley, A.; Henry, M. Simulating DrivAer with Structured Finite Difference Overset Grids; Technical Report, SAE Technical Paper; SAE: Warrendale, PA, USA, 2015. [Google Scholar]
  20. Guilmineau, E. Numerical simulations of flow around a realistic generic car model. SAE Int. J. Passeng. Cars Mech. Syst. 2014, 7, 646–653. [Google Scholar] [CrossRef]
  21. Jakirlic, S.; Kutej, L.; Basara, B.; Tropea, C. Computational study of the aerodynamics of a realistic car model by means of RANS and hybrid RANS/LES approaches. SAE Int. J. Passeng. Cars Mech. Syst. 2014, 7, 559–574. [Google Scholar] [CrossRef]
  22. Fu, C.; Uddin, M.; Robinson, C.; Guzman, A.; Bailey, D. Turbulence Models and Model Closure Coefficients Sensitivity of NASCAR Racecar RANS CFD Aerodynamic Predictions. SAE Int. J. Passeng. Cars Mech. Syst. 2017, 10, 330–344. [Google Scholar] [CrossRef]
  23. Argyropoulos, C.; Markatos, N. Recent advances on the numerical modelling of turbulent flows. Appl. Math. Model. 2015, 39, 693–732. [Google Scholar] [CrossRef]
  24. Shih, T.H.; Liou, W.W.; Shabbir, A.; Yang, Z.; Zhu, J. A new k-ϵ eddy viscosity model for high Reynolds number turbulent flows. Comput. Fluids 1995, 24, 227–238. [Google Scholar] [CrossRef]
  25. Abe, K.; Kondoh, T.; Nagano, Y. A new turbulence model for predicting fluid flow and heat transfer in separating and reattaching flows—I. Flow field calculations. Int. J. Heat Mass Transf. 1994, 37, 139–151. [Google Scholar] [CrossRef]
  26. Durbin, P.A. Near-wall turbulence closure modeling without “damping functions”. Theor. Comput. Fluid Dynam. 1991, 3, 1–13. [Google Scholar]
  27. Menter, F.R. Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  28. Menter, F.R. Review of the shear-stress transport turbulence model experience from an industrial perspective. Int. J. Comput. Fluid Dynam. 2009, 23, 305–316. [Google Scholar] [CrossRef]
  29. Gritskevich, M.S.; Garbaruk, A.V.; Schütze, J.; Menter, F.R. Development of DDES and IDDES formulations for the k-ω shear stress transport model. Flow Turbul. Combust. 2012, 88, 431–449. [Google Scholar] [CrossRef]
  30. Jones, W.P.; Launder, B.E. The prediction of laminarization with a two-equation model of turbulence. Int. J. Heat Mass Transf. 1972, 15, 301–314. [Google Scholar] [CrossRef]
  31. Launder, B.E.; Sharma, B.I. Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc. Lett. Heat Mass Transf. 1974, 1, 131–137. [Google Scholar] [CrossRef]
  32. Durbin, P. On the k-ϵ stagnation point anomaly. Int. J. Heat Fluid Flow 1996, 17, 89–90. [Google Scholar] [CrossRef]
  33. Launder, B.E.; Kato, M. Modeling flow-induced oscillations in turbulent flow around a square cylinder. In Proceedings of the ASME Fluid Engineering Conference, Washington, DC, USA, 20–24 June 1993; Volume 20. [Google Scholar]
  34. Wolfshtein, M. The velocity and temperature distribution in one-dimensional flow with turbulence augmentation and pressure gradient. Int. J. Heat Mass Transf. 1969, 12, 301–318. [Google Scholar] [CrossRef]
  35. Davidson, L.; Nielsen, P.V.; Sveningsson, A. Modifications of the V2 Model for Computing the Flow in a 3D Wall Jet. In Proceedings of the International Symposium on Turbulence, Heat and Mass Transfer, Antalya, Turkey, 12–17 October 2003. [Google Scholar]
  36. Giannuzzi, M. Low Reynolds turbulence model CFD simulation for complex electronic system: An industrial point of view. J. Phys. Conf. Ser. 2014, 525, 012002. [Google Scholar] [CrossRef]
  37. Durbin, P.A. Separated flow computations with the k-epsilon-v-squared model. AIAA J. 1995, 33, 659–664. [Google Scholar] [CrossRef]
  38. Wilcox, D.C. Reassessment of the scale-determining equation for advanced turbulence models. AIAA J. 1988, 26, 1299–1310. [Google Scholar] [CrossRef]
  39. Wilcox, D.C. Formulation of the kw turbulence model revisited. AIAA J. 2008, 46, 2823–2838. [Google Scholar] [CrossRef]
  40. Spalart, P.R. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In Proceedings of the First AFOSR International Conference on DNS/LES; Greyden Press: Dayton, OH, USA, 1997. [Google Scholar]
  41. Menter, F.R.; Kuntz, M.; Langtry, R. Ten years of industrial experience with the SST turbulence model. Turbul. Heat Mass Transf. 2003, 4, 625–632. [Google Scholar]
  42. Zhang, C.; Uddin, M.; Song, X.; Fu, C.; Foster, L. Simultaneous Improvement of Vehicle Under-Hood Airflow and Cooling Drag Using 3D CFD Simulation. In Proceedings of the SAE 2016 World Congress and Exhibition, Detroit, MI, USA, 12–14 April 2016; SAE International: Warrendale, PA, USA, 2016. [Google Scholar] [CrossRef]
  43. Howard, R.; Pourquie, M. Large eddy simulation of an Ahmed reference model. J. Turbul. 2002, 3, 1–18. [Google Scholar] [CrossRef]
Figure 1. Full vehicle model.
Figure 1. Full vehicle model.
Fluids 04 00148 g001
Figure 2. Inclined Condenser, Radiator, and Fan Module (CRFM) and air-duct in Case 2; (a) side view, (b) top view.
Figure 2. Inclined Condenser, Radiator, and Fan Module (CRFM) and air-duct in Case 2; (a) side view, (b) top view.
Fluids 04 00148 g002
Figure 3. Front grille opening; (a) front grille with no side grille sealing as used in Case 1, (b) the front grille with both side grilles sealed off as used in Case 2. The blue colored parts indicate how the side grilles were sealed off.
Figure 3. Front grille opening; (a) front grille with no side grille sealing as used in Case 1, (b) the front grille with both side grilles sealed off as used in Case 2. The blue colored parts indicate how the side grilles were sealed off.
Fluids 04 00148 g003
Figure 4. The computational domain.
Figure 4. The computational domain.
Fluids 04 00148 g004
Figure 5. The mesh on the central plane of the vehicle.
Figure 5. The mesh on the central plane of the vehicle.
Fluids 04 00148 g005
Figure 6. The wall y + values on the external surfaces of the vehicle.
Figure 6. The wall y + values on the external surfaces of the vehicle.
Fluids 04 00148 g006
Figure 7. C D for both cases from the predictions of all turbulence models and the experiment.
Figure 7. C D for both cases from the predictions of all turbulence models and the experiment.
Fluids 04 00148 g007
Figure 8. C L for both cases from the predictions of all turbulence models.
Figure 8. C L for both cases from the predictions of all turbulence models.
Fluids 04 00148 g008
Figure 9. C L / C D for both cases from the predictions of all turbulence models.
Figure 9. C L / C D for both cases from the predictions of all turbulence models.
Fluids 04 00148 g009
Figure 10. Accumulated drag coefficient along the streamwise direction for the baseline case obtained from all turbulence models (L represents the vehicle total length).
Figure 10. Accumulated drag coefficient along the streamwise direction for the baseline case obtained from all turbulence models (L represents the vehicle total length).
Fluids 04 00148 g010
Figure 11. Delta accumulated drag coefficient between Case 1 and 2 vehicle configurations along the streamwise direction obtained from all turbulence models (L represents the vehicle total length).
Figure 11. Delta accumulated drag coefficient between Case 1 and 2 vehicle configurations along the streamwise direction obtained from all turbulence models (L represents the vehicle total length).
Fluids 04 00148 g011
Figure 12. The pressure coefficient from the RKE model on the back of the vehicle for the baseline case.
Figure 12. The pressure coefficient from the RKE model on the back of the vehicle for the baseline case.
Fluids 04 00148 g012
Figure 13. Vector scene with constrained streamlines in the wake region (on the central plane of the vehicle) for the baseline case from all models studied.
Figure 13. Vector scene with constrained streamlines in the wake region (on the central plane of the vehicle) for the baseline case from all models studied.
Fluids 04 00148 g013
Figure 14. The delta pressure coefficient against the RKE model for all other models for the baseline case on the back of the vehicle.
Figure 14. The delta pressure coefficient against the RKE model for all other models for the baseline case on the back of the vehicle.
Fluids 04 00148 g014
Figure 15. Four selected lines (on the central plane of the vehicle) for quantitative analysis inside the boundary layer. Note that the four red lines only represent the location, not the actual length.
Figure 15. Four selected lines (on the central plane of the vehicle) for quantitative analysis inside the boundary layer. Note that the four red lines only represent the location, not the actual length.
Fluids 04 00148 g015
Figure 16. Streamwise velocity profiles for the baseline case at four selected locations on the central plane of the vehicle.
Figure 16. Streamwise velocity profiles for the baseline case at four selected locations on the central plane of the vehicle.
Fluids 04 00148 g016
Figure 17. Locations of selected plane sections.
Figure 17. Locations of selected plane sections.
Fluids 04 00148 g017
Figure 18. Streamwise velocity on selected plane sections from all models for the baseline case.
Figure 18. Streamwise velocity on selected plane sections from all models for the baseline case.
Fluids 04 00148 g018
Figure 19. Delta streamwise velocity on selected plane sections against the IDDES fine for all other models for the baseline case.
Figure 19. Delta streamwise velocity on selected plane sections against the IDDES fine for all other models for the baseline case.
Fluids 04 00148 g019
Figure 20. Delta pressure coefficient ( C P of Case 2 subtracted from C P of Case 1) distributions on the external surfaces of the vehicle.
Figure 20. Delta pressure coefficient ( C P of Case 2 subtracted from C P of Case 1) distributions on the external surfaces of the vehicle.
Fluids 04 00148 g020
Figure 21. Delta streamwise velocity (velocity field in Case 2 relative to that of Case 1) on selected plane sections from all models.
Figure 21. Delta streamwise velocity (velocity field in Case 2 relative to that of Case 1) on selected plane sections from all models.
Fluids 04 00148 g021
Table 1. Percent deviation of C D predictions from the experimental data for the baseline vehicle model. RKE, Realizable k ε ; AKN, Abe–Kondoh–Nagano; IDDES, Improved Delayed Detached Eddy Simulation.
Table 1. Percent deviation of C D predictions from the experimental data for the baseline vehicle model. RKE, Realizable k ε ; AKN, Abe–Kondoh–Nagano; IDDES, Improved Delayed Detached Eddy Simulation.
% C D DeviationRKEAKNV2FSSTDDESIDDESIDDES Fine
Case 11.3%13.8%10.0%12.9%11.5%10.9%10.8%
Table 2. Drag and lift coefficients for all models.
Table 2. Drag and lift coefficients for all models.
Case #Model C D C L C L / C D C LF C LR %Front
Case 1RKE0.3200.0980.305−0.0130.111−13.2%
AKN0.3600.0960.267−0.0350.132−36.5%
V2F0.3480.1290.371−0.0150.144−11.4%
SST0.3570.0960.269−0.0330.130−34.6%
DDES0.3520.0840.238−0.0410.126−49.5%
IDDES0.3500.0600.171−0.0450.106−75.4%
IDDES Fine0.3500.0720.206−0.0580.135−80.8%
Case 2RKE0.2900.0730.251−0.0210.122−29.1%
AKN0.3310.0500.150−0.0460.126−92.6%
V2F0.3220.0650.201−0.0420.136−65.7%
SST0.3370.0860.255−0.0370.152−43.0%
DDES0.3500.0800.226−0.0770.157−97.7%
IDDES0.3470.0880.253−0.0720.160−82.1%
Table 3. Delta drag and lift coefficients between Case 1 and Case 2 for all models.
Table 3. Delta drag and lift coefficients between Case 1 and Case 2 for all models.
Model Δ C D Δ C L Δ C L / C D Δ C LF Δ C LR Δ %Front
RKE0.0300.0250.0550.008−0.01115.9%
AKN0.0290.0460.1160.0110.00656.1%
V2F0.0260.0650.1710.0280.00754.3%
SST0.0200.0100.0140.004−0.0228.4%
DDES0.0020.0050.0120.036−0.03048.2%
IDDES0.003−0.028−0.0820.027−0.0536.7%

Share and Cite

MDPI and ACS Style

Zhang, C.; Bounds, C.P.; Foster, L.; Uddin, M. Turbulence Modeling Effects on the CFD Predictions of Flow over a Detailed Full-Scale Sedan Vehicle. Fluids 2019, 4, 148. https://doi.org/10.3390/fluids4030148

AMA Style

Zhang C, Bounds CP, Foster L, Uddin M. Turbulence Modeling Effects on the CFD Predictions of Flow over a Detailed Full-Scale Sedan Vehicle. Fluids. 2019; 4(3):148. https://doi.org/10.3390/fluids4030148

Chicago/Turabian Style

Zhang, Chunhui, Charles Patrick Bounds, Lee Foster, and Mesbah Uddin. 2019. "Turbulence Modeling Effects on the CFD Predictions of Flow over a Detailed Full-Scale Sedan Vehicle" Fluids 4, no. 3: 148. https://doi.org/10.3390/fluids4030148

APA Style

Zhang, C., Bounds, C. P., Foster, L., & Uddin, M. (2019). Turbulence Modeling Effects on the CFD Predictions of Flow over a Detailed Full-Scale Sedan Vehicle. Fluids, 4(3), 148. https://doi.org/10.3390/fluids4030148

Article Metrics

Back to TopTop