CFD Online Discussion Forums - chtMultiRegionSimpleFoam

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - chtMultiRegionSimpleFoam (https://www.cfd-online.com/Forums/openfoam-solving/123373-chtmultiregionsimplefoam.html)

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

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

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

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,

I will check it again, then.

Thanks a lot,
Samuele

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$

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

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

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.

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

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

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?

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

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

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?

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

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

Could you have a look?

Thanks a lot,
Samuele

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

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.

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