[philippconen]

Dear Foamers,

while doing a little benchmark with solids4Foam and ANSYS Mechanical I run some simulations of a Cantilever. The mesh is made in ANSYS and imported via fluentMeshToFoam.

Dealing with 150.000 cells made the model in OpenFOAM stuck. On my local computer I saw an increasing amount of used RAM until nothing more was available, so my computer crashed.
Also on the super computer of my university the simulation seems to stucks (without knowledge about the used resources - ?/64GB RAM).

My question would be, if the stucked simulation is related to the warning OpenFOAM gave me? What could be the reason for the stuck? How can I fix it?

Here are my pre logs:

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
| solids4Foam::convertCaseFormat start                                |
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

foam-extend-4.0 loaded: no changes made

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1-20200408                       |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-20200408
Exec   : fluentMeshToFoam mech0_005.msh -writeSets
Date   : Oct 05 2021
Time   : 10:54:04
Host   : "cc-v01"
PID    : 12758
CtrlDict : "/home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005/system/controlDict"
Case   : /home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Dimension of grid: 3
Number of points: 168861

number of faces: 468500
Number of cells: 150000
Reading points
Reading points
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Reading uniform faces
Other readCellGroupData: 2 1 249f0 1 4
Reading uniform cells
Read zone1:1 name:interior-solid patchTypeID:interior
Reading zone data
Read zone1:2 name:solid patchTypeID:solid
Reading zone data
Read zone1:5 name:support patchTypeID:wall
Reading zone data
Read zone1:6 name:tip patchTypeID:wall
Reading zone data
Read zone1:7 name:top patchTypeID:wall
Reading zone data
Read zone1:8 name:bot patchTypeID:wall
Reading zone data
Read zone1:9 name:right patchTypeID:wall
Reading zone data
Read zone1:10 name:left patchTypeID:wall
Reading zone data


FINISHED LEXING


dimension of grid: 3
Creating shapes for 3-D cells
Building patch-less mesh...--> FOAM Warning : 
    From function polyMesh::polyMesh(... construct from shapes...)
    in file meshes/polyMesh/polyMeshFromShapeMesh.C at line 627
    Found 37000 undefined faces in mesh; adding to default patch.
done.

Building boundary and internal patches.
maxZoneID: 10
Creating patch 0 for zone: 1 start: 1 end: 431500 type: interior name: interior-solid
Creating patch 1 for zone: 5 start: 431501 end: 432000 type: wall name: support
Creating patch 2 for zone: 6 start: 432001 end: 432500 type: wall name: tip
Creating patch 3 for zone: 7 start: 432501 end: 447500 type: wall name: top
Creating patch 4 for zone: 8 start: 447501 end: 462500 type: wall name: bot
Creating patch 5 for zone: 9 start: 462501 end: 465500 type: wall name: right
Creating patch 6 for zone: 10 start: 465501 end: 468500 type: wall name: left
Writing patch interior-solid of size 431500 to faceSet.
Patch interior-solid is internal to the mesh  and is not being added to the boundary.
Writing patch support of size 500 to faceSet.
Adding new patch support of type wall as patch 0
Writing patch tip of size 500 to faceSet.
Adding new patch tip of type wall as patch 1
Writing patch top of size 15000 to faceSet.
Adding new patch top of type wall as patch 2
Writing patch bot of size 15000 to faceSet.
Adding new patch bot of type wall as patch 3
Writing patch right of size 3000 to faceSet.
Adding new patch right of type wall as patch 4
Writing patch left of size 3000 to faceSet.
Adding new patch left of type wall as patch 5

Default patch type set to empty

Writing mesh... to "constant/polyMesh"  done.

Only one cell group: no set written

End

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1-20200408                       |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-20200408
Exec   : setsToZones -noFlipMap
Date   : Oct 05 2021
Time   : 10:54:07
Host   : "cc-v01"
PID    : 12759
CtrlDict : "/home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005/system/controlDict"
Case   : /home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Searched : "constant/polyMesh/sets"
Found    : 
7
(
interior-solid
tip
bot
right
top
left
support
)


