CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Error in thermophysical properties (chtMultiRegionFoam)

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree11Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 16, 2013, 00:35
Post Error in thermophysical properties (chtMultiRegionFoam)
  #1
Senior Member
 
Mominul MuKuT
Join Date: Mar 2009
Location: Bangladesh
Posts: 124
Rep Power: 17
mukut is on a distinguished road
Hello foamers!

I have modified multiRegionHeater Tutorial of chtMultiRegionFoam to following geometry, as I want to simulate plasma actuator induced flow:


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
mukut is offline   Reply With Quote

Old   October 21, 2013, 16:18
Default
  #2
otq
New Member
 
Christopher Hughes
Join Date: Oct 2012
Posts: 27
Rep Power: 13
otq is on a distinguished road
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.
otq is offline   Reply With Quote

Old   October 21, 2013, 21:26
Default
  #3
Senior Member
 
Mominul MuKuT
Join Date: Mar 2009
Location: Bangladesh
Posts: 124
Rep Power: 17
mukut is on a distinguished road
Quote:
Originally Posted by otq View Post
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..
mukut is offline   Reply With Quote

Old   January 30, 2014, 04:20
Default foam::error::printStack compressibleInterFoam parallel run
  #4
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 18
fabian_roesler is on a distinguished road
Hi mukut

How did you solve the problem with the foam::error:rintStack? 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 is offline   Reply With Quote

Old   February 6, 2014, 06:39
Default twoPhaseMixtureThermo bug?
  #5
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 18
fabian_roesler is on a distinguished road
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
fabian_roesler is offline   Reply With Quote

Old   February 6, 2014, 09:14
Default
  #6
Senior Member
 
Olivier
Join Date: Jun 2009
Location: France, grenoble
Posts: 272
Rep Power: 17
olivierG is on a distinguished road
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
vs1 likes this.
olivierG is offline   Reply With Quote

Old   May 2, 2014, 16:54
Red face
  #7
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 18
fabian_roesler is on a distinguished road
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
fabian_roesler is offline   Reply With Quote

Old   August 23, 2014, 17:40
Default
  #8
New Member
 
Cliff
Join Date: Aug 2014
Posts: 10
Rep Power: 11
cliffdub is on a distinguished road
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
cliffdub is offline   Reply With Quote

Old   August 23, 2014, 17:51
Default
  #9
otq
New Member
 
Christopher Hughes
Join Date: Oct 2012
Posts: 27
Rep Power: 13
otq is on a distinguished road
Quote:
Originally Posted by cliffdub View Post
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
otq is offline   Reply With Quote

Old   August 23, 2014, 19:13
Default
  #10
New Member
 
Cliff
Join Date: Aug 2014
Posts: 10
Rep Power: 11
cliffdub is on a distinguished road
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
cliffdub is offline   Reply With Quote

Old   August 26, 2014, 13:41
Default
  #11
Member
 
baran
Join Date: Aug 2014
Posts: 45
Rep Power: 11
baran_foam is on a distinguished road
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
Attached Images
File Type: jpg check1.jpg (31.8 KB, 121 views)

Last edited by baran_foam; August 26, 2014 at 13:44. Reason: picture is not good
baran_foam is offline   Reply With Quote

Old   August 27, 2014, 03:35
Default nonuniform value 0
  #12
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 18
fabian_roesler is on a distinguished road
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
thiagopl, vs1 and Zane like this.
fabian_roesler is offline   Reply With Quote

Old   August 27, 2014, 06:37
Default
  #13
Member
 
baran
Join Date: Aug 2014
Posts: 45
Rep Power: 11
baran_foam is on a distinguished road
Quote:
Originally Posted by fabian_roesler View Post
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

Last edited by wyldckat; August 30, 2014 at 09:41. Reason: [QUOTE] -> [CODE]
baran_foam is offline   Reply With Quote

Old   August 30, 2014, 09:50
Default
  #14
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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
vs1 likes this.
__________________
wyldckat is offline   Reply With Quote

Old   September 1, 2014, 00:55
Default
  #15
Member
 
baran
Join Date: Aug 2014
Posts: 45
Rep Power: 11
baran_foam is on a distinguished road
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
baran_foam is offline   Reply With Quote

Old   September 1, 2014, 04:25
Default
  #16
Senior Member
 
Fabian Roesler
Join Date: Mar 2009
Location: Germany
Posts: 213
Rep Power: 18
fabian_roesler is on a distinguished road
Quote:
Originally Posted by baran_foam View Post
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
fabian_roesler is offline   Reply With Quote

Old   September 2, 2014, 02:18
Default
  #17
Member
 
baran
Join Date: Aug 2014
Posts: 45
Rep Power: 11
baran_foam is on a distinguished road
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:rintStack(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

Last edited by baran_foam; September 2, 2014 at 02:50. Reason: to give full error details
baran_foam is offline   Reply With Quote

Old   September 3, 2014, 08:54
Default
  #18
Member
 
baran
Join Date: Aug 2014
Posts: 45
Rep Power: 11
baran_foam is on a distinguished road
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
wyldckat likes this.
baran_foam is offline   Reply With Quote

Old   October 28, 2015, 02:40
Default Error in chtMultiRegionFOAM
  #19
New Member
 
Federica
Join Date: Oct 2015
Posts: 13
Rep Power: 10
fkika is on a distinguished road
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!
fkika is offline   Reply With Quote

Old   October 31, 2015, 09:18
Default
  #20
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
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
vs1 likes this.
wyldckat is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, error, plasma actuator modeling, pre-processing

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
table properties for thermophysical properties romant OpenFOAM Running, Solving & CFD 1 August 12, 2014 09:41
thermophysical properties immortality OpenFOAM Running, Solving & CFD 0 December 2, 2012 07:49
polynomial thermophysical properties in a solid region (chtMultiRegionSimpleFoam) Koga OpenFOAM Programming & Development 0 November 15, 2012 05:14
how to incorporate temperature dependent thermophysical properties in fluent. CANDY Fluent UDF and Scheme Programming 4 October 22, 2012 05:19
thermophysical properties of ham Alex Ivancic Main CFD Forum 1 November 5, 1998 12:09


All times are GMT -4. The time now is 07:04.