COMSOL Forums: Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/
Most recent forum messagesWed, 12 Dec 2012 18:40:22 +0000COMSOL Forums: Hagen-Poiseuille Flowhttp://www.comsol.com/shared/images/logos/comsol_logo.gif
http://www.comsol.com/community/forums/general/thread/33518/
Re: Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/#p92235
Hi<br />
<br />
and have you checked the laminar inflow condition ?, note that the pressure inlet is counting also the extra length, I believe (check your results, as the plots will not show the fictuous inlet region, while the pressure difference should give you te drop in the laminar inflow addition length (at least is was so for some versions ago, this might have changed ;)<br />
<br />
--<br />
Good luck<br />
IvarWed, 12 Dec 2012 18:40:22 +000033518.1355337622.92235Re: Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/#p92226
Close to the theoretical solution, but you should be able to get even closer. Make sure you are using the correct value of the pressure gradient dp/dz. There is a relatively significant pressure loss at the entry which leads to an actual pressure gradient (throughout most of the pipe) that is smaller than what is estimated with P_inlet/ total length.<br />
<br />
Nagi Elabbasi<br />
Veryst EngineeringWed, 12 Dec 2012 17:26:24 +000033518.1355333184.92226Re: Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/#p92207
Thank you so much. I still can not get it right but you guys sure did help me on the way.<br />
<br />
I changed the model and now I have a 2D axissymmetric model with a radius of 0.05m and a length of 13 m.<br />
<br />
you are right I use 323.15K as temperature. <br />
<br />
I have an inlet pressure of 0.126024724<br />
and outlet still 0<br />
<br />
My maximum velocity calculated by the hagen poiseuille flow is 1.107E-2 m/s<br />
whereas COMSOL have it at 9.98E-3 m/s<br />
<br />
This is offcourse closer than before but still I would love it to be a little better.<br />
<br />
I have added the COMSOL file if anyone is interested<br />
Wed, 12 Dec 2012 11:16:51 +000033518.1355311011.92207Re: Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/#p92055
Hi<br />
<br />
i was wondering if there is not a factor 2 (/8 instead of /4 ?) missing in there. <br />
I did a quick run in 2D axi SPF aand this gives rather coherent values but is still mesh dependent, and the pressure inlet needs a bout a meter to stabilise fully.<br />
<br />
I'm not that astonished that the flow flattens in the middle with high resolution mesh. Also the temperature, you are not using 20°C are you ?, the eta/mu value doest really fit for me.<br />
<br />
--<br />
Good luck<br />
IvarMon, 10 Dec 2012 14:26:37 +000033518.1355149597.92055Re: Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/#p92052
Hi Line,<br />
<br />
Are you sure you are using a 2D axisymmetric model not a 2D planar model? If so you should only model a rectangle with width equal to the radius (not diameter) of the pipe, and COMSOL will automatically add the symmetry boundary condition on the left side of the rectangle.<br />
<br />
Nagi Elabbasi<br />
Veryst EngineeringMon, 10 Dec 2012 13:54:17 +000033518.1355147657.92052Hagen-Poiseuille Flow
http://www.comsol.com/community/forums/general/thread/33518/#p92040
I am having problems getting the same result when I try to model the Hagen poiseuille flow as when I calculate it.<br />
<br />
For more information on Hagen-Poiseuille <a href="http://en.wikipedia.org/wiki/Hagen%E2%80%93Poiseuille_equation" title="en.wikipedia.org/wiki/Hagen%E2%80%93Poiseuille_equation">en.wikipedia.org/wiki/Hagen%E2%80%93Poiseuille_equation</a><br />
<br />
I have a 2d model with a rectangel width 0.1m and length 13m.<br />
<br />
At the inlet I have a pressure of 2.52Pa using the pressure, no viscous stress option. And at the outlet 0Pa.<br />
<br />
When I calculate the velocity profile for the flow I use the viscosity provided by COMSOL<br />
<br />
eta=0.000547<br />
<br />
and I calculate the velocity from<br />
u(r)=(R^2-r^2)/(4*eta)*dp/dz<br />
<br />
in cylindrical coordinates. R is the radius ie the half width of the rectangel. r i the variable radius. <br />
dp/dz=1.9388Pa/m<br />
<br />
I have attached a picture with the development of the flow and the calculated velocity profile. What am I doing wrong here?<br />
Mon, 10 Dec 2012 09:05:34 +000033518.1355130334.92040