CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Verification & Validation

laminar pipe flow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 8, 2022, 16:52
Default laminar pipe flow
  #1
New Member
 
asarmakeeva's Avatar
 
Ana
Join Date: Jan 2018
Location: Washington, DC
Posts: 5
Rep Power: 8
asarmakeeva is on a distinguished road
Hi! I was running laminar pipe flow, when I sampled at the place around the entrance length ~0.04m, I was getting fine results. When I expand my pipe, I have some changes in maximum velocity(flow is still changing), and the error with the analytical solution gets bigger when I increased the number of mesh cells.

Initial settings:
length of the pipe = 0.2 meters
R = 0.00227
rho = 1.0
v = 0.5216m/s

Analytical maximum velocity = 1.043 m/s I am getting for the most refined case for 460k cells velocity 1.071 m/s and for the coarse mesh 57k cells max velocity is 1.059m/s. I know that to have some errors is fine, but this is a very simple case.

Maybe a problem in my boundary solution or this variation in velocity field is well known and fine?

Attaching the table to show how velocity is changing Screenshot from 2022-02-08 15-36-36.pngand link with the case and my postprocessing https://github.com/asarmakeeva/PipeFlowVerification


Thank you for your time and help.
asarmakeeva is offline   Reply With Quote

Old   February 9, 2022, 03:41
Default
  #2
Senior Member
 
Josh Williams
Join Date: Feb 2021
Location: Scotland
Posts: 112
Rep Power: 5
joshwilliams is on a distinguished road
Quote:
Originally Posted by asarmakeeva View Post
Hi! I was running laminar pipe flow, when I sampled at the place around the entrance length ~0.04m, I was getting fine results. When I expand my pipe, I have some changes in maximum velocity(flow is still changing), and the error with the analytical solution gets bigger when I increased the number of mesh cells.

Initial settings:
length of the pipe = 0.2 meters
R = 0.00227
rho = 1.0
v = 0.5216m/s

Analytical maximum velocity = 1.043 m/s I am getting for the most refined case for 460k cells velocity 1.071 m/s and for the coarse mesh 57k cells max velocity is 1.059m/s. I know that to have some errors is fine, but this is a very simple case.

Maybe a problem in my boundary solution or this variation in velocity field is well known and fine?

Attaching the table to show how velocity is changing Attachment 88304and link with the case and my postprocessing https://github.com/asarmakeeva/PipeFlowVerification


Thank you for your time and help.

Hi Ana,



The issue is that you have an ill-posed set of boundary conditions. Specifically, you define a velocity at the inlet and outlet in this case. You should change the velocity outlet to

Code:
 type         zeroGradient;

A good indicator of your BC issue is that the mesh is so fine -- way too fine in fact for the simple case of laminar flow in a pipe. As it is a simple solver (no multiphase or non-isothermal effects to introduce error), our first debugging check should usually be the boundary conditions.


Best,
Josh
joshwilliams is offline   Reply With Quote

Old   February 9, 2022, 16:49
Default
  #3
New Member
 
asarmakeeva's Avatar
 
Ana
Join Date: Jan 2018
Location: Washington, DC
Posts: 5
Rep Power: 8
asarmakeeva is on a distinguished road
Hi Josh, thank you for your answer,

because zero gradient is the boundary condition which extrapolates the quantity to the patch from the nearest cell value the issue is not in the the ill-posed bc at the outlet.

for the sanity check I run simulation with changes which you suggested and got same result.

Kind regards, Ana.

Last edited by asarmakeeva; February 9, 2022 at 18:23.
asarmakeeva is offline   Reply With Quote

Old   August 2, 2023, 11:15
Default
  #4
New Member
 
Karim Hassanieh
Join Date: Jul 2016
Posts: 19
Rep Power: 9
Karimh89 is on a distinguished road
Hi Anna
You also need to change the pressure at the outlet i agree the boundary condition is ill posed in summary this is what should be done
Inlet U-fixedValue
P-zeroGradient
Outlet U-zeroGradient
p-FixedValue
Also i noticed your running the case in transient state i beleive it should run fine under steady state condition

regards
Karim
Karimh89 is offline   Reply With Quote

Reply

Tags
pipe flow, verification

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Boundary conditions of laminar flow in pipe alireza.glz OpenFOAM 4 May 27, 2019 06:03
Is there a correlation for Nu number in laminar pipe flow with constant wall heat flu Shomaz ul Haq CFX 5 June 19, 2017 09:04
Radial velocity with a swirling core in laminar pipe flow MissCFD CFX 1 June 26, 2015 08:03
gravitational acceleration in laminar pipe flow atmcfd ANSYS 0 January 4, 2010 23:19
First steps - laminar flow in a pipe Maria Phoenics 8 November 27, 2001 12:26


All times are GMT -4. The time now is 18:11.