[bvieira]

Hello all!

I put together a solver combining rhoPimpleFoam and PimpleDyMFoam to be able to simulate a pitching airfoil by rotating the entire grid.
After being able to compile everything and run a case with no errors, I’m now facing a convergence problem.

I get an error of number of iterations exceeded:

 --> FOAM FATAL ERROR:
[0] Maximum number of iterations exceeded
[0]
[0]     From function specieThermo<Thermo>::T(scalar f, scalar T0, scalar (specieThermo<Thermo>::*F)(const scalar) const, scalar (specieThermo<Thermo>::*dFdT)(const scalar) const) const

I think my temperature is diverging and then dragging density with it.
My initial conditions are taken from a converged steady-state solution using rhoSimpleFoam.
I tried to play with the under-relaxation factors, as well as the number of inner correctors (nCorrectors), but that doesn’t seem to be helping much (runs longer, but still blows up).
Check mesh passes. This grid has also been used in Fluent before.
Also, I’m running with a pretty small timestep, with MaxCo being less than 1.5. I'm currently using 10 subiterations (nOuterCorrectors).

On my solver side, I think everything seems to be working. I have inspected the mesh during rotation, and it seems to be rotating as desired.
Just as a check, I run the same case, without pitching on rhoPimpleFoam, but still transient, and I was getting a similar problem.

Any ideas on what could be going wrong? Any suggestion on more stable schemes perhaps ? Or changes in the relaxations factors?

I'm attaching the fvSolution and the fvSchemes files.
Some of the underrelaxation factors that I've tried are put as comments in the fvSolution file.

I appreciate any inputs.

[karkar]

Hi bvieira,

did you manage to track down your problem?

Could you share your rhoPimpleDyMFoam solver? I currently need that as well and I guess I would be wasting time reinventing a wheel you already made round...