Error appears when using chtMultiRegionFoam

ruffy19 · May 23, 2016, 05:19

Hello Foam-Experts!

I'm really new to OpenFoam and the programming in it and so I'm not experienced with it. But for my studies I got the following task:
I have to simulate a solid cylinder hanging in air. The cylinder radiates heat and convection is on hand.
So I use "chtMultiRegionFoam" with the command

Code:

mpirun -np 3 chtMultiRegionFoam -parallel

to run the solution with its 2 declared regions and 3 generated processors. I used the same way to solve the case for a functional project called "planWall2D" for which I assume that it is programmed correctly, because I got this form the openfoamwiki.
Then I realized, that the error-messages for both cases are the same. This error-message appears after the thermophysical properties for region air were read and added and ends up like:

Code:

mpirun noticed that process rank 2 with PID 477 on node boot2docker exited on signal 8 <Floating point exception>

.
Because an error-message like that appears for my case as well as for the model-programm, I guess the error is not my fault. So I would like to ask you all, what the problem actually is.

By the way I use the windows-version of OpenFoam.

I really please you to help me and also if you need some more information about my problem, please answer me as soon as possible.

Thanks in advance
Raffael

ruffy19 · May 25, 2016, 03:35

I really need to be helped with this problem. I don't know what the error message means and what to do about that. It is really important for me.
I'm not experienced with OpenFOAM and maybe I did something wrong. But until now I found nothing...

I really please you to help me!

Thanks in advance
Raffael

Antimony · May 25, 2016, 04:09

Hi,

The error message that you have is a Floating Point Exception. The most likely culprit for this is a division by zero somewhere.

Look through the different boundary condition files and see if there is something that shouldn't be zero is zero.

What does the log of your run say apart from the mpirun run abort?

Cheers,
Antimony

ruffy19 · May 25, 2016, 05:41

Hello Antimony,

At first I have to say, that I'm very thankful for you're reply!

A division by zero? That's a good point to look for, thanks for that. I'll post about the results of my upcoming search later.

The whole error-message looks like this:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : May 25 2016
Time   : 09:27:07
Host   : "boot2docker"
PID    : 255
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.256"
"boot2docker.257"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
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 air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    Adding to thermoFluid

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

--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 256 on node boot2docker exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

Between the selected thermophysical package and the abortion-message of mpirun is a huge list of incomprehensible error-messages. Do you know what I means? Are they important? Unfortunately I don't know how to copy them.

Again thank you!
Raffael

Antimony · May 25, 2016, 21:38

Hi,

You can get the entire output to a log file by adding "2>&1" at the end of the function call.

Since it seems to stop after loading your thermophysical properties, is your temperature 0 at any point? Do take care that the units of temperature are in K...

Cheers,
Antimony

ruffy19 · May 31, 2016, 03:58

Hello Antimony,

I checked for zeros in all initial fields, espacially for the temperatures.
Thereby I realized, that after every execution of chtMultiRegionFoam all temperatures are set on default (293K), although I've set them before on whole other values. I corrected it, but still this error-message appears:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : May 31 2016
Time   : 07:44:39
Host   : "boot2docker"
PID    : 803
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.804"
"boot2docker.805"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
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 air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    Adding to thermoFluid

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

[1] #0  [0] #0  [2] #0Foam::error::printStack(Foam::Ostream&)  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[2] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[1] #1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[0] #2  ? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[1] #2  ? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[2] #2  ? in "/lib64/libc.so.6"
[0] #3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/lib64/libc.so.6"
[1] #3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/lib64/libc.so.6"
[2] #3  Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::calculate() in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[0] #4  Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[2] #4  Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[1] #4  Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[0]  in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[1] #5  Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&)[2] #5  Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&)#5  Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[2] #6  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&)[1] #6  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[0] #6  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[0] #7  [1] #7   in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so"
[2] #7  ??? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[0] # in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[2] #8  __libc_start_main8  __libc_start_main[1] #8  __libc_start_main in "/lib64/libc.so.6"
[0] #9   in "/lib64/libc.so.6"
[1] #9   in "/lib64/libc.so.6"
[2] #9  ??? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiR in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
 in "/opt/OpeegionFoam"
