ИСТИНА

In the report, we consider the Finite Element Method (FEM) solution of non-linear elliptic equation with Newman boundary conditions and with strong discontinuity in the solution. We apply dumped Newton method to resolve the nonlinearity and evaluate mixed FEM to resolve discontinuity on internal interface. The computational algorithm is implemented using package FreeFEM++. This problem is a part of mathematical model of Abricosov vortex in bi-layer superconductor structure. The elliptic equation of this model we solve is Usadel equation. Some results of simulation of bi-layer structure with vortex are presented.