CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (https://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Error in thermophysical properties (chtMultiRegionFoam) (https://www.cfd-online.com/Forums/openfoam-pre-processing/124938-error-thermophysical-properties-chtmultiregionfoam.html)

mukut October 15, 2013 23:35

Error in thermophysical properties (chtMultiRegionFoam)
 
Hello foamers!

I have modified multiRegionHeater Tutorial of chtMultiRegionFoam to following geometry, as I want to simulate plasma actuator induced flow:
http://i40.tinypic.com/2j61uag.jpg

In tutorials, following regions were created by toposetDict:

topAir
bottomAir
heater
leftSolid
rightSolid

Then I have modified the geometry of this tutorial like as above image in which I have defined following regions by topoSetDict

topAir (same as tutorial but differ in dimension)
bottomAir (same as tutorial but differ in dimension)
heater (as a dielectric material which is kapton film)
leftSolid (as a top electrode which is copper)
rightSolid (as a bottom electrode which is copper)
innerelec (as a inner electrode which is also copper)

I have set themophysical properties for each region according to their material properties except topAir and bottomAir as they are same air regions in my case like tutorials.

I have run following commands in the terminal:

Code:

blockMesh
topoSet
splitMeshRegions -cellZones -overwrite
decomposePar -allRegions
mpirun -np 4 chtMultiRegionFoam -parallel

But I found following error after executing
mpirun -np 4 chtMultiRegionFoam -parallel

Code:

mukut@mukut-Endeavor-MR3300:~/OpenFOAM/mukut-2.2.1/run/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater$  mpirun -np 4 chtMultiRegionFoam -parallel
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.2.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/

Build  : 2.2.1-57f3c3617a2d
Exec  : chtMultiRegionFoam -parallel
Date  : Oct 12 2013
Time  : 09:50:20
Host  : "mukut-Endeavor-MR3300"
PID    : 4482
Case  : /home/mukut/OpenFOAM/mukut-2.2.1/run/tutorials/heatTransfer/chtMultiRegionFoam/multiRegionHeater
nProcs : 4
Slaves :
3
(
"mukut-Endeavor-MR3300.4483"
"mukut-Endeavor-MR3300.4484"
"mukut-Endeavor-MR3300.4485"
)

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 : Disallowing user-supplied system call operations


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

Create time

Create fluid mesh for region bottomAir for time = 0

Create fluid mesh for region topAir for time = 0

Create solid mesh for region heater for time = 0

Create solid mesh for region leftSolid for time = 0

Create solid mesh for region rightSolid for time = 0

Create solid mesh for region innerelec for time = 0

*** Reading fluid mesh thermophysical properties for region bottomAir

    Adding to thermoFluid

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

    Adding to rhoFluid
[2] #0 
    Adding to UFluid

    Adding to phiFluid

Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[2] #2  in "/lib/x86_64-linux-gnu/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/openfoam221/platforms/linux64GccDPOpt/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/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
[2]  #5  Foam::autoPtr<Foam::rhoThermo>  Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&,  Foam::word const&) in  "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
[2]  #6  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word  const&) in  "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
[2] #7 
[2]  in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"
[2] #8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[2] #9 
[2]  in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"
[mukut-Endeavor-MR3300:04484] *** Process received signal ***
[mukut-Endeavor-MR3300:04484] Signal: Floating point exception (8)
[mukut-Endeavor-MR3300:04484] Signal code:  (-6)
[mukut-Endeavor-MR3300:04484] Failing at address: 0x3e800001184
[mukut-Endeavor-MR3300:04484] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f05aec694a0]
[mukut-Endeavor-MR3300:04484] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7f05aec69425]
[mukut-Endeavor-MR3300:04484] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7f05aec694a0]
[mukut-Endeavor-MR3300:04484]  [ 3]  /opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11heRhoThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEE9calculateEv+0x6cd)  [0x7f05b3ccbdad]
[mukut-Endeavor-MR3300:04484] [ 4]  /opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo27addfvMeshConstructorToTableINS_11heRhoThermoIS0_NS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_10perfectGasINS_6specieEEEEENS_16sensibleEnthalpyEEEEEEEEEE3NewERKNS_6fvMeshERKNS_4wordE+0x5c)  [0x7f05b3ceba6c]
[mukut-Endeavor-MR3300:04484] [ 5]  /opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11basicThermo3NewINS_9rhoThermoEEENS_7autoPtrIT_EERKNS_6fvMeshERKNS_4wordE+0x11b)  [0x7f05b3cabccb]
[mukut-Endeavor-MR3300:04484] [ 6]  /opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo3NewERKNS_6fvMeshERKNS_4wordE+0x9)  [0x7f05b3caa1d9]
[mukut-Endeavor-MR3300:04484] [ 7] chtMultiRegionFoam() [0x423ff8]
[mukut-Endeavor-MR3300:04484] [ 8] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7f05aec5476d]
[mukut-Endeavor-MR3300:04484] [ 9] chtMultiRegionFoam() [0x42c5ed]
[mukut-Endeavor-MR3300:04484] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 4484 on node mukut-Endeavor-MR3300 exited on signal 8 (Floating point exception).