[boot2docker:00803] *** Process received signal ***
[boot2docker:00803] Signal: Floating point exception (8)
[boot2docker:00803] Signal code:  (-6)
[boot2docker:00803] Failing at address: 0x1f400000323
[boot2docker:00803] [ 0] /lib64/libc.so.6(+0x326a0)[0x7f7016c636a0]
[boot2docker:00803] [ 1] /lib64/libc.so.6(gsignal+0x35)[0x7f7016c63625]
[boot2docker:00803] [ 2] /lib64/libc.so.6(+0x326a0)[0x7f7016c636a0]
[boot2docker:00803] [ 3] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11heRhoThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x646)[0x7f701cd2f3c6]
[boot2docker:00803] [ 4] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo27addfvMeshConstructorToTableINS_11heRhoThermoIS0_NS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEEEE3NewERKNS_6fvMeshERKNS_4wordE+0x5c)[0x7f701cd3b8dc]
[boot2docker:00803] [ 5] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11basicThermo3NewINS_9rhoThermoEEENS_7autoPtrIT_EERKNS_6fvMeshERKNS_4wordE+0x1ac)[0x7f701ccd4d1c]
[boot2docker:00803] [ 6] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo3NewERKNS_6fvMeshERKNS_4wordE+0x9)[0x7f701ccd3429]
[boot2docker:00803] [ 7] chtMultiRegionFoam[0x4a2024]
[boot2docker:00803] [ 8] /lib64/libc.so.6(__libc_start_main+0xfd)[0x7f7016c4fd5d]
[boot2docker:00803] [ 9] chtMultiRegionFoam[0x432619]
[boot2docker:00803] *** End of error message ***
[boot2docker:00805] *** Process received signal ***
[boot2docker:00805] Signal: Floating point exception (8)
[boot2docker:00805] Signal code:  (-6)
[boot2docker:00805] Failing at address: 0x1f400000325
[boot2docker:00805] [ 0] /lib64/libc.so.6(+0x326a0)[0x7efe48d2c6a0]
[boot2docker:00805] [ 1] /lib64/libc.so.6(gsignal+0x35)[0x7efe48d2c625]
[boot2docker:00805] [ 2] /lib64/libc.so.6(+0x326a0)[0x7efe48d2c6a0]
[boot2docker:00805] [ 3] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11heRhoThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x646)[0x7efe4edf83c6]
[boot2docker:00805] [ 4] nFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[boot2docker:00804] *** Process received signal ***
[boot2docker:00804] Signal: Floating point exception (8)
[boot2docker:00804] Signal code:  (-6)
[boot2docker:00804] Failing at address: 0x1f400000324
[boot2docker:00804] [ 0] /lib64/libc.so.6(+0x326a0)[0x7fc502d5f6a0]
[boot2docker:00804] [ 1] /lib64/libc.so.6(gsignal+0x35)[0x7fc502d5f625]
[boot2docker:00804] [ 2] /lib64/libc.so.6(+0x326a0)[0x7fc502d5f6a0]
[boot2docker:00804] [ 3] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11heRhoThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x646)[0x7fc508e2b3c6]
[boot2docker:00804] [ 4] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo27addfvMeshConstructorToTableINS_11heRhoThermoIS0_NS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEEEE3NewERKNS_6fvMeshERKNS_4wordE+0x5c)[0x7fc508e378dc]
[boot2docker:00804] [ 5] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11basicThermo3NewINS_9rhoThermoEEENS_7autoPtrIT_EERKNS_6fvMeshERKNS_4wordE+0x1ac)[0x7fc508dd0d1c]
[boot2docker:00804] [ 6] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo3NewERKNS_6fvMeshERKNS_4wordE+0x9)[0x7fc508dcf429]
[boot2docker:00804] [ 7] chtMultiRegionFoam[0x4a2024]
[boot2docker:00804] [ 8] /lib64/libc.so.6(__libc_start_main+0xfd)[0x7fc502d4bd5d]
[boot2docker:00804] [ 9] chtMultiRegionFoam[0x432619]
[boot2docker:00804] *** End of error message ***
/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo27addfvMeshConstructorToTableINS_11heRhoThermoIS0_NS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEEEE3NewERKNS_6fvMeshERKNS_4wordE+0x5c)[0x7efe4ee048dc]
[boot2docker:00805] [ 5] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam11basicThermo3NewINS_9rhoThermoEEENS_7autoPtrIT_EERKNS_6fvMeshERKNS_4wordE+0x1ac)[0x7efe4ed9dd1c]
[boot2docker:00805] [ 6] /opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo3NewERKNS_6fvMeshERKNS_4wordE+0x9)[0x7efe4ed9c429]
[boot2docker:00805] [ 7] chtMultiRegionFoam[0x4a2024]
[boot2docker:00805] [ 8] /lib64/libc.so.6(__libc_start_main+0xfd)[0x7efe48d18d5d]
[boot2docker:00805] [ 9] chtMultiRegionFoam[0x432619]
[boot2docker:00805] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 1 with PID 804 on node boot2docker exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

Could it be, that declarations as "empty" in the boundary conditions are considered as zero? Should I change them into specific values?

I have regions with different temperatures, so that a heat transfer takes place. Do I need a transferProperties-file in the constant-dictionary for that?

Thanks in advance

Raffael

P.S.: Thank you for the advice how to copy the whole error-message.

ruffy19 · June 2, 2016, 10:01

I used the command

Code:

unset FOAM_SIGFPE

