CFD Online Discussion Forums

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

samiam1000 September 11, 2013 07:22

chtMultiRegionSimpleFoam
 
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:

Code:

    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 September 11, 2013 17:22

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:
Code:

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 September 12, 2013 03:27

1 Attachment(s)
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 September 12, 2013 04:53

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 September 12, 2013 05:08

I will check it again, then.

Thanks a lot,
Samuele

samiam1000 September 12, 2013 06:07

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:
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 : 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 September 12, 2013 06:37

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 September 12, 2013 06:49

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

Ahmed Khattab September 12, 2013 06:59

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 September 12, 2013 08:19

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

wyldckat September 14, 2013 10:27

Greetings to all!

@Samuele:
Quote:

Originally Posted by samiam1000 (Post 451343)
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

Quote:

Originally Posted by samiam1000 (Post 451306)
I ran my case and I get a different error: could you help in solving this, too?

Code:

#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:
Code:

checkMesh -allGeometry -allTopology
Best regards,
Bruno

samiam1000 September 15, 2013 12:39

Dear Bruno,

thanks for answering and pardon for the late reply.

This is the output of the command you suggested:

Code:

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 September 15, 2013 14:58

Hi Samuele,

Quote:

Originally Posted by samiam1000 (Post 451828)
Code:

    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.


:eek: 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 September 16, 2013 03:04

1 Attachment(s)
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 September 18, 2013 06:41

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:
Code:

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 September 21, 2013 15:22

Hi Samuele,

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

Code:

#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 September 22, 2013 16:43

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

Could you have a look?

Thanks a lot,
Samuele

Tobi September 23, 2013 06:48

Code:

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
Code:

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 September 23, 2013 08:02

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 September 24, 2013 06:20

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


All times are GMT -4. The time now is 15:47.