[ARUN K RAJ]

Hi all,
I would like to know much about how condensation phenomenon could be done by using openFOAM solvers like twoPhaseEulerFoam and compressibleTwoPhaseEulerFOAM. Any detailed info abt tutorials related with any of these solvers or journals regarding multiphase flow would be most appreciated.
There are certain default tutorials like bubblecolumn and bed in these solvers but no detailed info abt what it does isn't provided.

Thank you,
Have a great day all

[wyldckat]

Greetings Arun,

I've made a very quick search online just now for "OpenFOAM condensation" and the following two links seem to be a good starting point:

• http://www.cfd-online.com/Forums/ope...-openfoam.html - on this thread you can find some vague concepts, namely that OpenFOAM might not have implemented exactly what you need.
• http://www.tfd.chalmers.se/~hani/kur...ChangeFoam.pdf - this tutorial is from here: http://www.tfd.chalmers.se/~hani/kurser/OS_CFD_2011/ - and seems to be a very good example to get started on this topic.

Best regards and good luck!
Bruno

[ARUN K RAJ]

Thank you Bruno,

I have started to solve condensation in a vertical plate geometry as a start using interPhaseChangeFoam. I have added the energy equation as per the Martin's tutorial related with interPhaseChangeFoam. But the problem i face now is i get many errors such as the ones below:

 #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/openfoam230/platforms/linuxGccDPOpt/lib/libfiniteVolume.so"
#7  
 at interPhaseChangeFoam.C:0
#8  
 in "/opt/openfoam230/platforms/linuxGccDPOpt/bin/interPhaseChangeFoam"
#9  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#10  
 in "/opt/openfoam230/platforms/linuxGccDPOpt/bin/interPhaseChangeFoam"
Floating point exception (core dumped)

It is most certainly related with my boundary conditions or may be the new addtions i have made to the solver using the tutorial. But i have tried changing the B.C's bt the same error shows up after few iterations. Could anyone suggest me hw to solve this problem...

[ron_OFuser]

Arjun,

This work is already available on Chalmers website as pointed by wyldckat. You might be messing with some terms carefully follow the instruction given in the report.

http://www.tfd.chalmers.se/~hani/kur...ChangeFoam.pdf

Cheers Chalmers

[wyldckat]

Quick answer:

Originally Posted by ARUN K RAJ

It is most certainly related with my boundary conditions or may be the new addtions i have made to the solver using the tutorial. But i have tried changing the B.C's bt the same error shows up after few iterations. Could anyone suggest me hw to solve this problem...

We're not psychic
This to say that if you have not yet solved this problem, then please follow the instructions given here: http://www.cfd-online.com/Forums/ope...-get-help.html

[ARUN K RAJ]

I would make my problem clear first of all. I am trying to study condensation phenomenon in a vertical flat plate geometry. The inlet is marked as top of plate. I am using steam-air mixture for analysis. I hve added the temperature to the solver " InterPhaseChangeFoam" using the tutorial by Martin andersen. http://www.pdfiz.net/k-7721458.html
Problems faced:
1. when i am trying to compile the solver i get these errors :

/opt/openfoam230/src/finiteVolume/lnInclude/alphaControls.H: In function ‘int main(int, char**)’:
/opt/openfoam230/src/finiteVolume/lnInclude/alphaControls.H:17:6: warning: unused variable ‘alphaApplyPrevCorr’ [-Wunused-variable]
/opt/openfoam230/src/finiteVolume/lnInclude/alphaControls.H:23:8: warning: unused variable ‘icAlpha’ [-Wunused-variable]
/opt/openfoam230/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]

2. Secondly if i proceed further with this solver, If i change my deltaT value to a higher value like 0.01 or 0.1 (usually deltaT was 1e-08 as in tutorial of cavitating bullet as exmaple of interphasechangefoam solver) the iterations stops since the courant number reaches a much higher value of the order of (e100 ++).

3. I am using cavitation model of Merkle

Any help would be most appreciated.

[wyldckat]

Hi Arun,

