TY - JOUR
T1 - Finite Element Methods for a Bi-Wave Equation Modeling D-Wave Superconductors
AU - Xiaobing Feng & Michael Neilan
JO - Journal of Computational Mathematics
VL - 3
SP - 331
EP - 353
PY - 2010
DA - 2010/06
SN - 28
DO - http://doi.org/10.4208/jcm.1001-m1001
UR - https://global-sci.org/intro/article_detail/jcm/8523.html
KW - Bi-wave operator, d-wave superconductors, Conforming finite elements, Error
estimates.
AB - <p style="text-align: justify;">In this paper we develop two conforming finite element methods for a fourth order bi-wave equation arising as a simplified Ginzburg-Landau-type model for $d$-wave superconductors in absence of applied magnetic field. Unlike the biharmonic operator $\Delta^2$, the bi-wave operator $\Box^2$ is not an elliptic operator, so the energy space for the bi-wave equation is much larger than the energy space for the biharmonic equation. This then makes it possible to construct low order conforming finite elements for the bi-wave equation. However, the existence and construction of such finite elements strongly depends on the mesh. In the paper, we first characterize mesh conditions which allow and not allow construction of low order conforming finite elements for approximating the bi-wave equation. We then construct a cubic and a quartic conforming finite element. It is proved that both elements have the desired approximation properties, and give optimal order error estimates in the energy norm, suboptimal (and optimal in some cases) order error estimates in the $H^1$ and $L^2$ norm. Finally, numerical experiments are presented to guage the efficiency of the proposed finite element methods and to validate the theoretical error bounds.</p>