Difference between revisions of "OpenFOAM/C2/2D-Laminar-Flow-in-a-channel/English"

From Script | Spoken-Tutorial
Jump to: navigation, search
Line 31: Line 31:
  
  
Validation with '''analytic result'''
+
and Validation with '''analytic result'''
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 3: System requirement
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 3 : 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  
 
| 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  
  
Line 53: Line 53:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 4: System Requirement  
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 4 : 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;"| The tutorials were recorded using the versions specified in previous slide  
+
* The tutorials were recorded using the versions specified in previous slide  
  
Subsequently the tutorials were edited to latest versions
+
* Subsequently the tutorials were edited to latest versions
  
To install latest system requirements go to Installation Sheet
+
* To install latest system requirements go to Installation Sheet
  
 +
| 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:none;padding:0.097cm;"| Slide 4 : Pre-requisite
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 5 : Pre-requisite
 
| 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 should know how to create '''geometry''' using '''OpenFOAM '''
 
| 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 should know how to create '''geometry''' using '''OpenFOAM '''
  
Line 69: Line 70:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 5 : About flow in a channel
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 6 : About flow in a channel
| 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;"| '''Channel flow''' problem description
+
| 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 simulate flow in a '''Channel''' to determine
 +
* '''Flow development length''' along the downstream
 +
|-
 +
|-
 +
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 7 : Channel Flow
 +
| 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;"| '''Channel flow''' problem description.
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 5 : boudnary patches
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 7 : Channel flow
 
| 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 '''boundary names''' and '''inlet conditions''' are shown in this figure
 
| 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 '''boundary names''' and '''inlet conditions''' are shown in this figure
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 6 :  
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 8 : Boundary Conditions
 
| 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 '''flow develpoment length''' is given by the formula  
 
| 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 '''flow develpoment length''' is given by the formula  
  
'''L= 0.05 * Re * D'''
+
'''L = 0.05 * Re * D'''
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 6 :
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 8 : Boundary Conditions
 
| 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;"| '''Re''' which is the '''Reynolds number'''
 
| 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;"| '''Re''' which is the '''Reynolds number'''
  
Line 93: Line 99:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 6 :
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 8 : Boundary Conditions
 
| 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;"| Using the '''formula''' length of the channel comes out to be '''5 meters '''
 
| 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;"| Using the '''formula''' length of the channel comes out to be '''5 meters '''
  
Line 106: Line 112:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide : file structure of the case
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 9 : File structure  
 
| 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 is a '''steady state 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;"| This is a '''steady state problem '''
  
  
 
Therefore we are using a '''steady state incompressible''' solver for this case
 
Therefore we are using a '''steady state incompressible''' solver for this case
 
+
This is the '''file structure''' of our case
 +
The '''folder''' should be created in the solver type that we choose
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Create a floder and name it as channel in simpleFoam folder
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Create a floder and name it as channel in simpleFoam folder
Line 139: Line 146:
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"|  
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;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;"| Paste it in the channel folder created and  
+
| 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;"| Paste it in the channel folder and  
  
  
Line 169: Line 176:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| incompresible 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;"| 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'''
 
| 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 217: Line 224:
  
  
'''kinematic viscosity (nu) is set at 0.001 m²/s'''.
+
'''kinematic viscosity (that is nu), in this case, is set at 0.001 m²/s (meter square per second)'''.
  
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| cd polyMesh 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;"| cd polyMesh 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 polyMesh 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 terminal type '''cd polyMesh and press enter'''
  
 
|-
 
|-
Line 257: Line 264:
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Boundary conditions and types  
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Boundary conditions and types  
| 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;"| Next we have the '''boundary conditions''' and their types which are set as '''inlet ,outlet,top,bottom '''
+
| 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;"| Next we have set up '''boundary conditions''' and their types which are set as '''inlet ,outlet,top and bottom '''
  
  
Line 264: Line 271:
 
|-
 
|-
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| FrontAndBack boundary type as empty
 
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| FrontAndBack boundary type as empty
| 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;"| Since it is a '''2D problem frontAndBack''' is kept as '''empty'''
+
| 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;"| As this is a '''2D geometry frontAndBack''' is kept as '''empty'''
  
 
|-
 
|-
Line 306: Line 313:
  
 
initial flow conditions: p, R and U
 
