[Elham]

Originally Posted by vonboett

one note though: interMixingFoam is either keeping the overall mass right but allows negative concentrations or volumetric concentrations greater than one, or you have the concentrations limited to the physical range but may loose mass...
I tried to optimize this a bit in debrisInterMixing-2.3, you can download the code with the corresponding paper if you want to do the same.

Thanks for the comprehensive reply. I am going to model evaporation of water in air. So I have two phases at first then there will be three and there are mass transfer between two of them. Based on what you described I cannot use interMixingFoam since I need to have three different continuity equations. Please correct me if I am wrong.

Regards,

Elham

[Elham]

Originally Posted by Elham

Thanks for the comprehensive reply. I am going to model evaporation of water in air. So I have two phases at first then there will be three and there are mass transfer between two of them. Based on what you described I cannot use interMixingFoam since I need to have three different continuity equations. Please correct me if I am wrong.

Regards,

Elham

Dear Albrecht

I have added source term on MULES of alpha2 and the fvScalarMatrix alpha2Eqn. The solver compiled but it crashed at the time of running.

#0  Foam::error::printStack(Foam::Ostream&) at ~/OpenFOAM/OpenFOAM-2.3.1/src/OSspecific/POSIX/printStack.C:221
#1  Foam::sigFpe::sigHandler(int) at ~/OpenFOAM/OpenFOAM-2.3.1/src/OSspecific/POSIX/signals/sigFpe.C:108
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#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 at ~/OpenFOAM/OpenFOAM-2.3.1/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGSolverScale.C:57
#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 at ~/OpenFOAM/OpenFOAM-2.3.1/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGSolverSolve.C:370 (discriminator 1)
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ~/OpenFOAM/OpenFOAM-2.3.1/src/OpenFOAM/matrices/lduMatrix/solvers/GAMG/GAMGSolverSolve.C:119
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ~/OpenFOAM/OpenFOAM-2.3.1/src/finiteVolume/fvMatrices/fvScalarMatrix/fvScalarMatrix.C:169 (discriminator 1)
#7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ~/OpenFOAM/OpenFOAM-2.3.1/src/finiteVolume/lnInclude/fvMatrixSolve.C:82
#8  
 at ~/OpenFOAM/elham-2.3.1/applications/solvers/multiphase/interFoam/interMixingHeatFoam/pEqn.H:58 (discriminator 5)
#9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10  
 at ??:?
Floating point exception (core dumped)
elham@elham-Precision-T7610:~/OpenFOAM/elham-2.3.1/tutorials/multiphase/interMixingHeatFoam/droplet$

Do you think the problem is in the code? Should I change the fluxes as well?

Regards,

Elham

[vonboett]

have you tried interPhaseChangeFoam?

[Elham]

Originally Posted by vonboett

have you tried interPhaseChangeFoam?

I could fix it. Since I have three phases I cannot use interPhaseChangeFoam. Now I am running three phases with source term. Thanks for the alphaEqns.H file that you provided.

Regards,

Elham

[Elham]

Originally Posted by vonboett

one note though: interMixingFoam is either keeping the overall mass right but allows negative concentrations or volumetric concentrations greater than one, or you have the concentrations limited to the physical range but may loose mass...
I tried to optimize this a bit in debrisInterMixing-2.3, you can download the code with the corresponding paper if you want to do the same.

Dear Albrecht,

In interMixingFoam there are diffusion terms in alpha equations,eg,

- fvm::laplacian(Dc23 + Dc32 + alphatab*turbulence->nut() , alpha3)

for alpha3, but they do not exist in the MULES equations. Would you please explain the reason? I tried to add diffusion term in the MULES but I got the error that there is not any off-diagonal term in the matrix.

Regards,

Elham

[vonboett]

Hi,

Interesting version you found there, where I experimented with eddy diffusivity... However, the code I ment is in the supplement of
https://www.geosci-model-dev.net/10/3963/2017/
and it is ment as "inspiration" and can't be taken as 1:1 copy to interMixingFoam because it has implemented other stuff to deal with debris flows. Phase1 stays separate, phase 2 and 3 can mix, but do not have momentum transfer.

[Elham]

Originally Posted by vonboett

Hi,

Interesting version you found there, where I experimented with eddy diffusivity... However, the code I ment is in the supplement of
https://www.geosci-model-dev.net/10/3963/2017/
and it is ment as "inspiration" and can't be taken as 1:1 copy to interMixingFoam because it has implemented other stuff to deal with debris flows. Phase1 stays separate, phase 2 and 3 can mix, but do not have momentum transfer.

Dear Albrecht,

I have checked your code. In the debrisInterMixingFoam.C you have added a resetting velocity and pressure feature in the PISO loop after outerCorrection loop to prevent melting or freezing of the material. It was a problem that I encountered and is really helpful. But I am after the origin of the error and is there any resource to give more insight?

Cheers,

Elham

[astroda]

Originally Posted by Elham

Thanks for the comprehensive reply. I am going to model evaporation of water in air. So I have two phases at first then there will be three and there are mass transfer between two of them. Based on what you described I cannot use interMixingFoam since I need to have three different continuity equations. Please correct me if I am wrong.

Regards,

Elham

Hi Elham,

How is this going, I am interested in what became of it! I am looking to model the distillation process, so from a miscible mixture to a pure vapour phase change.

Is there a working (or non working) version of your solver?

Kindly,

[Elham]

Originally Posted by astroda

Hi Elham,

How is this going, I am interested in what became of it! I am looking to model the distillation process, so from a miscible mixture to a pure vapour phase change.

Is there a working (or non working) version of your solver?

Kindly,

Dear Aaron,

My model based on interMixingFoam is working now, although not published yet. What else do you need to know now?