to disable floating point exception (or something like that, I actually don't really know the effect of this command..).
After this command, the error-message looks totally different:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam
Date   : Jun 02 2016
Time   : 13:55:33
Host   : "boot2docker"
PID    : 314
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 1
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam
Date   : Jun 02 2016
Time   : 13:55:34
Host   : "boot2docker"
PID    : 315
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 1
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam
Date   : Jun 02 2016
Time   : 13:55:34
Host   : "boot2docker"
PID    : 316
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 1
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create solid mesh for region cable for time = 0

Create time

Create fluid mesh for region air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    Adding to thermoFluid

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



--> FOAM FATAL IO ERROR: 
inconsistent patch and patchField types for 
    patch type empty and patchField type calculated

file: /home/ofuser/workingDir/FoamCases/KabelInLuft/0/air/T.boundaryField.FrontAndBack from line 30 to line 31.

    From function static Foam::tmp<Foam::fvPatchField<Type> > Foam::fvPatchField<Type>::New(const Foam::fvPatch&, const Foam::DimensionedField<Type, Foam::volMesh>&, const Foam::dictionary&) [with Type = double]
    in file /opt/OpenFOAM/OpenFOAM-v3.0+/src/finiteVolume/lnInclude/fvPatchFieldNew.C at line 163.

FOAM exiting

-------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:

  Process name: [[40566,1],0]
  Exit code:    1
--------------------------------------------------------------------------

Maybe you can help me to find out what is wrong now?

Thanks in advance
Raffael

jmdf · June 4, 2016, 18:52

From the error it appears that you are using a "calculated" boundary condition in a boundary with the type empty.

If you need to use the calculated boundary, check if that boundary is defined as patch or wall on the "boundary" file (constant/<region>/polyMesh/boundary).

I hope it helps you!

ruffy19 · June 6, 2016, 04:22

Hello João!

I thank you so much for your advice. This problem is solved, but now I'm standing in front of the next problem.... like I said, I'm very new to this and know almost nothing about OpenFOAM.

Now appears the following error:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : Jun 06 2016
Time   : 08:07:39
Host   : "boot2docker"
PID    : 250
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.251"
"boot2docker.252"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    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 laminar
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

[2] [1] 
[1] 
[1] --> FOAM FATAL IO ERROR: 
[1] keyword PIMPLE is undefined in dictionary "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor1/system/air/fvSolution"
[1] 
[1] file: /home/ofuser/workingDir/FoamCases/KabelInLuft/processor1/system/air/fvSolution from line 0 to line 0.

[2] 
[2] --> FOAM FATAL IO ERROR: 
[2] keyword PIMPLE is undefined in dictionary "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor2/system/air/fvSolution"
[2] 
[2] file: /home/ofuser/workingDir/FoamCases/KabelInLuft/processor2/system/air/fvSolution from line 0 to line 0.
[2] 
[2]     From function const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&) const
[2]     in file db/dictionary/dictionary.C at line 640.
[2] 
FOAM parallel run exiting
[2] 
[1] 
[1]     From function const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&) const
[1]     in file db/dictionary/dictionary.C at line 640.
[1] 
FOAM parallel run exiting
[1] 
[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: 
[0] keyword PIMPLE is undefined in dictionary "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor0/system/air/fvSolution"
[0] 
[0] file: /home/ofuser/workingDir/FoamCases/KabelInLuft/processor0/system/air/fvSolution from line 22 to line 78.
[0] 
[0]     From function const Foam::dictionary& Foam::dictionary::subDict(const Foam::word&) const
[0]     in file db/dictionary/dictionary.C at line 640.
[0] 
FOAM parallel run exiting
[0] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[boot2docker:00248] 2 more processes have sent help message help-mpi-api.txt / mpi-abort
[boot2docker:00248] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

I see here, that the keyword "PIMPLE" is undefined in fvSolutions in the processor-folders --> system.
But there is no system-folder in one of the processor-folders, so there is no fvSolutions. So I looked after the ordinary fvSolution to search for PIMPLE... but there is only SIMPLE. Here you can see my fvSolution:

Code:

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

solvers
{
    rho
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-7;
        relTol          0;
    }

    p_rgh
    {
        solver           GAMG;
        tolerance        1e-7;
        relTol           0.01;

        smoother         GaussSeidel;

        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator     faceAreaPair;
        mergeLevels      1;

        maxIter          100;
    }

    "(U|h|k|epsilon)"
    {
        solver           PBiCG;
        preconditioner   DILU;
        tolerance        1e-7;
        relTol           0.1;
    }
}

SIMPLE
{
    momentumPredictor on;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       100000;
    rhoMin          rhoMin [1 -3 0 0 0] 0.2;
    rhoMax          rhoMax [1 -3 0 0 0] 2;
}

