Numerical methods for time-fractional evolution equations with nonsmooth data: a concise overview

Over the past few decades, there has been substantial interest in evolution equations that involving a fractional-order derivative of order $\alpha\in(0,1)$ in time, due to their many successful applications in engineering, physics, biology and finance. Thus, it is of paramount importance to develop and to analyze efficient and accurate numerical methods for reliably simulating such models, and the literature on the topic is vast and fast growing. The present paper gives a concise overview on numerical schemes for the subdiffusion model with nonsmooth problem data, which are important for the numerical analysis of many problems arising in optimal control, inverse problems and stochastic analysis. We focus on the following aspects of the subdiffusion model: regularity theory, Galerkin finite element discretization in space, time-stepping schemes (including convolution quadrature and L1 type schemes), and space-time variational formulations, and compare the results with that for standard parabolic problems. Further, these aspects are showcased with illustrative numerical experiments and complemented with perspectives and pointers to relevant literature.Comment: 24 pages, 3 figure

A nested hybridizable discontinuous Galerkin method for computing second-harmonic generation in three-dimensional metallic nanostructures

In this paper, we develop a nested hybridizable discontinuous Galerkin (HDG) method to numerically solve the Maxwell's equations coupled with the hydrodynamic model for the conduction-band electrons in metals. By means of a static condensation to eliminate the degrees of freedom of the approximate solution defined in the elements, the HDG method yields a linear system in terms of the degrees of freedom of the approximate trace defined on the element boundaries. Furthermore, we propose to reorder these degrees of freedom so that the linear system accommodates a second static condensation to eliminate a large portion of the degrees of freedom of the approximate trace, thereby yielding a much smaller linear system. For the particular metallic structures considered in this paper, the resulting linear system obtained by means of nested static condensations is a block tridiagonal system, which can be solved efficiently. We apply the nested HDG method to compute the second harmonic generation (SHG) on a triangular coaxial periodic nanogap structure. This nonlinear optics phenomenon features rapid field variations and extreme boundary-layer structures that span multiple length scales. Numerical results show that the ability to identify structures which exhibit resonances at $\omega$ and $2\omega$ is paramount to excite the second harmonic response.Comment: 31 pages, 7 figure

Space-time adaptive finite elements for nonlocal parabolic variational inequalities

This article considers the error analysis of finite element discretizations and adaptive mesh refinement procedures for nonlocal dynamic contact and friction, both in the domain and on the boundary. For a large class of parabolic variational inequalities associated to the fractional Laplacian we obtain a priori and a posteriori error estimates and study the resulting space-time adaptive mesh-refinement procedures. Particular emphasis is placed on mixed formulations, which include the contact forces as a Lagrange multiplier. Corresponding results are presented for elliptic problems. Our numerical experiments for $2$-dimensional model problems confirm the theoretical results: They indicate the efficiency of the a posteriori error estimates and illustrate the convergence properties of space-time adaptive, as well as uniform and graded discretizations.Comment: 47 pages, 20 figure

Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity

Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft