Difference between revisions of "OpenFOAM/C3/Turbulent-Flow-in-a-Lid-driven-Cavity/English"
Nancyvarkey (Talk | contribs) |
Nancyvarkey (Talk | contribs) |
||
(13 intermediate revisions by one other user not shown) | |||
Line 28: | Line 28: | ||
| 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 tutorial I will show you | | 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 tutorial I will show you | ||
− | * Solving '''turbulent | + | * Solving '''turbulent case''' in '''OpenFOAM''' and |
* Plotting '''streamlines''' in '''ParaView''' | * Plotting '''streamlines''' in '''ParaView''' | ||
Line 37: | Line 37: | ||
System Requirement | System Requirement | ||
| 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 record this tutorial I am using | | 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 record this tutorial I am using | ||
− | * '''Linux operating system''' '''Ubuntu 12.04''' | + | * '''Linux operating system''' '''Ubuntu version 12.04''' |
− | * '''OpenFoam version 2.1.1''' | + | * '''OpenFoam version 2.1.1''' and |
* '''ParaView version 3.12.0''' | * '''ParaView version 3.12.0''' | ||
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 | + | | 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 | + | Please make a note this problem is already set up in '''pisoFoam solver''' in '''OpenFoam''' directory |
− | The '''boundary conditions''' are the '''Lid velocity U | + | The '''boundary conditions''' are the '''Lid velocity U is 1 m/s''' |
− | And we are solving this for a '''Reynolds number Re | + | And we are solving this for a '''Reynolds number Re equal to 10000''' |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 6: Solver | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 6: 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;"| We | + | | 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 using a '''Transient solver''' for '''incompressible turbulent flow of Newtonian fluids''' called as '''pisoFoam''' |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Steps in setting up the problem | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Steps in setting up the problem | ||
− | | 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 '''Terminal''' window by pressing '''Ctrl+Atl+T''' keys together | + | | 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 '''Terminal''' window by pressing '''Ctrl+Atl+T''' keys together |
|- | |- | ||
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;"| | + | | 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:none;padding:0.097cm;"| Type cd pisoFoam | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type cd 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;"| | + | | 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 pisoFoam''' (Note that '''F''' here is capital ) 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 | + | 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 |
− | + | 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 |
− | + | 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 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''' |
− | 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 transportProperties '''and press '''Enter''' | + | To do this, in the '''terminal''' type '''gedit space transportProperties '''and press '''Enter''' |
− | The '''transportModel '''we are using here is '''Newtonian '''and '''Viscosity''' is kept as '''1 e raise to -4''' | + | The '''transportModel '''we are using here is '''Newtonian '''and the '''Viscosity''' is kept as '''1 e raise to -4''' |
− | 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 | + | 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''' | |
− | 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 | ||
− | |||
− | + | 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 | + | 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 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 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 ''' | + | 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 373: | Line 374: | ||
− | | 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;"| | + | | 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 to visualise the '''stream lines''' |
− | *On | + | *On top of the menu bar of '''ParaView''' |
− | *Go to '''Filters > Common | + | *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 | + | 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 | + | 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 '''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 | + | | 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 ''' | + | To do this go to '''Filters > Data Analysis and Plot over line''' |
− | Save the data as''' | + | 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 ''' | + | | 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 results | + | 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. | ||
− | |||
− | |||
*'''Turbulent Flow in a Lid Driven Cavity''' | *'''Turbulent Flow in a Lid Driven Cavity''' | ||
*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 473: | Line 472: | ||
* Modify the '''grid size''' of the '''cavity''' | * Modify the '''grid size''' of the '''cavity''' | ||
− | * Change it to | + | * Change it to 100 100 1 |
− | * | + | * And visualise the results in '''ParaView''' using '''streamlines''' |
|- | |- | ||
Line 484: | 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 496: | Line 495: | ||
-Gives certificates to those who pass an online test | -Gives certificates to those who pass an online test | ||
− | -For more details | + | -For more details please write to |
− | contact | + | contact at the rate spoken hyphen tutorial dot org |
|- | |- | ||
Line 547: | 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 | + | | 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 | |
− | More information on the | + | 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.
|
|
Slide 1: | Hello and welcome to the spoken tutorial on modelling Turbulent flow in a Lid Driven Cavity using OpenFOAM |
Slide 2:
|
In this tutorial I will show you
|
Slide 3:
|
To record this tutorial I am using
|
Slide 4:
System Requirement
|
|
Slide 5 :
|
To practice this tutorial you should have some basic knowledge of
|
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
|
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
|
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
|
Now type ls and press Enter
|
Type cd ras
|
Our folder name is cavity
|
Type cd cavity
ls |
You can see the cavity folder. Let me clear this off
|
Point to the 3 folders
Boundary and Initial conditions
|
You can see three folders 0, constant and system
|
Inside the 0 folder
|
To do this, in the command terminal type cd space 0 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
|
Type cd.. | Now type cd space dot dot and press Enter
|
Type cd constant
Type ls |
Let us open the constant folder
|
PolyMesh folder and fluid property files | In this you will see the polyMesh folder containing
|
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
|
RASProperties | In the terminal type gedit (space) RASProperties and press Enter
|
turbulentProperties | Now in the command terminal, type gedit space turbulentproperties and press Enter
|
TransportModel
|
Now let us open the transportProperties model
|
Do not change the blockMeshDict file
|
We are not changing the geometry in this case
|
Type cd.. | In the terminal type cd space dot dot and press Enter
|
Meshing the geometry
|
Now, we are done with the setup
|
Running the solver : pisoFoam | Now we can run the solver
The iterations running can be seen in the terminal window
|
Post-processing the results in paraview | The iterations running will stop at the end of the time step
|
View the geometry
|
On the left hand side in the Object Inspector menu click on Apply
And from the drop down menu change from solid colour to capital U
|
Click on the Play button on VCR control for animation
|
Now on the top of the ParaView window you can see the VCR control
|
Visualise the streamlines
|
Now to visualise the stream lines
|
Streamlines on top | On the left hand side of the Object inspector menu you can see Apply. Click on it
|
Streamlines view | You can also change the orientation in which the stream lines are viewed
|
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 X axis Delete this |
You can see the X, Y and Z axes which are visible
|
Plot data over line
|
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
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.
This brings us to the end of the tutorial |
Slide 8: Assignment | As an assignment
|
http://spoken-tutorial.org/What_is_a_Spoken_Tutorial
| |
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
|
|
Slide 10:
Forum to answer questions
|
|
Slide 11:
Lab Migration Project
For more details visit this site: http://cfd.fossee.in/ |
|
Slide 12:
Case Study Project
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 |