relaxationFactors
{
    fields
    {
        rho             1;
        p_rgh           0.7;
    }
    equations
    {
        U               0.3;
        h               0.7;
        nuTilda         0.7;
        k               0.7;
        epsilon         0.7;
        omega           0.7;
        "ILambda.*"     0.7;
    }
}


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

What can I do now? What could this error mean?
I try my best, but I don't have much time for my project anymore, so I need your help!

Thanks in advance
Raffael

jmdf · June 6, 2016, 04:36

Since the fvSolution file is the same for all the processors, the solvers reads it from the system folder.

In this fvSolution file you need to have the PIMPLE keywords with the inputs you need/want. Did you construct the fvSolution yourself or did you copied it from another location? It appears that you copied it from a tutorial of a steady-state solver.

As you are using chtMultiRegionFoam look into the tutorials for this solver. For example: $FOAM_TUTORIALS/heatTransfer/chtMultiRegionFoam/multiRegionHeater
Check the fvSolution of this tutorial and update your own with the required keywords.

ruffy19 · June 7, 2016, 07:39

Hello João!

At first thank you for your very quick reply! Your advice took me further and I have the feeling, that my problems with OpenFOAM are almost solved.
But this time this Error appeared:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : Jun 07 2016
Time   : 11:33:59
Host   : "boot2docker"
PID    : 281
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.282"
"boot2docker.283"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    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 laminar
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding MRF

No MRF models present

    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region cable

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Selecting radiationModel none
    Adding fvOptions

No finite volume options present

Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.801 max: 1.56
Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.801 max: 1.56
Time = 300


Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
[2] 
[2] 
[2] --> FOAM FATAL ERROR: 
[2] cannot be called for a calculatedFvPatchField
    on patch Nature of field U in file "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor2/0/air/U"
    You are probably trying to solve for a field with a default boundary condition.
[2] 
[2]     From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientInternalCoeffs() const [with Type = Foam::Vector<double>]
[2]     in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.
[2] 
FOAM parallel run exiting
[2] 
[0] 
[0] 
[0] --> FOAM FATAL ERROR: 
[0] cannot be called for a calculatedFvPatchField
    on patch Nature of field U in file "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor0/0/air/U"
    You are probably trying to solve for a field with a default boundary condition.
[0] 
[0]     From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientInternalCoeffs() const [with Type = Foam::Vector<double>]
[0]     in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.
[0] 
FOAM parallel run exiting
[0] 
[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] cannot be called for a calculatedFvPatchField
    on patch Nature of field U in file "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor1/0/air/U"
    You are probably trying to solve for a field with a default boundary condition.
[1] 
[1]     From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientInternalCoeffs() const [with Type = Foam::Vector<double>]
[1]     in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.
[1] 
FOAM parallel run exiting
[1] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[boot2docker:00279] 2 more processes have sent help message help-mpi-api.txt / mpi-abort
[boot2docker:00279] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

I tried to find out that the problem actually is, but without success. Can you help me out? I know, that you know what to do at this point. Please help me.

Thanks in advance
Raffael

jmdf · June 7, 2016, 09:01

No problem

Let's see, from what i understand you are using a boundary condition named "calculated" in the U field (0/<region>/U). And the error says that you cannot use such type.

I'm not into your project so i can't give you a more detailed advice on what boundary type you should use. However, fixedValue or zeroGradient is normally used on the U field.

ruffy19 · June 8, 2016, 09:19

Hello João and all the other Foamers out there,

again thanks for your help! I changed the boundary conditions for "U" for both regions and by the way I realized, that "k" and "epsilon" are both zero when "U" is zero (I have no moving fluid in my system). So I changed them too.

But like always it doesn't work... This time the message looks like this:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : Jun 08 2016
Time   : 13:07:44
Host   : "boot2docker"
PID    : 320
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.321"
"boot2docker.322"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    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 laminar
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding MRF

No MRF models present

    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region cable

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Selecting radiationModel none
    Adding fvOptions

No finite volume options present

Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.801 max: 1.56
Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.801 max: 1.56
Time = 300


Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
[1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] cannot be called for a calculatedFvPatchField
    on patch Nature of field h in file "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor1/0/air/h"
    You are probably trying to solve for a field with a default boundary condition.
[1] 
[1]     From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientInternalCoeffs() const [with Type = double]
[1]     in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.
[1] 
FOAM parallel run exiting
[1] 
[2] 
[2] 
[2] --> FOAM FATAL ERROR: 
[2] cannot be called for a calculatedFvPatchField
    on patch Nature of field h in file "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor2/0/air/h"
    You are probably trying to solve for a field with a default boundary condition.
[2] 
[2]     From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientInternalCoeffs() const [with Type = double]
[2]     in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.
[2] 
FOAM parallel run exiting
[2] 
[0] 
[0] 
[0] --> FOAM FATAL ERROR: 
[0] cannot be called for a calculatedFvPatchField
    on patch Nature of field h in file "/home/ofuser/workingDir/FoamCases/KabelInLuft/processor0/0/air/h"
    You are probably trying to solve for a field with a default boundary condition.
