[viraj20feb]

Hello everyone,

I have three phases solid, liquid, and gas. I want to simulate heat transfer between all three phases, and mass transfer between liquid and gas. I am using multiphaseEulerFoam solver in OF-8. My simulation breaks down immediately in the first time step and gives me the following error:

#0 Foam::error:rintStack(Foam::Ostream&) at ~/OpenFOAM-8/OpenFOAM-8/src/OSspecific/POSIX/printStack.C:218
#1 Foam::sigFpe::sigHandler(int) at ~/OpenFOAM-8/OpenFOAM-8/src/OSspecific/POSIX/signals/sigFpe.C:104
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ~/OpenFOAM-8/OpenFOAM-8/src/OpenFOAM/fields/Fields/scalarField/scalarField.C:114 (discriminator 2)
#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/viraj/OpenFOAM-8/OpenFOAM-8/platforms/linux64GccDPInt32Debug/bin/multiphaseEulerFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ~/OpenFOAM-8/OpenFOAM-8/src/OpenFOAM/lnInclude/GeometricScalarField.C:107 (discriminator 8)
#6 Foam:hasePair::Pr() const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/phaseSystems/phasePair/phasePair/phasePair.C:135
#7 Foam::heatTransferModels::RanzMarshall::K(double) const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/interfacialModels/heatTransferModels/RanzMarshall/RanzMarshall.C:65
#8 Foam::heatTransferModel::K() const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/interfacialModels/heatTransferModels/heatTransferModel/heatTransferModel.C:78
#9 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::BlendedInterfacialModel<Foam::heatTransferMo del>::evaluate<double, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > (Foam::heatTransferModel::*)() const, Foam::word const&, Foam::dimensionSet const&, bool) const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseSystems/../../phaseSystems/lnInclude/BlendedInterfacialModel.C:175 (discriminator 6)
#10 Foam::BlendedInterfacialModel<Foam::heatTransferMo del>::K() const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseSystems/../../phaseSystems/lnInclude/BlendedInterfacialModel.C:397
#11 Foam::TwoResistanceHeatTransferPhaseSystem<Foam::M omentumTransferPhaseSystem<Foam:haseSystem> >::heatTransfer() const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseSystems/../../phaseSystems/lnInclude/TwoResistanceHeatTransferPhaseSystem.C:462 (discriminator 3)
#12 Foam::PhaseTransferPhaseSystem<Foam::TwoResistance HeatTransferPhaseSystem<Foam::MomentumTransferPhas eSystem<Foam:haseSystem> > >::heatTransfer() const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseSystems/../../phaseSystems/lnInclude/PhaseTransferPhaseSystem.C:332
#13 Foam::InterfaceCompositionPhaseChangePhaseSystem<F oam::PhaseTransferPhaseSystem<Foam::TwoResistanceH eatTransferPhaseSystem<Foam::MomentumTransferPhase System<Foam:haseSystem> > > >::heatTransfer() const at ~/OpenFOAM-8/OpenFOAM-8/applications/solvers/multiphase/multiphaseEulerFoam/multiphaseEulerFoam/multiphaseSystems/../../phaseSystems/lnInclude/InterfaceCompositionPhaseChangePhaseSystem.C:455
#14 ? in "/home/viraj/OpenFOAM-8/OpenFOAM-8/platforms/linux64GccDPInt32Debug/bin/multiphaseEulerFoam"
#15 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#16 ? in "/home/viraj/OpenFOAM-8/OpenFOAM-8/platforms/linux64GccDPInt32Debug/bin/multiphaseEulerFoam"

I believe that the error is in the phaseProperties file (attached). Can anyone please guide me in the right direction in order to get rid of this error?

Thank you!

[Jamessmp23]

Have you solved your case ?

How do you apply mass transfer equation in your case ? I'm trying to apply my mass transfer equation in K modelled but I found only a way to change mass transfer in empirical model.

Thanks in advanced.

[viraj20feb]

Originally Posted by Jamessmp23

Have you solved your case ?

How do you apply mass transfer equation in your case ? I'm trying to apply my mass transfer equation in K modelled but I found only a way to change mass transfer in empirical model.

Thanks in advanced.

Hi,

Yes, I solved my case. I had to use a blending model since I have multiple phases. I don't understand what your question is regarding the mass transfer equation. I am using the Frossling model to model mass transfer available in OpenFOAM. You may want to lookup interfacePhaseChangePhaseSystem to understand the solution procedure for mass transfer.

Best,
Viraj

[smoreno84]

I also want to simulate heat transfer between three phases, and mass transfer liquid and gas using multiphaseEulerFoam in OpenFoam V8.

I start executing and the solvers stops at the first iteration without saying anything.

log.multiphaseEulerFoam file and phaseProperties are as attach Files

Can anyone help me with this? I really apreciate your guides. thanks.


Reading field p_rgh

Courant Number mean: 0.01374 max: 0.01374

Starting time loop

Courant Number mean: 0.01374 max: 0.01374
Time = 0.001

PIMPLE: Iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.395966 0.336801 0.590764
water fraction, min, max = 3.05333e-06 0 0.05
Phase-sum volume fraction, min, max = 1.00002 0.909946 1.1632
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.395963 0.336801 0.590535
water fraction, min, max = 6.0958e-06 0 0.05
Phase-sum volume fraction, min, max = 1.00002 0.909946 1.1632
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.39596 0.336801 0.590307
water fraction, min, max = 9.12743e-06 0 0.05
Phase-sum volume fraction, min, max = 1.00002 0.909946 1.1632
smoothSolver: Solving for H2O.air, Initial residual = 0, Final residual = 0, No Iterations 1
Constructing momentum equations
Tf.airAndWater: min = 325, mean = 325, max = 325

[Kenmaster]

Originally Posted by smoreno84

I also want to simulate heat transfer between three phases, and mass transfer liquid and gas using multiphaseEulerFoam in OpenFoam V8.

I start executing and the solvers stops at the first iteration without saying anything.

log.multiphaseEulerFoam file and phaseProperties are as attach Files

Can anyone help me with this? I really apreciate your guides. thanks.


Reading field p_rgh

Courant Number mean: 0.01374 max: 0.01374

Starting time loop

Courant Number mean: 0.01374 max: 0.01374
Time = 0.001

PIMPLE: Iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.395966 0.336801 0.590764
water fraction, min, max = 3.05333e-06 0 0.05
Phase-sum volume fraction, min, max = 1.00002 0.909946 1.1632
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.395963 0.336801 0.590535
water fraction, min, max = 6.0958e-06 0 0.05
Phase-sum volume fraction, min, max = 1.00002 0.909946 1.1632
MULES: Solving for alpha.air
MULES: Solving for alpha.water
air fraction, min, max = 0.39596 0.336801 0.590307
water fraction, min, max = 9.12743e-06 0 0.05
Phase-sum volume fraction, min, max = 1.00002 0.909946 1.1632
smoothSolver: Solving for H2O.air, Initial residual = 0, Final residual = 0, No Iterations 1
Constructing momentum equations
Tf.airAndWater: min = 325, mean = 325, max = 325

Did you solve this?