The tensor-train mimetic finite difference method for three-dimensional Maxwell’s wave propagation equations

Abstract

Coupling the mimetic finite difference method with the tensor-train format results in a very effective method for low-rank numerical approximations of the solutions of the time-dependent Maxwell wave propagation equations in three dimensions. To this end, we discretize the curl operators on the primal/dual tensor product grid complex and we couple the space discretization with a staggered-in-time second-order accurate time-marching scheme. The resulting solver is accurate to the second order in time and space, and still compatible, so that the approximation of the magnetic flux field has zero discrete divergence with a discrepancy close to the machine precision level. Our approach is not limited to the second-order of accuracy. We can devise higher-order formulations in space through suitable extensions of the tensor-train stencil to compute the derivatives of the mimetic differential operators. Employing the tensor-train format improves the solver performance by orders of magnitude in terms of CPU time and memory storage. A final set of numerical experiments confirms this expectation.