Adding set right as a faceZone.
Adding set top as a faceZone.
Adding set left as a faceZone.
Adding set interior-solid as a faceZone.
Adding set tip as a faceZone.
Adding set bot as a faceZone.
Adding set support as a faceZone.
Writing mesh.

End

Here are my solver logs:

// ************************************************************************* //
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | foam-extend: Open Source CFD                    |
|  \\    /   O peration     | Version:     4.1-20200408                       |
|   \\  /    A nd           | Web:         http://www.foam-extend.org         |
|    \\/     M anipulation  | For copyright notice see file Copyright         |
\*---------------------------------------------------------------------------*/
Build  : 4.1-20200408
Exec   : solids4Foam
Date   : Oct 05 2021
Time   : 10:54:08
Host   : "cc-v01"
PID    : 12762
CtrlDict : "/home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005/system/controlDict"
Case   : /home/ms/p/pc196535/Schreibtisch/MasterthesisDLR/Simulationen/solids4Foam/Sensitivitaetsanalyse/Balken_Sensitivitaetsanalyse/Balken_tipLoad/Balken_tipLoad_0_005
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

/*---------------------------------------------------------------------------*\
|    For further information on the solids4Foam toolbox implementations,      |
|    please see the following publications:                                   |
|                                                                             |
|    P. Cardiff, A Karac, P. De Jaeger, H. Jasak, J. Nagy, A. Ivankovic,      |
|    Z. Tukovic: An open-source finite volume toolbox for solid mechanics and |
|    fluid-solid interaction simulations. arXiv:1808.10736v2, 2018, available |
|    at https://arxiv.org/abs/1808.10736.                                     |
|                                                                             |
|    Z. Tukovic, A. Karac, P. Cardiff, H. Jasak, A. Ivankovic: OpenFOAM       |
|    finite volume solver for fluid-solid interaction.  Transactions of       |
|    Famena, 42 (3), pp. 1-31, 2018, 10.21278/TOF.42301.                      |
\*---------------------------------------------------------------------------*/

Selecting physicsModel solid

Selecting solidModel coupledUnsLinearGeometryLinearElastic
Selecting dynamicFvMesh staticFvMesh
Patch tip    Traction boundary field 
Patch top    Traction boundary field 
Patch bot    Traction boundary field 
Patch right    Traction boundary field 
Patch left    Traction boundary field 
    under-relaxation method: fixed

/*---------------------------------------------------------------------------*\
|    For further information on this class, please see the following          |
|    publication:                                                             |
|                                                                             |
|    P. Cardiff, Z. Tukovic, H, Jasak, A. Ivankovic: A block-coupled finite   |
|    volume methodology for linear elasticity and unstructured meshes.        |
|    Computers and Structures, 2016, 175 100-122,                             |
|    DOI: 10.1016/-j.compstruc.2016.07.004.                                   |
\*---------------------------------------------------------------------------*/

Creating the mechanicalModel
Selecting mechanical law linearElastic
Time = 1

Evolving solid solver
Solving the momentum equation for D
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 1.5625e-06 -1.5625e-06 4.6875e-06 1.5625e-06 1.5625e-06 1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 1.5625e-06 -1.5625e-06 4.6875e-06 1.5625e-06 1.5625e-06 1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
--> FOAM Warning : 
    From function eigenValues(const tensor&)
    in file primitives/Tensor/tensor/tensor.C at line 180
    complex eigenvalues detected for tensor: (4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06 -1.5625e-06 -1.5625e-06 -1.5625e-06 4.6875e-06)
EigenSolver : writeMatlabFiles is false
EigenSparseLU Eigen linear solver

The results of the benchmark I can access until now are available in the attachments.

Greetings!