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?
$\begingroup$
$\endgroup$
2
-
4$\begingroup$ How are you stabilizing the equations in the high-Pe case? $\endgroup$– Wolfgang BangerthCommented 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$– EngDRCommented Jan 24, 2018 at 11:29
Add a comment
|