DriveAer

Download mesh
Download simulation
Link to youtube video

Introduction

intro

DrivAer is a full-scale vehicle model for exterior design and aerodynamic testing used in the field of automotive engineering to simulate and evaluate the exterior shape and aerodynamic characteristics of vehicles. The model was introduced to bridge the gap between simplified models and highly complex production vehicles, and CAD files and experimental results for many types of geometries are publicly available.

This example uses a fastback model with side mirrors and wheels.

https://www.epc.ed.tum.de/en/aer/research-groups/automotive/drivaer

Use steady-state incompressible flow conditions with moving ground and rotating wheel.

The experimental results in the paper show that drag coefficient(Cd) is 0.247 for ASME and 0.243 for SAE, and the simulation results show that Cd = 0.243, which is consistent with the SAE experimental results.

The simulation conditions are as follows

  • solver : buoyantSimpleNFoam
  • turbulence model : $Realizable$ $k-epsilon$ model
  • density : 1.205 $kg/m^3$
  • viscosity : 1.82e-5 $kg/ms$
  • flow condition : 30 $m/s$ at inlet

Start BaramFlow and load mesh

Run the program and select [New Case] from the launcher. In the launcher, select [Pressure-based] for [Solver Type] and [None] for [Multiphase Model].

Use the given polyMesh folder. In the top tab, click [File]-[Load Mesh]-[OpenFOAM] in that order and select the polyMesh folder.

General

For this example, we’ll use default conditions.

Models

For this example, we’ll use $Realizable$ $k-epsilon$ model for turbulence.

Materials

For this example, we will use the properties of air.

  • air
    • Density : 1.205 $kg/m^3$
    • Viscosity : 1.82e-5 $kg/ms$

Boundary Conditions

Each boundary condition is set as follows

  • minX : velocity Inlet
    • Velocity Specfication Method : Magnitude, Normal to Boundary
    • Profile Type : Constant
    • Velocity Magnitude : 30 (m/s)
    • Turbulent Intensity : 1 (%)
    • Turbulent Viscosity Ratio : 10

  • maxX : Pressure Outlet
    • Pressure : 0 (Pa)

  • minY, maxY, maxZ : Symemtry
  • minZ : Wall
    • Velocity Condition : Translational Moving Wall
    • Velocity : (30, 0, 0)

  • body_no_wheel : Wall
    • Velocity Condition : No Slip

  • Wheels_Front_Smooth : Wall
    • Velocity Condition : Rotational Moving Wall
    • Speed (RPM) : 898.8
    • Rotation-Axis Origin : (0.007, 0, 0)
    • Rotation-Axis Direction : (0, -1, 0)

  • Wheels_Rear_Smooth : Wall
    • Velocity Condition : Rotational Moving Wall
    • Speed (RPM) : 898.8
    • Rotation-Axis Origin : (2.7932, 0, 0)
    • Rotation-Axis Direction : (0, -1, 0)

Reference Values

Set the Reference Value for the aerodynamic coefficient calculation as follows.

  • Area : 1.08
  • Density : 1.205
  • Length : 4.6132
  • Pressure : 0
  • Velocity : 30

Numerical Conditions

In this example, we’ll change the settings as shown below.

  • Pressure-Velocity Coupling Scheme : SIMPLEC
  • Discretization Scheme
    • Pressure : Momentum Weighted Reconstruct
    • Momentum : Second Order Upwind
    • Turbulence : Second Order Upwind
  • Under-Relaxation Factors
    • Pressure : 0.9
    • Momentum : 0.9
    • Turbulence : 0.9
    • Density : 0.9
  • Convergence Criteria
    • Pressure : 0.00001
    • Momentum : 0.001
    • Turbulence : 0.001

Monitor

Monitor the force coefficients of car.

Select [Monitors]-[Add]-[Forces] and set values as shown below.

Initialization

Change the values as shown below

  • Velocity
    • X-Velocity : 30 (m/s)
    • Y-Velocity : 0 (m/s)
    • Z-Velocity : 0 (m/s)
  • Pressure
    • 0 (Pa)
  • Turbulence
    • Scale of Velocity : 30 (m/s)
    • Turbulent Intensity : 1 (%)
    • Turbulent Viscosity Ratio : 10

Enter the value and click the Initialize button at the bottom. Then click the [File]-[Save] menu to save the case file.

Run

Change the values as shown below, and click [Start Calculation] button.

  • Number of Iterations : 2,000
  • Save Interval : 100
  • Data Write Format : Binary
  • Selct [Parallel]-[Environment] in menu. Set Number of Cores as you want and select [Local Machine] for [Parallel Type].

When the calculation is started, you’ll see a graph of Residuals and Force monitor as shown below.

Post-processing

Scalar distribution at boundary

BARAM uses ParaView for post-processing. To start post-processing, click the ParaView button in [External Tools]. In this example, we will draw the pressure distribution and streamlines.

Change the Case Type to Decomposed Case.

  • In [Mesh Regions] select boundaries
    • Wheels_Front_Smooth, Wheels_Rear_Smooth, body_no_wheels

Utilize the [Extract block] function to extract the geometry of the vehicle and floor as shown below.

Change [Solid Color] to p_rgh.

Streamline

Draw the streamline of the flow around the vehicle.

Change p to [Solid Color].

In the Pipeline Browser on the left, click baram.foam once to activate it.

Click [Stream Tracer] icon.

Change settings as follows

  • Seed Type : Point Cloud
  • Center : (-1.5, 0.1, 0.1)
  • Radius : 0.2
  • Number of Points : 100
  • Coloring : U

You can see streamline as shown below.