CFD Online Discussion Forums
Page 1 of 2 1 2
Show 40 post(s) from this thread on one page

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   chtMultiRegionSimpleFoam: strange error (https://www.cfd-online.com/Forums/openfoam-solving/124620-chtmultiregionsimplefoam-strange-error.html)

samiam1000 October 9, 2013 10:22
chtMultiRegionSimpleFoam: strange error
 
Dear all,

trying to run my chtMultiRegionSimpleFoam simulation, I get this error:

Code:
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  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  : chtMultiRegionSimpleFoam
Date  : Oct 09 2013
Time  : 16:09:54
Host  : "lab-laptop"
PID    : 7471
Case  : /home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady
nProcs : 1
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 fluid for time = 0

Create solid mesh for region domain1 for time = 0

Create solid mesh for region domain2 for time = 0

Create solid mesh for region domain3 for time = 0



--> FOAM FATAL IO ERROR:
keyword type is undefined in dictionary ".T"

file: .T from line 49 to line 66.

    From function dictionary::lookupEntry(const word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 402.

FOAM exiting

lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady$
The point is that it does not suggest me which is the corrupted file. Any idea? Is the error supposed to be in the domain3 region? Or domain3 is ok and I have to look in the next ones?

Thanks a lot,
Samuele

wyldckat October 12, 2013 15:16
Hi Samuele,

Try this:
Code:
find . -name ".T"
It will tell you where the files named ".T" are located.

By the way, on Linux, files that start with "." are usually considered as hidden files. To see then, use the "-a" option with "ls" or "ll", for examples:
Code:
ls -a
ll -a
ls -la
The two last ones should be the same, because "ll" is an alias (when defined as such), as proved by this command:
Code:
alias ll
Best regards,
Bruno

samiam1000 October 14, 2013 04:04
I Bruno,

the point is that there was not such a file :-/..

..I've changed the folder case name and not I don't have this error.

But I have this new one:

Code:
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady$ chtMultiRegionSimpleFoam
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  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  : chtMultiRegionSimpleFoam
Date  : Oct 14 2013
Time  : 08:55:57
Host  : "lab-laptop"
PID    : 2667
Case  : /home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady
nProcs : 1
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 fluid for time = 0

Create solid mesh for region domain1 for time = 0

Create solid mesh for region domain2 for time = 0

Create solid mesh for region domain3 for time = 0

Create solid mesh for region domain5 for time = 0

Create solid mesh for region packs for time = 0

*** Reading fluid mesh thermophysical properties for region fluid

    Adding to thermoFluid

Selecting thermodynamics package heRhoThermo<pureMixture<const<hConst<perfectGas<specie>>,sensibleEnthalpy>>>
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/x86_64-linux-gnu/libc.so.6"
#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"
#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"
#5  Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#7 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Floating point exception (core dumped)
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady$
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady$
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady$
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady$
Could you help?

Thanks a lot,
Samuele

Tobi October 14, 2013 06:38
Hi Samuel,

the first error is nice - never seen that befor.

The last error "Floating Point exeption" occure if you are dividing by Zero.

Check all your BC and set always the "value uniform X" keyword.

But hmmm ... it seems that you get the error while updating the thermodynamics. Maybe when the constructor is selected.

But first check my first hint.
Your meshes are okay?
Boundary files okay too?

Regards
Tobi

samiam1000 October 14, 2013 07:49
Dear Tobi,

thanks for answering. Well, let me have a look at all the BC. Anyway, I am pretty confident it is ok! As far as the mesh is concerned, it should be alright.

I'll give you a feedback, soon.

Thanks a lot,
Samuele

samiam1000 October 21, 2013 05:39
Dear Tobi,

I have checked all the BC and I have done some modifications.

The point is that now I have done a small step ahead, but there still is something wrong:

Code:
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  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  : chtMultiRegionSimpleFoam
Date  : Oct 21 2013
Time  : 11:26:34
Host  : "lab-laptop"
PID    : 4069
Case  : /home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady
nProcs : 1
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 fluid for time = 0

Create solid mesh for region domain1 for time = 0

Create solid mesh for region domain2 for time = 0

Create solid mesh for region domain3 for time = 0

Create solid mesh for region domain5 for time = 0

Create solid mesh for region packs for time = 0

*** Reading fluid mesh thermophysical properties for region fluid

    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 turbulence

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
--> Upgrading k to employ run-time selectable wall functions
    Backup original k to k.old
    Writing updated k
--> Upgrading epsilon to employ run-time selectable wall functions
    Backup original epsilon to epsilon.old
    Writing updated epsilon
--> Creating mut to employ run-time selectable wall functions
    Writing new mut
--> Creating alphat to employ run-time selectable wall functions
    Writing new alphat
kEpsilonCoeffs
{
    Cmu            0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt            1;
}

    Adding to ghFluid

    Adding to ghfFluid

--> FOAM Warning :
    From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.1/src/OpenFOAM/lnInclude/Field.C at line 262
    Reading "/home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady/0/fluid/p_rgh" from line 18 to line 20
    expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
--> FOAM Warning :
    From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.1/src/OpenFOAM/lnInclude/Field.C at line 262
    Reading "/home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady/0/fluid/p_rgh.boundaryField.intake" from line 55 to line 20
    expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain1

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain2

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain3

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain5

    Adding to thermos

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

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigSegv::sigHandler(int) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::tmp<Foam::Field<double> > Foam::fvPatch::patchInternalField<double>(Foam::UList<double> const&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#4  Foam::fvPatchField<double>::patchInternalField() const in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#5  Foam::basicSymmetryFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6  Foam::symmetryFvPatchField<double>::symmetryFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7  Foam::fvPatchField<double>::adddictionaryConstructorToTable<Foam::symmetryFvPatchField<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8  Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#9  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::readField(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#10  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#11  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#12  at basicThermo.C:0
#13  Foam::basicThermo::lookupOrConstruct(Foam::fvMesh const&, char const*) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#14  Foam::basicThermo::basicThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#15  Foam::solidThermo::solidThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libsolidThermo.so"
#16  Foam::heThermo<Foam::solidThermo, Foam::pureMixture<Foam::constIsoSolidTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleEnthalpy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libsolidThermo.so"
#17  Foam::solidThermo::addfvMeshConstructorToTable<Foam::heSolidThermo<Foam::solidThermo, Foam::pureMixture<Foam::constIsoSolidTransport<Foam::species::thermo<Foam::hConstThermo<Foam::rhoConst<Foam::specie> >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libsolidThermo.so"
#18  Foam::autoPtr<Foam::solidThermo> Foam::basicThermo::New<Foam::solidThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libsolidThermo.so"
#19  Foam::solidThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libsolidThermo.so"
#20 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#21  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#22 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Segmentation fault (core dumped)
Any idea about this `new' (for me) error?

Tobi October 21, 2013 06:28
Hi,

first of all post the p_rgh file (0/fluid/p_rgh).
And the second error (Segmentation fault) is in your domain5. You have an error there.

samiam1000 October 21, 2013 06:47
Thanks Tobi.

Here is my p_rgh file:

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

dimensions      [ 1 -1 -2 0 0 0 0 ];

internalField  100000;

boundaryField
{
    wall-fluid
    {
        type            buoyantPressure;
        value          $internalField;
    }
    back
    {
        type            buoyantPressure;
        value          $internalField;
    }
    empty_1
    {
        type            buoyantPressure;
        value          $internalField;
    }
    empty_2-fluid
    {
        type            symmetryPlane;
    }
    glass
    {
        type            buoyantPressure;
        value          $internalField;
    }
    inlet_main
    {
        type            buoyantPressure;
        value          $internalField;
    }
    intake
    {
        type            fixedValue;
        value          $internalField;
    }
    pcm_surface
    {
        type            buoyantPressure;
        value          $internalField;
    }
    fluid_to_domain5
    {
        type            buoyantPressure;
        value          $internalField;
    }
    fluid_to_domain3
    {
        type            buoyantPressure;
        value          $internalField;
    }
    fluid_to_domain1
    {
        type            buoyantPressure;
        value          $internalField;
    }
    fluid_to_packs
    {
        type            buoyantPressure;
        value          $internalField;
    }
    fluid_to_domain2
    {
        type            buoyantPressure;
        value          $internalField;
    }
}


// ************************************************************************* //
and here the domain5 changeLogDict:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.1.0                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    object      changeDictionaryDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dictionaryReplacement
{
    boundary
    {
        empty_2-packs
        {
            type            symmetryPlane;
            nFaces          1342;
            startFace      47303;
        }
        wall_packs
        {
            type            wall;
            nFaces          705;
            startFace      48645;
        }
        domain5_to_fluid
        {
            type            mappedWall;
            nFaces          3997;
            startFace      49350;
            sampleMode      nearestPatchFace;
            sampleRegion    fluid;
            samplePatch    fluid_to_domain1;
            offsetMode      uniform;
            offset          (0 0 0);
        }
    }

    T
    {
        internalField  uniform 298;

        boundaryField
        {
            empty_2-packs
            {
                type            symmetryPlane;
            }
            wall_packs
            {
                type            zeroGradient;
            }
            domain5_to_fluid
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                value          uniform 298;
                neighbourFieldName T;
                kappa          solidThermo;
                kappaName      none;
            }
        }
    }

    p
    {
        internalField  uniform 100000;

        boundaryField
        {
            empty_2-packs
            {
                type            symmetryPlane;
            }
            wall_packs
            {
                type            calculated;
                value          uniform 0;
            }
            domain5_to_fluid
            {
                type            calculated;
                value          uniform 0;
            }
        }
    }
}
Thanks for answering and helping,
Samuele

samiam1000 October 21, 2013 09:26
Hi Tobi,

following your advice, I solved an error, but I now get this message:

Code:
Attempt to cast type calculated to type compressible::turbulentTemperatureCoupledBaffleMixed
However, I think that my BC are properly set. Where may the error be, according to you?

Thanks a lot,
Samuele

Tobi October 21, 2013 11:02
Hi

Code:
      domain5_to_fluid
        {
            type            mappedWall;
            nFaces          3997;
            startFace      49350;
            sampleMode      nearestPatchFace;
            sampleRegion    fluid;
            samplePatch    fluid_to_domain1;  // should be fluid_to_domain5
            offsetMode      uniform;
            offset          (0 0 0);
        }
The other error... can you give me more information?

samiam1000 October 21, 2013 11:12
Yeah: it's right.

About the `new error', here is the complete message:

Code:
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  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  : chtMultiRegionSimpleFoam
Date  : Oct 21 2013
Time  : 16:10:36
Host  : "lab-laptop"
PID    : 2401
Case  : /home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady
nProcs : 1
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 fluid for time = 0

Create solid mesh for region domain1 for time = 0

Create solid mesh for region domain2 for time = 0

Create solid mesh for region domain3 for time = 0

Create solid mesh for region domain5 for time = 0

Create solid mesh for region packs for time = 0

*** Reading fluid mesh thermophysical properties for region fluid

    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 turbulence

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu            0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt            1;
}

    Adding to ghFluid

    Adding to ghfFluid

--> FOAM Warning :
    From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.1/src/OpenFOAM/lnInclude/Field.C at line 262
    Reading "/home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady/0/fluid/p_rgh" from line 18 to line 20
    expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
--> FOAM Warning :
    From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.1/src/OpenFOAM/lnInclude/Field.C at line 262
    Reading "/home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady/0/fluid/p_rgh.boundaryField.intake" from line 55 to line 20
    expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain1

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain2

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain3

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region domain5

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

*** Reading solid mesh thermophysical properties for region packs

    Adding to thermos

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

    Adding to radiations

Radiation model not active: radiationProperties not found
Selecting radiationModel none
    Adding fvOptions

No finite volume options present

Time = 1


Solving for fluid region fluid
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00545361, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.018735, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.00880018, No Iterations 1


--> FOAM FATAL ERROR:
Attempt to cast type zeroGradient to type compressible::turbulentTemperatureCoupledBaffleMixed

    From function refCast<To>(From&)
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.1/src/OpenFOAM/lnInclude/typeInfo.H at line 114.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#3  Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4  Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5  Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#6 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#7 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 
 in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Aborted (core dumped)

Tobi October 21, 2013 12:02
T file of fluid please and the boundary file of your fluid domain.

samiam1000 October 23, 2013 02:56
Tobi,

pardon for the late answer. I was out of office and I did not have the case with me. Sorry,

Anyway, the file that you want is here:

Code:
    T
    {
        internalField  uniform 298;

        boundaryField
        {
            wall-fluid
            {
                type            zeroGradient;
                value          uniform 273;
            }
            back
            {
                type            zeroGradient;
//                type            fixedValue;
                value          uniform 273;
            }
            empty_1
            {
                type            zeroGradient;
                value          uniform 273;
            }
            empty_2-fluid
            {
                type            symmetryPlane;
            }
            glass
            {
                type            fixedValue;
                value          uniform 298;
            }
            inlet_main
            {
                type            fixedValue;
                value          uniform 273;
            }
            intake
            {
                type            zeroGradient;
                value          uniform 273;
            }
            pcm_surface
            {
                type            fixedValue;
                value          uniform 275;
            }
            fluid_to_domain5
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                neighbourFieldName T;
                kappa          fluidThermo;
                kappaName      none;
                value          uniform 278;
            }
            fluid_to_domain3
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                neighbourFieldName T;
                kappa          fluidThermo;
                kappaName      none;
                value          uniform 278;
            }
            fluid_to_domain1
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                neighbourFieldName T;
                kappa          fluidThermo;
                kappaName      none;
                value          uniform 278;
            }
            fluid_to_packs
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                neighbourFieldName T;
                kappa          fluidThermo;
                kappaName      none;
                value          uniform 278;
            }
            fluid_to_domain2
            {
                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                neighbourFieldName T;
                kappa          fluidThermo;
                kappaName      none;
                value          uniform 278;
            }
        }
    }