How can I get rid of it?

Thanks in advance.

Best regards,
Mukut

otq October 21, 2013 15:18

Did you make your own library?
libfluidThermophysicalModels.so?
The standard should just be basicThermophysicalModels.
If you did make your own library, you have to make sure that everything is connected and pointed to properly.

The error is referring to the calculate function in the sensible enthalpy section so it may also be a typo there.

If all you did was change the geometry by adding more solids, you shouldn't be getting any of those kinds of errors since they relate to the solver itself not the problem.
Granted that is from 2.1.1, but I doubt they would add fluid just like that.

mukut October 21, 2013 20:26

Quote:

Originally Posted by otq (Post 458146)
Did you make your own library?
libfluidThermophysicalModels.so?
The standard should just be basicThermophysicalModels.
If you did make your own library, you have to make sure that everything is connected and pointed to properly.

The error is referring to the calculate function in the sensible enthalpy section so it may also be a typo there.

If all you did was change the geometry by adding more solids, you shouldn't be getting any of those kinds of errors since they relate to the solver itself not the problem.
Granted that is from 2.1.1, but I doubt they would add fluid just like that.

Thank you so much. the problem was solved.. :)

fabian_roesler January 30, 2014 04:20

foam::error::printStack compressibleInterFoam parallel run
 
Hi mukut

How did you solve the problem with the foam::error::printStack? I am facing the same problem like you did in OpenFOAM 2.2.2. I perform simulations with compressibleInterFoam on a case with round about 2e6 cells. In serial the simulation runs perfect but would take weeks to finish. In parallel run I get the following error massage (I cleaned the full error output to show only errors from one calculation node):

Code:

/*---------------------------------------------------------------------------*\
 | =========                |                                                |
 | \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
 |  \\    /  O peration    | Version:  2.2.2                                |
 |  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
 |    \\/    M anipulation  |                                                |
 \*---------------------------------------------------------------------------*/
 Build  : 2.2.2-9739c53ec43f
 Exec  : compressibleInterFoam -parallel
 Date  : Jan 30 2014
 Time  : 08:40:54
 Host  : "nod00"
 PID    : 31247
 Case  : /home/me/OpenFOAM/me-2.2.2/run/test
 nProcs : 16
 Slaves : 
 15
 (
 "nod01.31248"
 "nod02.31249"
 "nod03.31250"
 "nod04.31251"
 "nod05.29274"
 "nod06.29275"
 "nod07.29276"
 "nod08.29277"
 "nod09.29278"
 "nod10.11981"
 "nod11.11982"
 "nod12.11983"
 "nod13.11984"
 "nod14.11985"
 "nod15.11986"
 )
 
 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 : Disallowing user-supplied system call operations
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 Create time
 
 Create mesh for time = 0
 
 
 Reading g
 
 PIMPLE: Operating solver in PISO mode
 
 Reading field p_rgh
 
 Reading field U
 
 Reading/calculating face flux field phi
 
 Constructing twoPhaseMixtureThermo
 
 Selecting thermodynamics package 
 {
    type            heRhoThermo;
    mixture        pureMixture;
    transport      const;
    thermo          hConst;
    equationOfState perfectFluid;
    specie          specie;
    energy          sensibleInternalEnergy;
 }
 
 [3] #0  Foam::error::printStack(Foam::Ostream&)
 [3] #1  Foam::sigFpe::sigHandler(int) in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
 [3] #2  in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
 [3] #3  Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::he(Foam::Field<double> const&, Foam::Field<double> const&, int) const[5]  at sigaction.c:0
 [3] #4  Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::init() in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so"
 [3] #5  Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so"
 [3] #6  Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectFluid<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so"
 [3] #7  Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so"
 [3] #8  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so"
 [3] #9  Foam::twoPhaseMixtureThermo::twoPhaseMixtureThermo(Foam::fvMesh const&) in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so"
 [3] #10  in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libtwoPhaseMixtureThermo.so"
 [3] #11  __libc_start_main in "/home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libtwoPhaseMixtureThermo.so"
 [3] #12  in "/lib64/libc.so.6"
 
 [nod03:31250] *** Process received signal ***
 [nod03:31250] Signal: Floating point exception (8)
 [nod03:31250] Signal code:  (-6)
 [nod03:31250] Failing at address: 0x2271800007a12
 [nod03:31250] [ 0] /lib64/libc.so.6() [0x3f64832920]
 [nod03:31250] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x3f648328a5]
 [nod03:31250] [ 2] /lib64/libc.so.6() [0x3f64832920]
 [nod03:31250] [ 3] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so(_ZNK4Foam8heThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_12perfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE2heERKNS_5FieldIdEESJ_i+0x9e) [0x7f35e50fcbae]
 [nod03:31250] [ 4] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam8heThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_12perfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEE4initEv+0x206) [0x7f35e5123c96]
 [nod03:31250] [ 5] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam8heThermoINS_9rhoThermoENS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_12perfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEEC2ERKNS_6fvMeshERKNS_4wordE+0x25a) [0x7f35e51306ba]
 [nod03:31250] [ 6] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo27addfvMeshConstructorToTableINS_11heRhoThermoIS0_NS_11pureMixtureINS_14constTransportINS_7species6thermoINS_12hConstThermoINS_12perfectFluidINS_6specieEEEEENS_22sensibleInternalEnergyEEEEEEEEEE3NewERKNS_6fvMeshERKNS_4wordE+0x3c) [0x7f35e51307cc]
 [nod03:31250] [ 7] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam11basicThermo3NewINS_9rhoThermoEEENS_7autoPtrIT_EERKNS_6fvMeshERKNS_4wordE+0x115) [0x7f35e50f2495]
 [nod03:31250] [ 8] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libfluidThermophysicalModels.so(_ZN4Foam9rhoThermo3NewERKNS_6fvMeshERKNS_4wordE+0x9) [0x7f35e50f0969]
 [nod03:31250] [ 9] /home/me/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64Gcc45DPOpt/lib/libtwoPhaseMixtureThermo.so(_ZN4Foam21twoPhaseMixtureThermoC2ERKNS_6fvMeshE+0x280) [0x7f35e53c64b0]
 [nod03:31250] [10] compressibleInterFoam() [0x42d908]
 [nod03:31250] [11] /lib64/libc.so.6(__libc_start_main+0xfd) [0x3f6481ecdd]
 [nod03:31250] [12] compressibleInterFoam() [0x425679]
 [nod03:31250] *** End of error message ***
 
 mpirun noticed that process rank 1 with PID 31248 on node nod03 exited on signal 8 (Floating point exception).
 --------------------------------------------------------------------------


The error occurs when OpenFOAM tries to initialize the thermophysical model of the air phase. I only changed the geometry from the tutorial case. Moreover, the error shows a problem of the thermophysical model, not the case itself. Thanks for your help.


Cheers


Fabian

fabian_roesler February 6, 2014 06:39

twoPhaseMixtureThermo bug?
 
Hi

I started from scratch to find the trigger of my error in compressibleInterFoam case. I took the tutorial case and ran it serial and parallel. Everything is fine. When I change the temperatures to the range I have to simulate my specific case (round about 900 C), then the serial case runs flawless. However the parallel run stops with printStack error after constructing the twoPhaseMixtureThermo. So the only problem is the increase in temperature? I also tried with janaf thermo model without success. Could this be a bug (OF 2.2.2)? Help is still welcome!

