Difference between revisions of "OpenFOAM/C3/Flow-over-a-flat-plate/English"
Line 11: | Line 11: | ||
|- | |- | ||
| Slide 1 | | Slide 1 | ||
− | | Hello and welcome to the spoken tutorial on Flow over a flat plate using '''OpenFOAM''' | + | | Hello and welcome to the spoken tutorial on Flow over a flat plate using '''OpenFOAM''' |
|- | |- | ||
Line 20: | Line 20: | ||
* Changing the grid spacing in '''meshing''' | * Changing the grid spacing in '''meshing''' | ||
* Post Processing results in '''ParaView''' and | * Post Processing results in '''ParaView''' and | ||
− | * Visualizing using '''Vector | + | * Visualizing using '''Vector Plots''' |
|- | |- | ||
Line 29: | Line 29: | ||
I am using | I am using | ||
− | * '''Linux''' Operating system '''Ubuntu''' version 12.04 | + | * '''Linux''' Operating system '''Ubuntu''' version 12.04 |
*'''OpenFOAM''' version 2.1.1 and | *'''OpenFOAM''' version 2.1.1 and | ||
*'''ParaView''' version 3.12.0 | *'''ParaView''' version 3.12.0 | ||
Line 53: | Line 53: | ||
Flow over a flat plate diagram | Flow over a flat plate diagram | ||
− | | We can visualise the growth of '''boundary layer''' | + | | We can visualise the growth of the '''boundary layer''' |
'''Boundary layer''' is a very thin region above the body | '''Boundary layer''' is a very thin region above the body | ||
− | where the velocity is 0.99 times the '''free stream velocity''' | + | where the velocity is 0.99 times the '''free stream velocity''' |
|- | |- | ||
| Slide 7: Diagram of boundary conditions. | | Slide 7: Diagram of boundary conditions. | ||
− | | This is a '''diagram''' of '''flow over the flat plate''' | + | | This is a '''diagram''' of the '''flow over the flat plate''' |
The''' boundary conditions''' are given as follows | The''' boundary conditions''' are given as follows | ||
Line 67: | Line 67: | ||
*You have the '''Inlet''' | *You have the '''Inlet''' | ||
*the '''Plate''' | *the '''Plate''' | ||
− | *'''Top''' | + | *'''Top''' which is the '''Farfield''' |
− | *and '''Outlet''' | + | *and '''Outlet''' which is the''' pressure outlet boundary''' |
|- | |- | ||
| Slide 8: Inlet parameters | | Slide 8: Inlet parameters | ||
| | | | ||
− | * The Free stream velocity '''U | + | * The Free stream velocity '''U is 1 m/s''', and |
− | * we are solving this for '''Reynolds | + | * we are solving this for a '''Reynolds number (Re) = 100''' |
|- | |- | ||
Line 82: | Line 82: | ||
Then go to the '''Run''' '''directory ''' | Then go to the '''Run''' '''directory ''' | ||
− | You will see '''Tutorials. '''Click on it | + | You will see '''Tutorials. '''Click on it |
− | '''Scroll down''' and then click on '''Incompressible | + | '''Scroll down''' and then click on '''Incompressible''' |
− | '''Scroll down''' | + | '''Scroll down''' |
− | You will see the '''simpleFoam '''folder | + | You will see the '''simpleFoam '''folder |
Click on it | Click on it | ||
− | This '''solver''' suits our case | + | This '''solver''' suits our case |
|- | |- | ||
| Right click >> Create new folder >> flatplate | | Right click >> Create new folder >> flatplate | ||
− | | In this, create a '''folder''' by the name '''flatplate | + | | In this, create a '''folder''' by the name '''flatplate''' |
− | Right click | + | Right click '''Create New Folder ''' '''flatplate''' |
|- | |- | ||
| Double-click '''pitzdaily '''folder | | Double-click '''pitzdaily '''folder | ||
− | | Now, let | + | | Now, let us open the '''pitzdaily case''' |
|- | |- | ||
Line 110: | Line 110: | ||
|- | |- | ||
| Select '''0''', '''constant '''and '''system '''folders | | Select '''0''', '''constant '''and '''system '''folders | ||
− | | Copy the '''three folders''' | + | | Copy the '''three folders''' '''0''', '''constant '''and '''system''' |
|- | |- | ||
Line 120: | Line 120: | ||
Paste them there | Paste them there | ||
− | | Now let us go one level back | + | | Now let us go one level back |
− | Paste these three '''folders''' inside the''' flatplate '''folder | + | Paste these three '''folders''' inside the''' flatplate '''folder |
|- | |- | ||
Line 130: | Line 130: | ||
|- | |- | ||
| | | | ||
− | | Change the geometry and boundary condition names in the '''blockMeshDict '''file | + | | Change the geometry and boundary condition names in the '''blockMeshDict '''file |
|- | |- | ||
Line 138: | Line 138: | ||
Let us open the '''blockMeshDict''' file . Scroll down | Let us open the '''blockMeshDict''' file . Scroll down | ||
− | The geometry is in '''meters''' | + | The geometry is in '''meters''' |
|- | |- | ||
Line 146: | Line 146: | ||
|- | |- | ||
| Simplegrading (1 3 1) | | Simplegrading (1 3 1) | ||
− | | | + | | You can see the '''simpleGrading ''' |
− | It is kept as (1 3 1) as we need a finer '''mesh''' near the plate | + | It is kept as (1 3 1) as we need a finer '''mesh''' near the plate |
|- | |- | ||
| | | | ||
− | | Now close this | + | | Now close this |
− | Go two '''levels''' back | + | Go two '''levels''' back |
|- | |- | ||
| | | | ||
− | | Similarly, make changes in the '''boundary condition''' names inside the '''files''' in the '''0''' folder | + | | Similarly, make changes in the '''boundary condition''' names inside the '''files''' in the '''0''' folder |
− | These '''files''' have '''pressure''' | + | These '''files''' have '''pressure''' '''velocity''' and '''wall''' functions |
|- | |- | ||
| | | | ||
− | | To '''calculate''' the values of '''wall''' functions | + | | To '''calculate''' the values of '''wall''' functions |
− | please refer to the earlier '''tutorials''' in the '''OpenFoam''' series | + | please refer to the earlier '''tutorials''' in the '''OpenFoam''' series |
− | Let us go one '''level''' back | + | Let us go one '''level''' back |
|- | |- | ||
Line 181: | Line 181: | ||
Press Ctrl+Alt+t keys simultaneously | Press Ctrl+Alt+t keys simultaneously | ||
− | | Now let us open the '''terminal window | + | | Now let us open the '''terminal window''' |
− | In the '''terminal window''', type '''run '''and press '''Enter | + | In the '''terminal window''', type '''run '''and press '''Enter''' |
|- | |- | ||
| Type cd tutorials | | Type cd tutorials | ||
− | | | + | | Now type '''cd space tutorials '' press''' Enter''' |
|- | |- | ||
| Type cd incompressible | | Type cd incompressible | ||
− | | | + | | Now type '''cd incompressible '''press''' Enter''' |
|- | |- | ||
| Type cd simpleFoam | | Type cd simpleFoam | ||
− | | | + | | Now type '''cd space simpleFoam '''press''' Enter''' |
|- | |- | ||
| Type ls | | Type ls | ||
− | | Now type '''ls''' and press '''Enter | + | | Now type '''ls''' and press '''Enter''' |
|- | |- | ||
Line 207: | Line 207: | ||
|- | |- | ||
| Type cd flatplate | | Type cd flatplate | ||
− | | Now type '''cd space flatplate '''and press '''Enter | + | | Now type '''cd space flatplate '''and press '''Enter''' |
|- | |- | ||
| Type ls | | Type ls | ||
− | | Now type '''ls''' and press '''Enter | + | | Now type '''ls''' and press '''Enter''' |
|- | |- | ||
| | | | ||
− | | You can see the three '''folders''' '''0 | + | | You can see the three '''folders''' '''0 constant''' and '''system''' |
|- | |- | ||
| Type blockMesh | | Type blockMesh | ||
− | | Now, we will mesh the geometry. | + | | Now, we will mesh the geometry. |
− | We are using a '''course mesh''' for this problem | + | We are using a '''course mesh''' for this problem |
− | '''Meshing''' can be done by typing '''blockMesh''' in the '''terminal''' | + | '''Meshing''' can be done by typing '''blockMesh''' in the '''terminal''' |
|- | |- | ||
| | | | ||
− | | Press '''Enter | + | | Press '''Enter''' |
− | '''Meshing''' has been done | + | '''Meshing''' has been done |
|- | |- | ||
| | | | ||
− | | Note that if there is some error in the '''blockMesh''' file | + | | Note that if there is some error in the '''blockMesh''' file |
− | it will be shown in the '''terminal''' window | + | it will be shown in the '''terminal''' window |
|- | |- | ||
Line 241: | Line 241: | ||
| To view the geometry, type “'''paraFoam'''” | | To view the geometry, type “'''paraFoam'''” | ||
− | press '''Enter | + | And press '''Enter''' |
|- | |- | ||
| Paraview window click on APPLY button on left hand side | | Paraview window click on APPLY button on left hand side | ||
− | | After the '''ParaView''' window opens, on the left hand side of the '''object inspector''' menu, click '''Apply''' | + | | After the '''ParaView''' window opens, on the left hand side of the '''object inspector''' menu, click '''Apply''' |
− | We can see the '''geometry''' | + | We can see the '''geometry''' |
− | Close the '''ParaView''' window | + | Close the '''ParaView''' window |
− | Let me switch back to the''' slides | + | Let me switch back to the''' slides''' |
|- | |- | ||
| Slide 9: solver | | Slide 9: solver | ||
− | | The solver we are using here is | + | | The solver we are using here is '''simpleFoam''' |
'''SimpleFoam''' is a '''steady state solver '''for '''incompressible ''' | '''SimpleFoam''' is a '''steady state solver '''for '''incompressible ''' | ||
Line 266: | Line 266: | ||
type simpleFoam | type simpleFoam | ||
− | | Let me switch back to the '''terminal '''window | + | | Let me switch back to the '''terminal '''window |
− | In the '''terminal window''' ,type '''simpleFoam''' and press '''Enter | + | In the '''terminal window''' ,type '''simpleFoam''' and press '''Enter''' |
− | You will see the '''iterations''' running in the '''terminal '''window | + | You will see the '''iterations''' running in the '''terminal '''window |
|- | |- | ||
| Type paraFoam | | Type paraFoam | ||
− | | Once the solving is done, type '''paraFoam''' to view the results | + | | Once the solving is done, type '''paraFoam''' to view the results |
|- | |- | ||
| In the Paraview window click on APPLY button on left hand side | | In the Paraview window click on APPLY button on left hand side | ||
− | | On the left hand side of the '''Object Inspector''' menu, click '''Apply to '''view the geometry | + | | On the left hand side of the '''Object Inspector''' menu, click '''Apply to '''view the geometry |
|- | |- | ||
Line 299: | Line 299: | ||
|- | |- | ||
| VCR control | | VCR control | ||
− | | Now on top of the '''ParaView''' window, you will see the '''VCR''' | + | | Now on top of the '''ParaView''' window, you will see the '''VCR''' control |
|- | |- | ||
| Click on Play button of VCR control | | Click on Play button of VCR control | ||
− | | Click on the '''Play''' button | + | | Click on the '''Play''' button |
|- | |- | ||
Line 347: | Line 347: | ||
|- | |- | ||
| | | | ||
− | | Click '''Apply''' on the left hand side of '''Object Inspector Menu''' | + | | Click '''Apply''' on the left hand side of '''Object Inspector Menu''' |
|- | |- | ||
| Changing vector size | | Changing vector size | ||
− | | You can change the number of '''vectors''' by changing their size at the bottom | + | | You can change the number of '''vectors''' by changing their size at the bottom |
|- | |- | ||
Line 357: | Line 357: | ||
set scale factor 0.1 | set scale factor 0.1 | ||
− | | Also, the size of the '''vectors''' can be changed by '''clicking''' on the '''Edit''' button | + | | Also, the size of the '''vectors''' can be changed by '''clicking''' on the '''Edit''' button |
The '''set scale''' '''factor''' can be changed to '''0.1''' | The '''set scale''' '''factor''' can be changed to '''0.1''' | ||
Line 363: | Line 363: | ||
|- | |- | ||
| Click the apply button | | Click the apply button | ||
− | | Again, click the '''Apply''' button | + | | Again, click the '''Apply''' button |
|- | |- | ||
Line 371: | Line 371: | ||
|- | |- | ||
| Click on ZoomToBox icon | | Click on ZoomToBox icon | ||
− | | To do this, in the '''Active Variable Control '''menu, click on '''zoomToBox '''option | + | | To do this, in the '''Active Variable Control '''menu, click on the '''zoomToBox '''option |
|- | |- | ||
Line 379: | Line 379: | ||
|- | |- | ||
| Parabolic variation of vector plot | | Parabolic variation of vector plot | ||
− | | We can see the '''parabolic variation''' of '''vector | + | | We can see the '''parabolic variation''' of '''vector plots '''as the '''flow''' moves over the''' plate''' |
|- | |- | ||
| Delete the vector plot | | Delete the vector plot | ||
− | | Delete this. Now delete the '''vector plot | + | | Delete this. Now delete the '''vector plot''' |
|- | |- | ||
Line 389: | Line 389: | ||
| Also, we can see that the color near to 1 | | Also, we can see that the color near to 1 | ||
− | corresponds to the '''velocity''' of '''0.99 '''times the free stream velocity | + | corresponds to the '''velocity''' of '''0.99 '''times the free stream velocity |
|- | |- | ||
| To plot the data along x and y axis | | To plot the data along x and y axis | ||
− | | | + | | We can also plot the '''variation''' of velocity along the '''x '''and '''y '''axes using the '''plot data over line''' |
|- | |- | ||
| Slide 10: Summary | | Slide 10: Summary | ||
− | | This brings us to the end of the tutorial | + | | This brings us to the end of the tutorial |
In this tutorial we learnt : | In this tutorial we learnt : | ||
Line 408: | Line 408: | ||
| As an Assignment, | | As an Assignment, | ||
− | Create a '''geometry''' of '''flow over | + | Create a '''geometry''' of '''flow over the flat plate''' |
'''Refine the grid spacing''' near the plate | '''Refine the grid spacing''' near the plate | ||
Line 433: | Line 433: | ||
-For more details, please write to | -For more details, please write to | ||
− | contact | + | contact(at the rate)spoken(hyphen)tutorial(dot)org |
|- | |- | ||
Line 476: | Line 476: | ||
| Spoken Tutorial project is a part of the Talk to a Teacher project, | | Spoken Tutorial project is a part of the Talk to a Teacher project, | ||
− | It is supported by the National Mission on Education through ICT, MHRD, Government of India | + | It is supported by the National Mission on Education through ICT, MHRD, Government of India |
More information on this mission is available at this URL http://spoken-tutorial.org/NMEICT-Intro | More information on this mission is available at this URL http://spoken-tutorial.org/NMEICT-Intro | ||
Line 482: | Line 482: | ||
|- | |- | ||
| About the contributor | | About the contributor | ||
− | | This is Rahul Joshi from IIT BOMBAY signing off | + | | This is Rahul Joshi from IIT BOMBAY signing off |
− | Thanks for joining | + | Thanks for joining |
|} | |} |
Revision as of 11:15, 30 July 2019
Tutorial: Flow over a flat plate using OpenFOAM.
Script and Narration: Rahul Joshi
Keywords: Video tutorial, CFD, Flat plate, Boundary layer, glyph (vector plotting).
Visual Cue | Narration |
---|---|
Slide 1 | Hello and welcome to the spoken tutorial on Flow over a flat plate using OpenFOAM |
Slide 2 : Learning Objectives | In this tutorial I will teach you about
|
Slide 3:
System Requirement |
To record this tutorial
I am using
|
Slide 4:
System Requirement
|
|
Slide 5: Flow over Flat Plate | Flow over flat plate is a fundamental problem in fluid mechanics |
Slide 6: Flow over Flat Plate
Flow over a flat plate diagram |
We can visualise the growth of the boundary layer
Boundary layer is a very thin region above the body where the velocity is 0.99 times the free stream velocity |
Slide 7: Diagram of boundary conditions. | This is a diagram of the flow over the flat plate
The boundary conditions are given as follows
|
Slide 8: Inlet parameters |
|
Click on home>> OpenFoam >> Run >> Tutorials >> Incompressible >> SimpleFoam | Now let us Go to the Home folder, In the Home folder, click on the OpenFoam folder
Then go to the Run directory You will see Tutorials. Click on it Scroll down and then click on Incompressible Scroll down You will see the simpleFoam folder Click on it This solver suits our case |
Right click >> Create new folder >> flatplate | In this, create a folder by the name flatplate
Right click Create New Folder flatplate |
Double-click pitzdaily folder | Now, let us open the pitzdaily case |
Zoom in | Let me zoom this. |
Select 0, constant and system folders | Copy the three folders 0, constant and system |
Copy them | Copy this |
Go back to SimpleFoam >> flatplate folder
Paste them there |
Now let us go one level back
Paste these three folders inside the flatplate folder |
Click on constant >> polyMesh | Open the constant folder and then the polyMesh folder |
Change the geometry and boundary condition names in the blockMeshDict file | |
I have already made the changes.
Let us open the blockMeshDict file . Scroll down The geometry is in meters | |
We have set the dimensions of the flatplate | |
Simplegrading (1 3 1) | You can see the simpleGrading
It is kept as (1 3 1) as we need a finer mesh near the plate |
Now close this
Go two levels back | |
Similarly, make changes in the boundary condition names inside the files in the 0 folder
These files have pressure velocity and wall functions | |
To calculate the values of wall functions
Let us go one level back | |
The system folder can be kept default
Let us close this | |
Let us Open the terminal window :
Press Ctrl+Alt+t keys simultaneously |
Now let us open the terminal window
In the terminal window, type run and press Enter |
Type cd tutorials | Now type cd space tutorials press' Enter |
Type cd incompressible | Now type cd incompressible press Enter |
Type cd simpleFoam | Now type cd space simpleFoam press Enter |
Type ls | Now type ls and press Enter |
We can see the flatplate folder. | |
Type cd flatplate | Now type cd space flatplate and press Enter |
Type ls | Now type ls and press Enter |
You can see the three folders 0 constant and system | |
Type blockMesh | Now, we will mesh the geometry.
We are using a course mesh for this problem Meshing can be done by typing blockMesh in the terminal |
Press Enter
Meshing has been done | |
Note that if there is some error in the blockMesh file
it will be shown in the terminal window | |
Type paraFoam | To view the geometry, type “paraFoam”
And press Enter |
Paraview window click on APPLY button on left hand side | After the ParaView window opens, on the left hand side of the object inspector menu, click Apply
We can see the geometry Close the ParaView window Let me switch back to the slides |
Slide 9: solver | The solver we are using here is simpleFoam
SimpleFoam is a steady state solver for incompressible and turbulent flows |
Demo :
type simpleFoam |
Let me switch back to the terminal window
In the terminal window ,type simpleFoam and press Enter
|
Type paraFoam | Once the solving is done, type paraFoam to view the results |
In the Paraview window click on APPLY button on left hand side | On the left hand side of the Object Inspector menu, click Apply to view the geometry |
Properties | Scroll down the properties panel of the Object Inspector menu for time step, regions and fields |
Change the drop down menu from Solid Color to U | To view the contours from the drop down menu,
in the Active Variable Control menu, change from solid color to capital U |
You can see the initial condition of the velocity | |
VCR control | Now on top of the ParaView window, you will see the VCR control |
Click on Play button of VCR control | Click on the Play button |
You will see the contour of Pressure or Velocity on the flat plate accordingly | |
Toggle on the Color legend | This is the velocity contour
Toggle on the Color legend |
Color legend left hand side top icon | To do this, click on the color legend icon on the Active Variable Control menu |
Click on APPLY button | Click Apply in the Object inspector menu |
Click on Display | In the Object inspector menu, click on Display |
Click on rescale to data range | Scroll down and click on Rescale to data range |
Shift color legend on top of the geometry | Let me shift this Color legend on top |
Top menu >> Filter > Common > glyph | To visualize the Vector Plot,
go to the Filters Menu > Common > glyph |
Go to the Properties in Object Inspector menu | |
Click Apply on the left hand side of Object Inspector Menu | |
Changing vector size | You can change the number of vectors by changing their size at the bottom |
Scroll down and click on edit button
set scale factor 0.1 |
Also, the size of the vectors can be changed by clicking on the Edit button
The set scale factor can be changed to 0.1 |
Click the apply button | Again, click the Apply button |
Now let me zoom this | |
Click on ZoomToBox icon | To do this, in the Active Variable Control menu, click on the zoomToBox option |
And zoom over any area that you desire | |
Parabolic variation of vector plot | We can see the parabolic variation of vector plots as the flow moves over the plate |
Delete the vector plot | Delete this. Now delete the vector plot |
Corresponding to color of 1 in color legend | Also, we can see that the color near to 1
corresponds to the velocity of 0.99 times the free stream velocity |
To plot the data along x and y axis | We can also plot the variation of velocity along the x and y axes using the plot data over line |
Slide 10: Summary | This brings us to the end of the tutorial
In this tutorial we learnt :
|
Slide 11: Assignment | As an Assignment,
Create a geometry of flow over the flat plate Refine the grid spacing near the plate |
Slide 12 : About Spoken tutorials |
http://spoken-tutorial.org/What_is_a_Spoken_Tutorial
|
Slide 13: Spoken Tutorial Worekshops | The Spoken Tutorial Project Team
-Conducts workshops using spoken tutorials -Gives certificates to those who pass an online test -For more details, please write to contact(at the rate)spoken(hyphen)tutorial(dot)org |
Slide 14: Forum to answer questions
|
|
Slide 15: Forum to answer questions
|
|
Slide 16: Lab Migration Project
For more details visit this site: http://cfd.fossee.in/ |
|
Slide 17: Case Study Project
For more details visit this site: http://cfd.fossee.in/ |
|
Slide 18:
Acknowledgement
|
Spoken Tutorial project is a part of the Talk to a Teacher project,
It is supported by the National Mission on Education through ICT, MHRD, Government of India More information on this mission is available at this URL http://spoken-tutorial.org/NMEICT-Intro |
About the contributor | This is Rahul Joshi from IIT BOMBAY signing off
Thanks for joining |