[samiam1000]

Dear All,

I am trying to run chtMultiRegionSimpleFoam.

I have prepared a case with 1 fluid region and 6 solid regions.

I have set my case and when I lauch it, I get this error:

    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

Time = 1


Solving for fluid region part_2-solid
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.00888973, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.00350072, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 1, Final residual = 0.00662661, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 1, Final residual = 0.000600295, No Iterations 1
Min/max T:270 300
GAMG:  Solving for p_rgh, Initial residual = 0.94945, Final residual = 0.00676794, No Iterations 7
time step continuity errors : sum local = 0.316082, global = -0.025282, cumulative = -0.025282
Min/max rho:1.15862 1.28736

Solving for solid region pcm


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

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

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField const& Foam::refCast<Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField const, Foam::fvPatchField<double> const>(Foam::fvPatchField<double> const&) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#3  Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#4  Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5  Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoeffs() in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#7  Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#8  
 in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#9  
 in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  
 in "/home/zampini/OpenFOAM/OpenFOAM-2.2.0/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Annullato
zampini@pc-zampini:~/Documenti/personali/Epta/SCC/steady$

I can not understand what it means. Could you help, please?

Thanks a lot,
Samuele

[wyldckat]

Hi Samuele,

Can you attach the boundary conditions for the "pcm" region, as well as any regions that are in touch with it?

Although, if I have to guess, I think that this sentence:

Attempt to cast type zeroGradient to type compressible::turbulentTemperatureCoupledBaffleMixed

means that for a certain patch, on one region you have "zeroGradient" and on the other region you have "turbulentTemperatureCoupledBaffleMixed". But the two seem to be incompatible.

Best regards,
Bruno

[samiam1000]

Dear Bruno,

thanks for answering, first. Also, about your suggestions, I completely agree. The point is that I chacked the boundary conditions twice and I can't find any error. I am attaching to this email all the changeLogDictionary I use. Could you kindly have a look and let me know what's wrong with it?

Thanks a lot,
Samuele

[Ahmed Khattab]

Dear Samuele,

it seems that you put this B.C for a patch in a region (compressible::turbulentTemperatureCoupledBaffleMi xed). then you put zeroGradient to the same patch in the other region which is not applicable. the B.C must be the same for same patch in different regions.

hope it helps.
BR,

[samiam1000]

I will check it again, then.

Thanks a lot,
Samuele

[samiam1000]

I have found a first error: I haven't set the right BC in the fluid region.

I ran my case and I get a different error: could you help in solving this, too?

Thanks a lot,
Samuele

The error message is:

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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 : Sep 12 2013
Time : 11:47:31
Host : "lab-laptop"
PID : 7539
Case : /home/lab/Documenti/Ethics/FRISBEE/CFD/SCC/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 part_2-solid for time = 0
Create solid mesh for region pcm for time = 0
Create solid mesh for region packs_1 for time = 0
Create solid mesh for region packs_2 for time = 0
Create solid mesh for region part_2-solid.1 for time = 0
Create solid mesh for region domain2 for time = 0
Create solid mesh for region domain5 for time = 0
*** Reading fluid mesh thermophysical properties for region part_2-solid
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 laminar
Adding to ghFluid
Adding to ghfFluid
Radiation model not active: radiationProperties not found
Selecting radiationModel none
Adding fvOptions
No finite volume options present
*** Reading solid mesh thermophysical properties for region pcm
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_1
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_2
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 part_2-solid.1
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 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
Time = 1
　
Solving for fluid region part_2-solid
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.00888973, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.00350072, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 0.00662661, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 0.00187789, No Iterations 1
Min/max T:270 300
GAMG: Solving for p_rgh, Initial residual = 0.949213, Final residual = 0.00675391, No Iterations 7
time step continuity errors : sum local = 0.316054, global = -0.02528, cumulative = -0.02528
Min/max rho:1.15862 1.28736
Solving for solid region pcm
DICPCG: Solving for h, Initial residual = 1, Final residual = 0.0743357, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 270 max(T) [0 0 0 1 0 0 0] 300
Solving for solid region packs_1
DICPCG: Solving for h, Initial residual = 1, Final residual = 0.0768139, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 288.897 max(T) [0 0 0 1 0 0 0] 300
Solving for solid region packs_2
#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::polyMeshTetDecomposition::findFaceBasePts(Foam::polyMesh const&, double, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::polyMesh::tetBasePtIs() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::mappedPatchBase::facePoints(Foam::polyPatch const&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#6 Foam::mappedPatchBase::calcMapping() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#7 void Foam::mappedPatchBase::distribute<double>(Foam::List<double>&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10 Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#11 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#12 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#13 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#15 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Segmentation fault (core dumped)
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/SCC/steady$

