Project Steady state simulation of flow over a throttle body
- Case setup : CONVERGE CFD STUDIO
- Commanad editor : Cygwin
- Post-processor : Paraview
GOVERNING EQUATIONS FOR THE FLOW THROUGH GEOMETRY (NAVIER-STOKES EQUATION)
The flow through the geometry is governed by three physical rules.
- Conservation of mass.
∂(ρu)/∂x = 0 - Newtons second law.
F = ma - Conservation of energy.
h1+ (u1^2)/2= h2+ (u2^2)/2
The air traveling past the throttle body is assumed to be an ideal gas, i.e. the ideal gas law can be utilized p = ρRT
THE CASE SETUP
- The case setup was done using the converge software. The geometry was imported to the converge and the was cleaned to check if there were any voids or gap in it.
- The diagnosis was done and all the errors were solved
- Boundary flagging was done and the geometry was divided into 4 different regions, namely inlet, outlet, elbow wall and the throttle-body.
- Boundary fencing technique was made used to group or flag the boundaries into different regions.
- The body has an inlet pressure of 150k Pa and an outlet pressure of 100k Pa.
- And further in the case setup the condtions such as law of wall were selected.
- A desired mesh of 2mm was assigned in all three dimensions of the body, ie,,, in other words around 5000 cells or elements were created in the body along each of their length in the respective axis (dx, dy, dz)
- The simulation that we are intrested here is a steady-state simulation where we are intrested only in the final results and not what is happening to the flow wrt time
- Turbulence model – RNG K-Ɛ
The Simulation parameters
- The simulation was ran for a maximum time step of 1sec, starting with a minimum time interval of 1e-9 for around 15000 cycles untill the flow in the throttle body reaches the steady state.
- The caluclation is done in such a way that, for every 150 cycles of the simulation running a 3d output file is getting dumped or printed to the output folder specified.
- Hence a total of 100 3d output files were generated to the output folder specified.
Post Processing the 3d output files
- So when all the 3d output files are getting dumped into the output files by the cygwin, it is important to convert it into format that the paraview can understand like .case format or .vtm format using cygwin.
- Once the post processed files are ready paraview can read those files and the desired contour of the different parameters can be obtained and the animation can be done.
Surface mesh with edges
Surface with elements
Fixed embedding layer around the throttle to capture data even more accurately
- The given geometry has been discritized into finite number of elements or cells inorder to visualize the behaviour of the flow through each and every element.
- If the meshing is done properly then it will resutl in the non-linear or too much deviation in the solution and the solution will not converge to the exact solution.
- Since the geometry is not too complicated the meshing is done in converge itself. There are about 35000 cells or elements created in the given geometric model.
- Around 3 layer of fixed embedding is added around the region of throttle body to make sure the mesh around it is generated quite uniformly, like atleast to have some cells around the throttle, as the throttle needs to rotate the clearence between the elbow (wall) and the throttle is very less.
PLOTS OF PRESSURE, VELOCITY, CELL-COUNT AND MASS-FLOW RATE
PRESSURE
The variation of the pressure is as shown in the plot. For a flow to occur there should be a pressure diffrence developed. Hence an inlet pressure of 150k Pa and an oultet pressure of 100k Pa is provided and it is this pressure diffrence the makes the flow to happen.
A higher pressure is experienced when the flow stream meets the valve surfaces.
When the flow starts to accelerate the pressure is low initally and starts increasing as the flow starts accelerating. The increasein pressure is observed till the flow reaches the throttle presesnt inside the body. Once the flow reaches the throttle high fluctuation of pressure is observed if the throttle body is rotating, else the pressure increases linearly and once it reaches the outlet of the body, the flow pressure decreases as at the exit the flow accelerates ie,,,it’s velocity increases, as pressure decreases. After some time the flow reaches thr steady state comdition and pressure remains almost a constant value.
TOTAL CELL-COUNT
The total cell count is a constant value throughout the simulation and it mainly depends on the mesh size that we have given during the time of setting up the case.
As shown in the above fig, there are four different ranks shown. Each rank corresponds to a number of processors used during running the simulation. Since the simulation is ran for 4 processors (to reduce the computational time involved and to lower the loads on only 1 processor), hence there are 4 rank indicated. The main intention of creatin the cells or elements is to visualize the flow through each cells and what is causing the flow to vary in that particular region of the cells created.
There are about 35000 cells created here for running the simulation.
VELOCITY
Due to the pressure difference at the inlet and the outlet the flow starts with zero velocity initially and starts accelerating. Once the flow starts accelerating it’s velocity increases gradually and reaches the maximum value at the exit of the throttle body.
The variation in the velocity of the flow in the above plot is due to it’s interaction with the throttle body as the flow encounters it. Once the flow encounters the throttle body as the body is exposed to a fluid flow stream it’s velocity decreases, because some amount of area of throttle body is exposed to the fluid flow which decelerates the flow.
A Maximum velocity of 200 m/s is achieved in this simulation at the throttle body outlet.
MASS-FLOW RATE
Mass-flow rate involves the terms like density, area and velocity of the flow. It is the quantity of fluid flowing through the given geometry in unit time interval in kg/s.
At the inlet the fluid is just accelerating the velocity of the flow will be lesser and hence the mass-flow rate will also be lesser. As the fluid moves towards the outlet the flow accelerated and it’s velocity increases thereby increasing the flow rate of the fluid passing through the cycles. At the throttle the mass-flow rate decreases and after certain amount of time the flow attains a steady state and the mass-flow rate is almost a constant value.
TKE (TOTAL KINETIC ENERGY)
The total kinetic energy is the amount of kinetic energy per unit mass. The tke increases as the flow accelerates from zero velocity. After
simulation reaches steady state, tke also attains steady state and remains a constant value.
THE CASE SETUP