Thursday 31 January 2019

Automotive Computational Fluid Dynamics (CFD) Analysis

     This post is about the numerical analysis of an Ahmed body. It was a new experience because of the area between the car's floor and the road which is different from most of the numerical analysis performed in open channel aeronautics and turbo-machinery.

     The numerical analysis was performed using the commercial software, SolidWorks Flow Simulation. The software employs κ − ε turbulence model with damping functions, SIMPLE-R (modified), as the numerical algorithm and second-order upwind and central approximations as the spatial discretization schemes for the convective fluxes and diffusive terms. The time derivatives are approximated with an implicit first-order Euler scheme. Flow simulation solves the Navier–Stokes equations, which are formulations of mass, momentum, and energy conservation laws for fluid flows. To predict turbulent flows, the Favre-averaged Navier–Stokes equations are used.

     The software generates Cartesian mesh using immersed boundary method. The mesh had a cell size of 0.035 m in the far field regions within the computational domain. A fine mesh was need between the road the the car's floor to make sure the interaction of the car's floor with the road was captured accurately. Therefore the mesh between the car's floor and the road was refined to have a cell size of 0.00875 m. Another mesh control was applied around the body to refine the mesh with a cells size of 0.0175 m to capture the trailing vortices. The resulting mesh had 209,580 total cells, among those cells, 31,783 cells were at the solid fluid boundary. The computational domain size was ~1L x 1.12L x 3L where L being the vehicle's length. The computational domain along with the computational mesh is shown in Fig. 1.



Fig. 1 Mesh, computational domain and the boundary conditions.

     The red arrows within the Fig. 1 represents the inlet boundary condition of ambient (free-stream) velocity and the blue arrows represent the outlet boundary condition of the ambient pressure. The green arrows represents the co-ordinates axes direction.

     The results from the numerical analysis were compared with [1-3]. The results are within 10% of the experimental results. The velocity (superimposed by the velocity streamlines) and pressure profiles around the car body at various free-stream velocities is shown in Fig. 2.


Fig. 3 Velocity and pressure plots. From the top, Row 1, L-R; ambient velocity of 30 and 40 m/s. Row 2, L-R; ambient velocity of 60 and 80 m/s. Row 3, ambient velocity of 105 m/s.

     It was a good experience learning about automotive CFD after spending a long time in aeronautic/turbo-machinery CFD. Thank you for reading. Please share my work. If you would like to collaborate on a project please reach out.


[1] F.J.Bello-Millán, T.Mäkelä, L.Parras, C.delPino, C.Ferrera, "Experimental study on Ahmed's body drag coefficient for different yaw angles", Journal of Wind Engineering and Industrial Aerodynamics, Volume 157, October 2016, Pages 140-144.

[2] Guilmineau E., Deng G.B., Queutey P., Visonneau M. (2018) Assessment of Hybrid LES Formulations for Flow Simulation Around the Ahmed Body. In: Deville M. et al. (eds) Turbulence and Interactions. TI 2015. Notes on Numerical Fluid Mechanics and Multidisciplinary Design, vol 135. Springer, Cham.

[3] A. Thacker, S.Aubrun, A.Leroy, P.Devinant, "Effects of suppressing the 3D separation on the rear slant on the flow structures around an Ahmed body", Journal of Wind Engineering and Industrial Aerodynamics, Volumes 107–108, August–September 2012, Pages 237-243.

Monday 7 January 2019

Vertical Axis Wind Turbine Computational Fluid Dynamics Analysis

     This post is be about the validation and verification of the computational fluid dynamics analysis of a three blade vertical axis wind turbine. The turbine had a diameter of 2 m with each blade being 1 m tall. The blades had an NACA-0018 airfoil cross section.

     The computational fluid dynamics analysis employed the κ-ε turbulence model with damping functions as the turbulence model, SIMPLE-R as the numerical algorithm. The spatial discretization schemes for the convective fluxes and diffusive terms used are the second order upwind and central approximations, respectively. An implicit first-order Euler scheme is employed to approximate the time derivatives.

     The Cartesian computational mesh with immersed boundary method had a total of 769,357 cells. Among those 769,357 cells, 166,188 cells were around the turbine blades. Mesh controls were employed to refine the mesh near the turbine blades. A time step of 3e-3 was employed. The computational domain inlet was 1.5 D away from the turbine and the outlet was 3D away. The computational domain walls on the sides were 1D x 1.5D, where D represents the turbine diameter. The mesh and the computational domain are shown in Fig. 1. The vertical teal arrow represents the force of gravity, the curved teal arrow represents the direction of turbine rotation. The dark blue arrow represents the direction of free stream velocity.

Fig. 1, Mesh and computational domain.

     The simulations ran at a tip-speed ratio of 1.87 at a wind speed of 4.03 m.s-1. The velocity distribution around the turbine after 4 revolutions is shown in Fig. 2. Validation of the numerical analysis was carried out using [1]. The results of power produced by the turbine were with in 4% of the experimental results [1]. An animation of the numerical analysis is also shown.

Fig. 1, Flow field around the turbine.

     Thank you for reading. If you would like to contribute to the research, both financially and scientifically, please feel free to reach out.





[1] Yi-Xin Peng, You-Lin Xu, Sheng Zhan and Kei-Man ShumHigh-solidity straight-bladed vertical axis wind turbine: Aerodynamic force measurements, Journal of Wind Engineering and Industrial Aerodynamics, January 2019.