[Kobi_Faust]

Hi all,

I'm having a multi region model, solving it with the solver chtMultiRegionFoam.
Everything runs well so far. Now I did some changes on the thermophysicalProperties of the fluid region. With this I'm getting a "Floating point exception" error. Any ideas where this comes from? Or what I am doing wrong?

Attached the old and new thermophysicalProperties file, and the terminal log.


Old thermophysicalProperties:

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

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectFluid;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        molWeight       18; //g per mol
		nMoles			1; 
    }
    equationOfState
    {
        R           3000;
        rho0        1027;
    }
    thermodynamics
    {
        Cp          4183; 
        Hf          0;
    }
    transport
    {
        mu              959e-6; 
        Pr              4.56; 
    }
}

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

New thermophysicalProperties:

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

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState icoPolynomial;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        molWeight       18; //g per mol
		nMoles			1; 
    }
    equationOfState
    {
        rhoCoeffs<8> ( 765.33 1.8142 -0.0035 0 0 0 0 0);
    }
    thermodynamics
    {
               Hf          0;
		Sf			0;
		CpCoeffs<8> ( 28070 -281.7 1.25 -0.00248 0.000001857 0 0 0);
    }
    transport
    {
        	muCoeffs<8> (0.0967 -0.0008207 0.000002344 -0.000000002244 0 0 0 0);
		kappaCoeffs<8> (-0.5752 0.006397 -0.000008151 0 0 0 0 0);
    }
}

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

log-file:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1812                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1812 OPENFOAM=1812
Arch   : "LSB;label=32;scalar=64"
Exec   : mychtMultiRegionFoam
Date   : May 31 2019
Time   : 13:51:07
Host   : hmc01.cluster.frm2.tum.de
PID    : 38852
I/O    : uncollated
Case   : /localhome/kfaust/Coupling_V6.5/OF_V6.5
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region Wasser for time = 0

Create solid mesh for region Meat_15 for time = 0

Create solid mesh for region Meat_30 for time = 0

Create solid mesh for region Plattenrahmen for time = 0

Create solid mesh for region Cladding_aussen for time = 0

Create solid mesh for region Cladding_innen for time = 0

Create solid mesh for region Kamm_oben for time = 0

Create solid mesh for region Kamm_unten for time = 0

Create solid mesh for region Rand_aussen for time = 0

Create solid mesh for region Rand_innen for time = 0

*** Reading fluid mesh thermophysical properties for region Wasser

    Adding to thermoFluid

Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       polynomial;
    thermo          hPolynomial;
    equationOfState icoPolynomial;
    specie          specie;
    energy          sensibleEnthalpy;
}

AMI: Creating addressing and weights between 657432 source faces and 656858 target faces
AMI: Patch source sum(weights) min:0.98604664 max:1.000062 average:0.99998044
AMI: Patch target sum(weights) min:0.98186748 max:1.0001693 average:0.99998884
    Adding to rhoFluid

    Adding to UFluid

    Adding to phiFluid

    Adding to hRefFluid

    Adding to ghFluid

    Adding to ghfFluid

    Adding to turbulenceFluid

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
    Adding to reactionFluid

Combustion model not active: combustionProperties not found
Selecting combustion model none
    Adding to radiationFluid

Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding to fieldsFluid

    Adding to QdotFluid

    Adding MRF

No MRF models present

    Adding fvOptions

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libc.so.6
#3  logf32x in /lib/x86_64-linux-gnu/libm.so.6
#4  Foam::compressible::alphatJayatillekeWallFunctionFvPatchScalarField::updateCoeffs() at ??:?
#5  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::evaluate() at ??:?
#6  Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > >::correctNut() at ??:?
#7  Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#8  ? at ??:?
#9  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#10  ? at ??:?
Floating point exception (core dumped)

Tank you for your help.
Kobi

[mwaqas]

Hello Kobi,


I never used polynomial function in chtMultiregionFOAM. But it seems like Prandtl number is missing in your new thermophysicalProperties. I am not really sure if it is not required in case of polynomial functions of physical properties.


Regards

[Kobi_Faust]

Hello Muhammad,

thanks for your answer. Thats true, that the Prandtl number is missing. But it didn't change when adding it.
If only one entry would be missing, the failure should be something like "keyword missing".

So it must be some different mistake.

I would appreciate further ideas.

Regards,
Kobi

[jherb]

The crash happens in

alphatJayatillekeWallFunctionFvPatchScalarField

You could try to find the exact location and the values of the variables/fields?

Also have you tried to decrease relaxation factors and/or time step size?