OPTIMAL ANALYSIS OF NON-UNIFORM GALERKIN-MIXED FINITE ELEMENT APPROXIMATIONS TO THE GINZBURG-LANDAU EQUATIONS IN SUPERCONDUCTIVITY.

This paper is concerned with new error analysis of a lowest-order backward Euler Galerkin-mixed finite element method for the time-dependent Ginzburg-Landau equations. The method is based on a commonly-used nonuniform approximation, in which a linear Lagrange element, the lowest-order Nédélec edge element, and the Raviart-Thomas face element are used for the order parameter ψ, the magnetic field curlA, and the magnetic potential A, respectively. This mixed method has been widely used in practical simulations due to its low cost and ease of implementation. In the Ginzburg-Landau model, the order parameter ψ is the most important variable, which indicates the state of the superconductor. An important feature of the method is the inconsistency of the approximation orders. A crucial question is how the first-order approximation of (curlA, A) influences the accuracy of ψ_h. The main purpose of this paper is to establish the second-order accuracy for the order parameter in a spatial direction, although the accuracy for (curlA, A) is first order only. Previous analysis only gave the first-order convergence for all three variables due to certain artificial pollution involved in the analysis. Our analysis is based on a nonstandard quasi-projection for ψ and the corresponding more precise estimates, including the H^-1 -norm. With the quasi-projection, we prove that the lower-order approximation to (curlA, A) does not pollute the accuracy of ψ_h. Our numerical experiments confirm the optimal convergence of ψ_h. The approach can be extended to many other multiphysics models.. [ABSTRACT FROM AUTHOR]