Originally Posted by ARUN K RAJ

Problems faced:
1. when i am trying to compile the solver i get these errors :

/opt/openfoam230/src/finiteVolume/lnInclude/alphaControls.H: In function ‘int main(int, char**)’:
/opt/openfoam230/src/finiteVolume/lnInclude/alphaControls.H:17:6: warning: unused variable ‘alphaApplyPrevCorr’ [-Wunused-variable]
/opt/openfoam230/src/finiteVolume/lnInclude/alphaControls.H:23:8: warning: unused variable ‘icAlpha’ [-Wunused-variable]
/opt/openfoam230/src/finiteVolume/lnInclude/readTimeControls.H:38:8: warning: unused variable ‘maxDeltaT’ [-Wunused-variable]

Those are warnings, not errors.
I'll quote myself from a post I made back in 2010 :

Originally Posted by wyldckat

Those warnings are actually normal. AFAIK, those solvers use header files that are common to various solvers and not all solvers use those variables, therefore it triggers warning messages on those that don't use them.

Originally Posted by ARUN K RAJ

2. Secondly if i proceed further with this solver, If i change my deltaT value to a higher value like 0.01 or 0.1 (usually deltaT was 1e-08 as in tutorial of cavitating bullet as exmaple of interphasechangefoam solver) the iterations stops since the courant number reaches a much higher value of the order of (e100 ++).

Uh... That's normal. Did you study the OpenFOAM User Guide? Because this is explained in the very first tutorial: http://www.openfoam.org/docs/user/cavity.php

Originally Posted by ARUN K RAJ

3. I am using cavitation model of Merkle

Who? Please be clearer regarding specific details, because even if I might be able to get you in the right direction, you should not assume that I know everything, because I do not know everything about everything .


Best regards,
Bruno

[ARUN K RAJ]

Thank you ALL for the valuable help and suggestions

I am trying to switch to twoPhaseEulerFoam as it seems to be the choice of most of the members here to perform condensation phenomenon in heat pipe. I have started studying abt the solver using the tutorial "bed". If anyone knws some papers relating to condensation & heat pipes using the solver pls let me know

[bibin.sme]

Originally Posted by ARUN K RAJ

Hi all,
I would like to know much about how condensation phenomenon could be done by using openFOAM solvers like twoPhaseEulerFoam and compressibleTwoPhaseEulerFOAM. Any detailed info abt tutorials related with any of these solvers or journals regarding multiphase flow would be most appreciated.
There are certain default tutorials like bubblecolumn and bed in these solvers but no detailed info abt what it does isn't provided.

Thank you,
Have a great day all

come to me I ll tell..

[Kanarya]

hi Arun,

why did you decided to change it to Eulerian model?

what is your application?
is it heat exchanger or?
thanks!

Originally Posted by ARUN K RAJ

Thank you ALL for the valuable help and suggestions

I am trying to switch to twoPhaseEulerFoam as it seems to be the choice of most of the members here to perform condensation phenomenon in heat pipe. I have started studying abt the solver using the tutorial "bed". If anyone knws some papers relating to condensation & heat pipes using the solver pls let me know

[gaza]

Hi Arun
I am also trying to make solver to model condensation. I followed Martin's tutorial and used Merkle model. However it only works for coarse mesh with Co = 0.0001 or lower and timestep of order 1e-8. Additionally results were unphysical.
I do not know why timestep is so small in interPhaseChangeFoam but in interFoam it is or order of 1e-3.

[Kanarya]

Hi Gaza,

It is better to share the code in order to see what is wrong with it..
general comments lie this it is very difficult to understand and fix it...
which version of OF are you using?
what can of application are you going to simulate?
thanks!

kanarya

Originally Posted by gaza

Hi Arun

I am also trying to make solver to model condensation. I followed Martin's tutorial and used Merkle model. However it only works for coarse mesh with Co = 0.0001 or lower and timestep of order 1e-8. Additionally results were unphysical.
I do not know why timestep is so small in interPhaseChangeFoam but in interFoam it is or order of 1e-3.