CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Laminar , steady state pipe flow

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 22, 2009, 04:36
Question Laminar , steady state pipe flow
  #1
andesameer
Guest
 
Posts: n/a
I would like to solve steady state, laminar flow of newtonian fluid in cylindrical pipe. Can you please suggest which solver would be more appropriate.

Thank a lot

Sam
  Reply With Quote

Old   June 22, 2009, 21:55
Thumbs up laminar solver
  #2
New Member
 
A B SUDHARSAN
Join Date: Jun 2009
Posts: 6
Rep Power: 9
sudhar is on a distinguished road
hi sam ,
the simple and the most basic solver for a laminar and incompressible flow is icoFoam.
sudhar is offline   Reply With Quote

Old   June 22, 2009, 23:16
Default
  #3
andesameer
Guest
 
Posts: n/a
Hi Sudharshan,

Thank you for your reply.

icofoam solver solves transient, laminar and incompressible newtonian flow. But I am looking to solve for steady state, laminar and incompressible newtonian flow.

Can I use simplefoam to solve laminar flow by changing code?
  Reply With Quote

Old   June 23, 2009, 04:30
Default
  #4
Senior Member
 
santos's Avatar
 
Jose Luis Santos
Join Date: Mar 2009
Location: Portugal
Posts: 215
Rep Power: 10
santos is on a distinguished road
Send a message via Skype™ to santos
Hi,

To use simpleFoam for laminar flow, you have to edit constant/RASProperties of your case and change RASModel to laminar.

Regards,
Jose Santos
santos is offline   Reply With Quote

Old   June 23, 2009, 13:36
Default
  #5
Member
 
Sven Winkler
Join Date: May 2009
Posts: 70
Rep Power: 9
sven is on a distinguished road
You can also use icoFoam. Although the solver is transient, you will get an steady-state solution after some time, since the flow itself becomes steady-state. as far as I know, a steady-state solver does exactly the same, it solves the equations with the time derivatives and the solution proceeds in time until a steady-state is reached. I did this for a rectangular channel flow and it worked well.
sven is offline   Reply With Quote

Old   June 24, 2009, 21:57
Post
  #6
New Member
 
A B SUDHARSAN
Join Date: Jun 2009
Posts: 6
Rep Power: 9
sudhar is on a distinguished road
hi sam,
you can use simpleFoam for your case but need to edit RAS properties before using. Good luck.

regards,
sudharsan
sudhar is offline   Reply With Quote

Old   June 24, 2009, 23:04
Default
  #7
andesameer
Guest
 
Posts: n/a
Thank you to all.

I will solve the my problem using simplefoam. I will let u know if I face any difficulty.

Thank you.
Sam
  Reply With Quote

Old   March 15, 2012, 10:57
Default
  #8
Senior Member
 
Goutam Saha
Join Date: Dec 2011
Location: UK
Posts: 131
Rep Power: 6
Goutam is on a distinguished road
Quote:
Originally Posted by andesameer View Post
Thank you to all.

I will solve the my problem using simplefoam. I will let u know if I face any difficulty.

Thank you.
Sam
Dear Sam

I am trying to solve a steady state incompressible laminar flow for a pipe with circular cross section, but I didn't get the fully developed parabolic shape for velocity profile.
I am using simpleFOAM (OF-2.1.0) where I have changed the RAS properties to laminar and set the turbulent off. Could you please suggest me, where is the problem? Do you able to solve this simple problem?

If you give me your email ID, I can send you the files.

Thanks

Goutam
Goutam is offline   Reply With Quote

Old   March 15, 2012, 18:24
Default
  #9
New Member
 
Lachlan Graham
Join Date: Dec 2009
Posts: 14
Rep Power: 8
Locky3827 is on a distinguished road
Hi Sam,
I have used icoFoam for this problem and nonNewtonianicoFoam for the corresponding non-Newtonian case. Both work very well, it is also easy to cross check with the analytical values.
Regards,
Lachlan
Locky3827 is offline   Reply With Quote

Old   March 16, 2012, 00:38
Default
  #10
Member
 
