[Elham]

Originally Posted by nimasam

Dear Elham
Please re-read my points, you will figure out how it is derived.
equation 18:

lets consider S=-mDot(1/rho1-alpha1(1/rho1-1/rho2)))
then source terms for boiling and condensations would be:
Sb=S*alpha1
Sc=S*(1-alpha1)
Also i mentioned in previous post that:

then right hand side of equation 18 becomes:

so alpha1*div(U) is appearing here
now consider boiling and condensation separately,
for boiling the source term is Sb=S*alpha1
so

for condensation Sb=S*(1.0-alpha1), then:

Sorry but I am still confused.

Base on

ddt(alpha1)+div(U,alpha1)+Sp*alpha1+Su=0

For boiling it should be :

Sp=- div(U) and Su= -S

[nimasam]

i said Sb=S*alpha1, if you consider Su=S, then you missed alpha1 in equation.
about negative signs, it is better to refer the code to see the exact formulation of MULES, i may missed the signs in post

[Elham]

Originally Posted by nimasam

i said Sb=S*alpha1, if you consider Su=S, then you missed alpha1 in equation.
about negative signs, it is better to refer the code to see the exact formulation of MULES, i may missed the signs in post

Dear Nima,

Thanks for all kind answers. I got the idea.You have separated the general source term into two parts; condensation and boiling so that S=Sb+Sc. Then for the boiling phase you just have Sb as source term and for condensation just Sc.

Cheers,

Elham

[claudiocor]

Originally Posted by nimasam

i have developed a solver based on interFOAM, to solve energy equation and besides consider mass transfer between two phases, this solver works correctly for one dimensional case (stephan phase change problem)
but for two case studies, temperature at interface behaves strangely
now any suggestion, cooperation or idea will be helpful

P.S
developed files and case studies are available in attachment
+
some descriptions can be found here:

http://www.4shared.com/document/-eBG...OF_method.html

Hi nimasam
I'm interested in this problem too. I have trying to solve the problem of water vapor condensation in vertical pipe. My problem is three-dimensional. Could you help me solve this problem? I'm used OpenFoam 4.1

[nimasam]

Dear Cláudio Corrêaread this post, i upload a solver based on OpenFOAM_2.2.0 which considers both evaporation and condensation and it is suitable for your work

[Elham]

Dear Nima,

I am wondering why you defined two different word classes, word& type and word& alpha1Name in phaseChangeHeatFoam?

In constructor of phaseChangeTwoPhaseMixture.C you have:

Foam:haseChangeTwoPhaseMixture:haseChangeTwoPh aseMixture
(
const word& type,
const volVectorField& U,
const surfaceScalarField& phi,
const word& alpha1Name
)

I would appreciate if you let me know.

Cheers,

Elham

[Elham]

Dear Niam,

Regarding to phaseChangeHeatFoam the energy equation (25) in your paper is :

ddT + div (UT) - lapplacian (k/rhoC, T) = -1/rhoC *mdotTriplePrime*Hlg

there is a term in TEqn.H in line 19

- fvm::Sp(fvc::div(phi),T)

which is not corresponding to any term in energy equation.


The whole left side of TEqn is as following:

    fvScalarMatrix TEqn
    (
          fvm::ddt(T)
        + fvm::div(phi, T)
        - fvm::Sp(fvc::div(phi),T)
        - Foam::fvm::laplacian( k/rhoC_ , T,"laplacian(alphaEff,T)")
    /*    + (
            fvc::div(fvc::absolute(phi, U), p)
          + fvc::ddt(rho, KE) + fvc::div(rhoPhi, KE)
          )/rhoC
    */
    );

Would you please explain a bout it?

Regards,

Elham

[wavefunction]

Hi Nima Sam.

Could you kindly send me your solver to the email ID saicharanb56@gmail.com ? I would really appreciate it.

I am trying to simulate sessile droplet evaporation and I believe your solver will be ideal for me. I am, however, new to OpenFOAM and do not possess enough skill to code my own solver.

Thanks in advance.

[nimasam]

Dear Saicharan

please look this post

[nimasam]

Originally Posted by Elham

Dear Niam,

Regarding to phaseChangeHeatFoam the energy equation (25) in your paper is :

