DESIGN AND CONVERGENCE ANALYSIS OF NUMERICAL METHODS FOR STOCHASTIC EVOLUTION EQUATIONS WITH LERAY–LIONS OPERATOR JER´ OMEˆ DRONIOU, BENIAMIN GOLDYS, AND KIM-NGAN LE Abstract. The gradient discretisation method (GDM) is a generic frame- work, covering many classical methods (Finite Elements, Finite Volumes, Dis- continuous Galerkin, etc.), for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for a general stochastic evo- lution problem based on a Leray–Lions type operator. The problem contains the stochastic p-Laplace equation as a particular case. The convergence of the Gradient Scheme (GS) solutions is proved by using Discrete Functional Anal- ysis techniques, Skorohod theorem and the Kolmogorov test. In particular, we provide an independent proof of the existence of weak martingale solutions for the problem. In this way, we lay foundations and provide techniques for proving convergence of the GS approximating stochastic partial differential equations. Keywords: p-Laplace equation, stochastic PDE, numerical methods, gradient discretisa- tion method, convergence analysis 1. Introduction The parabolic p-Laplacian problem occurs in many mathematical models of physical processes, such as nonlinear diffusion [1] and non-Newtonian flows [49]. However, in practical situations with large scales, rapid velocity and pressure fluctuations, the motion of flow becomes unsteady and it is described as being turbulent [51]. Turbulence is a combination of a slow oscillating (deterministic) component and a fast oscillating component that can be modelled as a white noise perturbation of regular fluid velocity field. Therefore, in order to investigate turbulence in the parabolic p-Laplacian problem, the first step is to develop the theory and numeri- cal algorithms for the stochastic parabolic p-Laplacian problem. Motivated by this arXiv:2004.12047v2 [math.NA] 29 Dec 2020 problem, we study in this paper a more general stochastic partial differential equa- tion based on a Leray–Lions type operator with homogeneous Dirichlet boundary condition. The model reads du ´ divpapu, ∇uqqdt “ fpuqdWt in p0,T qˆ Θ, up0, ¨q “ u0 in Θ, (1) u “ 0 on p0,T qˆBΘ, where T ą 0, Θ is an open bounded domain in Rd, d “ 1, 2, 3, and the initial 2 data u0 P L pΘq. Here, f is a continuous operator with linear growth acting between appropriate Banach spaces, see Section 2 for details. We assume that W “ tW ptq,t P r0,T su is a K-valued Wiener process with a trace class covariance Date: January 1, 2021. 1 2 JER´ OMEˆ DRONIOU, BENIAMIN GOLDYS, AND KIM-NGAN LE operator Q, for a certain Hilbert space K. Particular choices of a include the p- Laplace operator corresponding to apu, vq “ |v|p´1v for some p P p1, `8q (see also [45] for more general versions), nonlinear and nonlocal diffusion operators of the p´1 p 8 form apu, vq “ Λrus|v| v with Λrus : L pΘq Ñ L pΘ; SdpRqq uniformly elliptic (here, SdpRq denotes the space of symmetric dˆd matrices) – such a model appears in image smoothing with p “ 2 and Λrus “ gp|∇G ˚ u|q, where G is a Gaussian kernel, ˚ is the spatial convolution, and g : R` Ñ R` is smooth and decreasing [17]. If f “ 0 then the noise term vanishes, hence our stochastic model (1) includes deterministic equation as a special case. Some existence and uniqueness results for some particular forms of (1) can be found in the literature. In [20, 38, 39], a quasi-linear version is considered in which apu, ∇uq “ Apuq∇u, and an additional advective term divpBpuqq is added to the model; existence and uniqueness of suitable solutions are proved. [10] consid- ers a non-degenerate version of the p-Laplace model, in which apu, ∇uq “ p1 ` |∇u|qp´2∇u, and proves existence and regularity results. The analysis carried out in [50] only covers the straight p-Laplacian, and is restricted to p ě 2. Our as- sumption on the Leray–Lions operator are more general than in these references, in the sense that we accept models that are fully non-linear with respect to ∇u, and that may be non-monotone (a depending on both u and ∇u, in a non-linear way with respect to ∇u). Moreover, and contrary to these references, we propose an approach that has the double advantage of establishing the existence of a solution to (1), and of proving the convergence of a variety of numerical approximations of this model. Numerical methods of the deterministic version of model (1) (i.e. f “ 0) and their proofs of convergence are studied in [7,16,27,42] and the references cited therein. However, there is no numerical approximation of the stochastic model (1) due to difficulties arising in the nonlinear term and the infinite dimensional nature of the driving noise processes. There is an increasing number of numerical methods for the solution of stochas- tic evolution equations mentioned in the literature [43, 44, 53], where unique mild solutions are required and the approximate schemes are treated in terms of the semigroup approach. However, these assumptions are not applicable for a class of stochastic equations involving strongly nonlinear terms, such as Navier–Stokes, magnetohydrodynamics (MHD), Schr¨odinger, Landau–Lifshitz–Gilbert, Landau– Lifshitz–Bloch and nonlinear porous media equations. The stochastic Navier-Stokes equation [12, 14] and the stochastic Landau–Lifshitz–Gilbert equation [4, 5, 35, 36] are investigated by using the conforming finite element method to approximate their solutions. Furthermore, the convergence of the approximate solutions is also proved which implies the existence of weak martingale solutions. In the recent work [47], a general convergence theory for conforming finite element schemes of stochastic parabolic PDEs is developed by adapting ideas from [4, 5, 12, 14]. All these previous works, however, only deal with conforming approximations, which use for the spatial discretisation a subspace of the Sobolev space appearing in the weak formulation of the continuous problem. This usually imposes restrictions on the types of mesh that can be considered – typically, triangular/tetrahedral or quadrangular/hexahedral meshes. Moreover, conforming methods are known to be ill-suited in some applications, e.g. when mesh locking appears, when inf-sup stabil- ity is sought, or when some physical properties of the model must be respected (such NUMERICAL METHODS FOR STOCHASTIC PDES WITH LERAY–LIONS OPERATOR 3 as balance and conservativity of approximate fluxes). In such circumstances, non- conforming methods might be better suited; such methods include non-conforming finite elements and finite volume methods, and also recent high-order methods for polytopal meshes with cell and face unknowns – such as Hybrid-High Order schemes and Virtual Element Methods. We refer the reader to [2,8,18,21, 23] and reference therein for detailed presentations of these methods. In this work, we approximate (1) by using the Gradient Discretisation Method (GDM) [26] and an implicit Euler time stepping with uniform time steps. The GDM is a generic convergence analysis framework for a wide variety of methods (conforming or nonconforming) written in discrete variational formulation, and based on independent approximations of functions and gradients using the same degrees of freedom. Several well-known methods fall in the GDM framework, in particular: ‚ Galerkin methods, including the (standard or mass-lumped) conforming Finite Element methods [29], ‚ Nonconforming Finite Element methods, including the (standard or mass- lumped) nonconforming P1 scheme [29] and non-conforming Finite Element methods on polytopal meshes [28], ‚ Symmetric Interior Penalty Galerkin (SIPG) methods [32], ‚ Mixed Finite Element methods [26], ‚ Hybrid Mimetic Mixed methods and Mimetic Finite Difference methods [27], ‚ Hybrid High-Order and Virtual Elements Methods [22]. By writing numerical schemes for (1) and performing their analysis in the GDM framework, we provide a unified convergence result for all these methods. We refer to [22, 33] and to the complete monograph [26] for more details of the GDM and the methods it covers. Because the GDM encompasses non-conforming schemes, the functional spaces for the approximate solutions are not included in the classical (continuous) Sobolev spaces. Therefore, the usual Poincar´einequalities, Sobolev embeddings, Rellich or Aubin–Simon compactness theorems, or trace inequalities cannot be used. In the context of deterministic PDEs, a series of “Discrete Func- tional Analysis” results have been established to mimic these continuous functional analysis tools [26]. Our convergence analysis approach is based on the adaptation of these Discrete Functional Analysis techniques to the stochastic case, the Skorohod theorem and the Kolmogorov test; we show the convergence of the Gradient Scheme (GS) solutions to a weak martingale solution of (1). In this way, an independent proof of the existence of weak martingale solutions for the problem is provided. The paper is organised as follows. In Section 2 we recall the notations of the gradient discretisation method and propose the GS for approximating the stochastic model (1). Weak martingale solutions to (1) are defined and our main result is stated in this section. Section 3 provides priori estimates of approximated solutions and the noise term added at each step of the scheme in various norms. In Section 4, we first show the tightness of the sequence including the GS solutions and then prove the almost sure convergence in a certain norm, up to a change of probability space. The continuity of the limit and the martingale part are also proved in this section. Section 5 is devoted to the proof of the main theorem. Finally, in the Appendix we prove necessary results that are used in the course of the proof. 4 JER´ OMEˆ DRONIOU, BENIAMIN GOLDYS, AND KIM-NGAN LE 2. Gradient scheme and main results Before introducing the GS for approximation of (1), we introduce notations and assumptions used in the rest of the paper.