Analysis of linearized Galerkin-mixed FEMs for the time-dependent Ginzburg-Landau equations of superconductivity

A linearized backward Euler Galerkin-mixed finite element method is investigated for the time-dependent Ginzburg-Landau (TDGL) equations under the Lorentz gauge. By introducing the induced magnetic field σ = curl A as a new variable, the Galerkin-mixed FE scheme offers many advantages over conventional Lagrange type Galerkin FEMs. An optimal error estimate for the linearized Galerkin-mixed FE scheme is established unconditionally. Analysis is given under more general assumptions for the regularity of the solution of the TDGL equations, which includes the problem in two-dimensional nonconvex polygons and certain three dimensional polyhedrons, while the conventional Galerkin FEMs may not converge to a true solution in these cases. Numerical examples in both two and three dimensional spaces are presented to confirm our theoretical analysis. Numerical results show clearly the efficiency of the mixed method, particularly for problems on nonconvex domains.