[moblimostafa]

Dear Nima

Could you send me your your code?
mostafamobli@gmail.com

[yaoyina]

Hello, Nima Sam.

Would you be so kind to send your code to yao-yn14@mails.tsinghua.edu.cn?
Thanks in advance!

[chaolian]

Hi, Nima

Thanks for your contribution. I am interested in your solver.
Could you please send your solver to me: chaolian@mtu.edu

Thanks in advance.

[Kanarya]

Hi Nima,

I am also interested in your code
kah.recep@gmail.com

thanks in advance!!

[smraniaki]

Hello Nimasam

Thank you very much for your contribution to the world of OpenFoam. I am also interested to see your code. Would you please email it to rifel2002@yahoo.com
Smran

[gaza]

Hi Nima Sam,
I compiled your solver but the results are wrong.
I used the fluids properties the same as Bitan Shu. You can see here

http://www.cfd-online.com/Forums/ope...pvofhardt.html

I am also trying to do solver for evaporation/condensation but my results are also wrong.
Can anybody upload correct solver for evaporation/condensation purposes?

[Kanarya]

I do not think, there is correct one!
what is wrong with your results?

Originally Posted by gaza

Hi Nima Sam,
I compiled your solver but the results are wrong.
I used the fluids properties the same as Bitan Shu. You can see here

http://www.cfd-online.com/Forums/ope...pvofhardt.html

I am also trying to do solver for evaporation/condensation but my results are also wrong.
Can anybody upload correct solver for evaporation/condensation purposes?

[gaza]

Hi Kanarya,
In attachement

http://fluid.itcmp.pwr.wroc.pl/~pblasiak/download.html

there is stefan problem case with settings as in Bitan Shu conference paper:

http://pl.scribd.com/doc/28869783/Pa...0080214#scribd

Additionally there is file interfacePosStefProbl_analSol.dat where are analitycal results for this case for time 0-25 s. For 25 s interface position should be 3.8 mm but numerical value is much smaller (see figure in the attachement).

I created my own solver. Results are close to analytical but still incorrect.
I uploaded this with damBreak case here

http://www.cfd-online.com/Forums/ope...-modeling.html

In Nima Sam I also added additional term as in my solver in TEqn.H to avoid unboundedness in temperature field.

What do you think about Nima Sam solver and my?
Do you know what is wrong?
Did you validate Nima Sam solver on damBreak case?

[Kanarya]

Hi Gaza,

Which solver is yours?
Do you consider that your water or vapour is reached in saturation ?
If not, you need source and sink term for both phases...
For instance, if water is cool enough, it needs to gain energy in order to evaporate for this reason you need source term...

Best regards,

[gaza]

Hi Kanarya,
My solver is interEvapCondPhaseChangeFoam (first link). The second link is a test case for solver
evaPhaseChangeFoamU (solver created by Nima Sam). As you can see in the figure in attachement
numeirical results is not equal analytical, which should be 3.8mm for 25 s for fluid properties used (see Bitan Shu paper). Why do you think that Nima Sam solver works well?

I did not consider that water is in saturation condition because I do not know how to provide this or
temperature at interface equal to saturation temperature. I only used initial conditions for water close to saturation temperature.

Here you can see how to compile my solver

http://www.cfd-online.com/Forums/ope...-modeling.html

If you run stefan problem with my solver you can see that results are close to analytical but still not well. It is clear in damBreak case where you can see the difference in comparison to interFoam.
Do you know how can I improve my solver?

[gaza]

Originally Posted by Kanarya

Hi Gaza,

Which solver is yours?
Do you consider that your water or vapour is reached in saturation ?
If not, you need source and sink term for both phases...
For instance, if water is cool enough, it needs to gain energy in order to evaporate for this reason you need source term...

Best regards,

And one more question: what do you mean "you need source and sink term for both phases"?
My TEqn looks like

fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T)
+ fvm::div(rhoPhiCpf, T)
- fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoPhiCpf), T) // additional term to avoid uboundedness ot T
- fvm::laplacian(kappaf, T)
==
fvm::Sp(hEvap*(mCond + mVap)/T, T)
);

hEvap - enthalpy of vaporisation
mCond - mass due to condensation
mVap - mass due to evaporation

mCond and mVap are calculated using phase change model. In this case Ganapathy.

[Kanarya]

As I can see here you did not specified for which phase you are solving T equation right?so there is not volume fraction of alpha or sth?or I am wrong?

Originally Posted by gaza

And one more question: what do you mean "you need source and sink term for both phases"?
My TEqn looks like

fvScalarMatrix TEqn
(
fvm::ddt(rhoCp, T)
+ fvm::div(rhoPhiCpf, T)
- fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoPhiCpf), T) // additional term to avoid uboundedness ot T
- fvm::laplacian(kappaf, T)
==
fvm::Sp(hEvap*(mCond + mVap)/T, T)
);

hEvap - enthalpy of vaporisation
mCond - mass due to condensation
mVap - mass due to evaporation

mCond and mVap are calculated using phase change model. In this case Ganapathy.

[gaza]

Originally Posted by Kanarya

As I can see here you did not specified for which phase you are solving T equation right?so there is not volume fraction of alpha or sth?or I am wrong?

I created my solver based on interPhaseChangeFoam. In VOF energy equation is solved for both phases.

[Kanarya]

yes but if you look to the solver alpha.H it is only solving for one phase not like two fluid model, if you check the twoPhaseEulerFoam you will see...
what is you final gaol...what do you want to simulate with this solver ...probably two fluid model suits you better...

Originally Posted by gaza

I created my solver based on interPhaseChangeFoam. In VOF energy equation is solved for both phases.

[gaza]

Originally Posted by Kanarya

yes but if you look to the solver alpha.H it is only solving for one phase not like two fluid model, if you check the twoPhaseEulerFoam you will see...
what is you final gaol...what do you want to simulate with this solver ...probably two fluid model suits you better...

yes but I want to use one fluid approach because it is less computationally expensive.
I think that alpha equation in my solver is solved in wrong way but I do not know how to fix it.

Are you using Nima Sam solver? What do you think about results of it?

[Kanarya]

no I am not...I do not have it..could you please send it to me ?

Originally Posted by gaza

yes but I want to use one fluid approach because it is less computationally expensive.
I think that alpha equation in my solver is solved in wrong way but I do not know how to fix it.

Are you using Nima Sam solver? What do you think about results of it?

[gaza]

Originally Posted by Kanarya

no I am not...I do not have it..could you please send it to me ?

yes of course, here it is (third link):

http://fluid.itcmp.pwr.wroc.pl/~pblasiak/download.html

it is for OF1.6.x

[NightWing]

Can i get the final code???


email: arunkraj03@gmail.com

[thomas.]

Hi Nimasan,

could you please send me your code?

My email is thomas.schaub.h@gmail.com

Thank you very much in advance!

[LCC]

Hello, nimasam. I have some problems about two phase change progress. Can you send me your final code? My email is 18600156862@163.com. Thank you very much!