I a solving the 2d steady state convection-diffusion problem on the famous flow around a cylinder in rectangular domain benchmark. My numerical results show that with Q1 finite elements, the solver is converging for Pe=5000, whereas with Q2 discretization, the solver fails to converge. What can be the reason for poor performance in case of Q2 FEM discretization?

  • 4
    $\begingroup$ How are you stabilizing the equations in the high-Pe case? $\endgroup$ Jan 17, 2018 at 3:47
  • $\begingroup$ @ Wolfgang Bangerth I am using edge-oriented jump stabilization, the solver is now converging for the Q2 case, but still, performance of solver is not as good as that for Q1 case $\endgroup$
    – EngDR
    Jan 24, 2018 at 11:29


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Browse other questions tagged or ask your own question.