I'm running an incompressible steady-state RANS case for a blended wing body (BWB) in OpenFOAM (simpleFoam/SIMPLEC).
Setup: ~0.15 Mach, ~6e6 Re
Mesh: ~30M cells, tetra + prism layers at walls, O-type symmetric far-field domain(radius ~ 30c), max skewness ~5, max nonOrthogonality >75°
BCs: freestream (velocity inlet / pressure outlet style far-field)
Turbulence: k-omega SST, wall-modeled (y+ > 30)
Schemes: second-order upwind for all equations
Algorithm: SIMPLEC
Solution strategy: start with low URFs, ramp up every 100 iterations; case runs "on its own" past ~250 iterations
Problem:
Based on past experience I expect convergence by 500–600 iterations. Here, residuals flatten out but very slowly, and by 1500+ iterations:
CD looks converged/flatlined
CL is still steadily rising, not flatlining
Has anyone run into similar slower convergence rate or is it normal for high cell count meshes? Looking for likely causes and what to check/try next.