Approximating three-dimensional magnetohydrodynamics system forced by space-time white noise

The magnetohydrodynamics system consists of the Navier-Stokes and Maxwell's equations, coupled through multiples of nonlinear terms. Such a system forced by space-time white noise has been studied by physicists for decades, and the rigorous proof of its solution theory has been recently established in Yamazaki (2019, arXiv:1910.04820 [math.AP]) using the theory of paracontrolled distributions and a technique of coupled renormalizations. When an equation is well-posed, and it is approximated by replacing the differentiation operator by reasonable discretization schemes with a parameter, it is widely believed that a solution of the approximating equation should converge to the solution of the original equation as the parameter approaches zero. We prove otherwise in the case of the three-dimensional magnetohydrodynamics system forced by space-time white noise. Specifically, it is proven that the limit of the solution to the approximating system with an additional 32 drift terms solves the original system. These 32 drift terms depend on the choice of approximations, can be calculated explicitly in the process of renormalizations, and essentially represent a spatial version of It$\hat{\mathrm{o}}$-Stratonovich correction terms. In particular, the proof relies on the technique of coupled renormalizations again, as well as taking advantage of the special structure of the magnetohydrodynamics system on many occasions.