[stardust23]

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:

 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 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