Difference between revisions of "OpenFOAM/C3/Turbulent-Flow-in-a-Lid-driven-Cavity/English"

From Script | Spoken-Tutorial
Jump to: navigation, search
 
(6 intermediate revisions by one other user not shown)
Line 62: Line 62:
 
* '''Turbulence modelling'''
 
* '''Turbulence modelling'''
 
* Knowledge of how to solve '''flow''' in a '''Lid driven cavity'''
 
* Knowledge of how to solve '''flow''' in a '''Lid driven cavity'''
* If not so please refer to the relevant tutorial on our website.
+
* If not so please refer to the relevant tutorial on our website
  
 
|-
 
|-
Line 71: Line 71:
  
  
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| This problem is identical in '''geometry and boundary condition''' to the ''''Lid Driven Cavity'''' problem discussed in the basic level tutorial.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| This problem is identical in '''geometry and boundary condition''' to the ''''Lid Driven Cavity'''' problem discussed in the basic level tutorial
  
  
Please make a note this problem is already set up in '''pisoFoam solver''' in '''OpenFoam''' directory.
+
Please make a note this problem is already set up in '''pisoFoam solver''' in '''OpenFoam''' directory
  
  
Line 95: Line 95:
 
Type cd tutorials
 
Type cd tutorials
  
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In the '''terminal''' window type '''run''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In the '''terminal''' window type '''run''' and press '''Enter'''
  
  
Now type '''cd space tutorials''' and press '''Enter'''.
+
Now type '''cd space tutorials''' and press '''Enter'''
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd incompressible
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd incompressible
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now type '''cd space incompressible''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now type '''cd space incompressible''' and press '''Enter'''
  
 
|-
 
|-
Line 116: Line 116:
  
 
Cavity folder inside RAS
 
Cavity folder inside RAS
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now type''' ls''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now type''' ls''' and press '''Enter'''
  
  
Line 122: Line 122:
  
  
Our problem setup is inside '''ras''' folder which is called as '''reynolds averaged stress'''.
+
Our problem setup is inside '''ras''' folder which is called as '''reynolds averaged stress'''
  
 
|-
 
|-
Line 129: Line 129:
  
 
ls
 
ls
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Our folder name is '''cavity'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Our folder name is '''cavity'''
  
  
Now type '''cd space ras''' and press '''Enter'''.
+
Now type '''cd space ras''' and press '''Enter'''
  
  
Now type''' ls''' and press '''Enter'''.
+
Now type''' ls''' and press '''Enter'''
  
 
|-
 
|-
Line 141: Line 141:
  
 
ls
 
ls
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see the '''cavity''' folder. Let me clear this off.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see the '''cavity''' folder. Let me clear this off
  
  
Type '''cd space cavity''' and press '''Enter'''.
+
Now type '''cd space cavity''' and press '''Enter'''
  
  
Now type''' ls''' and press '''Enter'''.
+
Now type''' ls''' and press '''Enter'''
  
 
|-
 
|-
Line 156: Line 156:
  
 
0 folder
 
0 folder
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see three folders ''' 0, constant''' and '''system'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see three folders ''' 0, constant''' and '''system'''
  
  
The '''initial conditions''' are specified within the files in the ''''0'''' directory.
+
The '''initial conditions''' are specified within the files in the ''''0'''' directory
  
  
Let us take a look at the files in the ''''0'''' directory.
+
Now let us take a look at the files in the ''''0'''' directory  
  
 
|-
 
|-
Line 170: Line 170:
 
Type ls
 
Type ls
  
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| To do this, in the '''command terminal''' type '''cd space 0''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| To do this, in the '''command terminal''' type '''cd space 0''' and press '''Enter'''
  
  
Now type '''ls''' and press '''Enter'''.
+
Now type '''ls''' and press '''Enter'''
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"|Point to the files as per narration
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"|Point to the files as per narration
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see files named as '''epsilon, k, nut, nutilda, p, R''' and '''U'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see the files named as '''epsilon, k, nut, nutilda, p, R''' and '''U'''
  
  
These files are to be kept as default until the '''inlet parameters''' don't change.
+
These files are to be kept as default until the '''inlet parameters''' don't change
  
  
If any changes are to be done please refer to the tutorial on '''Simulating flow in a channel''' using OpenFoam, to calculate these values.
+
If any changes are to be done please refer to the tutorial on '''Simulating flow in a channel''' using OpenFoam, to calculate these values
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd..
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd..
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now type '''cd space dot dot''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now type '''cd space dot dot''' and press '''Enter'''
  
  
Let me clear this off.
+
Let me clear this off
  
 
|-
 
