• Edge-based nonlinear diffusion for finite element approximations of convection–diffusion equations and its relation to algebraic flux-correction schemes

      Barrenechea, Gabriel; Burman, Erik; Karakatsani, Fotini; University of Strathclyde; University College London; University of Chester (Springer, 2016-05-07)
      For the case of approximation of convection–diffusion equations using piecewise affine continuous finite elements a new edge-based nonlinear diffusion operator is proposed that makes the scheme satisfy a discrete maximum principle. The diffusion operator is shown to be Lipschitz continuous and linearity preserving. Using these properties we provide a full stability and error analysis, which, in the diffusion dominated regime, shows existence, uniqueness and optimal convergence. Then the algebraic flux correction method is recalled and we show that the present method can be interpreted as an algebraic flux correction method for a particular definition of the flux limiters. The performance of the method is illustrated on some numerical test cases in two space dimensions.
    • Numerical approximation of the Stochastic Cahn-Hilliard Equation near the Sharp Interface Limit

      Antonopoulou, Dimitra; Banas, Lubomir; Nurnberg, Robert; Prohl, Andreas; University of Chester; University of Bielefeld; Imperial College London; University of Tuebingen
      Abstract. We consider the stochastic Cahn-Hilliard equation with additive noise term that scales with the interfacial width parameter ε. We verify strong error estimates for a gradient flow structure-inheriting time-implicit discretization, where ε only enters polynomially; the proof is based on higher-moment estimates for iterates, and a (discrete) spectral estimate for its deterministic counterpart. For γ sufficiently large, convergence in probability of iterates towards the deterministic Hele-Shaw/Mullins-Sekerka problem in the sharp-interface limit ε → 0 is shown. These convergence results are partly generalized to a fully discrete finite element based discretization. We complement the theoretical results by computational studies to provide practical evidence concerning the effect of noise (depending on its ’strength’ γ) on the geometric evolution in the sharp-interface limit. For this purpose we compare the simulations with those from a fully discrete finite element numerical scheme for the (stochastic) Mullins-Sekerka problem. The computational results indicate that the limit for γ ≥ 1 is the deterministic problem, and for γ = 0 we obtain agreement with a (new) stochastic version of the Mullins-Sekerka problem.