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

Degassing boundary condition

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 14, 2009, 07:02
Default Hello all, Anyone please te
  #1
Member
 
Raghavendra
Join Date: Mar 2009
Location: Goteborg, Sweden
Posts: 94
Rep Power: 8
raagh77 is on a distinguished road
Send a message via Yahoo to raagh77 Send a message via Skype™ to raagh77
Hello all,

Anyone please tell me how "degassing boundary condition" (which is used in CFX) is implemented in OpenFOAM ?

Regards
Raghavendra
raagh77 is offline   Reply With Quote

Old   March 14, 2009, 07:10
Default Can you provide some informati
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Can you provide some information on this boundary condition:
- what is it used for
- how is it implemented

There may be some papers or information in publicly available sources as the name itself does not tell me what to do.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   March 14, 2009, 07:18
Default Hi..Hrvoje Jasak my task is t
  #3
Member
 
Raghavendra
Join Date: Mar 2009
Location: Goteborg, Sweden
Posts: 94
Rep Power: 8
raagh77 is on a distinguished road
Send a message via Yahoo to raagh77 Send a message via Skype™ to raagh77
Hi..Hrvoje Jasak
my task is to simulate air particle being injected (at the center) to the water as continuous medium. (Becker case)

(Domain being simple 2D rectangular region sides and bottom (except at the center) being walls. Velocity inlet at the center and pressure outlet at the top)


Boundary conditions

1) alpha 0 gradients at the wall
2) inlet uniform 1
3) 3/4 th water (continous phase) filled so alpha being eual to 0 and remaining 1/4 air(discrete phase) and alpha being 1
4) alpha being 1 at the "top" where pressure outlet boundary conditions are set

But I was suggested to implement degassing boundary condition at the top so that only gas (air) can pass through it but not the liquid (continuous phase - water )

I am just wondering how this condition can be implemented in OpenFOAM..

Regards
Raghavendra
raagh77 is offline   Reply With Quote

Old   March 14, 2009, 08:08
Default Hi, Degassing Boundary Cond
  #4
Senior Member
 
Holger Marschall
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 122
Rep Power: 10
holger_marschall is on a distinguished road
Send a message via Skype™ to holger_marschall
Hi,

Quote:
Degassing Boundary Condition is only available when one fluid is "Continuous" and the other "Dispersed Fluid". If the second fluid is
"Dispersed Solid", it is not available.

Degassing Boundary Conditions are used to model a free surface from which dispersed bubbles are permitted to escape, but the liquid phase
is not. They are useful for modelling flow in bubble columns. When Degassing Boundary Condition is selected as the Flow Specification of an outlet, the continuous phase sees this boundary as a free-slip wall and does not leave the domain. The dispersed phase sees this boundary as an outlet. However, the outlet pressure is not specified. Instead, a pressure distribution is computed on this fixedposition boundary, and can be interpreted as representing the weight of the surface height variations in the real flow.

Since the Degassing Boundary Condition does not specify a pressure value, a fixed pressure reference point will be set automatically for a domain with an incompressible flow with degassing boundaries, but no pressure boundaries.
from Physical Models Fluid Properties and Boundary Conditions in CFX-5

But in your case: why don't you simulate the free-surface in the bubble column as the free-surface it is The degassing b.c. in CFX is just a model for a quiscient free-surface being spatial fixed at the boundary. Just fill the bubble column to the level it was filled in the Becker experiment (compare to Henrik Rusche's PhD thesis).

Anyway I think you can "build" your boundary condition from the existing b.c.'s in OpenFOAM.

best regards
Holger
__________________
Holger Marschall
web: http://www.holger-marschall.info
mail: holgermarschall@yahoo.de
holger_marschall is offline   Reply With Quote

Old   March 14, 2009, 08:27
Default Hello Holger Marschall Than
  #5
Member
 
Raghavendra
Join Date: Mar 2009
Location: Goteborg, Sweden
Posts: 94
Rep Power: 8
raagh77 is on a distinguished road
Send a message via Yahoo to raagh77 Send a message via Skype™ to raagh77
Hello Holger Marschall

Thanks for your reply..
Yes, as you said am using the existing boundary conditions which is in twoPhaseEulerFoam tutorials.

But with slight modifications.
1.alpha 3/4 domain height 0 and top 1/4 is 1.
2. and K and epsilon being 1e-08 and 0.1 in the internal fields.
3. Small timestep 1e-04 (also with even less value!! )