B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 54
Rep Power: 7
skyinventorbt is on a distinguished road
switch off turbulence , use laminar in the RAS file of /const floder
skyinventorbt is offline   Reply With Quote

Old   March 27, 2013, 10:53
Default
  #11
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
Hi,

I would like to know, if it is possible to solve a laminar flow in a cylindrical pipe using 2d simulations in openfoam. If so, please can anyone let me know, how to construct the geometry. Further, should we change anything in the icoFoam solver.

Thanks

Regards
Vishal
nandiganavishal is offline   Reply With Quote

Old   March 27, 2013, 23:46
Default 2D simulations - Reg
  #12
Member
 
B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 54
Rep Power: 7
skyinventorbt is on a distinguished road
Dear Nandhigana Vishal,

You can create 2D simulations by considering a wedge shaped geometry as shown in Figure 5.3 in OpenFOAM manual 2.0.0.

For creating wedge refer Figure 5.7 in the same manual.
__________________
-------------------------------
B T KANNAN
PhD Research Scholar,
Indian Institute of Technology - Madras.
skyinventorbt is offline   Reply With Quote

Old   March 28, 2013, 02:12
Default
  #13
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
Quote:
Originally Posted by skyinventorbt View Post
Dear Nandhigana Vishal,

You can create 2D simulations by considering a wedge shaped geometry as shown in Figure 5.3 in OpenFOAM manual 2.0.0.

For creating wedge refer Figure 5.7 in the same manual.
Dear Kannan,

Thanks for your reply. I had a chance to go through the manual earlier, but could not exactly follow how the geometry and boundary conditions are incorporated using the wedge BC. For instance, I would like to create a cylinder of length 8 m and diameter 0.1 m. The inlet of the cylinder has a uniform velocity and the walls have no slip. How do we create this geometry and what angle should be specified for the wedge. Can you please illustrate this with an example along with the patch definitions.

Thanks for the help

Vishal
nandiganavishal is offline   Reply With Quote

Old   March 28, 2013, 07:15
Smile Simple solution
  #14
Member
 
B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 54
Rep Power: 7
skyinventorbt is on a distinguished road
Dear Vishal,

  1. Create an wedge (Angle of 5 Degree). Simply a hex collapsed to form a wedge as shown in manual.
  2. Axis is type empty
  3. Front and back face type wedge
  4. For velocity use type uniform and fixed value
Go through Cavity and pitzdaily examples in OpenFOAM tutorial once for understanding and kindly follow the above.
Of course you will find difficulty. But try once then you will get it.
__________________
-------------------------------
B T KANNAN
PhD Research Scholar,
Indian Institute of Technology - Madras.
skyinventorbt is offline   Reply With Quote

Old   April 1, 2013, 00:57
Default
  #15
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
Quote:
Originally Posted by skyinventorbt View Post
Dear Vishal,

  1. Create an wedge (Angle of 5 Degree). Simply a hex collapsed to form a wedge as shown in manual.
  2. Axis is type empty
  3. Front and back face type wedge
  4. For velocity use type uniform and fixed value
Go through Cavity and pitzdaily examples in OpenFOAM tutorial once for understanding and kindly follow the above.
Of course you will find difficulty. But try once then you will get it.
Hi Kannan,

Thanks for your reply. I got an understanding on the geometry construction for a simple cylinder. However, I would like to construct a concentric pipe type geometry. I would like to know if the same angle (=5 degree) should be specified for the inner cylinder too. In the manual, it says the angle should be less than 5 degrees to construct 2d axis symmetric cylinder type geometries. Could you explain why this should be the case.

Thanks

Regards
Vishal
nandiganavishal is offline   Reply With Quote

Old   April 1, 2013, 07:19
Default
  #16
Member
 
B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 54
Rep Power: 7
skyinventorbt is on a distinguished road
Dear Vishal,

For 2D simulations we must have one cell thickness for FVM. I guess this may be the reason (or) OpenFOAM defines 2D when the angle is less than 5 degrees with one cell in thickness.
skyinventorbt is offline   Reply With Quote

Old   April 1, 2013, 18:57
Default
  #17
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
Quote:
Originally Posted by skyinventorbt View Post
Dear Vishal,