[0] 
[0]     From function Foam::tmp<Foam::Field<Type> > Foam::calculatedFvPatchField<Type>::gradientInternalCoeffs() const [with Type = double]
[0]     in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 187.
[0] 
FOAM parallel run exiting
[0] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 2 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[boot2docker:00318] 2 more processes have sent help message help-mpi-api.txt / mpi-abort
[boot2docker:00318] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

The message says the same like before for "U" but this time for "h". The thing is... I have no field h - nowhere.

Where does it come from? What can I do now?

I'm sorry for annoying with my stupid questions. Every little problem is a huge unexplainable obstacle to me... I'm really a noob in OpenFOAM.

Thanks in advance
Raffael

jmdf · June 9, 2016, 03:13

The error is the same so you already know how to solve it, if needed

Regarding the h field, that's weird.. Do you have any file named "h" in your 0/ or 0/<region>/ folder?
The solver does not require such field although i recall having some problems for having "not needed" fields in the 0/ folder. Anyway, can you check that?

ruffy19 · June 9, 2016, 05:20

Hello João,

I know that the error-message is the same like before, but the problem is, that I have no h-file, neither in the 0-folder nor in the 0/<region>/ folder. I also checked for the processor-folders and there isn't h as well.
I just have these:

cellToRegion
epsilon
k
p
p_rgh
T
U

Why does OpenFOAM ask for "h"? Like you said... the solver doesn't need h.
Could it be, that the problem is in one of the other files? Maybe in p_rgh?
By the way, can I ask what the difference between p and p_rgh is? p_rgh is the only file with an h in it, so maybe this is meant.

Thanks in advance
Raffael

Antimony · June 9, 2016, 05:50

Hi,

"h" here should refer to enthalpy(?), basically it is associated with the temperature field, if my understanding is correct. Check your T boundary conditions...

For your other question: http://www.cfd-online.com/Forums/ope...tml#post276743

Cheers,
Antimony

ruffy19 · June 9, 2016, 09:27

Hello Antimony and João,

Thank you Antimony for the helpful advice. I checked my T-files and after I used chtMultiRegionFoam OpenFOAM complained about p_rgh with the same issue. So I checked all the files again and now a new error appeared. I looked through other threads, but there they didn't come to a solution to my problem now.
The error-message looks like this:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : Jun 09 2016
Time   : 12:25:12
Host   : "boot2docker"
PID    : 295
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.296"
"boot2docker.297"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    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 laminar
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding MRF

No MRF models present

    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region cable

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Selecting radiationModel none
    Adding fvOptions

No finite volume options present

Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.801 max: 1.56
Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.801 max: 1.56
Time = 300


Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 5.33e-08, No Iterations 12
Min/max T:265 356
GAMG:  Solving for p_rgh, Initial residual = 0.949, Final residual = 0.00943, No Iterations 10
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = 12.3, global = 1, cumulative = 1
GAMG:  Solving for p_rgh, Initial residual = 0.0532, Final residual = 9.32e-08, No Iterations 56
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = 3.68, global = -0.278, cumulative = 0.722

Solving for solid region cable
DICPCG:  Solving for h, Initial residual = 1, Final residual = 7.55e-07, No Iterations 18
Min/max T:332 352
ExecutionTime = 0.28 s  ClockTime = 0 s

Region: air Courant Number mean: 9.73e+06 max: 1.01e+11
Region: cable Diffusion Number mean: 0.801 max: 1.56
Time = 600


Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 3.03e-08, No Iterations 54
Min/max T:-2.85e+10 3.32e+09
GAMG:  Solving for p_rgh, Initial residual = 0.963, Final residual = nan, No Iterations 1000
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (air): sum local = nan, global = -nan, cumulative = -nan
GAMG:  Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 1000
[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: 
[0] wrong token type - expected Scalar, found on line 0 the word 'nan'
[0] 
[0] file: /home/ofuser/workingDir/FoamCases/KabelInLuft/processor0/system/air/data.solverPerformance.p_rgh at line 0.
[0] 
[0]     From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)
[0]     in file lnInclude/Scalar.C at line 93.[1] 
[1] 
[1] --> FOAM FATAL IO ERROR: 
[1] wrong token type - expected Scalar, found on line 0 the word 'nan'
[1] 
[1] file: /home/ofuser/workingDir/FoamCases/KabelInLuft/processor1/system/air/data.solverPerformance.p_rgh at line 0.
[1] 
[1]     From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)
[1]     in file lnInclude/Scalar.C at line 93.
[1] 
FOAM parallel run exiting
[1] 
[2] 
[2] 
[2] --> FOAM FATAL IO ERROR: 
[2] wrong token type - expected Scalar, found on line 0 the word 'nan'
[2] 
[2] file: /home/ofuser/workingDir/FoamCases/KabelInLuft/processor2/system/air/data.solverPerformance.p_rgh at line 0.
[2] 
[2]     From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)
[2]     in file lnInclude/Scalar.C at line 93.
[2] 
FOAM parallel run exiting
[2] 

