# buoyuantSimpleFoam & Boundaries

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

December 8, 2011, 04:38
buoyuantSimpleFoam & Boundaries
#1
Senior Member

Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 8
Hi everyone,

I'm stuck on a problem about pressure type boundaries with buoyantSimpleFoam.

My case is a simple tank with an inlet and an outlet. I try to force the convection movement by a difference of pressure between the inlet and the outlet :
- inlet : pressure : fixedValue = 150 kPa; U : pressureInletOutletVelocity
- outlet : pressure : fixedValue = 100 KPa; same for U.

The error message I got is :
Quote:
 Continuity error cannot be removed by adjusting the outflow. Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow. Total Flux : 0.00198728 Specified mass inflow : 3.98..e-06 Specified mass outflow : 1.75...e-05 Adjustable mass outflow : 0
I may have something wrong with p and p_rgh. Why these 2 are needed in the BC ?

If someone has any idea, I'd take it.

Aurélien

December 8, 2011, 05:24
#2
Member

Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 7
Quote:
 Originally Posted by Aurelien Thinat I try to force the convection movement by a difference of pressure between the inlet and the outlet : - inlet : pressure : fixedValue = 150 kPa; U : pressureInletOutletVelocity - outlet : pressure : fixedValue = 100 KPa; same for U.
I would try zeroGradient for U at the outlet.

pressureInletOutletVelocity calculates the velocity by the pressure you assigned.
Based on the mass inflow let's say with 150 kPa you get 3.98E-06 kg/s.
If you now choose a different pressure at the outlet and choose pressureInletOutletVelocity once again you'll get a mass flow according to your pressure. Since the inlet and outlet pressure differ, your mass flows are different as well and that is probably your problem because your "mass balance" is not equal to zero.

December 8, 2011, 05:37
#3
Senior Member

Aurelien Thinat
Join Date: Jul 2010
Posts: 165
Rep Power: 8
Hi Rob,

Yes, I guess you are right and I thought the same thing. So I changed my BCs into walls everywhere. Which is working.
And then, I changed the inlet to :
Quote:
 pressure p : { type fixedValue; value uniform 1e5; } pressure p_rgh { type buoyantPressure; value uniform 1e5; } U { type pressureInletVelocity; value uniform (0 0 0); }
and the outlet to :
Quote:
 pressure p : { type fixedValue; value uniform 1e5; } pressure p_rgh { type buoyantPressure; value uniform 1e5; } U { type zeroGradient; }
It's the same pressure value in inlet and outlet.
I got the same kind of error message :
Quote:
 Total flux : 0.00189... Specified mass inflow : 3.12...e-10 Specified mass outflow : 0 Adjustable mass outflow : 0

EDIT :
When I make a difference between the intlet pressure and outlet pressure (105000Pa and 95000 Pa). The solver iterates few times and then stop with the same error mesage :

Quote:
 Total flux : 1.23259 Specified mass inflow : 0.000621281 Specified mass outflow : 0.00104133 Adjustable mass outflow : 0

 December 8, 2011, 05:44 #4 Member   Rob Join Date: Sep 2011 Posts: 55 Rep Power: 7 I only gave buoyantSimpleFoam a quick try back then so I am not an expert there. But usually you only fix the pressure at the outlet and use zeroGradient for the inlet. I at least do it this way anytime. I do not know if you can use pressureInletVelocity then as a BC. What is the purpose or the goal of your simulation or what do you want to simulate?

 December 8, 2011, 05:53 #5 Senior Member   Aurelien Thinat Join Date: Jul 2010 Posts: 165 Rep Power: 8 It should be a pretty easy case (it is with other cfd codes at least) : I want to simulate the forced convection movement created by a difference of pressure between the inlet and the outlet. EDIT : I did what you suggested : a fixedValue of 90kPa at the outlet. zeroGradient pressure at the inlet. For U : pressureInletOutletVelocity at the inlet and zeeroGradient at the outlet. It computes. I'll keep you updated of the results. Thank you.

 December 9, 2011, 04:06 #6 Senior Member   Aurelien Thinat Join Date: Jul 2010 Posts: 165 Rep Power: 8 The case is now close to what it should be with fixedValue pressure and free velocity in inlet and outlet. I still have a problem with the wall boundaries in pressure. In the tutorials, they are using buoyantPressure for walls. When I use it, the case is running during something like 100 iterations. And when I use zeroGradient, it explodes after 3 iterations. If someone has any idea to solve this problem. Thank you. Aurelien

December 9, 2011, 04:30
#7
Member

Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 7
Quote:
 Originally Posted by Aurelien Thinat The case is now close to what it should be with fixedValue pressure and free velocity in inlet and outlet.
That's cool.

Quote:
 Originally Posted by Aurelien Thinat I still have a problem with the wall boundaries in pressure. In the tutorials, they are using buoyantPressure for walls. When I use it, the case is running during something like 100 iterations. And when I use zeroGradient, it explodes after 3 iterations.
Does this mean your case explodes with buoyantPressure for walls as well? Or does your case converge within the 100 iterations?

 December 9, 2011, 04:31 #8 Senior Member   Aurelien Thinat Join Date: Jul 2010 Posts: 165 Rep Power: 8 It explodes after 100 or 200 iterations (depends on the relaxation factor I used).

 December 9, 2011, 04:33 #9 Member   Rob Join Date: Sep 2011 Posts: 55 Rep Power: 7 Maybe you should post your schemes/solution files and of course your BC's. Maybe the reason for the simulation blowing up is something else.

 December 9, 2011, 11:27 #10 Senior Member   Aurelien Thinat Join Date: Jul 2010 Posts: 165 Rep Power: 8 It seems that buoyantSimpleFoam wasn't able to deal with this type of study. I switched to rhoSimplecFoam and it's starting. But now I have some problems to keep the computation stable. EDIT : It's now working well. If someone could list the limits of the different "buoyant*" solvers, or confirm they are not suitable for high pressure gradients. Thank you. Last edited by Aurelien Thinat; December 10, 2011 at 17:43.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Eran FloEFD, FloWorks & FloTHERM 3 August 11, 2009 04:23 PK FLUENT 0 July 12, 2007 11:58 swetha FLUENT 1 November 26, 2006 23:02 Jared CD-adapco 4 August 5, 2005 19:36

All times are GMT -4. The time now is 03:27.

 Contact Us - CFD Online - Top