I continued the solution a little more, the residuals fluctuate in numbers between 1e-02 and 1e-03 and repeatedly after each iteration the message (Stabilizing pressure coupled to enhance linear solver robustness) appears on the console and it is quite clear that the time between both iterations with the message Coupled Pressure, compared to Two iterations without the coupled pusher message is more.
I want to make the grid so small that the CL and CD values are almost constant. Does the increase in the number of cells in the grid lead to the fluctuation of the residuals?
According to similar studies, I think my geometry should have different results up to about 3.5 million cells.
I must do the mesh independence to ensure the correctness of the results.