[luitzor]

Hi Foamers!

I need some help with a propeller simulation (transient). As a quick note, I'm quite new to OpenFoam. A brief description of the project:

I created the propeller geometry with Solidworks and then created the enclosure and the mesh with Salome. There are two regions: the main enclosure (that doesn't rotate) which represents the freestream flow and the propeller enclosure (which should rotate) that contains the propeller. The meshes of both domains are created separately and then merged using mergeMeshes. The rotor cellZone region is then set using setSet. I also use topoSetDict to add all AMI faces to a faceSet. Finally, due to the size of the simulation, I run it in parallel using 20 procs (using the singleProcessorFaceSets to keep all AMI faces on the same processor). The simulation blows after a few iterations. I must note that the simulation blows when it's run in series too.

I've checked the mesh and it's okay (see attached log file).

This is the log of a failed parallel simulation.

| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : pimpleDyMFoam -parallel
Date   : Apr 19 2021
Time   : 15:56:34
Host   : "93a46648863c"
PID    : 165
Case   : /home/openfoam
nProcs : 20
Slaves : 
19
(
"93a46648863c.166"
"93a46648863c.167"
"93a46648863c.168"
"93a46648863c.169"
"93a46648863c.170"
"93a46648863c.171"
"93a46648863c.172"
"93a46648863c.173"
"93a46648863c.174"
"93a46648863c.175"
"93a46648863c.176"
"93a46648863c.177"
"93a46648863c.178"
"93a46648863c.179"
"93a46648863c.180"
"93a46648863c.181"
"93a46648863c.182"
"93a46648863c.183"
"93a46648863c.184"
)

Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function rotatingMotion
Applying solid body motion to cellZone rotor

PIMPLE: no residual control data found. Calculations will employ 2 corrector loops

Reading field p

Reading field U

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 165774 source faces and 170168 target faces
AMI: Patch source sum(weights) min/max/average = 0.998918, 1.0014, 1
AMI: Patch target sum(weights) min/max/average = 0.998764, 1.00002, 0.999999
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

No MRF models present

Reading/calculating face velocity Uf

No finite volume options present

Courant Number mean: 3.12403e-06 max: 0.000570603

Starting time loop

Courant Number mean: 3.0931e-06 max: 0.000564953
deltaT = 1.17647e-05
Time = 1.17647e-05

AMI: Creating addressing and weights between 165774 source faces and 170168 target faces
AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1
AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 0.00635745, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 0.369898, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 0.36776, No Iterations 1000
GAMG:  Solving for p, Initial residual = 1, Final residual = 0.0071022, No Iterations 11
time step continuity errors : sum local = 5.33589e-08, global = -3.80218e-08, cumulative = -3.80218e-08
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.47015, Final residual = 0.0836701, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.265834, Final residual = 0.0436782, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 0.10333, Final residual = 0.00688327, No Iterations 1000
GAMG:  Solving for p, Initial residual = 0.0941766, Final residual = 8.35714e-07, No Iterations 22
time step continuity errors : sum local = 5.98351e-09, global = 3.17688e-09, cumulative = -3.48449e-08
smoothSolver:  Solving for epsilon, Initial residual = 0.994799, Final residual = 0.00754865, No Iterations 1000
smoothSolver:  Solving for k, Initial residual = 1, Final residual = 2.46881e-06, No Iterations 1000
ExecutionTime = 857.17 s  ClockTime = 866 s

Courant Number mean: 0.00943062 max: 83392.3
deltaT = 2.82153e-10
Time = 1.1765e-05

AMI: Creating addressing and weights between 165774 source faces and 170168 target faces
AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1
AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.7144, Final residual = 0.0433931, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.775929, Final residual = 0.0502612, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.82386, Final residual = 0.0560173, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.910335, Final residual = 0.00880902, No Iterations 4
time step continuity errors : sum local = 3.13285e-10, global = 2.682e-12, cumulative = -3.48422e-08
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.530458, Final residual = 0.472351, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.437057, Final residual = 0.265858, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 0.371612, Final residual = 0.0227143, No Iterations 1000
GAMG:  Solving for p, Initial residual = 1.58002e-08, Final residual = 1.58002e-08, No Iterations 0
time step continuity errors : sum local = 6.73282e-06, global = 3.84411e-12, cumulative = -3.48384e-08
smoothSolver:  Solving for epsilon, Initial residual = 0.192127, Final residual = 7.48896e-07, No Iterations 4
smoothSolver:  Solving for k, Initial residual = 0.827074, Final residual = 0.106168, No Iterations 1000
ExecutionTime = 1320.87 s  ClockTime = 1331 s

Courant Number mean: 0.000144388 max: 1942.94
deltaT = 2.9044e-13
Time = 1.1765e-05

AMI: Creating addressing and weights between 165774 source faces and 170168 target faces
AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1
AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.768038, Final residual = 0.0766915, No Iterations 118
smoothSolver:  Solving for Uy, Initial residual = 0.81215, Final residual = 0.0808922, No Iterations 127
smoothSolver:  Solving for Uz, Initial residual = 0.902063, Final residual = 0.0894033, No Iterations 22
GAMG:  Solving for p, Initial residual = 0.880168, Final residual = 0.00271608, No Iterations 2
time step continuity errors : sum local = 0.000258367, global = 4.23745e-07, cumulative = 3.88907e-07
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.580994, Final residual = 0.0598996, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.720866, Final residual = 0.1011, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 0.615497, Final residual = 0.137957, No Iterations 1000
GAMG:  Solving for p, Initial residual = 0.000621124, Final residual = 9.73629e-07, No Iterations 22
time step continuity errors : sum local = 4.05813e-06, global = -1.78319e-07, cumulative = 2.10588e-07
smoothSolver:  Solving for epsilon, Initial residual = 1, Final residual = 0.106181, No Iterations 1000
smoothSolver:  Solving for k, Initial residual = 0.999862, Final residual = 0.154266, No Iterations 1000
ExecutionTime = 1944.37 s  ClockTime = 1959 s

Courant Number mean: 0.00645368 max: 6.77553e+07
deltaT = 8.5732e-21
Time = 1.1765e-05

AMI: Creating addressing and weights between 165774 source faces and 170168 target faces
AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1
AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.739661, Final residual = 0.0312613, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.873737, Final residual = 0.05126, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.800713, Final residual = 0.0290762, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.917227, Final residual = 0.00712336, No Iterations 5
time step continuity errors : sum local = 6.51386e-13, global = -6.0799e-15, cumulative = 2.10588e-07
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.58764, Final residual = 0.390951, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.510518, Final residual = 0.0398422, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 0.612049, Final residual = 0.0284529, No Iterations 1000
GAMG:  Solving for p, Initial residual = 1.86722e-09, Final residual = 1.86722e-09, No Iterations 0
time step continuity errors : sum local = 8.88265e-09, global = -1.68866e-15, cumulative = 2.10588e-07
smoothSolver:  Solving for epsilon, Initial residual = 0.203351, Final residual = 1.17043e-07, No Iterations 1
bounding epsilon, min: -2.63513e+29 max: 3.13976e+42 average: 2.11151e+37
smoothSolver:  Solving for k, Initial residual = 0.996435, Final residual = 0.967348, No Iterations 1000
ExecutionTime = 2438.95 s  ClockTime = 2456 s

Courant Number mean: 3.76917e-08 max: 706.294
deltaT = 2.42766e-23
Time = 1.1765e-05

AMI: Creating addressing and weights between 165774 source faces and 170168 target faces
AMI: Patch source sum(weights) min/max/average = 0.998859, 1.00151, 1
AMI: Patch target sum(weights) min/max/average = 0.998615, 1.00002, 0.999998
PIMPLE: iteration 1
smoothSolver:  Solving for Ux, Initial residual = 0.927617, Final residual = 0.0460816, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.484562, Final residual = 0.028569, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.983595, Final residual = 0.0748445, No Iterations 3
GAMG:  Solving for p, Initial residual = 0.796851, Final residual = 0.00460583, No Iterations 1
time step continuity errors : sum local = 1.00549e-06, global = -1.28245e-10, cumulative = 2.1046e-07
PIMPLE: iteration 2
smoothSolver:  Solving for Ux, Initial residual = 0.473033, Final residual = 0.458308, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 0.609904, Final residual = 0.380841, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 0.565501, Final residual = 0.00530204, No Iterations 1000
GAMG:  Solving for p, Initial residual = 1.06618e-10, Final residual = 1.06618e-10, No Iterations 0
time step continuity errors : sum local = 1.57834e-05, global = -1.45424e-10, cumulative = 2.10314e-07
[14] #0  Foam::error::printStack(Foam::Ostream&)[9] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
[14] #1  Foam::sigFpe::sigHandler(int) at ??:?
[9] #1  Foam::sigFpe::sigHandler(int) at ??:?
[9] #2  ? at ??:?
[14] #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #3  Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in "/lib/x86_64-linux-gnu/libc.so.6"
[14] #3  Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
[9] #4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[9] #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[14] #4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
[9] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[14] #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[14] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[9] #7   at ??:?
[14] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&)Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
[9] #8   at ??:?
[14] #8  Foam::fvMatrix<double>::solve()Foam::fvMatrix<double>::solve() at ??:?
[9] #9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
[14] #9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
[9] #10  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
[14] #10  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
[9] #11  ? at ??:?
[14] #11  ? at ??:?
[9] #12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[9] #13   at ??:?
[14] #12  __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6"
[14] #13   at ??:?
[93a46648863c:00174] *** Process received signal ***
[93a46648863c:00174] Signal: Floating point exception (8)
[93a46648863c:00174] Signal code:  (-6)
[93a46648863c:00174] Failing at address: 0x3e9000000ae
[93a46648863c:00174] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f3c98ee14b0]
[93a46648863c:00174] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7f3c98ee1428]
[93a46648863c:00174] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7f3c98ee14b0]
[93a46648863c:00174] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x744)[0x7f3c9a134634]
[93a46648863c:00174] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x28)[0x7f3c9a134888]
[93a46648863c:00174] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4f6)[0x7f3c9a12b8f6]
[93a46648863c:00174] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7f3c9ceef95b]
[93a46648863c:00174] [ 7] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x470f41]
[93a46648863c:00174] [ 8] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xd4)[0x471194]
[93a46648863c:00174] [ 9] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solveIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMatrixIS2_EEEE+0x2f)[0x7f3c9e601f3f]
[93a46648863c:00174] [10] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x9a9)[0x7f3c9e66e099]
[93a46648863c:00174] [11] pimpleDyMFoam[0x425849]
[93a46648863c:00174] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7f3c98ecc830]
[93a46648863c:00174] [13] pimpleDyMFoam[0x4286c9]
[93a46648863c:00174] *** End of error message ***
? at ??:?
[93a46648863c:00179] *** Process received signal ***
[93a46648863c:00179] Signal: Floating point exception (8)
[93a46648863c:00179] Signal code:  (-6)
[93a46648863c:00179] Failing at address: 0x3e9000000b3
[93a46648863c:00179] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7ff09b5174b0]
[93a46648863c:00179] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x38)[0x7ff09b517428]
[93a46648863c:00179] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x354b0)[0x7ff09b5174b0]
[93a46648863c:00179] [ 3] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam22symGaussSeidelSmoother6smoothERKNS_4wordERNS_5FieldIdEERKNS_9lduMatrixERKS5_RKNS_10FieldFieldIS4_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEEhi+0x744)[0x7ff09c76a634]
[93a46648863c:00179] [ 4] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam22symGaussSeidelSmoother6smoothERNS_5FieldIdEERKS2_hi+0x28)[0x7ff09c76a888]
[93a46648863c:00179] [ 5] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam12smoothSolver5solveERNS_5FieldIdEERKS2_h+0x4f6)[0x7ff09c7618f6]
[93a46648863c:00179] [ 6] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x15b)[0x7ff09f52595b]
[93a46648863c:00179] [ 7] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x191)[0x470f41]
[93a46648863c:00179] [ 8] pimpleDyMFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xd4)[0x471194]
[93a46648863c:00179] [ 9] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam5solveIdEENS_17SolverPerformanceIT_EERKNS_3tmpINS_8fvMatrixIS2_EEEE+0x2f)[0x7ff0a0c37f3f]
[93a46648863c:00179] [10] /opt/openfoam4/platforms/linux64GccDPInt32Opt/lib/libincompressibleTurbulenceModels.so(_ZN4Foam9RASModels8kEpsilonINS_29IncompressibleTurbulenceModelINS_14transportModelEEEE7correctEv+0x9a9)[0x7ff0a0ca4099]
[93a46648863c:00179] [11] pimpleDyMFoam[0x425849]
[93a46648863c:00179] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xf0)[0x7ff09b502830]
[93a46648863c:00179] [13] pimpleDyMFoam[0x4286c9]
[93a46648863c:00179] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 14 with PID 179 on node 93a46648863c exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------
OpenFOAM-4.1(31)

I've done some research, but I cannot find the exact problem I have. It could have something to do with the boundary conditions or the mesh, but I've literally run out of ideas. The boundary conditions are quite standard (I've attached the P and U boundary conditions. The other boundary conditions files can be found using the link below).

Group_1 is the outlet and Wall_noSlip is the wall of the outer domain (the name is misleading. It's actually a slip wall condition). I'm aware that this names are confusing and I'll fix it now.

This is the link to the project files (I deleted /constant/faces to reduce it's size).

https://drive.google.com/drive/folde...bE?usp=sharing

Any help would be appreciated!