[0] 
FOAM parallel run exiting
[0] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[boot2docker:00293] 1 more process has sent help message help-mpi-api.txt / mpi-abort
[boot2docker:00293] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

As far as this I never came before. My presumption is, that I configured something wrong in the system-folders. I attached my system-folders, so you know better, what is going on. If you need more information, I'll attach the rest too.

Please check with me, where the mistake can be.

Thanks in advance
Raffael

jmdf · June 10, 2016, 02:36

First question, are you using a time-step of 300?
If so, you need to lower that value, otherwise you won't be able to obtain any results. Anyway, if changing that the same problem appears:

Code:

Min/max T:-2.85e+10 3.32e+09

The temperature values are blowing up. To solve this or to try to stabilise the simulation, you need to play with relaxation factor and tolerances of your linear solvers (all of these are in the fvSolution file).
You can also try another schemes (fvSchemes). Are you using first or second order schemes?

The mesh could be also the reason of your problem/divergence. What you can do is write the results on the previous time-step and check where are located the cells with such T values. Is it clear?

ruffy19 · June 10, 2016, 09:39

Hello João,

thank you for your help! But this time I realized, what a noob I am.

To your first questeion: Yes, I wanted to create a kind of fast motion. But I followed your advice and took the tenth part of it - so it's 30 now.

What do you mean by your code? I tried to use it after the first error but OpenFOAM didn't know what to do with that. Can you tell me more about the code and the purpose behind that?

Your advice to play with the relaxation factors and tolerances I followed as good as I can. The matter is... I don't know how the change of a relaxation factor affects on the simulation or the output.

What I also didn't understand was, what do you mean by first or second order scheme? Another scheme? I really don't know anything about the fvscheme-file... I just copied it from the tutorial "multiRegionHeater".

Quote:

What you can do is write the results on the previous time-step and check where are located the cells with such T values.

How can I do that? Do you mean by creating an own time-folder for 30s?

Quote:

Is it clear?

Unfortunately no...

My mesh could be a problem? I put my system-folder and my blockMeshDict in the attachment, so you know better what this is all about.