[Ahmed Khattab]

Dear Samuele,

this error always appears when there is a problem with three factors: cell size, time step, velocity value. you must compromise the three factors to get smooth run.

hint: adjust dimensions and cell size such that patches boundaries lays on cells boundaries not in middle of it.

hope it helps,

Best Regards,

Ahmed

[samiam1000]

Do you mean that it could be necessary to remesh my geometry?

[Ahmed Khattab]

no i don't mean so, i only drag your attention to revise it, but you can change your velocity value, or time step. you can go back to the user manual to know how to set time step.

[samiam1000]

I think that being a steady simulation, the time-step does not influence the solution: is this right?

[wyldckat]

Greetings to all!

@Samuele:

Originally Posted by samiam1000

I think that being a steady simulation, the time-step does not influence the solution: is this right?

AFAIK, that is correct. The time step doesn't matter in steady-state simulations... at least not usually.
The unusual example might be LTSInterFoam: http://www.openfoam.org/version2.0.0/steady-vof.php

Originally Posted by samiam1000

I ran my case and I get a different error: could you help in solving this, too?

#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::polyMeshTetDecomposition::findFaceBasePts(Foam::polyMesh const&, double, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::polyMesh::tetBasePtIs() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"

It gave a SIGSEGV signal (#1) when trying to perform "polyMeshTetDecomposition :: findFaceBasePts" (#3). I'm guessing that this means that there is something very wrong with your mesh. Do a complete check mesh by running:

checkMesh -allGeometry -allTopology

Best regards,
Bruno

[samiam1000]

Dear Bruno,

thanks for answering and pardon for the late reply.

This is the output of the command you suggested:

zampini@pc-zampini:~/Documenti/personali/Epta/SCC/steady$ checkMesh -allGeometry -allTopology
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.0
Exec   : checkMesh -allGeometry -allTopology
Date   : Sep 15 2013
Time   : 18:37:51
Host   : "pc-zampini"
PID    : 1944
Case   : /home/zampini/Documenti/personali/Epta/SCC/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 polyMesh for time = 0

Enabling all (cell, face, edge, point) topology checks.

Enabling all geometry checks.

Time = 0

Mesh stats
    points:           621150
    faces:            1816951
    internal faces:   1761449
    cells:            596440
    faces per cell:   5.9996
    boundary patches: 13
    point zones:      0
    face zones:       14
    cell zones:       5

Overall number of cells of each type:
    hexahedra:     596200
    prisms:        240
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Topological cell zip-up check OK.
    Number of identical duplicate faces (baffle faces): 3200
    Face-face connectivity OK.
  <<Writing 6400 faces with non-standard edge connectivity to set edgeFaces
  <<Writing 4 cells with two non-boundary faces to set twoInternalFacesCells
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                   Bounding box
    deflettori          6400     3403     multiply connected (shared edge)   (-0.1 -0.05 0.45) (1.3 0.3 0.5)
    foam-pcm            4200     4402     ok (non-closed singly connected)   (-0.05 0 -0.05) (1.25 0.3 -0.05)
    foam-part_2-solid   14391    14810    ok (non-closed singly connected)   (-0.1 -0.05 -0.05) (1.3 0.3 0.75)
    foam-part_2-solid.1 8400     8662     ok (non-closed singly connected)   (-0.05 -0.05 -0.05) (1.25 0 0.45)
    glass               5600     5781     ok (non-closed singly connected)   (-0.1 -0.05 0.75) (1.3 0.3 0.75)
    inlet_1             400      451      ok (non-closed singly connected)   (0.55 -0.05 -0.05) (0.6 0.3 -0.05)
    inlet_2             400      451      ok (non-closed singly connected)   (0.6 -0.05 -0.05) (0.65 0.3 -0.05)
    intake_1            400      451      ok (non-closed singly connected)   (-0.1 -0.05 -0.05) (-0.05 0.3 -0.05)
    intake_2            400      451      ok (non-closed singly connected)   (1.25 -0.05 -0.05) (1.3 0.3 -0.05)
    symmetry-packs_1    2250     2346     ok (non-closed singly connected)   (0 0.3 0) (0.5 0.3 0.45)
    symmetry-pcm        2500     2832     ok (non-closed singly connected)   (-0.05 0.3 -0.05) (1.25 0.3 0.45)
    symmetry-part_2-solid7911     8250     ok (non-closed singly connected)   (-0.1 0.3 -0.05) (1.3 0.3 0.75)
    symmetry-packs_2    2250     2346     ok (non-closed singly connected)   (0.7 0.3 0) (1.2 0.3 0.45)
  <<Writing 3391 conflicting points to set nonManifoldPoints

Checking geometry...
    Overall domain bounding box (-0.1 -0.05 -0.05) (1.3 0.3 0.75)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (3.24011e-16 3.56288e-17 1.47307e-15) OK.
    Max cell openness = 3.10576e-16 OK.
    Max aspect ratio = 13.8392 OK.
    Minimum face area = 7.67269e-06. Maximum face area = 0.000250165.  Face area magnitudes OK.
    Min volume = 3.06526e-08. Max volume = 2.41892e-06.  Total volume = 0.392.  Cell volumes OK.
    Mesh non-orthogonality Max: 54.2303 average: 11.4654
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 1.05523 OK.
    Coupled point location match (average 0) OK.
    Face tets OK.
    Min/max edge length = 0.00290397 0.0199751 OK.
    All angles in faces OK.
    Face flatness (1 = flat, 0 = butterfly) : average = 1  min = 1
    All face flatness OK.
    Cell determinant (wellposedness) : minimum: 0 average: 5.679
 ***Cells with small determinant found, number of cells: 80
  <<Writing 80 under-determined cells to set underdeterminedCells
    Concave cell check OK.

Failed 1 mesh checks.

End

Could you help?

[wyldckat]

Hi Samuele,

Originally Posted by samiam1000

    Cell determinant (wellposedness) : minimum: 0 average: 5.679
 ***Cells with small determinant found, number of cells: 80
  <<Writing 80 under-determined cells to set underdeterminedCells
    Concave cell check OK.

Oh, this is bad, very bad! Cell determinant values of "0" basically means that there is either a contorted cell that crosses over itself or that it's a cell without volume. This is probably what's triggering the crash!

paraFoam provides you with the ability to also see the sets. Turn on that option and choose to see "underdeterminedCells" that should appear in the same list as the patches. Then try to see where exactly where the problem cells are and try to re-do your mesh.

Another possibility is to follow the example shown here: http://openfoamwiki.net/index.php/SetSet#Usage_example - more specifically, to only remove the cells associated to "underdeterminedCells". But keep in mind that this kind of cell removal strategy has certain limitations, such as possibly and wrongly removing some important cells.

Best regards,
Bruno

[samiam1000]

Dear Bruno,

I am attaching a picture of the whole volume where it is evident where the underdeterminedCells are.

I can't understand what's wrong with them. Do you have any idea?

First of all, I will try your suggestions.

Thanks a lot,

Samuele

[samiam1000]

Dear Bruno, Dear All,

after having re-meshed my geometry, I get a very strange result.

First of all, all the mesh checks are ok! Hance I thought that my simulation would have started immediately, but..
..but I got this error:

lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/SCC/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 : Sep 18 2013
Time : 12:36:14
Host : "lab-laptop"
PID : 4908
Case : /home/lab/Documenti/Ethics/FRISBEE/CFD/SCC/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 part_2-solid for time = 0
Create solid mesh for region pcm for time = 0
Create solid mesh for region packs_1 for time = 0
Create solid mesh for region packs_2 for time = 0
Create solid mesh for region part_2-solid.1 for time = 0
Create solid mesh for region domain2 for time = 0
Create solid mesh for region domain5 for time = 0
*** Reading fluid mesh thermophysical properties for region part_2-solid
Adding to thermoFluid
Selecting thermodynamics package 
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
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::fluidThermo::fluidThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#16 Foam::rhoThermo::rhoThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#17 Foam::heThermo<Foam::rhoThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#18 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"
#19 Foam::autoPtr<Foam::rhoThermo> Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#20 Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#21 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#22 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#23 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Segmentation fault (core dumped)
lab@lab-laptop:~/Documenti/Ethics/FRISBEE/CFD/SCC/steady$

And trying to open the geometry with paraview, it crashes when I want to view the patches and not the internal mesh. I get a segmentation fault.
And this happens for each region.

Do you have any idea? Could you help?

[wyldckat]

Hi Samuele,

There seems to be a problem with a patch that is defined to be a symmetry plane:

#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"

But without access to the files, I cannot see the actual problem for myself.

Best regards,
Bruno

[samiam1000]

Here is the case: https://www.dropbox.com/sh/tgdwuqkfodgdffk/zFNdkvUgjH

Could you have a look?

Thanks a lot,
Samuele

[Tobi]

Solving for solid region packs_2
#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::polyMeshTetDecomposition::findFaceBasePts(Foam::polyMesh const&, double, bool) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::polyMesh::tetBasePtIs() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::mappedPatchBase::facePoints(Foam::polyPatch const&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#6 Foam::mappedPatchBase::calcMapping() const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libmeshTools.so"
#7 void Foam::mappedPatchBase::distribute<double>(Foam::List<double>&) const in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8 Foam::compressible::turbulentTemperatureCoupledBaffleMixedFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so"
#9 Foam::mixedFvPatchField<double>::evaluate(Foam::UPstream::commsTypes) in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#10 Foam::mixedEnergyFvPatchScalarField::updateCoeffs() in "/opt/openfoam221/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#11 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#12 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#13 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
#14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#15 
in "/opt/openfoam221/platforms/linux64GccDPOpt/bin/chtMultiRegionSimpleFoam"
Segmentation fault (core dumped)

This is not a problem of

this error always appears when there is a problem with three factors: cell size, time step, velocity value. you must compromise the three factors to get smooth run.

... if I am not wrong I know this error and I searched long time to find the trivial Problem.
If I am right you have a Problem in the Boundary Conditions of T in packs_2 or in your fluid Region. Maybe your boundary file is wrong (patch type).

Maybe you have no value set.?

Your mesh seems okay.

Regards
Tobi

[samiam1000]

Dear Tobi,

thanks for answering. Actually, I have solved this very problem (it was due to a bad definition of the boudary conditions) and the simulation's running.
However, the temperature seems to be meaningless: I do have a max temperature of about 740000 K. Too much, I say.

I am going to check this problem, too. Any idea to begin to investigate the issue?

Thanks a lot,
Samuele.

[Tobi]

Hi,

write out the first 10 or 20 integrations and have a look at your Domain.
You will be able to see the regions where you get the high temperature values. It could be possible that this Problem occure due to a mesh Problem. otherwise you see if your BC are incorrect or your Settings are wrong.

Good luck