[Kanarya]

I did compile it as well but doesnt work...
could you send me your case only not the code if you do not mind...
I can try my code in your working case
thanks in advance!!!

Originally Posted by mihaipruna

As I posted earlier in the thread, to get it to work for me, I had to recompile
phaseChangeTwoPhaseMixture

[mihaipruna]

why not test it like they did in the tutorial, by removing the bullet from the cavitatingBullet case and setting a temperature BC?

I don't have access to my case files right now.

[Kanarya]

Hi,
I did the test with their tutorial...when I set the alpha=0, then nothing is changing and alpha remains always zero which is wrong right?
But when I put some value to alpha like 0.00001, it seems it does something...
did you try to give alpha=0 for initial condition?

Thanks!

Originally Posted by mihaipruna

why not test it like they did in the tutorial, by removing the bullet from the cavitatingBullet case and setting a temperature BC?

I don't have access to my case files right now.

[mihaipruna]

no, I had initial alpha=1.I was modeling boiling.
How many time steps did you run it for?

[Kanarya]

I run for 0.05s...
did you read the master thesis Roman Thiele?
he modified it with different models for condensation...
what do you think?
are you still interested in this topic?

thanks!

Originally Posted by mihaipruna

no, I had initial alpha=1.I was modeling boiling.
How many time steps did you run it for?

[Kanarya]

HI,

I modified the solver so called interPhaseChangeFoam by following the instructions form Andersen tutorial...Is is work fine with cavitaingBullet case but not in simple pipe flow..it gives me following message:

Selecting phaseChange model SchnerrSauer
Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading transportProperties

Calculating field g.h


PIMPLE: Operating solver in PISO mode

time step continuity errors : sum local = 1.15855e-05, global = -1.19144e-16, cumulative = -1.19144e-16
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 3.48266e+19, No Iterations 50
GAMG:  Solving for pcorr, Initial residual = 0.965555, Final residual = 1.08509e+21, No Iterations 50
time step continuity errors : sum local = 4.51124e+35, global = -1.41141e+28, cumulative = -1.41141e+28
Courant Number mean: 3.09991e+35 max: 5.41686e+38

Starting time loop

Courant Number mean: 0.00057227 max: 1
deltaT = 1.84609e-47
Time = 1.84609e-47

#0  Foam::error::printStack(Foam::Ostream&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/usr/local/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  
 at sigaction.c:0
#3  Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/usr/local/OpenFOAM/OpenFOAM-2.1.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4  void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#5  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#6  
 in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#7  
 in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#8  
 in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#9  
 in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
#10  __libc_start_main in "/lib64/libc.so.6"
#11  
 in "/home/links/rek209/OpenFOAM/rek209-2.1.0/platforms/linux64GccDPOpt/bin/myInterPhaseChangeFoam"
Floating point exception

I think, there is problem with the BC for p_rgh but I am not sure
And BC I have used :

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 1e5;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

   
     outlet
    {
        type            fixedValue;
        value           $internalField;
    }

    hot_wall
    {
       type            zeroGradient;
        
    }
    ocold_wall
    {
        type            zeroGradient;
        
    }
    icold_wall
    {
      type            zeroGradient;
       
    }

}

// ************************************************************************* //

and for alpha:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      alpha1;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform  0.00001;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
    }

    walls
    {
        type            symmetry;
    }

}

// ************************************************************************* //

Can you help me with this problem or any suggestions?

Thanks in advance!

[Kanarya]

Hi,

is there anybody interested in my problem?

[wyldckat]

Originally Posted by Kanarya

is there anybody interested in my problem?

Please follow the instructions given on this thread: http://www.cfd-online.com/Forums/ope...-get-help.html
Otherwise, most people will not even invest time in asking you for more information

[MiladBagheri]

Hello everyone,

I know this is an old thread, nevertheless I would like to share the solver here. I hope it is useful for someone.

You can download it from my GitHub:

https://github.com/MiladBagheri/Nume...ng-Regimes.git

You may also check some of the videos I made using this solver to simulate different boiling regimes.

https://www.youtube.com/channel/UCif...AINBMI_rFifgVA