[bFOAMer]

Hi all!
I'm a Foamer beginner and I hope you can help me. I have been trying to solve a heat transfer case with chtMultiRegionSimpleFoam and kOmegaSST turbulence model but it doesn't work. The error message is :

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.0.1-d8a290b55d28
Exec : chtMultiRegionSimpleFoam
Date : Jul 11 2016
Time : 18:14:06
Host : "debvbox"
PID : 3085
Case : /home/pepito/OpenFOAM/pepito-3.0.1/run/scambiatore002
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create fluid mesh for region aria for time = 0

Create fluid mesh for region co2 for time = 0

Create solid mesh for region alluminio for time = 0

*** Reading fluid mesh thermophysical properties for region aria

Adding to thermoFluid

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}

Adding to rhoFluid

Adding to UFluid

Adding to phiFluid

Adding to gFluid

Adding to hRefFluid

Adding to ghFluid

Adding to ghfFluid

Adding to turbulence

Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSST
Selecting patchDistMethod meshWave
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#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 ??:?
#4 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::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5 Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::F2() const at ??:?
#6 Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::F23() const at ??:?
#7 Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::correctNut(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#8 Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#9 Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::kOmegaSST(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&, Foam::word const&) at ??:?
#10 Foam::RASModel<Foam::EddyDiffusivity<Foam::Thermal Diffusivity<Foam::CompressibleTurbulenceModel<Foam ::fluidThermo> > > >::adddictionaryConstructorToTable<Foam::RASModels ::kOmegaSST<Foam::EddyDiffusivity<Foam::ThermalDif fusivity<Foam::CompressibleTurbulenceModel<Foam::f luidThermo> > > > >::New(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#11 Foam::RASModel<Foam::EddyDiffusivity<Foam::Thermal Diffusivity<Foam::CompressibleTurbulenceModel<Foam ::fluidThermo> > > >::New(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#12 Foam::TurbulenceModel<Foam::geometricOneField, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>, Foam::compressibleTurbulenceModel, Foam::fluidThermo>::adddictionaryConstructorToTabl e<Foam::RASModel<Foam::EddyDiffusivity<Foam::Therm alDiffusivity<Foam::CompressibleTurbulenceModel<Fo am::fluidThermo> > > > >::NewTurbulenceModel(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::fluidThermo const&, Foam::word const&) at ??:?
#13 ? at ??:?
#14 ? at ??:?
#15 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#16 ? at ??:?
Floating point exception


I don't understand where the error is. Can you say me how I can solve this problem?

Thank you for your help!

[alexeym]

Hi,

To trace the origin of the error, you need to look at line #3 and line #5.

#3

Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

The error is in division, usually this happens when denominator is 0.

#5

Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::F2() const at ??:?

This is the function where error happened (it you compile OF in debug mode, backtrace will also contain file names and line numbers instead of ?:??).

F2 function could be found at https://github.com/OpenFOAM/OpenFOAM...OmegaSST.C#L68

template<class BasicTurbulenceModel>
tmp<volScalarField> kOmegaSST<BasicTurbulenceModel>::kOmegaSST::F2() const
{
    tmp<volScalarField> arg2 = min
    (
        max
        (
            (scalar(2)/betaStar_)*sqrt(k_)/(omega_*y_),
            scalar(500)*(this->mu()/this->rho_)/(sqr(y_)*omega_)
        ),
        scalar(100)
    );

    return tanh(sqr(arg2));
}

There are betaStar_, omega_ and y_ in denominators. So I guess, you initial or boundary conditions contain 0 in omega somewhere, as betaStart_ is 0.09 and y_ is in general positive.

[bFOAMer]

Originally Posted by alexeym

Hi,

To trace the origin of the error, you need to look at line #3 and line #5.

#3

Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

The error is in division, usually this happens when denominator is 0.

#5

Foam::RASModels::kOmegaSST<Foam::EddyDiffusivity<F oam::ThermalDiffusivity<Foam::CompressibleTurbulen ceModel<Foam::fluidThermo> > > >::F2() const at ??:?

This is the function where error happened (it you compile OF in debug mode, backtrace will also contain file names and line numbers instead of ?:??).

F2 function could be found at https://github.com/OpenFOAM/OpenFOAM...OmegaSST.C#L68

template<class BasicTurbulenceModel>
tmp<volScalarField> kOmegaSST<BasicTurbulenceModel>::kOmegaSST::F2() const
{
    tmp<volScalarField> arg2 = min
    (
        max
        (
            (scalar(2)/betaStar_)*sqrt(k_)/(omega_*y_),
            scalar(500)*(this->mu()/this->rho_)/(sqr(y_)*omega_)
        ),
        scalar(100)
    );

    return tanh(sqr(arg2));
}

There are betaStar_, omega_ and y_ in denominators. So I guess, you initial or boundary conditions contain 0 in omega somewhere, as betaStart_ is 0.09 and y_ is in general positive.

Hi alexeym, thank you very much indeed for your precious help!
Now I have understood what the message means. Unfortunately I don't know the way to compile OpenFOAM in debug mode because I'm still learning to use OF but I will be grateful if you help me. I have already looked for some numerical error in omega files but i haven't found anything. I can upload the casefolder.zip if you prefer.

I'm sorry for the trouble and thank you again and again!

[alexeym]

Upload your case, otherwise it become rather "guess what is wrong" game.