Cheers

Fabian

olivierG February 6, 2014 09:14

hello,

Did you check if after decomposePar, you did not get in the proc../0/T files some "value (non)uniform 0;" ?
Because some tools (merge / stitch patch, maybe other ...) seem to do that (i don"t know why).

regards,
olivier

fabian_roesler May 2, 2014 15:54

Hi

Sorry for this huge delay. I think you have a good point here. I used snappy to generate the mesh. To not have to reconstruct I copied the 0 time folder into the processor directories after meshing. Then I ran setFields. The important point know: I use

Code:

"$WM_PROJECT_DIR/etc/caseDicts/setConstraintTypes"
in my fields. So I do not set the processor boundaries myself. This could lead to the error. I will check this an confirm the solution here asap.

Cheers

Fabian

cliffdub August 23, 2014 16:40

Dear Mukut or anyone else who might be able to help,

I am trying to use chtMultiRegionFoam to model transient conduction between two solids (no fluid). I am having a very difficult time and I hope you don't mind that I trouble you with such a basic problem. Is there an possibility that you could have a look at my two-solid heat transfer case and see where I have gone wrong?

https://www.dropbox.com/sh/1fxb42g1b...92kKyZgha?dl=0

I currently do not have any fluids in the region properties. I suppose this could cause trouble for chtMultiRegion.

Thank you,
Cliff

otq August 23, 2014 16:51

Quote:

Originally Posted by cliffdub (Post 507465)
Dear Mukut or anyone else who might be able to help,

I am trying to use chtMultiRegionFoam to model transient conduction between two solids (no fluid). I am having a very difficult time and I hope you don't mind that I trouble you with such a basic problem. Is there an possibility that you could have a look at my two-solid heat transfer case and see where I have gone wrong?

https://www.dropbox.com/sh/1fxb42g1b...92kKyZgha?dl=0

I currently do not have any fluids in the region properties. I suppose this could cause trouble for chtMultiRegion.

Thank you,
Cliff

Cliff I would suggest changing the region properties file. You haven't changed it. You are telling it that there are fluid regions involved, but then not providing them.

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

fluidRegionNames ( );

solidRegionNames ( solid1 solid2);

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

is a good place to start. Then you have to make sure that you have the appropriate folders in your system folder for the regions, including the changeDictionaryDict files in their respective system/regionName folder

cliffdub August 23, 2014 18:13

Thank you very much. It's running now. Much joy here.

You wouldn't happen to know why in the controlDict I am not able to successfully run with

adjustTimeStep yes;

when I remove this, it runs, albeit slowly. The adjustable time step works well for the multiRegionHeater tutorial case, but not in my simple case.

Cliff

baran_foam August 26, 2014 12:41

1 Attachment(s)
Hello everyone,
I am trying to solve conjugate heat transfer problem in chtMultiRegionFoam with geometry of cavity with some insulation. That means two region...solid and air. But it is showing some problem when it read thermophysical property..please suggest something..

https://www.dropbox.com/s/kdt1prx29c...heck1.PNG?dl=0

-Baran

fabian_roesler August 27, 2014 02:35

nonuniform value 0
 
Hi

I usually get those errors when the initial fields and boundary conditions of pressure p and temperature T have nonphysical values. Sometimes, for example when you decompose a case for a parallel run, the fields get initialized with nonuniform values. So when you find a pressure or temperature boundary with value 0 this causes the error. Most obvious, the thermodynamic model can't handle a fluid with 0 K temperature and 0 Pa pressure. Best way to get around this is using changeDictionary and overwrite all fields before running.

Cheers

Fabian

baran_foam August 27, 2014 05:37

Quote:

Originally Posted by fabian_roesler (Post 507943)
Hi

I usually get those errors when the initial fields and boundary conditions of pressure p and temperature T have nonphysical values. Sometimes, for example when you decompose a case for a parallel run, the fields get initialized with nonuniform values. So when you find a pressure or temperature boundary with value 0 this causes the error. Most obvious, the thermodynamic model can't handle a fluid with 0 K temperature and 0 Pa pressure. Best way to get around this is using changeDictionary and overwrite all fields before running.

Cheers

Fabian

Very very thanks for your guidance...Its working for this now. But for my case solver is not reading "changeDictionary"... I manually change bc in 0 directory...but it is showing problem regarding boundary condition of T..
https://www.dropbox.com/s/s8faeevvgd...heck3.PNG?dl=0

After that I change bc manually in 0 directory to
Code:

            "cavityAir_to_.*"
            {
                type                         
                compressible::turbulentTemperatureCoupledBaffleMixed;
                neighbourFieldName T;
                K              basicThermo;
                KName          none;
                value          uniform 300;
            }

After that it is showing following error..
https://www.dropbox.com/s/3vs6c1h64y...heck2.PNG?dl=0

Maybe this is very basic thing..but...as I am new to OF...so struggling to understand this thing...please suggest ...

-baran

wyldckat August 30, 2014 08:50

Greetings to all!

@baran: I saw the PM you sent me and came looking for public posts you had. I'm glad the first issue is solved.

But the new problem is because you're using the boundary condition type "calculated" where it should not be used.
Keep in mind that the "calculated" type of boundary condition means that the result on that patch for that specific field, is a result of the flow on the rest of the domain and that it can be calculated based on the other fields or from the flow itself. In other words, the error message you're getting is implying that you cannot use the "calculated" type for your patch.

Beyond this, I suggest that you take a few steps back and test ideas using a simpler test case. For example, you can use the case given here: http://openfoamwiki.net/index.php/Ge..._-_planeWall2D - since it's a simpler case than yours, and since it's one where it's easier to calculate the analytical solution for comparison, you can then make small changes to make it closer and closer to your problem at hand. For example, this way you can test various thermodynamic properties with this test case, to assess what you can and cannot do.
In other words: isolate and conquer one small problem at a time ;)

