2D CurlCurl

Hello everybody,

I still try to use a error estimator as presented by Bürg for an hp-refinement method for a Maxwell problem.
Therefore, I need to compute the term
[tex]
$$\begin{array}{r}\Vert \mathrm{\nabla }\times ({\mu }^{-1}\mathrm{\nabla }\times u_h)-\epsilon {\omega }^2{u}_h{\Vert }_{L^2(K)}^2\end{array}$$
[/tex]
on every cell K. I tried to set a new HCurl gridfunction for the curl of u_h, but since curl(u_h) returns the scalar curl
[tex]
$$\begin{array}{r}\mathrm{curl}(u)=(\frac{\partial u_2}{\partial x_1}-\frac{\partial u_1}{\partial x_2})\end{array}$$
[/tex]
this does not work (different dimensions).

Has anybody a suggestion how fix this issue?

Best regards,
Phil

Hi Phil,

you can set the curl, which is a scalar, into a L2 GridFunction. Then, you can use the gradient of this GridFunction to compute the curl going from scalar to vector to compute your norm.

Best,
Michael