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 acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Browse other questions tagged or ask your own question.