initial flow conditions: 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;"| It should contain various files named
+
| 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;"| It should contain various files such as
  
  
Line 319: Line 326:
  
 
'''p , R and capital U''' which the are '''initial conditions of the flow'''
 
'''p , R and capital U''' which the are '''initial conditions of the flow'''
 +
 +
Now let me switch back to the slides
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 7 : formula to calculate epsilon and k.
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 10 : Calculate K.
 
| 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;"| '''Calculate k''' which is the '''turbulent kinetic energy'''
 
| 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;"| '''Calculate k''' which is the '''turbulent kinetic energy'''
  
Line 337: Line 346:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 8: calculate epsilon
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 11 : Calculate epsilon
  
  
epsilon
+
epsilon - rate of disspiation turbulent kinetic energy
  
  
Line 380: Line 389:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| In the terminal window
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| In the terminal window type '''cd ..''' 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 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;"| Now in the terminal Type '''cd (space) ..(dot dot) and press enter'''
  
Line 392: Line 401:
  
 
terminal window type : blockMesh
 
terminal window type : blockMesh
| 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 need to '''mesh the geometry'''
+
| 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 need to '''mesh the geometry'''
  
  
Line 407: Line 416:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 8 : Solver type
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 12 : 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;"| The '''type of solver''' we are using here is '''simpleFoam'''
 
| 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 '''type of solver''' we are using here is '''simpleFoam'''
  
Line 476: Line 485:
  
  
you can see the final value of '''velocity magnitude'''
+
you can see the final value of the '''velocity magnitude'''
  
 
|-
 
|-
Line 515: Line 524:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 8:Validation
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 13 :Validation
 
| 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 results obtained can be '''validated '''with
 
| 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 results obtained can be '''validated '''with
  
Line 531: Line 540:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 9 : Summary
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 14 : 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;"| In this tutorial we learnt
 
| 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 we learnt
  
Line 547: Line 556:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 10 :
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 15 :
  
 
Assignment
 
Assignment
Line 557: Line 566:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 11 : Abour spoken tutorials
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 16 : About the Spoken Tutorial 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;"| The video available at this URL:
+
| 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;"| Watch the video available at this URL:
  
 
http://spoken-tutorial.org/What_is_a_Spoken_Tutorial  
 
http://spoken-tutorial.org/What_is_a_Spoken_Tutorial  
Line 567: Line 576:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 12 :
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 17 : Spoken Tutorial Workshops
 
+
About 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;"| The Spoken Tutorial Project Team
 
| 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 Spoken Tutorial Project Team
  
Line 579: Line 586:
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide:  
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 18 :  
  
 
Forum to answer questions
 
Forum to answer questions
| 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;"| Do you have questions on THIS Spoken Tutorial?
+
* Do you have questions on THIS Spoken Tutorial?
Choose the minute and second where you have the question
+
* Choose the minute and second where you have the question
Explain your question briefly
+
* Explain your question briefly
Someone from the FOSSEE team will answer them. Please visit
+
* Someone from the FOSSEE team will answer them. Please visit
 
http://forums.spoken-tutorial.org/
 
http://forums.spoken-tutorial.org/
 +
| 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:none;padding:0.097cm;"| Slide:  
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 19 :  
  
 
Forum to answer questions
 
Forum to answer questions
| 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;"| Questions not related to the Spoken Tutorial?
+
* Questions not related to the Spoken Tutorial?
Do you have general/technical questions on the Software?
+
* Do you have general/technical questions on the Software?
Please visit the FOSSEE forum
+
* Please visit the FOSSEE forum
 
http://forums.fossee.in/
 
http://forums.fossee.in/
Choose the Software and post your question
+
* Choose the Software and post your question
 +
| 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:none;padding:0.097cm;"| Slide:  
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 20 :  
  
 
Lab Migration Project
 
Lab Migration 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;"| We coordinate migration from commercial CFD software like ANSYS to OpenFOAM
+
* We coordinate migration from commercial CFD software like ANSYS to OpenFOAM
We conduct free Workshops and provide solutions to CFD Problem Statements in OpenFOAM
+
* We conduct free Workshops and provide solutions to CFD Problem Statements in OpenFOAM
 
For more details visit this site:
 
For more details visit this site:
 
http://cfd.fossee.in/
 
http://cfd.fossee.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;"|
  
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide:  
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 21:  
  
 
Case Study Project
 
Case Study 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;"|  We invite students to solve a feasible CFD problem statement of reasonable complexity using OpenFOAM
+
* 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
+
* We give honorarium and certificate to those who do this
 
