Hi,
I’m trying to compute the order of convergence for a fourth-order variational inequality using the penalty method, in FreeFEM++. But I got unexpected results.
I computed the L2 and the H1 errors against an exact solution, but the convergence rate seems to be inconsistent with the theoretical results. Note that the exact solution has H3 global regularity.
Rate-VI-4.edp (4.1 KB)
Attached here is the code that I have used.
Best regards.