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

Problems with simpleFoam with flow in pipe

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 19, 2010, 14:08
Default Problems with simpleFoam with flow in pipe
  #1
New Member
 
Fernando
Join Date: Feb 2010
Posts: 28
Rep Power: 7
fsalvucci is on a distinguished road
Hi everyone!! I need a little help. I'm performing a 3D simulation of laminar, incompressible flow in a straight pipe. When analizing my results, i found them a little bit confusing. The thing is, that, instead of finding that the flow is parabolic, it develops into a parabolic-like profile but "translated" from the origin, that is, a parabolic manner but that in the walls is not zero. The detalis of the results can be seen in the attached images.

As boundary conditions for U, i set a parabolic profile with max_vel = 0.46 m/s at inlet, zeroGradient at outlet, and zero for the wall.
For p, i set it to zero at the outlet. (see the o files for U and p attached)

I also attach the fvSolution and fvSchemes files.

In addition, i check the residuals, and they seem to converge nicely (see Ux residual plot attached)

The tubes dimensions are diamater = 0.003 m and lenght = 0.01 m, the cinematic viscosity (nu) is 3e-6, and turbulence is off.

Anyone may have an idea of what is happening?
Attached Files
File Type: gz 3D-pipe.tar.gz (68.7 KB, 23 views)
fsalvucci is offline   Reply With Quote

Old   July 19, 2010, 14:11
Default
  #2
New Member
 
Fernando
Join Date: Feb 2010
Posts: 28
Rep Power: 7
fsalvucci is on a distinguished road
Extra-comment: I also performed the simulations with a surfaceNomalFixedValue for inlet velocity, and happens the same, the flows develope into this "shifted" parabolic profile.

Helpppp!!! jajajaja
fsalvucci is offline   Reply With Quote

Old   July 19, 2010, 15:43
Default
  #3
Senior Member
 
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 11
MartinB is on a distinguished road
Hi Fernando,
here is a case with surfaceNormalFixedValue as boundary condition and your specifications.
My velocity profile looks ok...
Run "blockMesh" first...
Martin
Attached Files
File Type: gz 3-D_pipe_eigen.tar.gz (2.4 KB, 36 views)
MartinB is offline   Reply With Quote

Old   July 19, 2010, 16:24
Default
  #4
New Member
 
Fernando
Join Date: Feb 2010
Posts: 28
Rep Power: 7
fsalvucci is on a distinguished road
Thanks a lot!!! Apparently, although i was setting turbulence off, in RASmodel i was setting kepsilonmodel, instead of laminar.... . I think that would correct my results!

Graciassssssssss
fsalvucci is offline   Reply With Quote

Old   April 17, 2012, 04:20
Default
  #5
New Member
 
Mak
Join Date: Jul 2010
Location: United States
Posts: 10
Rep Power: 7
trinath2rao is on a distinguished road
Quote:
Originally Posted by MartinB View Post
Hi Fernando,
here is a case with surfaceNormalFixedValue as boundary condition and your specifications.
My velocity profile looks ok...
Run "blockMesh" first...
Martin
Hi Martin ,

I am entry level Open Foam User,

When I saw in the system/fvSchemes and system/fvSolution files..I found k and epsilon equation are set to solve..
==================================================
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) Gauss upwind;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}
===========================================


Since the problem is laminar, why are you solving for k and e equation..?

can you help me ???
trinath2rao is offline   Reply With Quote

Old   April 17, 2012, 04:25
Default
  #6
Senior Member
 
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 11
MartinB is on a distinguished road
Hi,

these settings are ignored if simulation is set to be laminar in constant/turbulenceProperties or constant/RASProperties.

Martin
MartinB is offline   Reply With Quote

Old   April 17, 2012, 08:31
Default
  #7
New Member
 
Mak
Join Date: Jul 2010
Location: United States
Posts: 10
Rep Power: 7
trinath2rao is on a distinguished road
Dear Martin,

Thank You for quick reply...
trinath2rao 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
Disturbed flow field at outlet boundary (Multiphase flow through pipe) Michiel CFX 17 April 21, 2010 10:14
Pipe flow with obstacle - HELP Min FLUENT 5 April 10, 2007 17:23
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11
flow in a vibrating pipe xu CFX 0 October 30, 2001 09:35
Hydrostatic pressure in 2-phase flow modeling (CFX4.2) HB &DS CFX 0 January 9, 2000 14:19


All times are GMT -4. The time now is 17:44.