Best regards,
Bruno

baran_foam August 31, 2014 23:55

Greeting everyone,

@bruno.....thanks for your guidance....I went through the tutorial of multiRegionHeater of chtMultiRegionFoam. According to this I modified my case as well as boundary conditions which was implemented in wrong way initially. But still there is some problem in turbulence model which is observed in log.chtMultiRegionFoam though I used same thing like in tutorial .

Quote:

Pstream initialized with:
floatTransfer : 0
nProcsSimpleSum : 0
commsType : nonBlocking
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create fluid mesh for region cavityAir for time = 0

Create solid mesh for region consT for time = 0

Create solid mesh for region front for time = 0

Create solid mesh for region back for time = 0

*** Reading fluid mesh thermophysical properties for region cavityAir

Adding to thermoFluid

Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Adding to rhoFluid

Adding to kappaFluid

Adding to UFluid

Adding to phiFluid

Adding to gFluid

Adding to turbulence

Selecting turbulence model type laminar
Adding to ghFluid

Adding to ghfFluid

Selecting radiationModel none
Adding to KFluid

Adding to dpdtFluid

*** Reading solid mesh thermophysical properties for region consT

Adding to thermos

Constructed constSolidThermo with
rho : rho [1 -3 0 0 0 0 0] 200
Cp : Cp [0 2 -2 -1 0 0 0] 1050
K : K [1 1 -3 -1 0 0 0] 0.2
Hf : Hf [0 2 -2 0 0 0 0] 0
emissivity : emissivity [0 0 0 0 0 0 0] 0
kappa : kappa [0 -1 0 0 0 0 0] 0
sigmaS : sigmaS [0 -1 0 0 0 0 0] 0

*** Reading solid mesh thermophysical properties for region front

Adding to thermos

Constructed constSolidThermo with
rho : rho [1 -3 0 0 0 0 0] 200
Cp : Cp [0 2 -2 -1 0 0 0] 1050
K : K [1 1 -3 -1 0 0 0] 0.2
Hf : Hf [0 2 -2 0 0 0 0] 0
emissivity : emissivity [0 0 0 0 0 0 0] 0
kappa : kappa [0 -1 0 0 0 0 0] 0
sigmaS : sigmaS [0 -1 0 0 0 0 0] 0

