Leszek Demkowicz

C. Carstensen, L. Demkowicz, J. Gopalakrishnan

Discontinuous Petrov Galerkin (DPG) methods are made easily implementable using `broken' test spaces, i.e., spaces of functions with no continuity constraints across mesh element interfaces. Broken spaces derivable from a standard exact sequence of first order (unbroken) Sobolev spaces are of particular interest. A characterization of interface spaces that connect the broken spaces to their unbroken counterparts is provided. Stability of certain formulations using the broken spaces can be derived from the stability of analogues that use unbroken spaces. This technique is used to provide a complete error analysis of DPG methods for Maxwell equations with perfect electric boundary conditions. The technique also permits considerable simplifications of previous analyses of DPG methods for other equations. Reliability and efficiency estimates for an error indicator also follow. Finally, the equivalence of stability for various formulations of the same Maxwell problem is proved, including the strong form, the ultraweak form, and a spectrum of forms in between.

Weifeng Qiu, Leszek Demkowicz

The paper presents a generalization of Arnold-Falk-Winther elements for three dimensional linear elasticity, to meshes with elements of variable order. The generalization is straightforward but the stability analysis involves a non-trivial modification of involved interpolation operators. The analysis addresses only the h-convergence.

Jamie Bramwell, Leszek Demkowicz, Jay Gopalakrishnan et al.

We present two new methods for linear elasticity with simultaneously yield stress and displacement approximations of optimal accuracy in both the mesh size h and polynomial degree p. This is achieved within the recently developed discontinuous Petrov-Galerkin (DPG) framework. In this framework, both the stress and the displacement approximations are discontinuous across element interfaces. We study locking-free convergence properties and the interrelationships between the two DPG methods.

S. Nagaraj, J. Grosek, S. Petrides et al.

We propose a three dimensional Discontinuous Petrov-Galerkin Maxwell approach for modeling Raman gain in fiber laser amplifiers. In contrast with popular beam propagation models, we are interested in a truly full vectorial approach. We apply the ultraweak DPG formulation, which is known to carry desirable properties for high-frequency wave propagation problems, to the coupled Maxwell signal/pump system and use a nonlinear iterative scheme to account for the Raman gain. This paper also introduces a novel and practical full-vectorial formulation of the electric polarization term for Raman gain that emphasizes the fact that the computer modeler is only given a measured bulk Raman gain coefficient. Our results provide promising qualitative corroboration of the model and methodology used.

Jaime Mora, Leszek Demkowicz

Numerical integration of the stiffness matrix in higher order finite element (FE) methods is recognized as one of the heaviest computational tasks in a FE solver. The problem becomes even more relevant when computing the Gram matrix in the algorithm of the Discontinuous Petrov Galerkin (DPG) FE methodology. Making use of 3D tensor-product shape functions, and the concept of sum factorization, known from standard high order FE and spectral methods, here we take advantage of this idea for the entire exact sequence of FE spaces defined on the hexahedron. The key piece to the presented algorithms is the exact sequence for the one-dimensional element, and use of hierarchical shape functions. Consistent with existing results, the presented algorithms for the integration of $H^1$, $H(\text{curl})$, $H(\text{div})$, and $L^2$ inner products, have the $O(p^7)$ computational complexity. Additionally, a modified version of the algorithms is proposed when the element map can be simplified, resulting in the reduced $O(p^6)$ complexity. Use of Legendre polynomials for shape functions is critical in this implementation. Computational experiments performed with $H^1$, $H(\text{div})$ and $H(\text{curl})$ test shape functions show good correspondence with the expected rates.

Brendan Keith, Leszek Demkowicz, Jay Gopalakrishnan

We introduce a cousin of the DPG method - the DPG* method - discuss their relationship and compare the two methods through numerical experiments.

Weifeng Qiu, Leszek Demkowicz

We continue our study on variable order Arnold-Falk-Winther elements for 2D elasticity in context of both affine and parametric curvilinear elements. We present an $h$-stability result for affine elements, and an asymptotic stability result for curvilinear elements. Both theoretical results are confirmed with numerical experiments.

Jaime Mora-Paz, Stefan Henneking, Leszek Demkowicz et al.

With the goal of accurately extracting the optical field losses in a three-dimensional (3D), circularly coiled waveguide (e.g., bent optical fiber), this effort presents the numerical methodologies that are implemented for an envelope Maxwell model that propagates electromagnetic fields as an entirely boundary value problem. Our unique modeling approach includes an ultraweak variational formulation of the envelope Maxwell model in the curved geometry of the bending, which is discretized by the discontinuous Petrov-Galerkin (DPG) method, which permits residual-driven mesh and polynomial-order adaptivity. This also, then, requires a unique approach for constructing perfectly matched layers (PMLs) as absorbing boundary conditions in both the direction of optical field propagation and in the tangential directions, where unguided energy escapes the waveguide. Our coiled waveguide modeling technology extracts the mode confinement losses from the propagation of the coherent optical field through the bent waveguide. We verify our simulations against the semi-analytical results from the analogous bent slab waveguide problem, and we successfully demonstrate stable convergence to loss values for the 3D coiled optical fiber problem, which has never been done previously for our specific modeling approach.

Kamil Doległo, Anna Paszyńska, Maciej Paszyński et al.

This paper deals with the following important research question. Traditionally, the neural network employs non-linear activation functions concatenated with linear operators to approximate a given physical phenomenon. They "fill the space" with the concatenations of the activation functions and linear operators and adjust their coefficients to approximate the physical phenomena. We claim that it is better to "fill the space" with linear combinations of smooth higher-order B-splines base functions as employed by isogeometric analysis and utilize the neural networks to adjust the coefficients of linear combinations. In other words, the possibilities of using neural networks for approximating the B-spline base functions' coefficients and by approximating the solution directly are evaluated. Solving differential equations with neural networks has been proposed by Maziar Raissi et al. in 2017 by introducing Physics-informed Neural Networks (PINN), which naturally encode underlying physical laws as prior information. Approximation of coefficients using a function as an input leverages the well-known capability of neural networks being universal function approximators. In essence, in the PINN approach the network approximates the value of the given field at a given point. We present an alternative approach, where the physcial quantity is approximated as a linear combination of smooth B-spline basis functions, and the neural network approximates the coefficients of B-splines. This research compares results from the DNN approximating the coefficients of the linear combination of B-spline basis functions, with the DNN approximating the solution directly. We show that our approach is cheaper and more accurate when approximating smooth physical fields.

Leszek Demkowicz, Jay Gopalakrishnan, Sriram Nagaraj et al.

A spacetime Discontinuous Petrov Galerkin (DPG) method for the linear time-dependent Schrodinger equation is proposed. The spacetime approach is particularly attractive for capturing irregular solutions. Motivated by the fact that some irregular Schrodinger solutions cannot be solutions of certain first order reformulations, the proposed spacetime method uses the second order Schrodinger operator. Two variational formulations are proved to be well-posed: a strong formulation (with no relaxation of the original equation) and a weak formulation (also called the ultraweak formulation, that transfers all derivatives onto test functions). The convergence of the DPG method based on the ultraweak formulation is investigated using an interpolation operator. A standalone appendix analyzes the ultraweak formulation for general differential operators. Reports of numerical experiments motivated by pulse propagation in dispersive optical fibers are also included.

Daniele Boffi, Martin Costabel, Monique Dauge et al.

In this paper we prove the discrete compactness property for a wide class of p-version finite element approximations of non-elliptic variational eigenvalue problems in two and three space dimensions. In a very general framework, we find sufficient conditions for the p-version of a generalized discrete compactness property, which is formulated in the setting of discrete differential forms of any order on a d-dimensional polyhedral domain. One of the main tools for the analysis is a recently introduced smoothed Poincaré lifting operator [M. Costabel and A. McIntosh, On Bogovskii and regularized Poincaré integral operators for de Rham complexes on Lipschitz domains, Math. Z., (2010)]. For forms of order 1 our analysis shows that several widely used families of edge finite elements satisfy the discrete compactness property in p-version and hence provide convergent solutions to the Maxwell eigenvalue problem. In particular, Nédélec elements on triangles and tetrahedra (first and second kind) and on parallelograms and parallelepipeds (first kind) are covered by our theory.