For more details visit this site:
 
For more details visit this site:
 
http://cfd.fossee.in/
 
http://cfd.fossee.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;"| 
  
  
 
|-
 
|-
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 13
+
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 22 :
  
 
Acknowledgement
 
Acknowledgement
Line 630: Line 641:
 
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 same is available at the following URL link http://spoken-tutorial.org/NMEICT-Intro
+
More information on the same is available at the following URL http://spoken-tutorial.org/NMEICT-Intro
  
 
|-
 
|-

Revision as of 11:32, 30 May 2019

Tutorial: Simulating Laminar flow in a channel.


Script and Narration : Rahul Joshi


Keywords: Video tutorial,CFD,laminar flow,simpleFoam,channel


Visual Cue
Narration
Slide 1 Hello and welcome to the spoken tutorial on Simulating Laminar Flow in a Channel using openfoam
Slide 2 : Learning Objectives In this tutorial I will show you


2D geometry of channel


Meshing the Geometry


Solving and Post Processing results in Paraview


and Validation with analytic result

Slide 3 : System requirement To record this tutorial


I am using Linux Operating system Ubuntu 12.04 .


OpenFOAM version 2.1.1


ParaView version 3.12.0


Note that OpenFOAM version 2.1.1 is supported on ubuntu version 12.04


Hence forth all the tutorials will be covered using OpenFOAM version 2.1.1 and ubuntu version 12.04

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 : Pre-requisite You should know how to create geometry using OpenFOAM


If not, please refer to the relevant tutorials on our website

Slide 6 : About flow in a channel
  • We simulate flow in a Channel to determine
  • Flow development length along the downstream
Slide 7 : Channel Flow Channel flow problem description.
Slide 7 : Channel flow The boundary names and inlet conditions are shown in this figure
Slide 8 : Boundary Conditions The flow develpoment length is given by the formula

L = 0.05 * Re * D

Slide 8 : Boundary Conditions Re which is the Reynolds number


and D which is the channel height



Slide 8 : Boundary Conditions Using the formula length of the channel comes out to be 5 meters


and height is kept as 1 meters.


The Inlet velocity is 1 meters per second


And we are solving this for a Reynolds number ( Re ) equal 100

Slide 9 : File structure This is a steady state problem


Therefore we are using a steady state incompressible solver for this case This is the file structure of our case The folder should be created in the solver type that we choose

Create a floder and name it as channel in simpleFoam folder


right click >> new folder

I have already created a folder in


simpleFoam folder of incompressible flow solvers

Nmae the folder channel The folder is named as channel


Now let me switch to the folder

copy 0,constant and system from pitzDaily to this channel folder Copy 0,Constant and System folders of any other case file in the simpleFoam directory
I have copied the file structure of case of pitzDaily
Paste it in the channel folder and


make the necessary changes in the geometry,


boundary faces and boundary condition

Now let me open the command terminal
Press Ctrl+Alt+t keys simultaneously To do this press Ctrl+Alt +t keys simultaneously on your keyboard



run and press enter In the terminal

Type run and press enter

now type cd space tutorials and press enter
incompressible and press enter now type cd space incompressible and press enter
simpleFoam and press enter type cd space simpleFoam and press enter
channel and press enter now type cd channel and press enter
ls and press enter now type ls and press enter
You will see three folders 0 , Constant and System
Type in command terminal

cd constant and press enter

now type cd constant and press enter
ls and press enter now type ls and press enter
polymesh folder and 2 other files In this you will see files containing


properties of fluid and a folder named polymesh

gedit RASProperties and press enter RASProperties contains Reynolds-averaged stress model
Transportproperties transportProperties contain the transport model


and


kinematic viscosity (that is nu), in this case, is set at 0.001 m²/s (meter square per second).

cd polyMesh and press enter Now in the terminal type cd polyMesh and press enter
ls and press enter Now type ls and press enter
you will see the blockMeshDict file here
gedit blockMeshDict and press enter To open the blockMeshDict file in the


terminal type gedit blockMeshDict and press enter


Scroll down

In the blockMeshDict file covertTometers is set to 1


Set the vertices for channel

The geomery is in meters so the convertTometers is set to 1


next we have defined the vertices of the channel