*** Reading solid mesh thermophysical properties for region back

Adding to thermos

Constructed constSolidThermo with
rho : rho [1 -3 0 0 0 0 0] 200
Cp : Cp [0 2 -2 -1 0 0 0] 1050
K : K [1 1 -3 -1 0 0 0] 0.2
Hf : Hf [0 2 -2 0 0 0 0] 0
emissivity : emissivity [0 0 0 0 0 0 0] 0
kappa : kappa [0 -1 0 0 0 0 0] 0
sigmaS : sigmaS [0 -1 0 0 0 0 0] 0

Region: cavityAir Courant Number mean: 9.5725798e-06 max: 1.0265911e-05
Region: consT Diffusion Number mean: 8.5923526e-11 max: 1.6744986e-10
Region: front Diffusion Number mean: 9.4968853e-11 max: 1.1397229e-10
Region: back Diffusion Number mean: 9.4968853e-11 max: 1.1397229e-10
deltaT = 10
Region: cavityAir Courant Number mean: 0.095725798 max: 0.10265911
Region: consT Diffusion Number mean: 8.5923526e-07 max: 1.6744986e-06
Region: front Diffusion Number mean: 9.4968853e-07 max: 1.1397229e-06
Region: back Diffusion Number mean: 9.4968853e-07 max: 1.1397229e-06
deltaT = 10
Time = 10


Solving for fluid region cavityAir
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 8.8229963e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.1851732e-11, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.9674932e-12, No Iterations 3
[1]
[1]
[1] --> FOAM FATAL ERROR:
[1]
request for turbulenceModel turbulenceModel from objectRegistry back failed
available objects of type turbulenceModel are

0
(
)
[1]
[1]
[1] From function objectRegistry::lookupObject<Type>(const word&) const
[1] in file /opt/apps/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 131.
[1]
FOAM parallel run aborting
I am unable to sort it out this issue.....!!! Please suggest something...

Regards,
Baran

fabian_roesler September 1, 2014 03:25

Quote:

Originally Posted by baran_foam (Post 507986)
Very very thanks for your guidance...Its working for this now. But for my case solver is not reading "changeDictionary"... I manually change bc in 0 directory...

Hi baren

probably you found out already that changeDictionary has to be called before starting the solution process with the solver. You have to prepare a changeDictionaryDict for all your regions. And then execute the change with the command

Code:

changeDictionary
in the shell. The error you get from the solver concerning you turbulence model could be due to a wrong thermoType. Check this thread:

http://www.cfd-online.com/Forums/ope...tml#post381177

