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$ – Wolfgang Bangerth Jan 17 '18 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 '18 at 11:29

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.