CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

chtMultiRegionSimpleFoam

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

Like Tree1Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   September 11, 2013, 07:22
Default chtMultiRegionSimpleFoam
  #1
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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

Last edited by samiam1000; September 11, 2013 at 09:43.
samiam1000 is offline   Reply With Quote

Old   September 11, 2013, 17:22
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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
wyldckat is offline   Reply With Quote

Old   September 12, 2013, 03:27
Default
  #3
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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
Attached Files
File Type: zip archivio.zip (7.2 KB, 10 views)
samiam1000 is offline   Reply With Quote

Old   September 12, 2013, 04:53
Default
  #4
Senior Member
 
Ahmed Khattab's Avatar
 
ahmed
Join Date: Feb 2010
Posts: 161
Blog Entries: 1
Rep Power: 7
Ahmed Khattab is on a distinguished road
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,
Ahmed Khattab is offline   Reply With Quote

Old   September 12, 2013, 05:08
Default
  #5
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
I will check it again, then.

Thanks a lot,
Samuele
samiam1000 is offline   Reply With Quote

Old   September 12, 2013, 06:07
Default
  #6
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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$
samiam1000 is offline   Reply With Quote

Old   September 12, 2013, 06:37
Default
  #7
Senior Member
 
Ahmed Khattab's Avatar
 
ahmed
Join Date: Feb 2010
Posts: 161
Blog Entries: 1
Rep Power: 7
Ahmed Khattab is on a distinguished road
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
Ahmed Khattab is offline   Reply With Quote

Old   September 12, 2013, 06:49
Default
  #8
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
Do you mean that it could be necessary to remesh my geometry?
samiam1000 is offline   Reply With Quote

Old   September 12, 2013, 06:59
Default
  #9
Senior Member
 
Ahmed Khattab's Avatar
 
ahmed
Join Date: Feb 2010
Posts: 161
Blog Entries: 1
Rep Power: 7
Ahmed Khattab is on a distinguished road
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.
Ahmed Khattab is offline   Reply With Quote

Old   September 12, 2013, 08:19
Default
  #10
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
I think that being a steady simulation, the time-step does not influence the solution: is this right?
samiam1000 is offline   Reply With Quote

Old   September 14, 2013, 10:27
Default
  #11
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

@Samuele:
Quote:
Originally Posted by samiam1000 View Post
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 View Post
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
wyldckat is offline   Reply With Quote

Old   September 15, 2013, 12:39
Default
  #12
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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?
samiam1000 is offline   Reply With Quote

Old   September 15, 2013, 14:58
Default
  #13
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Samuele,

Quote:
Originally Posted by samiam1000 View Post
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.
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
wyldckat is offline   Reply With Quote

Old   September 16, 2013, 03:04
Default
  #14
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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
Attached Images
File Type: jpg Schermata del 2013-09-16 09:01:04.jpg (86.8 KB, 40 views)
samiam1000 is offline   Reply With Quote

Old   September 18, 2013, 06:41
Default
  #15
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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?
samiam1000 is offline   Reply With Quote

Old   September 21, 2013, 15:22
Default
  #16
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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
wyldckat is offline   Reply With Quote

Old   September 22, 2013, 16:43
Default
  #17
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
Here is the case: https://www.dropbox.com/sh/tgdwuqkfodgdffk/zFNdkvUgjH

Could you have a look?

Thanks a lot,
Samuele
samiam1000 is offline   Reply With Quote

Old   September 23, 2013, 06:48
Default
  #18
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,087
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
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
Tobi is offline   Reply With Quote

Old   September 23, 2013, 08:02
Default
  #19
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 490
Rep Power: 9
samiam1000 is on a distinguished road
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.
samiam1000 is offline   Reply With Quote

Old   September 24, 2013, 06:20
Default
  #20
Senior Member
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Leoben (Austria)
Posts: 1,087
Blog Entries: 6
Rep Power: 19
Tobi will become famous soon enough
Send a message via ICQ to Tobi Send a message via Skype™ to 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
Tobi is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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

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



All times are GMT -4. The time now is 22:23.