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

BCs for pressure driven flow using rhoCentralFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 31, 2011, 16:41
Default BCs for pressure driven flow using rhoCentralFoam
  #1
Member
 
Eric M. Tridas
Join Date: May 2011
Location: Tampa, Florida
Posts: 48
Rep Power: 6
eric.m.tridas is on a distinguished road
Hey Foamers,

I'm attempting to model the airflow through a capillary into a vacuum system. The pump flow rate is 33 cubic meters per hour (which corresponds to an average velocity of 18.67 m/s through a KF-25 flange). The capillary is 65 mm long and 0.5 mm in diameter. I am trying to calculate the pressure in the vacuum chamber based on this pumping rate, therefore I do not want to impose any pressure boundary conditions on the vacuum side of the capillary. The inlet conditions are atmospheric.

The problem is that the simulation doesn't do anything when started. I believe the problem lies in my inlet boundary conditions. My BCs are as follows:

p:
inlet - fixedValue, value 101325
outlet - zeroGradient

U:
inlet - inletOutlet, inletValue uniform (0 0 0), value uniform (0 0 0)
outlet - fixedValue, value uniform (0 -18.674 0)

I've attached an image of my geometry for reference. The pump (outlet) boundary is on the middle of the face starting at (0.075, -0.05). If anyone is working on something similar or has encountered problems of this type your help would be greatly appreciated. Also, I can post my case if anyone is interested. Thanks in advance.

-Eric
Attached Images
File Type: png geometry.png (4.2 KB, 61 views)
eric.m.tridas is offline   Reply With Quote

Old   September 1, 2011, 02:40
Default
  #2
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 214
Rep Power: 10
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Dear Eric,

I think you should change your inlet velocity boundary condition. To a pressureInletOutlet boundary condition. This gives a similar behavior as the inletOutlet, but the inletValue is calculated from the pressure field. You may want to give a small initial velocity into the domain using the value as below. Or you may want to leave it at uniform (0 0 0)

U:
inlet -
{
type pressureInletOutlet;
value uniform (0.1 0 0);
}

Hope this helps,

Tom
tomf is offline   Reply With Quote

Old   September 1, 2011, 13:19
Default
  #3
Member
 
Eric M. Tridas
Join Date: May 2011
Location: Tampa, Florida
Posts: 48
Rep Power: 6
eric.m.tridas is on a distinguished road
Thanks for your reply Tom. Unfortunately the performance of the simulation remains the same while using your suggested boundary conditions. Do you have any other potential solutions?

-Eric
eric.m.tridas is offline   Reply With Quote

Old   September 2, 2011, 04:01
Default
  #4
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Delft, Netherlands
Posts: 214
Rep Power: 10
tomf is on a distinguished road
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hello Eric,

Unfortunately I do not have any other suggestion regarding the boundary conditions. Did you check if the simulation would run if you apply a pressure at the outlet and some inflow at the inlet? If it does not run then either, something else may be wrong.

Tom
tomf is offline   Reply With Quote

Old   September 2, 2011, 09:54
Default
  #5
New Member
 
Muhammad Umer Ijaz Chaudrey
Join Date: Aug 2011
Location: Eindhoven, The Netherlands
Posts: 26
Rep Power: 5
umer.chaudrey is on a distinguished road
Send a message via Skype™ to umer.chaudrey
Hello Eric,

Try giving a pressure boundary condition at the outlet but lower value than the inlets atmospheric. Or even for a test case, try giving "0" pressure at the outlet and try running the case. It should run that way, I have been working on a slightly similar case and I used these boundary conditions.

Regards,

Umer
umer.chaudrey is offline   Reply With Quote

Old   September 6, 2011, 13:38
Default
  #6
Member
 
Eric M. Tridas
Join Date: May 2011
Location: Tampa, Florida
Posts: 48
Rep Power: 6
eric.m.tridas is on a distinguished road
Hey guys,

Thanks for your replies. I've actually figured out part of the problem. I was expecting some transonic behavior from the system as the air flowed from atmosphere through the narrow section shown in the geometry. Because of this I had set my writeInterval as 1e-8. Even though I had put my writeControl as adjustableRunTime wtih a maxCo of 0.5, the flow wasn't moving as fast as I had anticipated. This forced the deltaT to be the value of my writeInterval. My simulation wasn't actually stalling it was just being forced into a very small time step which caused the appearance no motion.

Unfortunately I get some strange/unrealistic velocities and pressures as the simulation progresses in the narrow region of the model. As the simulation is started, the imposed velocity on the outlet begins to "pull" on the surrounding fluid, and this influence propagates from the outlet through the chamber. These strange velocities and pressures observed occur even before the fluid near the exit of the capillary is affected. I've attached a few pictures to show how the simulation looks before and after the strange behavior occurs as well as a close up image of the narrow region.
Attached Images
File Type: png before.png (19.1 KB, 37 views)
File Type: png after.png (26.1 KB, 37 views)
File Type: png closeup.png (26.4 KB, 44 views)
eric.m.tridas is offline   Reply With Quote

Old   October 28, 2011, 10:08
Default
  #7
Senior Member
 
Join Date: Mar 2011
Posts: 174
Rep Power: 6
anon_a is on a distinguished road
Hello Eric

A very important note: You must calculate your Knudsen number. If Kn > 0.01 (or maybe 0.1) your solution will be wrong, the Newton-Fourier constitutive equations are not valid and you need to use kinetic theory. You might want to check F. Sharipov & V. Seleznev, Data on internal rarefied gas flows. J. Phys. Chem. Ref. Data. 27(3), 657-706 (1998). Not even the dsmcFoam solver will work due to the very low velocities (your channel is extemely long compared to the radius so you will not reach transonic velocities).

It would also help to plot a legend along with your results to have an idea on the magnitude of the velocities.

Regards,
S.P.
anon_a is offline   Reply With Quote

Old   January 10, 2015, 10:31
Default
  #8
New Member
 
M. Sabouri
Join Date: Nov 2011
Posts: 21
Rep Power: 5
Moslem is on a distinguished road
Dear Eric,
I see that this is an old post. I wish you can help me on almost the same problem. I'm working on a similar problem. I want to simulate a nozzle flow expanding into the vacuum. I've implemented different pressure boundary conditions at the outlet. unfortunately non reproduced the correct flow condition. Would you please suggest the correct outlet boundary condition?
Regards.
Moslem 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
help with compressible flow BC's (need subsonic flow) meangreen Main CFD Forum 5 July 24, 2010 13:16
2D Lid Driven Cavity Flow simulation using MATLAB josephlm Main CFD Forum 3 June 24, 2010 01:58
steam flow in a pipe driven by a pressure gradient between inlet and outlet SalvoCalvo Main CFD Forum 0 March 11, 2010 07:52
Convergence Problem with body force driven flow Jinfeng FLUENT 1 December 9, 2009 04:54
Turbulence model in flow driven by surface tension Z. ZENG Main CFD Forum 7 April 28, 1999 07:18


All times are GMT -4. The time now is 16:07.