For 2D simulations we must have one cell thickness for FVM. I guess this may be the reason (or) OpenFOAM defines 2D when the angle is less than 5 degrees with one cell in thickness.
Hi Kannan,

Thanks for the reply. I have constructed the wedge geometry for a simple laminar pipe flow.

I wanted to simulate a laminar flow with inlet velocty = 1m/s in a cylindrical pipe of Diameter = 0.2 m and length = 8 m. I have considered kinematic viscosity = 2e-3 Pa.s and density = 1kg/m^3.

I have tried using both SimpleFoam (switching off turbulence) and icoFoam. However, the code does not converge. I am attaching the files. Please let me know, if the geometry construction is accurate, considered (theta = 5 degree).

Thanks

Regards
Vishal
Attached Files
File Type: zip simpleFoam.zip (8.6 KB, 39 views)
nandiganavishal is offline   Reply With Quote

Old   April 2, 2013, 00:03
Default
  #18
Member
 
B T KANNAN
Join Date: Jul 2011
Location: CHENNAI (MADRAS), INDIA
Posts: 54
Rep Power: 7
skyinventorbt is on a distinguished road
Dear Vishal,
The geometry is fine but there are errors in grid. Kindly check the mesh using ""checkMesh"" and correct it. What I found is
***Number of non-orthogonality errors: 19800.
***Error in face pyramids: 59900 faces are incorrectly oriented.
Failed 2 mesh checks.


Kindly correct these errors and submit the grid for solution.
__________________
-------------------------------
B T KANNAN
PhD Research Scholar,
Indian Institute of Technology - Madras.
skyinventorbt is offline   Reply With Quote

Old   April 2, 2013, 01:12
Default
  #19
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
Quote:
Originally Posted by skyinventorbt View Post
Dear Vishal,
The geometry is fine but there are errors in grid. Kindly check the mesh using ""checkMesh"" and correct it. What I found is
***Number of non-orthogonality errors: 19800.
***Error in face pyramids: 59900 faces are incorrectly oriented.
Failed 2 mesh checks.


Kindly correct these errors and submit the grid for solution.

Thanks, very much Kannan. I found the bug. I had changed the direction of z axis in my new blockMesh file and that did the trick. Now the solver works. I would now like to check the solution and would then move on to making a concentric pipe geometry. Is there anything I should make a note of while constructing the concentric pipe geometry. Let me know.

Thanks

Vishal
nandiganavishal is offline   Reply With Quote

Old   April 2, 2013, 13:22
Default
  #20
Senior Member
 
Vishal Nandigana
Join Date: Mar 2009
Location: Champaign, Illinois, U.S.A
Posts: 208
Rep Power: 10
nandiganavishal is on a distinguished road
Quote:
Originally Posted by nandiganavishal View Post
Thanks, very much Kannan. I found the bug. I had changed the direction of z axis in my new blockMesh file and that did the trick. Now the solver works. I would now like to check the solution and would then move on to making a concentric pipe geometry. Is there anything I should make a note of while constructing the concentric pipe geometry. Let me know.

Thanks

Vishal
Hi Kannan,

As we have created a wedge whose faces are not aligned to the cooridnate plane. So the maximum height in the radial direction now corresponds to (r*cos(theta/2)) and not exactly "r". Should we rescale the values obtained for the velocity and pressure in the radial direction ? Please let me know.

Thanks

Regards
Vishal
nandiganavishal is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SimpleFoam unable to reach steady state francois OpenFOAM Running, Solving & CFD 8 November 19, 2009 14:33
plz help,urgent, vof model steady state Garima Chaudhary FLUENT 2 May 30, 2007 04:38
steady state, laminar vof_model Garima Chaudhary FLUENT 0 May 24, 2007 03:11
Flow laminar and stationary of water in a pipe manuel OpenFOAM Running, Solving & CFD 6 March 24, 2007 19:23
buoyancy driven flow in steady state in CFX4.3 raymondyin CFX 11 May 7, 2001 06:15


All times are GMT -4. The time now is 22:39.