An efficient implementation of an implicit FEM scheme for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a modelling tool for processes with anomalous diffusion or spatial heterogeneity. However, the presence of a fractional differential operator causes memory (time fractional) or nonlocality (space fractional) issues, which impose a number of computational constraints. In this paper we develop efficient, scalable techniques for solving fractional-in-space reaction diffusion equations using the finite element method on both structured and unstructured grids, and robust techniques for computing the fractional power of a matrix times a vector. Our approach is show-cased by solving the fractional Fisher and fractional Allen-Cahn reaction-diffusion equations in two and three spatial dimensions, and analysing the speed of the travelling wave and size of the interface in terms of the fractional power of the underlying Laplacian operator

A comprehensive study of non-adaptive and residual-based adaptive sampling for physics-informed neural networks

Physics-informed neural networks (PINNs) have shown to be an effective tool for solving forward and inverse problems of partial differential equations (PDEs). PINNs embed the PDEs into the loss of the neural network, and this PDE loss is evaluated at a set of scattered residual points. The distribution of these points are highly important to the performance of PINNs. However, in the existing studies on PINNs, only a few simple residual point sampling methods have mainly been used. Here, we present a comprehensive study of two categories of sampling: non-adaptive uniform sampling and adaptive nonuniform sampling. We consider six uniform sampling, including (1) equispaced uniform grid, (2) uniformly random sampling, (3) Latin hypercube sampling, (4) Halton sequence, (5) Hammersley sequence, and (6) Sobol sequence. We also consider a resampling strategy for uniform sampling. To improve the sampling efficiency and the accuracy of PINNs, we propose two new residual-based adaptive sampling methods: residual-based adaptive distribution (RAD) and residual-based adaptive refinement with distribution (RAR-D), which dynamically improve the distribution of residual points based on the PDE residuals during training. Hence, we have considered a total of 10 different sampling methods, including six non-adaptive uniform sampling, uniform sampling with resampling, two proposed adaptive sampling, and an existing adaptive sampling. We extensively tested the performance of these sampling methods for four forward problems and two inverse problems in many setups. Our numerical results presented in this study are summarized from more than 6000 simulations of PINNs. We show that the proposed adaptive sampling methods of RAD and RAR-D significantly improve the accuracy of PINNs with fewer residual points. The results obtained in this study can also be used as a practical guideline in choosing sampling methods

NURBS-SEM: A hybrid spectral element method on NURBS maps for the solution of elliptic PDEs on surfaces

Non Uniform Rational B-spline (NURBS) patches are a standard way to describe complex geometries in Computer Aided Design tools, and have gained a lot of popularity in recent years also for the approximation of partial differential equations, via the Isogeometric Analysis (IGA) paradigm. However, spectral accuracy in IGA is limited to relatively small NURBS patch degrees (roughly p 648), since local condition numbers grow very rapidly for higher degrees. On the other hand, traditional Spectral Element Methods (SEM) guarantee spectral accuracy but often require complex and expensive meshing techniques, like transfinite mapping, that result anyway in inexact geometries. In this work we propose a hybrid NURBS-SEM approximation method that achieves spectral accuracy and maintains exact geometry representation by combining the advantages of IGA and SEM. As a prototypical problem on non trivial geometries, we consider the Laplace\u2013Beltrami and Allen\u2013Cahn equations on a surface. On these problems, we present a comparison of several instances of NURBS-SEM with the standard Galerkin and Collocation Isogeometric Analysis (IGA)

Schnelle Löser für Partielle Differentialgleichungen

This workshop was well attended by 52 participants with broad geographic representation from 11 countries and 3 continents. It was a nice blend of researchers with various backgrounds

Energy preserving model order reduction of the nonlinear Schr\"odinger equation

An energy preserving reduced order model is developed for two dimensional nonlinear Schr\"odinger equation (NLSE) with plane wave solutions and with an external potential. The NLSE is discretized in space by the symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting system of Hamiltonian ordinary differential equations are integrated in time by the energy preserving average vector field (AVF) method. The mass and energy preserving reduced order model (ROM) is constructed by proper orthogonal decomposition (POD) Galerkin projection. The nonlinearities are computed for the ROM efficiently by discrete empirical interpolation method (DEIM) and dynamic mode decomposition (DMD). Preservation of the semi-discrete energy and mass are shown for the full order model (FOM) and for the ROM which ensures the long term stability of the solutions. Numerical simulations illustrate the preservation of the energy and mass in the reduced order model for the two dimensional NLSE with and without the external potential. The POD-DMD makes a remarkable improvement in computational speed-up over the POD-DEIM. Both methods approximate accurately the FOM, whereas POD-DEIM is more accurate than the POD-DMD