CFD Online Discussion Forums - Pump impeller steady simulation crashing on cyclicAMI boundaries

Dear Foamers,
I’m trying to run the steady-state case of a pump impeller with SRFSimpleFoam (OpenFoam_v7). To reduce the computational cost, I’d like to simulate only one blade with cyclicAMI patches on the periodicity boundaries.
Long story short, the simulation crashes due to these cyclic interfaces. I tried to modify the mesh and it seems that the problem is related to the thin prismatic cells used for hub and tip BLs on the AMI patches (removing these cells the simulation starts normally).
I’m sharing the complete case that is not working hoping that someone could help checking the setup and possibly finding a way to solve the problem: https://drive.google.com/drive/folde..._B?usp=sharing

This is the solver output that I get when running SRFSimpleFoam:

Code:

 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

 Create time

 

 

 Create mesh for time = 0

 

 

 

 

 SIMPLE: No convergence criteria found

 

 

 Reading field p

 

 

 AMI: Creating addressing and weights between 13169 source faces and 13291 target faces

 AMI: Patch source sum(weights) min/max/average = 0, 1.01349369, 0.999697566

 AMI: Patch target sum(weights) min/max/average = 0, 1.01662743, 0.999685469

 Reading field Urel

 

 

 Reading/calculating face flux field phi

 

 

 Creating SRF model

 

 

 Selecting SRFModel rpm

 Selecting incompressible transport model Newtonian

 Selecting turbulence model type RAS

 Selecting RAS turbulence model kOmegaSST

 Selecting patchDistMethod meshWave

 bounding k, min: 0 max: 0.006935 average: 0.006935

 bounding omega, min: 0 max: 83.28 average: 83.28

 RAS

 {

     RASModel        kOmegaSST;

     turbulence      on;

     printCoeffs     on;

     alphaK1         0.85;

     alphaK2         1;

     alphaOmega1     0.5;

     alphaOmega2     0.856;

     gamma1          0.555555556;

     gamma2          0.44;

     beta1           0.075;

     beta2           0.0828;

     betaStar        0.09;

     a1              0.31;

     b1              1;

     c1              10;

     F3              false;

 }

 

 

 Creating finite volume options from "system/fvOptions"

 

 

 Selecting finite volume options model type limitVelocity

     Source: limitU

     - selecting all cells

     - selected 1202804 cell(s) with volume 0.000712983466

 

 

 Starting time loop

 

 

 Time = 1

 

 

 DILUPBiCGStab:  Solving for Urelx, Initial residual = 1, Final residual = 4.28508452e-08, No Iterations 5

 DILUPBiCGStab:  Solving for Urely, Initial residual = 1, Final residual = 7.93255556e-08, No Iterations 4

 DILUPBiCGStab:  Solving for Urelz, Initial residual = 1, Final residual = 3.60965556e-08, No Iterations 4

 #0  Foam::error::printStack(Foam::Ostream&) at ??:?

 #1  Foam::sigFpe::sigHandler(int) at ??:?

 #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"

 #3  Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?

 #4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/user/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/SRFSimpleFoam"

 #5  ? in "/home/user/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/SRFSimpleFoam"

 #6  ? in "/home/user/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/SRFSimpleFoam"

 #7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

 #8  ? in "/home/user/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt32Opt/bin/SRFSimpleFoam"

 Floating point exception