|-
Line 195: Line 195:
  
 
Type ls
 
Type ls
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Let us open the '''constant''' folder.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Let us open the '''constant''' folder
  
  
To do this type '''cd space constant''' and press '''Enter'''.
+
To do this type '''cd space constant''' and press '''Enter'''
  
  
Now type '''ls''' and press '''Enter'''.
+
Now type '''ls''' and press '''Enter'''
  
 
|-
 
|-
Line 208: Line 208:
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In this you will see the '''polyMesh''' folder containing  
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In this you will see the '''polyMesh''' folder containing  
 
*the geometry of the '''case''' inside '''blockMeshDict'''
 
*the geometry of the '''case''' inside '''blockMeshDict'''
*and the '''fluid properties'''.
+
*and the '''fluid properties'''
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Point to the two files as per narration
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Point to the two files as per narration
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In this '''case''' you will see two more files other than '''transportProperties''' named as '''RASProperties''' and''' turbulenceProperties'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In this '''case''' you will see two more files other than '''transportProperties''' named as '''RASProperties''' and''' turbulenceProperties'''
  
  
Let us open these two '''files'''.
+
Let us open these two '''files'''
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| RASProperties
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| RASProperties
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In the '''terminal''' type''' gedit (space) RASProperties''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In the '''terminal''' type''' gedit (space) RASProperties''' and press '''Enter'''
  
  
Scroll down.
+
Scroll down
  
  
Line 228: Line 228:
  
  
Close this.
+
Close this
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| turbulentProperties
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| turbulentProperties
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now in the '''command terminal''', type '''gedit space turbulentproperties''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now in the '''command terminal''', type '''gedit space turbulentproperties''' and press '''Enter'''
  
  
Scroll down.
+
Scroll down
  
  
'''The simulation type model''' for this '''case''' is kept as '''RASModel'''.
+
'''The simulation type model''' for this '''case''' is kept as '''RASModel'''
  
  
Close this.
+
Close this
  
 
|-
 
|-
Line 248: Line 248:
  
 
Change the value of viscosity
 
Change the value of viscosity
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now let us open the '''transportProperties model'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now let us open the '''transportProperties model'''
  
  
To do this, in the '''terminal''' type '''gedit space transportProperties '''and press '''Enter'''.
+
To do this, in the '''terminal''' type '''gedit space transportProperties '''and press '''Enter'''
  
  
Line 257: Line 257:
  
  
Close this.
+
Close this
  
 
|-
 
|-
Line 264: Line 264:
  
 
The system folder is to be kept default
 
The system folder is to be kept default
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| We are not changing the '''geometry''' in this '''case'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| We are not changing the '''geometry''' in this '''case'''
  
  
So we need not go inside the '''polyMesh '''folder and look at the '''blockMeshDict''' file.
+
So we need not go inside the '''polyMesh '''folder and look at the '''blockMeshDict''' file
  
  
It can be kept as it is.
+
It can be kept as it is
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd..
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd..
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In the '''terminal''' type '''cd space dot dot''' and press '''Enter'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| In the '''terminal''' type '''cd space dot dot''' and press '''Enter'''
  
  
We will keep the '''system''' folder default as there are no changes inside it.
+
We will keep the '''system''' folder default as there are no changes inside it
  
 
|-
 
|-
Line 287: Line 287:
  
 
Meshing is done
 
Meshing is done
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now, we are done with the setup.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now, we are done with the setup
  
  
We can '''mesh''' the '''geometry'''.
+
We can '''mesh''' the '''geometry'''
  
  
To do this in the '''terminal''' window '''type blockMesh''' and press '''Enter'''.
+
To do this in the '''terminal''' window '''type blockMesh''' and press '''Enter'''
  
  
'''Meshing '''has been done.
+
'''Meshing '''has been done
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Running the solver : pisoFoam
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Running the solver : pisoFoam
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now we can '''run''' the '''solver'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now we can '''run''' the '''solver'''
  
  
To do this in the '''terminal''' type '''pisoFoam''' and press '''Enter'''.
+
To do this in the '''terminal''' type '''pisoFoam''' and press '''Enter'''
  
 +
The '''iterations''' running can be seen in the '''terminal''' window
  
The '''iterations''' running can be seen in the '''terminal''' window.
 
  
 
+
It may take some time for the '''iterations''' to stop
It may take some time for the '''iterations''' to stop.
+
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Post-processing the results in paraview
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Post-processing the results in paraview
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| The '''iterations''' running will stop at the end of the '''time step'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| The '''iterations''' running will stop at the end of the '''time step'''
  
  
To visualize the results let us open the '''ParaView''' window.
+
To visualize the results let us open the '''ParaView''' window
  
  
To do this in the '''terminal''' type '''paraFoam''' and press '''Enter'''.
+
To do this in the '''terminal''' type '''paraFoam''' and press '''Enter'''
  
  
This will open the '''ParaView''' window.
+
This will open the '''ParaView''' window
  
 
|-
 
|-
Line 332: Line 331:
  
 
Change the drop down menu from solid color to U
 
Change the drop down menu from solid color to U
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| On the left hand side in the '''Object Inspector menu''' click on '''Apply'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| On the left hand side in the '''Object Inspector menu''' click on '''Apply'''
 +
 
  
 +
You can see the '''lid driven cavity geometry'''
  
You can see the '''lid driven cavity geometry'''.
 
  
 +
A common visualization is '''surface plots'''
  
A common visualization is '''surface plots'''.
 
  
 +
Change the display to '''Surface''' in the column.
  
Change the display to '''Surface''' in the column and from the drop down menu change from '''solid colour''' to '''U'''
+
And from the drop down menu change from '''solid colour''' to capital '''U'''
  
  
You can see the '''initial condition''' of '''velocity'''.
+
You can see the '''initial condition''' of '''velocity'''
  
 
|-
 
|-
Line 351: Line 352:
  
 
Toggle on the color legend  
 
Toggle on the color legend  
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now on top of the '''ParaView''' window you can see the '''VCR control'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Now on the top of the '''ParaView''' window you can see the '''VCR control'''
  
  
Click on the '''Play''' button.
+
Click on the '''Play''' button
  
  
You can see the motion of the fluid inside the '''cavity'''.
+
You can see the motion of the fluid inside the '''cavity'''
  
  
You can also '''toggle''' on the '''colour legend''' on the left hand side top of '''ParaView active variable control menu'''.
+
You can also '''toggle''' on the '''colour legend''' from the left hand side top of '''ParaView active variable control menu'''
  
  
Click on it. You can see the colour '''legend'''.
+
Click on it. You can see the colour '''legend'''
  
 
|-
 
|-
Line 376: Line 377:
  
 
*On top of the menu bar of '''ParaView'''
 
*On top of the menu bar of '''ParaView'''
*Go to '''Filters > Common and  Stream Tracers'''
+
*Go to '''Filters > Common and  Stream Tracer'''
 
*Click on it
 
*Click on it
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Streamlines on top
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Streamlines on top
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| On the left hand side of the '''Object inspector menu''' you can see '''Apply'''.  Click on it.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| On the left hand side of the '''Object inspector menu''' you can see '''Apply'''.  Click on it
  
  
You can see the '''stream lines''' at the centre of the '''lid driven cavity'''.
+
You can see the '''stream lines''' at the centre of the '''lid driven cavity'''
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Streamlines view
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Streamlines view
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can also change the '''orientation''' in which the '''stream lines''' are viewed.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can also change the '''orientation''' in which the '''stream lines''' are viewed
  
  
To do this , scroll down.
+
To do this , scroll down
  
  
You can see the''' seed type.'''
+
You can see the''' seed type'''
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Shift to right >> change point source to line source.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Shift to right >> change point source to line source
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Let me shift this to the right and change from''' point source''' to '''line source'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Let me shift this to the right and change from''' point source''' to '''line source'''
  
 
|-
 
|-
Line 408: Line 409:
 
Click on the X axis  
 
Click on the X axis  
  
Delete this.
+
Delete this
  
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see the '''X, Y''' and '''Z axes''' which are visible.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can see the '''X, Y''' and '''Z axes''' which are visible
  
  
Select any one of these axis in which you would like to view the '''stream lines'''.
+
Select any of these axis in which you would like to view the '''stream lines'''
  
  
I will select the '''Y axis''' and click '''Apply'''.
+
I will select the '''Y axis''' and click '''Apply'''
  
  
You can see the '''streamlines''' along the Y axis.
+
You can see the '''streamlines''' along the Y axis
  
  
Similarly you can select the '''X axis''' and plot the '''streamlines''' along the '''X axis'''.
+
Similarly you can select the '''X axis''' and plot the '''streamlines''' along the '''X axis'''
  
  
Now delete this.
+
Now delete this
  
 
|-
 
|-
Line 432: Line 433:
  
 
Save as .csv format
 
Save as .csv format
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can also plot the '''velocity''' along '''X''' and '''Y''' axis using '''plot over line'''.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can also plot the '''velocity''' along '''X''' and '''Y''' axis using '''plot over line'''
  
To do this go to '''Filter > Data Analysis and Plot over line'''.
+
To do this go to '''Filters > Data Analysis and Plot over line'''
  
Save the data as''' dot csv''' file from the file menu.
+
Save the data as''' dot csv''' from the '''File''' menu
  
Click on '''Save Data'''.
+
Click on '''Save Data'''
  
 
|-
 
|-
Line 447: Line 448:
  
 
For Re= 10000
 
For Re= 10000
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can plot this data in '''libre office spreadsheet '''or any other '''plotting''' software of your choice.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| You can plot this data in '''LibreOffice spreadsheet '''or any other '''plotting''' software of your choice
  
Now let me switch back to the slides.
+
Now let me switch back to the slides
  
The results obtained can be '''validated''' by using results of '''Ghia et.al for Reynolds Number , Re is equal to 10000'''
+
The results obtained can be '''validated''' by using results of '''Ghia et.al for Reynolds Number, Re equal to 10000'''
  
 
|-
 
|-
Line 457: Line 458:
  
 
Summary
 
Summary
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| That's all we have in this tutorial.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| That's all we have in this tutorial
  
 
Let us summarise.
 
Let us summarise.
Line 463: Line 464:
 
*and plotting '''stream lines '''in''' ParaView'''
 
*and plotting '''stream lines '''in''' ParaView'''
  
This brings us to the end of the tutorial.
+
This brings us to the end of the tutorial
  
 
|-
 
|-
Line 471: Line 472:
 
* Modify the '''grid size''' of the '''cavity'''
 
* Modify the '''grid size''' of the '''cavity'''
  
* Change it to (100 100 1)
+
* Change it to 100 100 1
  
 
* And visualise the results in '''ParaView''' using '''streamlines'''
 
* And visualise the results in '''ParaView''' using '''streamlines'''
Line 482: Line 483:
 
http://spoken-tutorial.org/What_is_a_Spoken_Tutorial  
 
http://spoken-tutorial.org/What_is_a_Spoken_Tutorial  
  
* It summarizes the Spoken Tutorial project.
+
* It summarizes the Spoken Tutorial project  
  
* If you do not have good bandwidth, you can download and watch it.
+
* If you do not have good bandwidth, you can download and watch it
  
 
|-
 
|-
Line 494: Line 495:
 
-Gives certificates to those who pass an online test  
 
-Gives certificates to those who pass an online test  
  
-For more details, please write to  
+
-For more details please write to  
  
contact@spoken-tutorial.org  
+
contact at the rate spoken hyphen tutorial dot org  
  
 
|-
 
|-
Line 545: Line 546:
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 13: Acknowledgements
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 13: Acknowledgements
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Spoken Tutorials project is a part of the Talk to a Teacher project.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:0.05pt solid #000000;padding:0.097cm;"| Spoken Tutorials project is a part of the Talk to a Teacher project
  
It is supported by the National Mission on Education through ICT, MHRD, Government of India.
+
It is supported by the National Mission on Education through ICT, MHRD, Government of India  
  
 
More information on the this mission is available at this URL http://spoken-tutorial.org/NMEICT-Intro
 
More information on the this mission is available at this URL http://spoken-tutorial.org/NMEICT-Intro

Latest revision as of 14:41, 5 August 2019

Tutorial: Turbulent flow in a lid driven cavity


Script : Chaitanya Talnikar, Shekhar Mishra , Rahul Joshi


Narration : Rahul Joshi


Keywords: Video tutorial ,CFD,Turbulent Flow in Lid driven cavity,OpenFOAM.


Visual Cue
Narration
Slide 1: Hello and welcome to the spoken tutorial on modelling Turbulent flow in a Lid Driven Cavity using OpenFOAM
Slide 2:


Learning Objectives

In this tutorial I will show you
  • Solving turbulent case in OpenFOAM and
  • Plotting streamlines in ParaView
Slide 3:


System Requirement

To record this tutorial I am using
  • Linux operating system Ubuntu version 12.04
  • OpenFoam version 2.1.1 and
  • ParaView version 3.12.0
Slide 4:

System Requirement

  • The tutorials were recorded using the versions specified in previous slide
  • Subsequently the tutorials were edited to latest versions
  • To install latest system requirements go to Installation Sheet
Slide 5 :


Prerequisites

To practice this tutorial you should have some basic knowledge of
  • Turbulence modelling
  • Knowledge of how to solve flow in a Lid driven cavity
  • If not so please refer to the relevant tutorial on our website
Demo:

Set up working Directory


This problem is identical in geometry and boundary condition to the 'Lid Driven Cavity' problem discussed in the basic level tutorial


Please make a note this problem is already set up in pisoFoam solver in OpenFoam directory


The boundary conditions are the Lid velocity U is 1 m/s


And we are solving this for a Reynolds number Re equal to 10000

Slide 6: Solver We are using a Transient solver for incompressible turbulent flow of Newtonian fluids called as pisoFoam
Steps in setting up the problem Now let us open the Terminal window by pressing Ctrl+Atl+T keys together
In the terminal window

Type cd tutorials

In the terminal window type run and press Enter


Now type cd space tutorials and press Enter

Type cd incompressible Now type cd space incompressible and press Enter
Type cd pisoFoam Now type cd space pisoFoam (Note that F here is capital ) and press Enter
Type ls


Two Folders les and ras


Cavity folder inside RAS

Now type ls and press Enter


In this you will see two folders les and ras


Our problem setup is inside ras folder which is called as reynolds averaged stress

Type cd ras


ls

Our folder name is cavity


Now type cd space ras and press Enter


Now type ls and press Enter

Type cd cavity

ls

You can see the cavity folder. Let me clear this off


Now type cd space cavity and press Enter


Now type ls and press Enter

Point to the 3 folders

Boundary and Initial conditions


0 folder

You can see three folders 0, constant and system


The initial conditions are specified within the files in the '0' directory


Now let us take a look at the files in the '0' directory

Inside the 0 folder


Type ls

To do this, in the command terminal type cd space 0 and press Enter


Now type ls and press Enter

Point to the files as per narration You can see the files named as epsilon, k, nut, nutilda, p, R and U


These files are to be kept as default until the inlet parameters don't change


If any changes are to be done please refer to the tutorial on Simulating flow in a channel using OpenFoam, to calculate these values

Type cd.. Now type cd space dot dot and press Enter


Let me clear this off

Type cd constant

Type ls

Let us open the constant folder


To do this type cd space constant and press Enter


Now type ls and press Enter

PolyMesh folder and fluid property files In this you will see the polyMesh folder containing
  • the geometry of the case inside blockMeshDict
  • and the fluid properties
Point to the two files as per narration In this case you will see two more files other than transportProperties named as RASProperties and turbulenceProperties


Let us open these two files

RASProperties In the terminal type gedit (space) RASProperties and press Enter


Scroll down


RASProperties contain the Reynolds average stress model for this case, which is kept as kEpsilon


Close this

turbulentProperties Now in the command terminal, type gedit space turbulentproperties and press Enter


Scroll down


The simulation type model for this case is kept as RASModel


Close this

TransportModel


Change the value of viscosity

Now let us open the transportProperties model


To do this, in the terminal type gedit space transportProperties and press Enter


The transportModel we are using here is Newtonian and the Viscosity is kept as 1 e raise to -4


Close this

Do not change the blockMeshDict file


The system folder is to be kept default

We are not changing the geometry in this case


So we need not go inside the polyMesh folder and look at the blockMeshDict file


It can be kept as it is

Type cd.. In the terminal type cd space dot dot and press Enter


We will keep the system folder default as there are no changes inside it

Meshing the geometry


blockMesh


Meshing is done

Now, we are done with the setup


We can mesh the geometry


To do this in the terminal window type blockMesh and press Enter


Meshing has been done

Running the solver : pisoFoam Now we can run the solver


To do this in the terminal type pisoFoam and press Enter

The iterations running can be seen in the terminal window


It may take some time for the iterations to stop

Post-processing the results in paraview The iterations running will stop at the end of the time step


To visualize the results let us open the ParaView window


To do this in the terminal type paraFoam and press Enter


This will open the ParaView window

View the geometry


Lid driven cavity geometry


Change the drop down menu from solid color to U

On the left hand side in the Object Inspector menu click on Apply


You can see the lid driven cavity geometry


A common visualization is surface plots


Change the display to Surface in the column.

And from the drop down menu change from solid colour to capital U


You can see the initial condition of velocity

Click on the Play button on VCR control for animation


Toggle on the color legend

Now on the top of the ParaView window you can see the VCR control


Click on the Play button


You can see the motion of the fluid inside the cavity


You can also toggle on the colour legend from the left hand side top of ParaView active variable control menu


Click on it. You can see the colour legend

Visualise the streamlines


Filters > Common > Stream Tracers


Now to visualise the stream lines
  • On top of the menu bar of ParaView
  • Go to Filters > Common and Stream Tracer
  • Click on it
Streamlines on top On the left hand side of the Object inspector menu you can see Apply. Click on it


You can see the stream lines at the centre of the lid driven cavity

Streamlines view You can also change the orientation in which the stream lines are viewed


To do this , scroll down


You can see the seed type

Shift to right >> change point source to line source Let me shift this to the right and change from point source to line source
Plot streamlines about X, Y and Z axis


Click on the Y axis

Click on the X axis

Delete this

You can see the X, Y and Z axes which are visible


Select any of these axis in which you would like to view the stream lines


I will select the Y axis and click Apply


You can see the streamlines along the Y axis


Similarly you can select the X axis and plot the streamlines along the X axis


Now delete this

Plot data over line


Save as .csv format

You can also plot the velocity along X and Y axis using plot over line

To do this go to Filters > Data Analysis and Plot over line

Save the data as dot csv from the File menu

Click on Save Data

Plot the results


Validate the results with Ghia et.al.

For Re= 10000

You can plot this data in LibreOffice spreadsheet or any other plotting software of your choice

Now let me switch back to the slides

The results obtained can be validated by using results of Ghia et.al for Reynolds Number, Re equal to 10000

Slide 7:

Summary

That's all we have in this tutorial

Let us summarise.

  • Turbulent Flow in a Lid Driven Cavity
  • and plotting stream lines in ParaView

This brings us to the end of the tutorial

Slide 8: Assignment As an assignment
  • Modify the grid size of the cavity
  • Change it to 100 100 1
  • And visualise the results in ParaView using streamlines
  • Watch the video available at this URL:

http://spoken-tutorial.org/What_is_a_Spoken_Tutorial

  • It summarizes the Spoken Tutorial project
  • If you do not have good bandwidth, you can download and watch it
The Spoken Tutorial Project Team

-Conducts workshops using spoken tutorials

-Gives certificates to those who pass an online test

-For more details please write to

contact at the rate spoken hyphen tutorial dot org

Slide 9:

Forum to answer questions

  • Do you have questions on THIS Spoken Tutorial?
  • Choose the minute and second where you have the question
  • Explain your question briefly
  • Someone from the FOSSEE team will answer them. Please visit

http://forums.spoken-tutorial.org/

Slide 10:

Forum to answer questions

  • Questions not related to the Spoken Tutorial?
  • Do you have general/technical questions on the Software?
  • Please visit the FOSSEE forum

http://forums.fossee.in/

  • Choose the Software and post your question


Slide 11:

Lab Migration Project

  • We coordinate migration from commercial CFD software like ANSYS to OpenFOAM
  • We conduct free Workshops and provide solutions to CFD Problem Statements in OpenFOAM

For more details visit this site: http://cfd.fossee.in/


Slide 12:

Case Study Project

  • We invite students to solve a feasible CFD problem statement of reasonable complexity using OpenFOAM
  • We give honorarium and certificate to those who do this

For more details visit this site: http://cfd.fossee.in/


Slide 13: Acknowledgements Spoken Tutorials project is a part of the Talk to a Teacher project

It is supported by the National Mission on Education through ICT, MHRD, Government of India

More information on the this mission is available at this URL http://spoken-tutorial.org/NMEICT-Intro

About the contributor The script is contributed by Shekhar Mishra and Chaitanya talnikar

This is Rahul Joshi from IIT BOMBAY signing off.

Thanks for joining

Contributors and Content Editors

DeepaVedartham, Nancyvarkey, Pravin1389, Rahuljoshi