Is there something strange?

samiam1000 October 23, 2013 04:37
Tobi,

I have solved that problem, changing the BC. But I still have troubles: now I get this message:

Code:
Solving for solid region domain1


--> FOAM FATAL ERROR:

    gradientInternalCoeffs cannot be called for a calculatedFvPatchField
    on patch domain1_to_fluid of field h in file "/home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady/0/domain1/h"
    You are probably trying to solve for a field with a default boundary condition.

    From function calculatedFvPatchField<Type>::gradientInternalCoeffs() const
    in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 199.

FOAM exiting
The point is that I don't have any h file, anywhere.

Any idea about this?

Tobi October 23, 2013 07:03
Hi,

the T-file Looks good.
The Problem in the cht solver is always that the Domains are related to each other and the Problem could be in both. I worked a Long time with the cht but at the Moment I am not working anymore with it.

Therefor I am getting rusty in the error Messages :D

Is it possible to upload you 0 Folder and the constant Folder, best way would be the whole case.

Maybe its not possible for you (due to secret) it would be very good to have the 0 Folder and the constant/***/polyMesh/boundary files to have a look at the whole case. It would be the best way to give you a correct answer.

Otherwise you should always get more Information :)

So if you solving "domain1_to_fluid" it is interessting to have the related files of both Domains.



