NAJan 20, 2015
Discontinuous Galerkin approximation of linear parabolic problems with dynamic boundary conditionsPaola F. Antonietti, Maurizio Grasselli, Simone Stangalino et al.
In this paper we propose and analyze a Discontinuous Galerkin method for a linear parabolic problem with dynamic boundary conditions. We present the formulation and prove stability and optimal a priori error estimates for the fully discrete scheme. More precisely, using polynomials of degree $p\geq 1$ on meshes with granularity $h$ along with a backward Euler time-stepping scheme with time-step $Δt$, we prove that the fully-discrete solution is bounded by the data and it converges, in a suitable (mesh-dependent) energy norm, to the exact solution with optimal order $h^p + Δt$. The sharpness of the theoretical estimates are verified through several numerical experiments.
42.2NAApr 21
Singularities in phase separation models: a spectral element approach for the nonlocal Cahn-Hilliard equationAndrés Miniguano-Trujillo, Andrea Poiatti, Maurizio Grasselli et al.
The nonlocal Cahn-Hilliard equation provides a natural extension of the classical model for phase separation by incorporating long-range interactions through a singular convolution kernel. While this formulation admits a rich existence and regularity theory, its numerical approximation remains challenging: discretisation of the nonlocal term leads to dense operators, and the singularity of the kernel requires special treatment in collocation-based schemes. In this work, we develop an efficient and error-controlled numerical framework for the nonlocal Cahn-Hilliard system with constant mobility, logarithmic potential, Newtonian interaction kernel, and no-flux boundary conditions. Our approach is based on a pseudospectral multishape method that accurately approximates the action of singular convolution operators. We present high-resolution numerical solutions for this nonlocal system of equations that can be achieved with limited computational resources.