Posts

F1 CAR SIMULATION

Image
In one of my previous posts, I covered a full body automotive steady state simulation with resolved boundary layers, whose non-dimensionless wall distance y+ was less than 1. In the following simulation, I have used wall function to run a kOmegaSST turbulence model on a boundary layer-less grid. Again, a steady state solver was used to compute the force Coefficients of a formula 1 car. The car model was downloaded from Grab cad:  https://grabcad.com/library   The average value of the drag coefficient was found to be 0.72 . The converged drag plot can be accessed below. Here is a worthwhile observation, while using snappyHexMesh tool to create a mesh which employs wall function, sometimes it is required that at least one boundary layer be extruded in order to fall within the range 30 < y+ < 300. Again, this depends on the nature of the background mesh and the level of refinement used inside the fluid domain. Kindly write to me if you require assistance in the case set-up. Th

WIND TURBINE CFD SIMULATION

Image
I decided to run a wind turbine simulation on OpenFOAM. This is an Arbitrary Mesh Interface setup, run using pimpleDyMFoam. The turbine stem/post has been neglected in the study, in order to ease the simulation process. Wall function was operated on a coarse mesh, which was created using the snappyHexMesh tool and the total number of cells worked up to 3,05,294.  A close eye has to be kept on the usage of the wall function while performing the grid study. Because drastically refining the cells, without creating boundary layers, can cause the Courant number to shoot up while performing the simulation. Three cylindrical geometries were used in the mesh. A small cylinder to house the wind turbine and to form the AMI interface with the intermediate cylinder. A large cylinder to envelop the entire fluid domain. Figure No.1: Pressure Side Figure No.2: Suction Side Below are some of the animations compiled from the simulation.  Velocity distribution along the blade:

AUTOMOTIVE CFD FOR DETERMINING AERODYNAMIC DESIGN COEFFICIENTS

Image
A rather unorthodox concept car geometry was chosen for automotive CFD analysis and the force Coefficients were calculated.  Figure No.1: STL geometry Turbulence properties were calculated from the wheel base length of the car and the Aerodynamic coefficients were calculated from the frontal area of the car. Here are the steps adhered to calculate the frontal area of the car using paraview. Here is a brief video tutorial on how to calculate the projection Area or the frontal area of a car.   The projection area and the wheelbase length values are called inside the force-coefficient file, within the system folder.  Figure No.2: lRef (base length) & Aref (Projection Area) The simulation was performed on a grid whose dimensionless wall-distance value was less than 1.0. A kOmegaSST turbulence model was adopted for a steady-state solver.  Figure No.3: Meshed STL geometry Along with the car, the ground patch or the lowerWall is also considered a

AIRCRAFT OpenFOAM SIMULATION (BOEING 747)

Image
Before I begin to furnish the results, a quick look at the computational resources that were utilized for the simulation. I possess a minimal capability laptop manufactured by Hewlett-Packard, which runs on 2 core(s) and a 4 GB RAM. The objective of the simulation is to set up a 3D mesh for an Aircraft flow simulation with a resolved boundary layer of Y+ < 1, and to study the Residuals . The kOmegaSST turbulence model was adopted. The STL file downloaded online is a scaled Boeing 747 model, with the elevator tabs deflected downwards. Figure No.1: Rear View The cell size in the blockMesh is 0.05, owing to the fact that a scaled 747 model was selected, otherwise a cell size of 0.5 can be executed for an un-scaled model. Regardless, any of the following cell size can be applied (0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0). Four boundary layers were generated with an expansion ratio of 1.2,  adding to the total cell count of over 3 million cells. A steady-state solver

MATLAB SCRIPT TO CALCULATE THE TURBULENCE PROPERTIES AND TO ESTIMATE THE WALL DISTANCE

This script is an adaptation of the original script written by Dr. Nguyen Vinh Tan, Agency for Science, Technology and Research, Singapore.  A matlab script with (.m) extension can be adopted from the following formula and the same can be run using mathwork or octave o n the linux platform. The turbulence properties and mesh layer thicknesses can be predicted using the same.  rho=1000; %This is density, a value 1.225 for air can be assigned nu=1.0e-6; %kinematic viscosity for air/water L=0.5; %L is the wheelbase length of a car/length of a flatplate/diameter for a cylinder simulation % If the reynolds number is assigned and the flow velocity has to be found out use: ReL = 60000; %global reynolds number Va = ReL*nu/L %If the flow velocity is assigned and the Reynolds number has to be calculated use: Va = 30; ReL = Va*L/nu %Turulence intensity can either be assigned or calculated I = 0.02; % 2 percentage(assigned) I = 0.05; % 5 percentage(assigned

snappyHexMesh Tutorial For a Complex Geometry and External Aerodynamics

Image
The geometry chosen for external fluid dynamics simulation is that of an offshore structure. In that, a deep-draft semi-submersible structure  is scaled scaled down to a value  of 0.01 (scaled down by a factor of 100). The semi-submersible structures are employed at locations where the oil extraction from the sea bed is considerably deep. These are movable structures which are half submerged during application. A working deck comes right on top of the structure where the cranes and oil rigs are fitted. The semi-submersible consists of four hull appendages which connect the working deck with a pontoon structure. The pontoon structure gives buoyancy to the semi-submersible in deep waters but does not come directly in contact with the ocean bed.  Figure No.1: DDS The effects of flow induced vibration on the above four columns during high flow currents is indeed the case of interest. Once a good mesh is achieved, the platform is set for the study on the effects of Vortex Induced

CFD ANALYSIS OF FLOW AROUND A CYLINDER (RESULTS)

Image
Here are few results compiled from the Cylinder simulation. I have chosen the following Reynolds numbers for the study: 30, 200, 10000, 1*10^6.  Reynolds Number 30: Figure No.1: Flow Visualization  At a low Reynolds number of 30, the vortex shedding has not begun and so the force Coefficient of lift could not be calculated, as shedding is reason for the lift to exist. Reynolds Number 200: Figure No. 2: Force Coefficient & Strouhal number For a Reynolds number of 200, the vortex shedding takes place and thus producing lift and drag. The force Coefficients are calculated in the above plot as Clrms = 0.3357 Cdmean = 1.2201 St = 0.19. Reynolds Number 3500: Figure No. 3: Force Coefficient & Strouhal number The force coefficients were calculated to be the following: Clrms = 1.014 Cdmean = 0.18 St = 0.19. Reynoldes Number 10000: For Reynolds number of 10^4, a DES based turbulence model was employed and the force coefficients were calculat