Week 3 - External flow simulation over an Ahmed body.

AIM: To study the aerodynamics of the Ahmed's body and Run the simulation for the velocity of 25 m/sec with the default air properties in fluent.

OBJECTIVE:

 1. To Describe Ahmed's body and its importance.

 2. To Explain the reason for the negative pressure in the wake region. 

 3. To Explain the significance of the point of separation.

 

1.Ahmed's body and its importance:

The Ahmed body represents a simplified, ground vehicle geometry of a bluff body type. Its shape is simple enough to allow for accurate flow simulation but retains some important practical features relevant to automobile bodies. This model describes how to calculate the turbulent flow field around a simple car-like geometry.

The airflow around the Ahmed body captures the essential flow features around an automobile and was first defined and characterized in the experimental work of Ahmed. Although it has a very simple shape, the Ahmed body allows us to capture characteristic features that are relevant to bodies in the automobile industry.

ANSYS is a general-purpose FEA and CFD simulation software. It gives a numerical prediction of the behaviour of a system. These numerical predictions need to be validated by experimental data.

This model is also used to describe the turbulent flow field around a car-like geometry. People working in vehicle aerodynamics use Ahmed body to validate their numerical model because the experimental data from wind tunnel testing is available for Ahmed body. Once the numerical model is validated, it is used to design new models of the car.

The Ahmed body was described originally by S. R. Ahmed in 1984. Three main features were seen in the wake:

1. The A recirculation region that is formed as the flow separates at the top of the vertical back surface of the model

2. The B recirculation region that is formed due to the separation at the base of the model.

3. The c-pillar vortices that form as the vorticity in the side boundary layers roll up over the slant edges.

The wake was shown to be highly dependent on slant angle. For slant angles less than 12°, the flow remains attached over the slant. The flow is essentially two-dimensional and has low drag. Between 12° and 30° the flow becomes much more three-dimensional as the c-pillar vortices form. These reach maximum strength at 30°. The drag increases significantly as the low pressure cores act on the rear surfaces. Past 30° the flow separates fully off the slant. This results in a sudden decrease in drag and weaker c-pillar vortices.

The Geometry of Ahmed's Body is as shown in the below figure.

 

 

 2. Explain the reason for the negative pressure in the wake region.

The negative pressure is due to the velocity difference applied at a certain Reynold number and the rear angle, which at a low Reynolds number when this difference in pressure happens the wake obtained increases the drag due to the pressure drop behind the Ahmed body, which at a high Reynolds number due to the vortices generated it produces small amounts of pressure thrusting the body.

The air becomes almost stagnant as it strikes Ahmed body which results in air exerting very high pressure on the front face of Ahmed body. The airflow then gets divided between the upper and lower surface of the Ahmed body. The higher pressure air on the front surface accelerates as it travels over the curved nose surface of Ahmed body, causing the pressure to drop. This lower pressure creates lifts over the roof surface as the air passes over it. As the air continues to flow and make its way to the rear, a notch is created by the rear slant owing to flow separation, leaving a vacuum or low-pressure space which the air is not able to fill properly. The resulting lower pressure creates lift and the drag that then acts upon the surface area of the rear slant.

The following figure shows the pressure thrusting.

 

 

3. The significance of the point of separation.

Flow separation or boundary layer separation is the detachment of a boundary layer from a surface into a wake. Separation occurs in flow that is slowing down, with pressure increasing, after passing the thickest part of a streamline body or passing through a widening passage. 

The point of separation defined by the Ahmed body and the flow, the alpha angle depending on which value is at it causes a change in the flow direction, forces required, lift on the body, induced drag due to “tip vortex”, at a critical value flow in slant stalls and reduced lift-reducing as well drag, all depending on the value of alfa and the Reynolds number.

 Procedure.

• Geometry: 
• First, we need to set up the geometry of Ahmed's Body for our simulation.

 

 

 

                 

                                                Isometric view of the three-dimensional model of Ahmed body.

 

 

                                          

                           

                                                                       Front view of the Ahmed body.

 

 

                            

                                                                           Side view of Ahmed body.

 

                  

                                                                     Rear view of Ahmed body.

 

Creating the enclosures around the Ahmed Body for baseline mesh generation.

                  

Outer enclosure

• The outer enclosure is at a distance of 2m from the Front side of the Ahmed body in the X-axis.
• The outer enclosure is at a distance of 5m from the Rear side of the Ahmed body in the X-axis.
• From the top, it is at a distance of 1m in the Y-axis.
• At the bottom, it is merged means it is at zero distance.
• From the left and right side of the Ahmed body, it is 0.5m along Z-axis.

Inner enclosure

• The inner enclosure is at a distance of 0.5m from the Front side of the Ahmed body in the X-axis.
• The inner enclosure is at a distance of 1m from the Rear side of the Ahmed body in the X-axis.
• From the top, it is at a distance of 0.5m in the Y-axis.
• At the bottom, it is merged means it is at zero distance.
• From the left and right side of the Ahmed body, it is 0.26m along Z-axis.

 

