Difference between revisions of "OpenFOAM/C2/2D-Laminar-Flow-in-a-channel/English"
Rahuljoshi (Talk | contribs) |
Nancyvarkey (Talk | contribs) |
||
(4 intermediate revisions by one other user not shown) | |||
Line 15: | Line 15: | ||
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 1 | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Slide 1 | ||
− | | 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;"| Hello and welcome to the spoken tutorial on '''Simulating Laminar Flow in a Channel 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;"| Hello and welcome to the spoken tutorial on '''Simulating 2D Laminar Flow in a Channel using openfoam''' |
|- | |- | ||
Line 21: | Line 21: | ||
| 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 | ||
+ | *'''2D geometry of channel ''' | ||
+ | *'''Meshing the Geometry''' | ||
+ | *'''Solving''' and '''Post Processing results''' in '''Paraview''' | ||
+ | *and Validation using '''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:0.05pt solid #000000;padding:0.097cm;"| 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''' | |
− | + | ||
|- | |- | ||
− | | 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 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 | |
− | + | | 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 5 : Prerequisite |
− | | 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;"| As a prerequisite for this tutorial, you should know how to create '''geometry''' using '''OpenFOAM.''' |
− | If not, please refer to the | + | If not, please refer to the relevant tutorials on our website. |
|- | |- | ||
− | | 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 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;"| ''' | + | | 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 | + | | 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;"| '''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 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: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 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 | + | | 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''' | ||
and '''D''' which is the '''channel height''' | and '''D''' which is the '''channel height''' | ||
− | |||
− | |||
− | |||
|- | |- | ||
− | | 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 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''' the '''length''' of the '''channel''' comes out to be '''5 meters ''' |
Line 94: | Line 93: | ||
− | And we are solving this for a '''Reynolds number | + | And we are solving this for a '''Reynolds number ( Re ) equal 100''' |
|- | |- | ||
− | | 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 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 | + | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Create a '''folder''' and name it as channel in '''simpleFoam folder'''. |
Line 111: | Line 111: | ||
− | '''simpleFoam folder''' of''' incompressible flow solvers''' | + | '''simpleFoam folder''' of''' incompressible flow solvers'''. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Nmae the folder channel | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Nmae the folder 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;"| The '''folder''' is named as''' 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;"| The '''folder''' is named as''' channel'''. |
− | Now let me switch to the '''folder''' | + | Now let me switch to the '''folder'''. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| copy 0,constant and system from pitzDaily to this channel folder | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| copy 0,constant and system from pitzDaily to this channel 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;"| ''' | + | | 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;"| Copy '''0, Constant''' and '''System folders''' of any other '''case'''. |
|- | |- | ||
| 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;"| I have copied the''' file structure''' of case of''' pitzDaily''' | + | | 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 '''simpleFoam directory''', I have copied the''' file structure''' of '''case''' of''' pitzDaily'''. |
|- | |- | ||
| 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 | + | | 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'''. |
− | make the necessary changes in the '''geometry | + | And make the necessary changes in the '''geometry, boundary faces''' and '''boundary condition'''. |
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
| 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;"| Now let me ''' | + | | 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 me open the '''command terminal'''. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Press Ctrl+Alt+t keys simultaneously | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Press Ctrl+Alt+t keys simultaneously | ||
− | | 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 press '''Ctrl+Alt +t''' keys simultaneously on your '''keyboard''' | + | | 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 press '''Ctrl+Alt +t''' keys simultaneously on your '''keyboard'''. |
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| run 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;"| 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 | + | | 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 '''run''' 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;"| | | 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;"| | + | | 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 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;"| | + | | 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;"| | + | | 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;"| simpleFoam 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;"| simpleFoam 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;"| | + | | 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;"| Type '''cd space simpleFoam''' 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;"| channel 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;"| channel 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;"| | + | | 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 channel''' 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;"| 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;"| 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;"| | + | | 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: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;"| You will see three | + | | 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 will see three '''folders 0 , Constant and System''' |
|- | |- | ||
Line 183: | Line 175: | ||
cd constant and press enter | cd constant 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;"| | + | | 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 constant''' 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;"| 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;"| 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;"| | + | | 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:none;padding:0.097cm;"| polymesh folder and 2 other files | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| polymesh folder and 2 other files | ||
− | | 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 '''files''' 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 '''files''' containing '''properties of fluid''' and a''' folder''' named '''polymesh'''. |
− | + | ||
− | + | ||
− | '''properties of fluid''' and a''' folder named polymesh''' | + | |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| gedit RASProperties 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;"| gedit 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;"| '''RASProperties''' contains '''Reynolds-averaged stress ''' | + | | 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;"| '''RASProperties''' contains '''Reynolds-averaged stress model'''. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Transportproperties | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Transportproperties | ||
− | | 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;"| '''transportProperties''' contain the '''transport 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;"| '''transportProperties''' contain the '''transport model'''. |
− | + | ||
− | + | ||
− | + | ||
− | '''kinematic viscosity (nu) is set at 0.001 m²/s'''. | + | And '''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 | + | | 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'''. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| 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;"| 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 | + | | 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: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;"| | + | | 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 will see the '''blockMeshDict''' file here. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| gedit blockMeshDict 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;"| gedit blockMeshDict 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 open the '''blockMeshDict''' file in the | + | | 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 open the '''blockMeshDict''' file in the '''terminal '''type '''gedit blockMeshDict and press '''Enter'''. |
− | + | Scroll down. | |
− | + | ||
− | + | ||
− | Scroll down | + | |
|- | |- | ||
Line 237: | Line 220: | ||
Set the vertices for channel | Set the vertices for 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;"| The | + | | 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 '''geometry''' is in '''meters '''so the '''convertTometers''' is set to 1. |
− | + | Next we have defined the '''vertices of the channel'''. | |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Mesh size for channel | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Mesh size for 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;"| We have used a '''100 X 100''' | + | | 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 have used a '''100 X 100 mesh size''' here and '''cell spacing''' is kept as '''( 1 1 1 )'''. |
|- | |- | ||
| 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 | + | | 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 '''inlet ,outlet,top''' and '''bottom'''. |
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
| 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;"| | + | | 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'''. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Due to a 2 dimensional geometry | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Due to a 2 dimensional 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;"| Also this being a '''simple 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;"| Also this being a '''simple geometry, mergePatchPair''' and '''edges''' are to be kept '''empty'''. |
− | + | Close the '''blockMeshDict''' file. | |
− | + | ||
− | + | ||
− | Close the '''blockMeshDict''' file | + | |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Terminal window | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Terminal window | ||
− | | 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 '''command terminal''' | + | | 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 '''command 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:none;padding:0.097cm;"| Terminal window | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Terminal window | ||
− | | 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;"| Again type '''cd space .. (dot dot)''' and press | + | | 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;"| Again 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:none;padding:0.097cm;"| In the | + | | 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''' type '''cd 0''' |
type ls and press enter | 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 in the | + | | 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 0 (Zero)''' and press '''Enter'''. |
− | Now | + | 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;"| In terminal you can see these slides | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| In terminal you can see these slides | ||
− | | 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 contains the ''' | + | | 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 contains the '''initial boundary conditions ''' and''' wall functions''' for the '''channel case'''. |
− | + | ||
− | + | ||
− | and''' wall functions''' for the '''channel case''' | + | |
|- | |- | ||
Line 297: | Line 271: | ||
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 | + | | 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 '''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. | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
− | | 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 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''' 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. | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | And '''U' ( dash ) = 0.05''' times '''u''' actual | + | |
|- | |- | ||
− | | 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 11 : Calculate epsilon |
− | epsilon | + | epsilon - rate of disspiation turbulent kinetic energy |
Line 338: | Line 299: | ||
l – length of channel | l – length of 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;"| Calculate '''epsilon''' from the '''formula''' given | + | | 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 '''epsilon''' from the '''formula''' given: |
− | Where '''epsilon''' is the''' rate of disspiation turbulent kinetic energy''' | + | Where '''epsilon''' is the''' rate of disspiation turbulent kinetic energy'''. |
Line 347: | Line 308: | ||
− | And''' l''' is the '''length of the channel''' | + | And''' l''' is the '''length of the channel'''. |
− | Let me minimise this | + | Let me minimise this. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Change only the boundary names | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Change only the boundary names | ||
− | | 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;"| Change only the '''boundary names''' in each of the above | + | | 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;"| Change only the '''boundary names''' in each of the above files. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Do not chnge the values of nu,nuTilde and R | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Do not chnge the values of nu,nuTilde and R | ||
− | | 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;"| Note that the values of '''nut, nuTilda and R ''' | + | | 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;"| Note that the values of '''nut, nuTilda''' and '''R ''' are kept to be default. |
− | + | ||
− | + | ||
− | are | + | |
|- | |- | ||
| 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;"| Rest of the | + | | 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;"| Rest of the files should contain '''initial value''' for each of the '''boundary faces'''. |
− | + | ||
− | + | ||
− | for each of the '''boundary faces''' | + | |
|- | |- | ||
− | | 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 | + | | 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:none;padding:0.097cm;"| No change in system folder | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| No change in system 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;"| There | + | | 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;"| There are no changes to be done in the '''system folder'''. |
|- | |- | ||
Line 383: | Line 338: | ||
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;"| | + | | 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'''. |
− | To do this | + | To do this, in the '''command terminal''' type '''blockMesh''' and press '''Enter'''. |
− | + | The '''Meshing''' is done. | |
− | + | Let me switch back to the slide. | |
− | + | ||
− | + | ||
− | Let me switch back to the slide | + | |
|- | |- | ||
− | | 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 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'''. |
− | It is a '''Steady-state''' | + | It is a '''Steady-state solver''' for '''incompressible''' and '''turbulent flows'''. |
− | + | Let me minimise this. | |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Type simpleFoam 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 simpleFoam 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 '''command terminal''' ''' | + | | 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 '''command terminal''' type '''simpleFoam''' 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;"| Iterations in terminal window | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Iterations in terminal window | ||
− | | 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;"| '''Iterations | + | | 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;"| '''Iterations''' running will be seen in the '''command terminal'''. |
|- | |- | ||
| 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;"| '''Iterations '''running may take some time | + | | 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;"| '''Iterations '''running may take some time. |
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Iterations converge or stop at end of time step | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Iterations converge or stop at end of 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''' will stop once the solution is | + | | 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''' will stop once the '''solution''' is '''converged''' or it reaches its '''end time value'''. |
− | + | ||
− | + | ||
− | '''converged''' or it reaches its '''end time value''' | + | |
|- | |- | ||
Line 430: | Line 379: | ||
type: paraFoam and press enter | type: paraFoam 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 view the results 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;"| To view the results in '''ParaView''' in the '''terminal,''' type '''paraFoam and press '''Enter'''. |
− | + | ||
− | + | ||
− | type '''paraFoam and press | + | |
− | This will open up the ''' | + | This will open up the '''ParaView''' window. |
|- | |- | ||
Line 443: | Line 389: | ||
View the geometry | View 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;"| On | + | | 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 left hand side of the '''ParaView''' window click '''Apply'''. |
Line 450: | Line 396: | ||
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Change from solid color to U | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Change 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 top of '''active variable control menu | + | | 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 top of '''active variable control menu''' change the drop-down menu from '''solid color''' to capital '''U'''. |
− | + | ||
− | + | ||
− | the | + | |
|- | |- | ||
Line 461: | Line 404: | ||
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| VCR control click PLAY button | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| VCR control click PLAY button | ||
− | | 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 top of the ''' | + | | 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 top of the '''ParaView''' window, click on the '''play''' button of the '''VCR control'''. |
− | + | You can see the final value of the '''velocity magnitude'''. | |
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
| style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Color legend from top left | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Color legend from top left | ||
− | | 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;"| Also | + | | 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;"| Also toggle on the '''color legend''' from the left hand side top of the '''active variable control menu'''. |
− | + | Click '''Apply''' again. | |
− | + | ||
− | + | ||
− | + | ||
|- | |- | ||
Line 484: | Line 421: | ||
click on rescale to data | click on rescale to data | ||
− | | 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 go to ''' | + | | 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 go to '''Display'''. |
− | Scroll down | + | Scroll down. You can see '''Rescale''', click on it. |
− | + | ||
− | + | ||
− | You can see '''Rescale''' , 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;"| Check the color legend for this | | style="border-top:none;border-bottom:0.05pt solid #000000;border-left:0.05pt solid #000000;border-right:none;padding:0.097cm;"| Check the color legend for 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;"| We can see that once the '''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;"| We can see that once the '''flow''' has fully developed, it attains a '''maximum uniform velocity''' at the '''center'''. |
− | + | ||
− | + | ||
− | it attains a '''maximum uniform velocity''' at the '''center''' | + | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
+ | 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 | + | | 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 the '''analytical solution''' for '''laminar flow''' in a |
+ | '''channel''' which is '''u(max)=1.5 Uavg''' | ||
− | + | Using '''OpenFOAM''' we obtain a result of '''u(max)''' equal to '''1.48 meters per second''' which is a good match. | |
− | + | This brings us to the end of the tutorial. | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | This brings us to the end of the tutorial | + | |
|- | |- | ||
− | | 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 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 | ||
− | + | *The '''file structure of channel''' | |
− | The ''' | + | *Obtained '''solution''' using '''steady state solver''' |
− | + | *Viewed the '''geometry''' in '''ParaView ''' | |
− | + | *'''Validation''' with '''analytic results''' | |
− | ''' | + | |
− | + | ||
− | + | ||
− | Viewed the geometry in ''' | + | |
− | + | ||
− | + | ||
− | '''Validation''' with '''analytic results''' | + | |
|- | |- | ||
− | | 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 15 : |
Assignment | Assignment | ||
| 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 an assignment: | | 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 an assignment: | ||
− | Solve the problem for '''Reynold Number 1500''' and | + | *Solve the problem for '''Reynold Number 1500''' and |
− | + | *'''Validate''' it with the '''analytical result''' | |
− | ''' | + | |
|- | |- | ||
− | | 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 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;"| | + | | 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. |
− | + | ||
− | + | ||
It summarizes the Spoken Tutorial project. | It summarizes the Spoken Tutorial project. | ||
Line 558: | Line 471: | ||
|- | |- | ||
− | | 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 17 : Spoken Tutorial Workshops |
− | + | ||
− | + | ||
| 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 569: | Line 480: | ||
-For more details, please write to [mailto:contact@spoken-tutorial.org contact@spoken-tutorial.org] | -For more details, please write to [mailto:contact@spoken-tutorial.org contact@spoken-tutorial.org] | ||
+ | |- | ||
+ | | 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 | ||
+ | * 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/ | ||
+ | | 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 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 | ||
+ | | 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 | |
+ | * 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/ | ||
+ | | 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 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/ | ||
+ | | 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 22: | ||
+ | |||
+ | Acknowledgement | ||
Line 583: | Line 532: | ||
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 | + | More information on the same is available at the following URL http://spoken-tutorial.org/NMEICT-Intro |
|- | |- |
Latest revision as of 07:44, 13 June 2019
Tutorial: Simulating Laminar flow in a channel.
Script and Narration : Rahul Joshi
Keywords: Video tutorial,CFD,laminar flow,simpleFoam,channel
|
|
Slide 1 | Hello and welcome to the spoken tutorial on Simulating 2D Laminar Flow in a Channel using openfoam |
Slide 2 : Learning Objectives | In this tutorial I will show you
|
Slide 3 : System requirement | To record this tutorial I am using
|
Slide 4 : System Requirement
|
|
Slide 5 : Prerequisite | As a prerequisite for this tutorial, you should know how to create geometry using OpenFOAM.
|
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
|
Slide 8 : Boundary Conditions | Using the formula the length of the channel comes out to be 5 meters
|
Slide 9 : File structure | This is a steady state problem
|
Create a folder and name it as channel in simpleFoam folder.
|
I have already created a folder in
|
Nmae the folder channel | The folder is named as channel.
|
copy 0,constant and system from pitzDaily to this channel folder | Copy 0, Constant and System folders of any other case. |
In the simpleFoam directory, I have copied the file structure of case of pitzDaily. | |
Paste it in the channel folder.
| |
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 space 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 space 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.
|
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.
|
In the blockMeshDict file covertTometers is set to 1
|
The geometry is in meters so the convertTometers is set to 1.
|
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 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 be kept empty.
|
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 terminal type cd 0
|
Now in the terminal type cd space 0 (Zero) and press Enter.
|
In terminal you can see these slides | This contains the initial boundary conditions and wall functions for the channel case. |
Wall functions : epsilon, k, nut, nutilda
|
It should contain various files such as epsilon ,k, nut,nuTilda which are the wall functions.
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.
|
Slide 11 : Calculate epsilon
|
Calculate epsilon from the formula given:
|
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 kept to be 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 are no changes to be done in the system folder. |
Mesh the geometry
|
Now, we need to mesh the geometry.
|
Slide 12 : Solver | The type of solver we are using here is simpleFoam.
|
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.
|
In paraview window
|
On left hand side of the ParaView window click Apply.
|
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.
|
Color legend from top left | Also toggle on the color legend from the left hand side top of the active variable control menu.
|
In object inspector menu
|
Now go to Display.
|
Check the color legend for this | We can see that once the flow has fully developed, it attains a maximum uniform velocity at the center.
|
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
|
Slide 14 : Summary | In this tutorial we learnt
|
Slide 15 :
Assignment |
As an assignment:
|
Slide 16 : About the Spoken Tutorial Project | Watch the video available at this URL.
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
|
|
Slide 19:
Forum to answer questions
|
|
Slide 20:
Lab Migration Project
For more details visit this site: http://cfd.fossee.in/ |
|
Slide 21:
Case Study Project
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 |