Menu

Discussion Forum

New Discussion
Filter

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Hagen-Poiseuille Flow

Posted Dec 10, 2012, 5:05 a.m. EST Computational Fluid Dynamics (CFD) 5 Replies

Line Bergmann
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

I am having problems getting the same result when I try to model the Hagen poiseuille flow as when I calculate it.

For more information on Hagen-Poiseuille en.wikipedia.org/wiki/Hagen%E2%80%93Poiseuille_equation

I have a 2d model with a rectangel width 0.1m and length 13m.

At the inlet I have a pressure of 2.52Pa using the pressure, no viscous stress option. And at the outlet 0Pa.

When I calculate the velocity profile for the flow I use the viscosity provided by COMSOL

eta=0.000547

and I calculate the velocity from
u(r)=(R^2-r^2)/(4*eta)*dp/dz

in cylindrical coordinates. R is the radius ie the half width of the rectangel. r i the variable radius.
dp/dz=1.9388Pa/m

I have attached a picture with the development of the flow and the calculated velocity profile. What am I doing wrong here?
5 Replies Last Post Dec 12, 2012, 2:40 p.m. EST
Nagi Elabbasi Facebook Reality Labs
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 10, 2012, 9:54 a.m. EST
Hi Line,

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.

Nagi Elabbasi
Veryst Engineering
Hi Line, 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. Nagi Elabbasi Veryst Engineering
Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 10, 2012, 10:26 a.m. EST
Hi

i was wondering if there is not a factor 2 (/8 instead of /4 ?) missing in there.
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.

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.

--
Good luck
Ivar
Hi i was wondering if there is not a factor 2 (/8 instead of /4 ?) missing in there. 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. 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. -- Good luck Ivar
Line Bergmann
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 12, 2012, 7:16 a.m. EST
Thank you so much. I still can not get it right but you guys sure did help me on the way.

I changed the model and now I have a 2D axissymmetric model with a radius of 0.05m and a length of 13 m.

you are right I use 323.15K as temperature.

I have an inlet pressure of 0.126024724
and outlet still 0

My maximum velocity calculated by the hagen poiseuille flow is 1.107E-2 m/s
whereas COMSOL have it at 9.98E-3 m/s

This is offcourse closer than before but still I would love it to be a little better.

I have added the COMSOL file if anyone is interested
Thank you so much. I still can not get it right but you guys sure did help me on the way. I changed the model and now I have a 2D axissymmetric model with a radius of 0.05m and a length of 13 m. you are right I use 323.15K as temperature. I have an inlet pressure of 0.126024724 and outlet still 0 My maximum velocity calculated by the hagen poiseuille flow is 1.107E-2 m/s whereas COMSOL have it at 9.98E-3 m/s This is offcourse closer than before but still I would love it to be a little better. I have added the COMSOL file if anyone is interested
Nagi Elabbasi Facebook Reality Labs
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 12, 2012, 1:26 p.m. EST
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.

Nagi Elabbasi
Veryst Engineering
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. Nagi Elabbasi Veryst Engineering
Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago Dec 12, 2012, 2:40 p.m. EST
Hi

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 ;)

--
Good luck
Ivar
Hi 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 ;) -- Good luck Ivar

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.

Suggested Content