Interference region.

 

Initially, there was a huge interference between the Outer and Inner encloser so that the mesh gets overlapped which gets not acceptable so to avoid this we made share topology and get deleted the Interference region that by using the interference command after this we get a good interaction between the Inner and Outer enclosers so that we can get the mesh properly.

 

MESH:

• Generating the mesh on the enclosures around Ahmed Body, the outer enclosure has Hexahedral grids to capture the flow structure and the inner enclosure has the Tetrahedral grids and inflation layers has created along the walls of the Ahmed body to capture the point of separation at the wake region clearly.
• By using the named selection command we are assigning the boundaries that are front side of the Ahmed body as inlet, Rear side as outlet and also the top, bottom, left and right sides as symmetry condition and Ahmed body as car-wall.

    SETUP:

      In the setup, we are going to set up the conditions for the simulation.

• in the setup>physics>steady state and density-based solver.
• setup>physics>viscous model>k-epsilon>ok.
• setup>physics>material, select air as a flow material>ok
• setup>physics>boundaries>inlet  -velocity as 25 m/s, outlet gauge pressure as 0, other conditions will be the default values we are considering>click on ok
• setup>solution>report>residuals, uncheck all the convergence criteria>ok
• setup>solution>report>definitions>new>force report>drag>drag coefficients>car-wall>select report file and report plot >save/display
• setup>solution>method>Hybrid>initialize
• Result>create>plane>point and normal>ok
• setup>solution>run calculation

   The Boundary conditions.

    The inlet velocity will be 25 m/s, the density of the fluid will be 1.225 kg/m^3, dynamic viscosity will be 1.7894x10^-5 kg/m-s.

SOLUTION:

• In the solution only we will be taking the solved simulation results.
• Residuals plot, drag and lift plots and the reports we are going to get it in the solutions.

Post processing:

• here, we will get the vector plots, and the the velocity and pressure contour planes.

 

RESULTS.

CASE 1.

In the first case the simulation was done with the following cell count.

• Y-plus = 50
• Viscous model: K-epsilon
• 1st box domain =100 mm size of hexahedral mesh
• 2nd box domain = 50 mm size of tetrahedral mesh
• Inflation layer thickness = 0.7 mm
• Elements = 77,052
• Nodes = 25,415

Velocity contour plane.

Pressure contour plane.

Residuals.

 

Drag coefficient plot.

 

Vector plot.

 

 

CASE 2.

In the second case the simulation was done for the following cell count.

• Y-plus = 125
• Viscous model: K-epsilon
• 1st box domain =75 mm size of hexahedral mesh
• 2nd box domain = 35 mm size of tetrahedral mesh
• Inflation layer thickness = 1.79 mm
• Elements = 1,95,799
• Nodes = 56,779

 

Velocity contour plane.

Pressure Contour plane.

Residuals plot.

Drag coefficient plot.

Vector plot.

 

CASE 3.

In the third case the simulation was done for the following cell count.

• Y-plus = 150
• Viscous model: K-epsilon
• 1st box domain =75 mm size of hexahedral mesh
• 2nd box domain = 35 mm size of tetrahedral mesh
• Inflation layer thickness = 1.43 mm
• Elements = 1,96,895
• Nodes = 56,975

Velocity contour plane.

Pressure contour plane.

Residuals plot.

Lift coefficient plot.

 

Drag coefficient plot.

 

CASE 4.

In the third case the simulation was done for the following cell count.

• Y-plus = 5
• Viscous model: K-Omega
• 1st box domain =75 mm size of hexahedral mesh
• 2nd box domain = 35 mm size of tetrahedral mesh
• Inflation layer thickness = 0.0716 mm
• Elements = 2,01,422
• Nodes = 57,752

Velocity contour plane.

Pressure contour plane.

Residuals plot.

Drag coefficient plot.

Lift coefficient plot.

Vector plot.

 

Comparision.

K-epsilon	ELEMENTS	Cd	Cl
CASE 1	77052	0.3622	0.25149
CASE 2	1,95,799	0.3349	0.2780
CASE 3	1,96,895	0.3327	0.2831
K-Omega
CASE 4	201422	0.3006	0.2029

 

Conclusion.

There is a variation of pressure across the Ahmed body. The pressure in front of the body is positive because that's the location where the flow first comes in contact with the body and comes to rest. At the corner, the flow accelerates, thereby decreasing the pressure.

The reason behind negative pressure in the wake region is due to flow separation, recirculating flow. It's a phenomenon that generally happens as the air moving past the body isn't able to pass through the immediate backside of the body creating a region, where there is no movement of air and thereby creating a pressure gradient.

As the wake flow is the main contributor to the drag force(pressure drag), the point of separation is an important aspect of the flow. The more attached the flow is to the body the recirculation zone will be smaller and hence reduce the adverse effect of pressure drag, improving the fuel efficiency of the body.