But not to miss, my current Error-message:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : chtMultiRegionFoam -parallel
Date   : Jun 10 2016
Time   : 13:14:04
Host   : "boot2docker"
PID    : 479
Case   : /home/ofuser/workingDir/FoamCases/KabelInLuft
nProcs : 3
Slaves : 
2
(
"boot2docker.480"
"boot2docker.481"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create fluid mesh for region air for time = 0

Create solid mesh for region cable for time = 0

*** Reading fluid mesh thermophysical properties for region air

    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 laminar
Selecting radiationModel none
    Adding to KFluid

    Adding to dpdtFluid

    Adding MRF

No MRF models present

    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region cable

    Adding to thermos

Selecting thermodynamics package 
{
    type            heSolidThermo;
    mixture         pureMixture;
    transport       constIso;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleEnthalpy;
}

    Adding to radiations

Selecting radiationModel none
    Adding fvOptions

No finite volume options present

Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.0801 max: 0.156
Region: air Courant Number mean: 0 max: 0
Region: cable Diffusion Number mean: 0.0801 max: 0.156
Time = 30


Solving for fluid region air
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG:  Solving for h, Initial residual = 1, Final residual = 9.65e-06, No Iterations 16
Min/max T:258 316
GAMG:  Solving for p_rgh, Initial residual = 0.972, Final residual = 0.0936, No Iterations 3
[0] 
[0] 
[0] --> FOAM FATAL ERROR: 
[0] temporary of type N4Foam8fvMatrixIdEE deallocated
[0] 
[0]     From function const T& Foam::tmp<T>::operator()() const [with T = Foam::fvMatrix<double>]
[0]     in file /opt/OpenFOAM/OpenFOAM-v3.0+/src/OpenFOAM/lnInclude/tmpI.H at line [1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] temporary of type N4Foam8fvMatrixIdEE deallocated
[1] 
[1]     From function const T& Foam::tmp<T>::operator()() const [with T = Foam::fvMatrix<double>]
[1]     in file /opt/OpenFOAM/OpenFOAM-v3.0+/src/OpenFOAM/lnInclude/tmpI.H at line 220.
[1] [2] 
[2] 
[2] --> FOAM FATAL ERROR: 
[2] temporary of type N4Foam8fvMatrixIdEE deallocated
[2] 
[2]     From function const T& Foam::tmp<T>::operator()() const [with T = Foam::fvMatrix<double>]
[2]     in file /opt/OpenFOAM/OpenFOAM-v3.0+/src/OpenFOAM/lnInclude/tmpI.H at line 220.
[0] 
FOAM parallel run aborting
[0] 

FOAM parallel run aborting
[1] 
220.
[2] 
FOAM parallel run aborting
[2] 
[2] #0  [0] #0  Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&)[1] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[1] #1  Foam::error::abort()[2] #1  Foam::error::abort() in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[0] #1  Foam::error::abort() in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[2] #2   in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[1] #2   in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/lib/libOpenFOAM.so"
[0] #2  ??? at chtMultiRegionFoam.C:0
[0] #3   at chtMultiRegionFoam.C:0
[1] #3   at chtMultiRegionFoam.C:0
[2] #3  ??? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[0] #4   in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam""

[2] #4  [1] #4  ??? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[0] #5  __libc_start_main in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[2] #5  __libc_start_main in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[1] #5  __libc_start_main in "/lib64/libc.so.6"
[0] #6   in "/lib64/libc.so.6"
[2] #6   in "/lib64/libc.so.6"
[1] #6  ??? in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD 
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
 in "/opt/OpenFOAM/OpenFOAM-v3.0+/platforms/linux64Gcc48DPInt32Opt/bin/chtMultiRegionFoam"
[boot2docker:00477] 2 more processes have sent help message help-mpi-api.txt / mpi-abort
[boot2docker:00477] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

This time I'm really not able to read something clear out of this. It is not possible to locate the mistake anymore. But I know you know what to do.

Thanks in advance
Raffael

jmdf · June 13, 2016, 08:29

I advise you to study how simulations work which will help you a lot on how to stabilise your case. Check this link regarding the first/second order schemes. The mesh could obviously be a problem, it affects the way that the equations are solved. Again, you should look for more information regarding this subject.

To be honest I have no idea about the new error you have. Don't know if any user of the forum have experience with this.

To write the results, you need to adjust the writeInterval variable (system/controlDict) to match with the time-step when the error occurs. Anyway the time-step still looks too high for me, however I don't know if it's part of your setup.
Do you have any information about the error you have at the moment?

June 7, 2016, 09:01		#12
jmdf Member João Ferreira Join Date: Nov 2014 Location: Braga, Portugal Posts: 53 Rep Power: 11	No problem Let's see, from what i understand you are using a boundary condition named "calculated" in the U field (0/<region>/U). And the error says that you cannot use such type. I'm not into your project so i can't give you a more detailed advice on what boundary type you should use. However, fixedValue or zeroGradient is normally used on the U field. ruffy19 likes this.

June 13, 2016, 08:29		#20
jmdf Member João Ferreira Join Date: Nov 2014 Location: Braga, Portugal Posts: 53 Rep Power: 11	I advise you to study how simulations work which will help you a lot on how to stabilise your case. Check this link regarding the first/second order schemes. The mesh could obviously be a problem, it affects the way that the equations are solved. Again, you should look for more information regarding this subject. To be honest I have no idea about the new error you have. Don't know if any user of the forum have experience with this. To write the results, you need to adjust the writeInterval variable (system/controlDict) to match with the time-step when the error occurs. Anyway the time-step still looks too high for me, however I don't know if it's part of your setup. Do you have any information about the error you have at the moment? ruffy19 likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Error in thermophysical properties (chtMultiRegionFoam)	mukut	OpenFOAM Pre-Processing	28	November 23, 2021 06:34
Error in chtMultiRegionFoam	kirankarki	OpenFOAM	6	August 21, 2018 08:00
What does this error when use solver chtMultiRegionFoam mean	sajad6	OpenFOAM Running, Solving & CFD	7	October 6, 2014 07:38
chtmultiregionFoam error	oilsok	OpenFOAM Running, Solving & CFD	1	June 12, 2014 11:19
chtmultiregionFoam	alvora	OpenFOAM	9	February 23, 2011 03:06

May 23, 2016, 05:19	Error appears when using chtMultiRegionFoam	#1
ruffy19 New Member Raffael H. Join Date: May 2016 Location: Ludwigshafen Posts: 11 Rep Power: 9	Hello Foam-Experts! I'm really new to OpenFoam and the programming in it and so I'm not experienced with it. But for my studies I got the following task: I have to simulate a solid cylinder hanging in air. The cylinder radiates heat and convection is on hand. So I use "chtMultiRegionFoam" with the command Code: mpirun -np 3 chtMultiRegionFoam -parallel to run the solution with its 2 declared regions and 3 generated processors. I used the same way to solve the case for a functional project called "planWall2D" for which I assume that it is programmed correctly, because I got this form the openfoamwiki. Then I realized, that the error-messages for both cases are the same. This error-message appears after the thermophysical properties for region air were read and added and ends up like: Code: mpirun noticed that process rank 2 with PID 477 on node boot2docker exited on signal 8 <Floating point exception> . Because an error-message like that appears for my case as well as for the model-programm, I guess the error is not my fault. So I would like to ask you all, what the problem actually is. By the way I use the windows-version of OpenFoam. I really please you to help me and also if you need some more information about my problem, please answer me as soon as possible. Thanks in advance Raffael

May 25, 2016, 03:35		#2
ruffy19 New Member Raffael H. Join Date: May 2016 Location: Ludwigshafen Posts: 11 Rep Power: 9	I really need to be helped with this problem. I don't know what the error message means and what to do about that. It is really important for me. I'm not experienced with OpenFOAM and maybe I did something wrong. But until now I found nothing... I really please you to help me! Thanks in advance Raffael

May 25, 2016, 04:09		#3
Antimony Senior Member Join Date: Aug 2013 Posts: 407 Rep Power: 15	Hi, The error message that you have is a Floating Point Exception. The most likely culprit for this is a division by zero somewhere. Look through the different boundary condition files and see if there is something that shouldn't be zero is zero. What does the log of your run say apart from the mpirun run abort? Cheers, Antimony

May 25, 2016, 21:38		#5
Antimony Senior Member Join Date: Aug 2013 Posts: 407 Rep Power: 15	Hi, You can get the entire output to a log file by adding "2>&1" at the end of the function call. Since it seems to stop after loading your thermophysical properties, is your temperature 0 at any point? Do take care that the units of temperature are in K... Cheers, Antimony ruffy19 likes this.

June 4, 2016, 18:52		#8
jmdf Member João Ferreira Join Date: Nov 2014 Location: Braga, Portugal Posts: 53 Rep Power: 11	From the error it appears that you are using a "calculated" boundary condition in a boundary with the type empty. If you need to use the calculated boundary, check if that boundary is defined as patch or wall on the "boundary" file (constant/<region>/polyMesh/boundary). I hope it helps you! dokeun and ruffy19 like this.

June 6, 2016, 04:36		#10
jmdf Member João Ferreira Join Date: Nov 2014 Location: Braga, Portugal Posts: 53 Rep Power: 11	Since the fvSolution file is the same for all the processors, the solvers reads it from the system folder. In this fvSolution file you need to have the PIMPLE keywords with the inputs you need/want. Did you construct the fvSolution yourself or did you copied it from another location? It appears that you copied it from a tutorial of a steady-state solver. As you are using chtMultiRegionFoam look into the tutorials for this solver. For example: $FOAM_TUTORIALS/heatTransfer/chtMultiRegionFoam/multiRegionHeater Check the fvSolution of this tutorial and update your own with the required keywords.

June 9, 2016, 03:13		#14
jmdf Member João Ferreira Join Date: Nov 2014 Location: Braga, Portugal Posts: 53 Rep Power: 11	The error is the same so you already know how to solve it, if needed Regarding the h field, that's weird.. Do you have any file named "h" in your 0/ or 0/<region>/ folder? The solver does not require such field although i recall having some problems for having "not needed" fields in the 0/ folder. Anyway, can you check that?

June 9, 2016, 05:20		#15
ruffy19 New Member Raffael H. Join Date: May 2016 Location: Ludwigshafen Posts: 11 Rep Power: 9	Hello João, I know that the error-message is the same like before, but the problem is, that I have no h-file, neither in the 0-folder nor in the 0/<region>/ folder. I also checked for the processor-folders and there isn't h as well. I just have these: cellToRegion epsilon k p p_rgh T U Why does OpenFOAM ask for "h"? Like you said... the solver doesn't need h. Could it be, that the problem is in one of the other files? Maybe in p_rgh? By the way, can I ask what the difference between p and p_rgh is? p_rgh is the only file with an h in it, so maybe this is meant. Thanks in advance Raffael

June 9, 2016, 05:50		#16
Antimony Senior Member Join Date: Aug 2013 Posts: 407 Rep Power: 15	Hi, "h" here should refer to enthalpy(?), basically it is associated with the temperature field, if my understanding is correct. Check your T boundary conditions... For your other question: http://www.cfd-online.com/Forums/ope...tml#post276743 Cheers, Antimony

June 10, 2016, 02:36		#18
jmdf Member João Ferreira Join Date: Nov 2014 Location: Braga, Portugal Posts: 53 Rep Power: 11	First question, are you using a time-step of 300? If so, you need to lower that value, otherwise you won't be able to obtain any results. Anyway, if changing that the same problem appears: Code: Min/max T:-2.85e+10 3.32e+09 The temperature values are blowing up. To solve this or to try to stabilise the simulation, you need to play with relaxation factor and tolerances of your linear solvers (all of these are in the fvSolution file). You can also try another schemes (fvSchemes). Are you using first or second order schemes? The mesh could be also the reason of your problem/divergence. What you can do is write the results on the previous time-step and check where are located the cells with such T values. Is it clear?