Unfortunately I don't know the way to compile OpenFOAM in debug mode because I'm still learning to use OF but I will be grateful if you help me.

Well, in this case file names and line numbers are of no use.

[bFOAMer]

Originally Posted by alexeym

Upload your case, otherwise it become rather "guess what is wrong" game.



Well, in this case file names and line numbers are of no use.

you can find my case at:
https://www.dropbox.com/s/h2iqrnit2q...er.tar.gz?dl=0

[alexeym]

Unfortunately I do not have 3.0.1 installed and with 3.0.x I can not reproduce your error. After corrections to boundary conditions (since Feb. 24 value field is required for calculated BC, co2/k/contact_region_3-trg type should be kqRWallFunction not kqrWallFunction, co2/alphat/contact_region_3-trg type should be compressible::alphatJayatillekeWallFunction) I get the following output:

...
Build  : 3.0.x-1c517d6c8485
Exec   : chtMultiRegionSimpleFoam
...
Time = 1


Solving for fluid region aria
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00979818, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.006705962, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.0005758433, No Iterations 1


--> FOAM FATAL ERROR: 
Supply either a patchName or a coupleGroup for patch contact_region-src in region aria

[bFOAMer]

Originally Posted by alexeym

Unfortunately I do not have 3.0.1 installed and with 3.0.x I can not reproduce your error. After corrections to boundary conditions (since Feb. 24 value field is required for calculated BC, co2/k/contact_region_3-trg type should be kqRWallFunction not kqrWallFunction, co2/alphat/contact_region_3-trg type should be compressible::alphatJayatillekeWallFunction) I get the following output:

...
Build  : 3.0.x-1c517d6c8485
Exec   : chtMultiRegionSimpleFoam
...
Time = 1


Solving for fluid region aria
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00979818, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.006705962, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.0005758433, No Iterations 1


--> FOAM FATAL ERROR: 
Supply either a patchName or a coupleGroup for patch contact_region-src in region aria

Thank you so much Alexey, you are a genius I believe! However after correction I got the same output so I think I've also made a mistake in the setting of thermal coupling... I'm desperate... I don't know how to solve this new problem too...

[bFOAMer]

ok, I think I've just found the solution to my problem...
thanks for your precious help Alexey!!!

[yangzhuan]

Originally Posted by alexeym

Hi,

To trace the origin of the error, you need to look at line #3 and line #5.

#3

Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?

The error is in division, usually this happens when denominator is 0.

Hi, I meet a question, and I have confused for a long time. I really don't know how to solve it,Could you help me analyze the mistake?

Courant Number mean: 3.40147e-05 max: 0.502766
Interface Courant Number mean: 1.48331e-05 max: 0.0517021
deltaT = 4.061e-68
Time = 4.71186
PIMPLE: iteration 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03688e-06, Final residual = 2.93968e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03667e-06, Final residual = 2.93849e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03639e-06, Final residual = 2.93744e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = -2.21102e-28  Max(alpha.water) = 1
smoothSolver:  Solving for alpha.water, Initial residual = 2.03613e-06, Final residual = 2.93615e-11, No Iterations 1
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.956421  Min(alpha.water) = 0  Max(alpha.water) = 1
DICPCG:  Solving for p_rgh, Initial residual = 0.00709543, Final residual = 0.000345669, No Iterations 8
time step continuity errors : sum local = 5.43246e-10, global = 1.84738e-10, cumulative = 2.56315e-06
DICPCG:  Solving for p_rgh, Initial residual = 0.00538723, Final residual = 6.52414e-08, No Iterations 224
time step continuity errors : sum local = 1.02756e-13, global = -6.67559e-16, cumulative = 2.56315e-06
[3] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #1  Foam::sigFpe::sigHandler(int) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #2  ? in "/lib64/libc.so.6"
[3] #3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #4  double Foam::gSumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&, int) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #5  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libOpenFOAM.so"
[3] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libfiniteVolume.so"
[3] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"
[3] #8  Foam::fvMatrix<double>::solve() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"
[3] #9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so"
[3] #10  Foam::RASModels::realizableKE<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/lib/libincompressibleTurbulenceModels.so"
[3] #11  ? in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"
[3] #12  __libc_start_main in "/lib64/libc.so.6"
[3] #13  ? in "/opt/software/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64GccDPInt32Prof/bin/interFoam"

[alexvaleije]

Originally Posted by bFOAMer

ok, I think I've just found the solution to my problem...
thanks for your precious help Alexey!!!

Hi!!

I know it's been a long time since this happened, but do you remember how did you solve this problem? I'm having the same issue and I haven't found any solution yet.

Regards,
Alex

[batteryFoamer]

Originally Posted by alexvaleije

Hi!!

I know it's been a long time since this happened, but do you remember how did you solve this problem? I'm having the same issue and I haven't found any solution yet.

Regards,
Alex

I also met the same question, could you share how to solve this problem?

[batteryFoamer]

Originally Posted by bFOAMer

ok, I think I've just found the solution to my problem...
thanks for your precious help Alexey!!!

I also met the same question, could you share how to solve this problem?

[Ary]

Originally Posted by batteryFoamer

I also met the same question, could you share how to solve this problem?

It is a long time again...

My solver compilation is OK, but

I met similar errors and the computation was crushed finally.

The problem was solved when I fixed UEqn where the pressure term shall be a minus term.

So, maybe you shall check your equations carefully.

May your computation converge!

Ary