For your Problem it could be possible that you use "fluidThermo" in the T file for your solid Domain. there should be "solidThermo".

Regards Tobi

samiam1000 October 23, 2013 08:20
I can share the case with you via Skype. I have tried to contact you there.

If you are online, please let me know. Otherwise, please write me at samuele.zampini@polimi.com and I'll send you tha case.

Thank you very much,
Samuele

Tobi October 23, 2013 08:39
Hi Samuel,

i am online today evening (in 2 hours).
Otherwise you can send me the case Tobias.Holzmann@Holzmann-cfd.de

Regards

Tobi

Ahmed Khattab October 23, 2013 11:10
Quote:
Originally Posted by samiam1000 (Post 458454)
Tobi,

I have solved that problem, changing the BC. But I still have troubles: now I get this message:

Code:
Solving for solid region domain1


--> FOAM FATAL ERROR:

    gradientInternalCoeffs cannot be called for a calculatedFvPatchField
    on patch domain1_to_fluid of field h in file "/home/lab/Documenti/Ethics/FRISBEE/CFD/testVerticalTN_PCM/PCMPacks/steady/0/domain1/h"
    You are probably trying to solve for a field with a default boundary condition.

    From function calculatedFvPatchField<Type>::gradientInternalCoeffs() const
    in file fields/fvPatchFields/basic/calculated/calculatedFvPatchField.C at line 199.

