OpenFOAM/C3/Simulating-Hagen-Poiseuille-flow/English-timed
From Script | Spoken-Tutorial
Revision as of 15:26, 1 December 2015 by Pratik kamble (Talk | contribs)
Time | Narration |
00:02 | Hello and welcome to the spoken tutorial on simulating Hagen-Poiseuille flow in OpenFOAM. |
00:09 | In this tutorial we will see
To create and mesh 3D cylindrical pipe. To simulate the Hagen-Poiseuille flow having fixed pressure ratio across boundaries. and To visualize the velocity contour in ParaView. |
00:25 | To record this tutorial, I am using
Linux Operating system Ubuntu 12.04 OpenFOAM version 2.1.1 and ParaView version 3.12.0 |
00:38 | To practice this tutorial learner should have the knowledge of Basic Fluid Dynamicsand Hagen-Poiseuille flow |
00:46 | Here is, Hagen-Poiseuille Flow Diagram. We can see the dimensions and boundaries of the pipe. |
00:52 | Viscosity of the fluid used, that is, water is given.Pressure at the inlet is 20 Pascals and at the outlet is 0 Pascals. |
01:04 | As it is an incompressible flow, only the pressure difference is of importance. |
01:10 | Formulas and Analytical Solution: For Hagen-Poiseuille flow,'Pressure drop along the pipe is: ('P1 minus P2) equals (32 mew Uaverage L) upon (D square) |
01:25 | By substituting the values from the previous diagram, we get,'Uaverage equals to 0.208 m/s'Maximum Velocity is given as,Two times the average velocity, which would be, 0.416 m/s |
01:44 | Reynolds Number for the flow is,Uaverage into D upon nu, that comes out to be, 2080Hence, the flow is transient. |
01:56 | Type of solver used here is,IcoFOAM |
02:01 | It is a Transient Solver It is used for incompressible, laminar flow of Newtonian fluids. |
02:08 | Pressure Boundary Conditions used,At Inlet: fixedPressureAt Outlet: fixedPressure
At Walls: ZeroGradient |
02:19 | Velocity Boundary Conditions used,At Inlet: pressureInletVelocity At Outlet: zeroGradient At Walls: fixedValue |
02:28 | For executing this case,First, Let's create the case directory in the 'icoFoam' folder.
And Give it some name.I have named it as '3dpipe'. |
02:41 | To know the location of this folder, go through the tutorial on lid driven cavity.Copy this '0' (zero), 'constant' and 'system' folders of lid driven cavity problem in the newly created folder. |
02:54 | Let's go inside the '3dpipe' folder. |
02:58 | I have already copied the folders into my 3dpipe folder and modified the files in it. |
03:05 | Now, let's go into the '0' folder.And open the 'P' file.This is the pressure boundary condition file. |
03:14 | Note that the dimensions are in (meter square) per (second square) (m2/s2). |
03:20 | Hence the pressure value in pascals is divided by the density, that is, 1000 Kg/m3 (Kg per meter cube), and written here. |
03:29 | Let's close the file. |
03:32 | File containing the velocity boundary condition is as seen:lets open the file we can see the velocity boundary condition for inlet, outlet and fixed walls |
03:43 | Let's close the file and come out of the '0' folder. |
03:48 | To see the blocking strategy, let me switch back to the slides. |
03:54 | To create a 3D geometry of a pipe I have made a 2D circular geometry and extruded the length in z direction. |
04:03 | Numbering Pattern is as shown. You can also see the dimension of the mesh. |
04:11 | To see the blockMeshDict file, let's minimize the slides. |
04:16 | Let's go into the folder 'constant', and then 'polyMesh'. lets open theblockMeshDict ' file. You can see the vertices, logs,edges and boundaries for inlet, outlet and fixed wall. |
04:37 | Let's close the file and lets come out of the polyMesh folder. |
04:42 | We see the 'transportProperties' file. Lets open the file Note the dynamic viscosity value, here, is 1e-06. |
04:53 | Let's close the file and come out of the 'folder ''constant' . |
04:59 | Let's go into the 'system' folder.Now, let's have a look at the 'controlDict' file. |
05:07 | The solution converges after 18 seconds therefore the final time step is kept 19.The time step has been set to 1e-03. |
05:20 | Let's close the file.Let's close the 'Home' folder. |
05:26 | Now to execute the case, we will, first, go inside the '3dpipe' folder through terminal.Let's open the terminal by pressing 'control', 'alt' and 't' key, altogether. |
05:40 | Type run and press Enter |
05:44 | Type cd (space) tutorials and press Enter |
05:50 | Type cd (space) incompressible and press Enter |
05:55 | Type cd (space) icoFoam and press Enter |
05:59 | Type cd (space) 3Dpipe and press Enter |
06:05 | Now to create the mesh, type blockMesh and press Enter. Meshing has been done. |
06:16 | To start the iterations type icoFoam and press Enter. We see the iterations are running. |
06:27 | Iterations has been done.After the iterations end type paraFoam for postprocessing the results and press Enter. It will open the" paraview". This is " paraview" |
06:41 | Let's click on Apply on the left hand side of the Object inspector menu to see the geometry. |
06:49 | Let's rotate the geometry for a better view. |
06:52 | Click on the active variable control menu and select U in the drop-down menu. |
07:01 | At the top, in VCR toolbar, click on Play button. |
07:06 | Go to Object Inspector menu, go to Display, click on Rescale to data range. |
07:16 | To view the half section, go to the toolbar named common, click on Clips go to object inspector menu properties and press Apply.Lets Zoom in |
07:35 | Let's open the color legend. |
07:38 | We can see the maximum velocity is near to the actual maximum velocityi.e 0.4 metersp/s |
07:46 | To view the graph Go to Filtersat the top Data Analysisand press Plot Over Lines. |
07:56 | Press Y axis and press Apply. |
08:00 | We can see the parabolic profile for Hagen-Poiseuille flow. |
08:05 | Let's close the graph. lets Close ParaView.And switch to the slides. |
08:12 | In this tutorial we have learned:
|
08:30 | As an assignment,
Change the geometry parameters such as length and diameter. Change the corresponding pressure ratio. and Use the fluid of different viscosity. |
08:43 | Watch the video available at the following link * It summarises the Spoken Tutorial project.If you do not have good bandwidth, you can download and watch it |
08:54 | The Spoken Tutorial Project Team Conducts workshops using spoken tutorials Gives certificates for those who pass an online test For more details, please write to contact at spoken hyphen tutorial dot org |
09:11 | 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 spoken hyphen tutorial dot org slash NMEICT hyphen Intro |