but after few iterations
simulation crashes
1.sometimes with high courant number
2.sometimes with high continuity error and high cournat number.

Do I need to do some modifications for the solver and recompile it again
(also I tried with bubbleFoam Solver but without success )

awaiting for your reply..
raagh77 is offline   Reply With Quote

Old   March 14, 2009, 08:36
Default Hi Raghavendra, what do you
  #6
Senior Member
 
Holger Marschall
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 122
Rep Power: 10
holger_marschall is on a distinguished road
Send a message via Skype™ to holger_marschall
Hi Raghavendra,

what do you mean by "sometimes"? What else did you change? Can you post the error message?

However, it would be much more interesting to code that b.c. ... Do you have any informations on how the outlet pressure is computed? Are there some papers/reports?

best regards
Holger
__________________
Holger Marschall
web: http://www.holger-marschall.info
mail: holgermarschall@yahoo.de
holger_marschall is offline   Reply With Quote

Old   March 14, 2009, 09:00
Default "sometimes" I meant was changi
  #7
Member
 
Raghavendra
Join Date: Mar 2009
Location: Goteborg, Sweden
Posts: 94
Rep Power: 8
raagh77 is on a distinguished road
Send a message via Yahoo to raagh77 Send a message via Skype™ to raagh77
"sometimes" I meant was changing timestep to smaller values..but didn't work

This is the error message i get at the last iterations

#0 Foam::error::printStack(Foam:stream&) in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Uninterpreted: [0xb7fda400]
#3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<foam::field,> const&, Foam::UPtrList<foam::lduinterfacefield> const&, Foam::Field<double> const&, unsigned char) const in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<foam::lduma trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<foam::field<double> >&, Foam::PtrList<foam::field<double> >&, unsigned char) const in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#7 Foam::fvMatrix<double>::solve(Foam::Istream&) in "/home/ragh/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libfiniteVolume.so"
#8 main in "/home/ragh/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/twoPhaseEulerFo am"
#9 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#10 Foam::regIOobject::readIfModified() in "/home/ragh/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/twoPhaseEulerFo am"
Floating point exception


at this timestep
Courant Number mean: 0.00349919 max: 67.8872


can I mail you the case setup files?
It is just 8mb..

Regards
Raghavendra
raagh77 is offline   Reply With Quote

Old   March 14, 2009, 09:14
Default it is just 1.6mb ..
  #8
Member
 
Raghavendra
Join Date: Mar 2009
Location: Goteborg, Sweden
Posts: 94
Rep Power: 8
raagh77 is on a distinguished road
Send a message via Yahoo to raagh77 Send a message via Skype™ to raagh77
it is just 1.6mb ..
raagh77 is offline   Reply With Quote

Old   March 14, 2009, 09:25
Default Sure. Please finde my e-mail i
  #9
Senior Member
 
Holger Marschall
Join Date: Mar 2009
Location: Darmstadt, Germany
Posts: 122
Rep Power: 10
holger_marschall is on a distinguished road
Send a message via Skype™ to holger_marschall
Sure. Please finde my e-mail in my profile!

regards
Holger
__________________
Holger Marschall
web: http://www.holger-marschall.info
mail: holgermarschall@yahoo.de
holger_marschall is offline   Reply With Quote

Old   July 27, 2013, 04:39
Default Degass BC
  #10
Member
 
Jeong Kim
Join Date: Feb 2010
Posts: 41
Rep Power: 7
enoch is on a distinguished road
To. Raghavendra

Have you ever implemented this BC?

Thanks.
enoch is offline   Reply With Quote

Old   March 27, 2015, 12:19
Default
  #11
New Member
 
Join Date: May 2014
Location: Germany
Posts: 23
Rep Power: 3
hester is on a distinguished road
Hello Raghavendra,

I am also curious about how to implement the degassing boundary condition. Have you ever figured out how to calculate the outlet pressure?

Regards,
hester
hester 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
Degassing boundary condition Ayman Dohdoh FLUENT 0 October 28, 2007 05:41
Degassing Condition Siqui FLUENT 0 June 13, 2005 17:49
Degassing Condition.......Please Help Paresh Jain CFX 1 February 19, 2004 13:19
Degassing Condition in CFX 5.6 Paresh Jain CFX 2 February 14, 2004 01:27
Degassing Boundary Condition ???? thomas CFX 2 February 12, 2004 03:59


All times are GMT -4. The time now is 17:50.