FOAM exiting
The point is that I don't have any h file, anywhere.

Any idea about this?

Hi,

i experienced such a problem before, already today. i found that my mistake was in changeDirectoryDict where i wrongly wrote than name of interface patch, i simply change uppercase letter to lowercase. that's it.

i hope it will be helpful for you. :)

samiam1000 October 23, 2013 11:24
I'm pretty sure about the patch names, but I'll check them again!

Thanks!

Tobi October 23, 2013 13:09
Hi,


1. your changeDict has a wrong entry like ahmed mentioned in your domain2
Code:
            domain2_to_fluid
            {
                type            mappedWall;
                nFaces          3997;
                startFace      49350;
                sampleMode      nearestPatchFace;
                sampleRegion    fluid;
                samplePatch    fluid_to_domain1;  //domain2
                offsetMode      uniform;
                offset          (0 0 0);
            }
        }
2. wrong schemes and solution for your solids

3. wrong changeDict in domain1
Code:
    T
    {
        internalField  uniform 298;

        boundaryField
        {
                    empty_2-packs
                    {
                        type            symmetryPlane;
                    }
                    wall_packs
                    {
                        type            zeroGradient        // forget ";"
                    }
                    domain1_to_fluid
                    {
                                type            compressible::turbulentTemperatureCoupledBaffleMixed;
                                value          uniform 298;
                                neighbourFieldName T;
                                kappa          solidThermo;
                                kappaName      none;
                    }
        }
4. after changeDict on domain2 there are still one error in your "0/domain2/T" file:
Code:
  empty_2-packs
    {
        type            symmetryPlane;
    }
    wall_packs
    {
        type            zeroGradient
    }

    domain1_to_fluid
    {
      type compressible::turbulentTemperatureCoupledBaffleMixed;
        value          uniform 298;
        neighbourFieldName T;
        kappa          solidThermo;
        kappaName      none;
    }
    domain1_to_fluid                      // wrong - not needed
    {
        type            calculated;
        value          uniform 0;
        neighbourFieldName T;
        kappa          solidThermo;
        kappaName      none;
    }
After the changes you will get the following output:
Code:
Time = 1


Solving for fluid region fluid
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0024031, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.0099621, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.0132661, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.0987243, No Iterations 4
Min/max T:269.145 301.75
DICPCG:  Solving for p_rgh, Initial residual = 0.999968, Final residual = 0.0094627, No Iterations 19
time step continuity errors : sum local = 22257.7, global = -0.00312494, cumulative = -0.00312494
Min/max rho:0.001 1.29413

Solving for solid region domain1
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0184498, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 297.999 max(T) [0 0 0 1 0 0 0] 298.001

Solving for solid region domain2
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0187101, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 276.999 max(T) [0 0 0 1 0 0 0] 277.003

Solving for solid region domain3
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0906404, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 297.999 max(T) [0 0 0 1 0 0 0] 298.001

Solving for solid region domain5
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0785481, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 297.999 max(T) [0 0 0 1 0 0 0] 298.002

Solving for solid region packs
DICPCG:  Solving for h, Initial residual = 1, Final residual = 0.0262329, No Iterations 2
Min/max T:min(T) [0 0 0 1 0 0 0] 297.999 max(T) [0 0 0 1 0 0 0] 298.001
ExecutionTime = 5.9 s  ClockTime = 6 s
- the second iteration is bad but I switched off the turbulence modelling and do not check the BC.

Have fun and good luck.

Regard
Tobi


All times are GMT -4. The time now is 22:01.
Page 1 of 2 1 2
Show 40 post(s) from this thread on one page