Took me 1 minute to Google and read :(

Regards Fabian

baran_foam September 2, 2014 01:18

Hello everyone,
@ fabian , thanks...for your suggestion...i sort it out the problem......of thermophysical properties...

but there is some mpi running error is coming for my case...its not running..
Quote:

Constructed constSolidThermo with
rho : rho [1 -3 0 0 0 0 0] 200
Cp : Cp [0 2 -2 -1 0 0 0] 1050
K : K [1 1 -3 -1 0 0 0] 0.2
Hf : Hf [0 2 -2 0 0 0 0] 0
emissivity : emissivity [0 0 0 0 0 0 0] 0
kappa : kappa [0 -1 0 0 0 0 0] 0
sigmaS : sigmaS [0 -1 0 0 0 0 0] 0

Region: cavityAir Courant Number mean: 9.5725798e-06 max: 1.0265911e-05
Region: consT Diffusion Number mean: 8.5923526e-11 max: 1.6744986e-10
Region: front Diffusion Number mean: 9.4968853e-11 max: 1.1397229e-10
Region: back Diffusion Number mean: 9.4968853e-11 max: 1.1397229e-10
deltaT = 10
Region: cavityAir Courant Number mean: 0.095725798 max: 0.10265911
Region: consT Diffusion Number mean: 8.5923526e-07 max: 1.6744986e-06
Region: front Diffusion Number mean: 9.4968853e-07 max: 1.1397229e-06
Region: back Diffusion Number mean: 9.4968853e-07 max: 1.1397229e-06
deltaT = 10
Time = 10


Solving for fluid region cavityAir
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 8.8229963e-09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.1851732e-11, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.9674932e-12, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 2.7586618e-08, No Iterations 2
Min/max T:-15.871339 1000
[0] #0 Foam::error::printStack(Foam::Ostream&)--------------------------------------------------------------------------
An MPI process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your MPI job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

Local host: INDUX08.com (PID 19427)
MPI_COMM_WORLD rank: 0

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1 Foam::sigFpe::sigHandler(int) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2
[0] at sigaction.c:0
[0] #3 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::Field<double> const&, Foam::Field<double> const&, Foam::Field<double> const&) const in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4 Foam::GAMGSolver::scalingFactor(Foam::Field<double >&, Foam::lduMatrix const&, Foam::Field<double>&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #8
[0] in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"
[0] #9 __libc_start_main in "/lib64/libc.so.6"
[0] #10
[0] in "/opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/chtMultiRegionFoam"
[INDUX08:19427] *** Process received signal ***
[INDUX08:19427] Signal: Floating point exception (8)
[INDUX08:19427] Signal code: (-6)
[INDUX08:19427] Failing at address: 0x22400004be3
[INDUX08:19427] [ 0] /lib64/libc.so.6() [0x3b17a32900]
[INDUX08:19427] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x3b17a32885]
[INDUX08:19427] [ 2] /lib64/libc.so.6() [0x3b17a32900]
[INDUX08:19427] [ 3] /opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver13scalingFacto rERNS_5FieldIdEERKS2_S5_S5_+0x52) [0x7ff2d4855b32]
[INDUX08:19427] [ 4] /opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver13scalingFacto rERNS_5FieldIdEERKNS_9lduMatrixES3_RKNS_10FieldFie ldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEE RKS2_h+0xa6) [0x7ff2d4855cc6]
[INDUX08:19427] [ 5] /opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7 PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS 8_S9_S9_S9_RNS1_IS8_EESD_h+0x12da) [0x7ff2d485892a]
[INDUX08:19427] [ 6] /opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5Fi eldIdEERKS2_h+0x492) [0x7ff2d4859ac2]
[INDUX08:19427] [ 7] /opt/apps/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERKNS _10dictionaryE+0x151) [0x7ff2d5722b61]
[INDUX08:19427] [ 8] chtMultiRegionFoam() [0x444e3a]
[INDUX08:19427] [ 9] /lib64/libc.so.6(__libc_start_main+0xfd) [0x3b17a1ecdd]
[INDUX08:19427] [10] chtMultiRegionFoam() [0x4233b9]
[INDUX08:19427] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 19427 on node
INDUX08.com exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------


then i go through this...
http://webcache.googleusercontent.co...&ct=clnk&gl=in
and change my deltaT to 10^-8 ....but still it showing same error....please suggest something...

Regards,
Baran

baran_foam September 3, 2014 07:54

Hello everyone,
my problem is solved....I did the following mistakes...

1. I convert one gambit mesh file into openfoam ,so its scaling is forgot to change..
2. Adjustable run time is on...so deltaT is changed to some different value...which create the floating point exception issue...

Regards,
Baran

fkika October 28, 2015 02:40

Error in chtMultiRegionFOAM
 
Hi everyone!

I have this error during the execution of my simulation:

Solving for fluid region water
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 1.167835e-009, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 5.9374103e-008, No Iterations 4


--> FOAM FATAL ERROR:
Kappa defined to employ fluidThermo method, but thermo package not available

From function temperatureCoupledBase::kappa(const scalarField&) const
in file derivedFvPatchFields/temperatureCoupledBase/temperatureCoupledBase.C at line 116.

FOAM exiting

The solver is ChtmultiRegionFOAM.

Anyone can help me?

Thank you so much!

wyldckat October 31, 2015 09:18

Quick questions @fkika: Please provide more details, such as:
  1. Which OpenFOAM version are you using?
  2. Which installation instructions did you follow for installing OpenFOAM?
  3. Can you reproduce this error with one of OpenFOAM's tutorials?
  4. How is the content of the files "thermophysicalProperties" in your case? If you don't know where they are, run:
    Code:

    find . -name thermophysicalProperties


All times are GMT -4. The time now is 14:26.