ddT + div (UT) - lapplacian (k/rhoC, T) = -1/rhoC *mdotTriplePrime*Hlg

there is a term in TEqn.H in line 19

- fvm::Sp(fvc::div(phi),T)

which is not corresponding to any term in energy equation.


The whole left side of TEqn is as following:

    fvScalarMatrix TEqn
    (
          fvm::ddt(T)
        + fvm::div(phi, T)
        - fvm::Sp(fvc::div(phi),T)
        - Foam::fvm::laplacian( k/rhoC_ , T,"laplacian(alphaEff,T)")
    /*    + (
            fvc::div(fvc::absolute(phi, U), p)
          + fvc::ddt(rho, KE) + fvc::div(rhoPhi, KE)
          )/rhoC
    */
    );

Would you please explain a bout it?

Regards,

Elham

fvc::div(phi)=0, so theoretically it is zero, but numerically it is not, and it help numerical stability

[Elham]

Thanks.

Elham

[gaza]

Originally Posted by Elham

Dear Niam,

Regarding to phaseChangeHeatFoam the energy equation (25) in your paper is :

ddT + div (UT) - lapplacian (k/rhoC, T) = -1/rhoC *mdotTriplePrime*Hlg

there is a term in TEqn.H in line 19

- fvm::Sp(fvc::div(phi),T)

which is not corresponding to any term in energy equation.


The whole left side of TEqn is as following:

    fvScalarMatrix TEqn
    (
          fvm::ddt(T)
        + fvm::div(phi, T)
        - fvm::Sp(fvc::div(phi),T)
        - Foam::fvm::laplacian( k/rhoC_ , T,"laplacian(alphaEff,T)")
    /*    + (
            fvc::div(fvc::absolute(phi, U), p)
          + fvc::ddt(rho, KE) + fvc::div(rhoPhi, KE)
          )/rhoC
    */
    );

Would you please explain a bout it?

Regards,

Elham

Dear Elham,
It is also written in the documentation
https://openfoam.org/release/2-2-0/n...s-boundedness/


According to this link the same can be achieved with
setting divergence scheme as bounded in fvScheme file, e.g.

 div(phi,T)      bounded Gauss upwind;

[saicharan662000@gmail.com]

Hi nimasam,
Could you tell me how you solved your problem if not possible to send .If you can send you can mail me to this mail : saicharan662000@gmail.com
Thanks and regards

[mikethe1wheelnut]

Originally Posted by nimasam

Dear foamer

The latest version of my solver (phaseChangeHeatFoam), several test cases and published papers are available in github.

https://github.com/NimaSam/phaseChangeHeatFoam/


please inform me about possible bugs.


Best Regards

To anybody like me finding this post all these years later, it uses openfoam 2.2.0, for which installation instructions exist only to up only to ubuntu 13.10, which is, in this case, unsupported and obsolete. If you follow those instructions, they will not succeed. Because the files they depend on won't be there. I have contacted the author about this and been directed to.. I can't recall where, (..I could track it down..) but it was one of the documentation pages you'll find if you search for evaporation/concensation/openfoam. So this is merely a confirmation of what you surely already expect. I can vouch that the author is very friendly if you are reasonable with whatever request you might make :-)


Best regards, Mike :-)

[BUAA LSM]

Please clicHello, Nima Sam.

Would you be so kind to send your code to zy2104210@buaa.edu.cn?
Thanks in advance!k one of the Quick Reply icons in the posts above to activate Quick Reply.

[nimasam]

Originally Posted by BUAA LSM

Please clicHello, Nima Sam.

Would you be so kind to send your code to zy2104210@buaa.edu.cn?
Thanks in advance!k one of the Quick Reply icons in the posts above to activate Quick Reply.

Hi dear foamer
it is interesting that after 6-7 years still i can see some people interested in phaseChangeHeatFoam. The solver is validating for film boiling alot with different authors. However, It was written for OpenFOAM-2.3 and I never tried to upgrade it higher version of OpenFOAM.


https://github.com/NimaSam/phaseChangeHeatFoam





but you can find lots of attempt in literature for developing phase change solver in OpenFOAM. e.g:


1- OpenFOAM standard solver (interCondensatingEvaporatingFoam from esi version)


2- boilingFoam (https://github.com/fmuni/boilingFoam-PUBLIC)