Mesh size for channel We have used a 100 X 100 mesh size here and cell spacing is kept as ( 1 1 1 )
Boundary conditions and types Next we have set up boundary conditions and their types which are set as inlet ,outlet,top and bottom



FrontAndBack boundary type as empty As this is a 2D geometry frontAndBack is kept as empty
Due to a 2 dimensional geometry Also this being a simple geometry


mergePatchPair and edges are to kept empty.


Close the blockMeshDict file

Terminal window In the command terminal Type cd space ..(dot dot) and press enter
Terminal window Again type cd space .. (dot dot) and press enter
In the terninal window type cd 0


type ls and press enter

Now in the terminal type cd space 0 (Zero) and press enter


Now Type ls and press enter

In terminal you can see these slides This contains the intial boundary conditions


and wall functions for the channel case

Wall functions : epsilon, k, nut, nutilda


initial flow conditions: p, R and U

It should contain various files such as


epsilon ,k, nut,nuTilda


which are the wall functions


and


p , R and capital U which the are initial conditions of the flow

Now let me switch back to the slides

Slide 10 : Calculate K. Calculate k which is the turbulent kinetic energy


from the formula given in the slide


Where Ux, Uy and Uz are the velocity


components in the x, y and z direction


And U' ( dash ) = 0.05 times u actual

Slide 11 : Calculate epsilon


epsilon - rate of disspiation turbulent kinetic energy


C mu - constant


l – length of channel

Calculate epsilon from the formula given


Where epsilon is the rate of disspiation turbulent kinetic energy


C mu is a constant and its value is 0.09


And l is the length of the channel


Let me minimise this

Change only the boundary names Change only the boundary names in each of the above files
Do not chnge the values of nu,nuTilde and R Note that the values of nut, nuTilda and R


are to kept default

Rest of the files should contain initial value


for each of the boundary faces

In the terminal window type cd .. and press enter Now in the terminal Type cd (space) ..(dot dot) and press enter
No change in system folder There is no change to be done in the system folder
Mesh the geometry


terminal window type : blockMesh

Now, We need to mesh the geometry


To do this


In the command terminal type blockMesh and press enter


Meshing is done


Let me switch back to the slide

Slide 12 : Solver The type of solver we are using here is simpleFoam


It is a Steady-state solver for incompressible, turbulent flow


let me minimise this

Type simpleFoam and press enter In the command terminal type simpleFoam and Press enter
Iterations in terminal window Iterations running will be seen in the command terminal
Iterations running may take some time
Iterations converge or stop at end of time step The iterations will stop once the solution is


converged or it reaches its end time value

In terminal

type: paraFoam and press enter

To view the results in paraview in the terminal


type paraFoam and press enter


This will open up the paraview window

In paraview window


View the geometry

On let hand side of the paraview window click Apply


The geometry can be seen here.

Change from solid color to U On top of active variable control menu change


the drop down menu from solid color to capital U

Look at left side of the channel geometry You can see the initial state of velocity magnitude at inlet.
VCR control click PLAY button On top of the paraview window click on the


play button of the VCR control


you can see the final value of the velocity magnitude

Color legend from top left Also toggle on the color legend from the


left hand side top of active variable control menu


Click APPLY again

In object inspector menu


click on rescale to data

Now go to display


Scroll down


You can see Rescale , click on it

Check the color legend for this We can see that once the flow has fully devloped


it attains a maximum uniform velocity at the center


Now let me switch back to the slides



Slide 13 :Validation The results obtained can be validated with


the analytical solution for laminar flow in a


channel which is u(max)=1.5 Uavg


Using openfoam we obatain a velocity of 1.48 meters per second which is a good match


This brings us to the end of the tutorial

Slide 14 : Summary In this tutorial we learnt


The File structure of channel


Obtained solution using steady state solver


Viewed the geometry in paraview


Validation with analytic results

Slide 15 :

Assignment

As an assignment:

Solve the problem for Reynold Number 1500 and

validate it with the analytical result

Slide 16 : About the Spoken Tutorial Project 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.

Slide 17 : Spoken Tutorial Workshops 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@spoken-tutorial.org

Slide 18 :

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 19 :

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 20 :

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 21:

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 22 :

Acknowledgement


Spoken Tutorials are part of Talk to a Teacher project,

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

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

About the contributor This is Rahul Joshi from IIT BOMBAY signing off.

Thanks for joining

Contributors and Content Editors

DeepaVedartham, Nancyvarkey, Pravin1389, Rahuljoshi