Partial differential equations are ubiquitous in mathematical modeling, and so is the need of their numerical solution. When discretizing partial differential equations, their structure, the regularity of their solutions, and the robustness with respect to certain parameters play prominent roles.
The workshop intends to assess these roles, looking also at their interplay.
Participation is free. If you plan to attend, please send an email to Pietro Zanotti not later than February 3, 2020.
09:0009:05 Opening
09:0509:50 R. Hiptmair  Discretizing the Advection of Differential Forms
09:5510:40 E. Georgoulis  Hypocoercivitycompatible Finite Element Methods for Kinetic Equations
10:4011:00 Coffee break in Aula C
11:0011:45 L. Beirão da Veiga  An Introduction to Curved Virtual Elements
11:5012:35 G. Rozza  Structure Preserving Reduced Order Methods for Fluidstructure Interaction Parametric Problems: State of the Art and Perspectives
12:3514:30 Lunch break
14:3015:15 R. Abgrall  Some Comments on Structure Preserving High Order Schemes: the Example of Entropy Preserving Schemes, and Other Constraints
15:2016:05 A. Linke  Gradientrobustness: a New Concept Assuring Accurate Spatial Discretizations for Vectorvalued PDEs
16:0516:25 Coffee break in Aula C
16:2517:10 G. Sangalli  A Fast Solver for Highdegree IGA
17:1518:00 Ch. Kreuzer  Oscillation in A Posteriori Error Analysis
The workshop is supported by

Some Comments on Structure Preserving High Order Schemes: the Example of Entropy Preserving Schemes, and Other Constraints
Rémi Abgrall (University of Zurich, Switzerland)
In this talk, I will consider the following question. Given a scheme designed for approximating an hyperbolic system of equation, how can it be modified so that an additional constraint is satisfied without destroying the accuracy. A typical example is fluid mechanics and the entropy. Another example, again in fluid mechanics, is, starting from a non conservative formulation of the problem, how can we build an approximation that will guaranty that the converged solution  if they exist  are weak solution of the problem ? If time permits, I will also consider the case of differential constraints.
An Introduction to Curved Virtual Elements
Lourenco Beirão da Veiga (University of MilanoBicocca, Italy)
The Virtual Element Method (VEM) was introduced in [1,2] as a generalization of the Finite Element Method that allows for general polygonal and polyhedral meshes. Polytopal meshes can be very useful for a wide range of reasons, including meshing of the domain (such as cracks) and data (such as inclusions) features, automatic use of hanging nodes, moving meshes, adaptivity. By avoiding the explicit construction of the local basis functions, Virtual Elements can easily handle general polygons/polyhedrons without the need of an overly complex construction.
The scope of the present talk is to present Virtual Elements with curved faces, introduced in [3] and further developed in [4]. Indeed, all the VEM papers in the literature make use of polygonal and polyhedral meshes, i.e. with straight edges and faces. On the other hand, as acknowledged in the finite element (FEM) literature, especially for high order methods the approximation of the domain by facets introduces an error that can dominate the analysis. This issue has lead, for example, to the development of non affine isoparametric FEM elements and to Isogeometric Analysis.
In the context of Virtual Elements, one can exploit the peculiar construction of the method that (1) does not need an explicit expression of the basis functions and (2) is directly defined in physical space, i.e. no reference element is used. This allows to define discrete spaces also on elements that are curved in such a way to exactly represent the domain of interest. The needed ingredient is a (piecewise regular) parametrization of the boundary of the domain.
In the first part of the talk we present the curved VEM on a model elliptic problem, including the definition of the involved discrete space, the theoretical aspects and some numerical tests. In the second part, we introduce a variant of the discrete space that is more suitable for solid mechanics problems, and develop the method in the realm of small deformation elasticity and inelasticity, including a final set of numerical tests.
[1] L. Beirao da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. D. Marini, A. Russo, Basic principles of virtual element methods, Math. Models Methods Appl. Sci., 23:199214 (2013)Hypocoercivitycompatible Finite Element Methods for Kinetic Equations
Emmanuil Georgoulis (University of Leicester, UK, and Technical University of Athens, Greece)
We shall discuss some recent developments on a family of Galerkin finite element methods for kinetic equations with degenerate diffusion, such as the classical Kolmogorov and FokkerPlanck equations. In particular, diffusion is not present in all spatial directions. Nonetheless, the solutions to these problems admit typically decay properties to some long time equilibrium, depending on closure by suitable boundary/decayatinfinity conditions. A key attribute of the proposed family of methods is that they also admit similar decay properties for very general families of triangulations. The method construction uses ideas by the general theory of hypocoercivity developed by Villani, along with judicious choice of numerical flux functions. These developments turn out to be sufficient to imply that the proposed finite element methods admit a priori error bounds with constants independent of the final time, despite these equations' degenerate diffusion nature. Thus, the new methods provably allow for robust error analysis for final times tending to infinity.
Discretizing the Advection of Differential Forms
Ralf Hiptmair (ETH Zurich, Switzerland)
Oscillation in A Posteriori Error Analysis
Christian Kreuzer (Technical University of Dortmund, Germany)
A posteriori error estimators are a key tool for the quality assessment of given finite element approximations to an unknown PDE solution as well as for the application of adaptive techniques.
Typically, the estimators are equivalent to the error up to an additive term, the so called oscillation. It is a common believe that this is the price for the `computability' of the estimator and that the oscillation is of higher order than the error. Cohen, DeVore, and Nochetto [CDN:2012], however, presented an example, where the error vanishes with the generic optimal rate, but the oscillation does not. Interestingly, in this example, the local $H^{1}$norms are assumed to be computed exactly and thus the computability of the estimator cannot be the reason for the asymptotic overestimation. In particular, this proves both believes wrong in general.
In this talk, we present a new approach to posteriori error analysis, where the oscillation is dominated by the error. The crucial step is a new splitting of the data into oscillation and oscillation free data. Moreover, the estimator is computable if the discrete linear system can essentially be assembled exactly.
[CDN:2012] A. Cohen, R. DeVore, and R. H. Nochetto, Convergence Rates of AFEM with H^{1} Data, Found. Comput. Math. 12 (2012):671718This is a joint work with A. Veeser.
Gradientrobustness: a New Concept Assuring Accurate Spatial Discretizations for Vectorvalued PDEs
Alexander Linke (WIAS, Berlin, Germany)
Vectorvalued PDEs like the incompressible NavierStokes equations (in primitive variables velocity and pressure) describe the time evolution of a vectorvalued quantity like the momentum density. For vectorvalued PDEs it is quite natural to derive formally on the continuous level a derived time evolution of the vorticity and the divergence. Accordingly, any space discretization for a vectorvalued PDE (implicitly) delivers a discrete vorticity and discrete divergence equation.
While the celebrated infsup stability will be shown to assure an accurate discrete divergence equation, the talk will actually focus on the question, which structural properties allow for an implicitly defined accurate discrete vorticity equation. The key observation is that the L^2orthogonality of divergencefree vector fields and gradient fields is the weak equivalent to the vector calculus identity $\nabla \times \nabla \psi = \mathbf{0}$ for arbitrary smooth scalar fields $\psi$.
In the context of the incompressible NavierStokes equations, the concept of pressurerobustness was introduced in 2016, in order to discriminate between space discretizations with accurate and inaccurate discrete vorticity equations. H(div)conforming finite element spaces have been found as an important means to realize the L^2orthogonality between discretely divergencefree vector fields and (arbitrary) gradient fields. Further, it is shown that spatial discretizations that are not pressurerobust may suffer from i) a degradation of the (preasymptotic) convergence rate, and ii) large, parameterdependent constants in error estimates. Typical flows that benefit from pressurerobustness are quasihydrostatic flows, quasigeostrophic flows and vortexdominated high Reynolds number flows.
Last but not least, the talk shows how to extend the concept of pressurerobustness to more general vectorvalued PDEs, leading to the concept of gradientrobustness. Gradient robustness assures that dominant and extreme gradient fields in a vectorvalued PDE will not lead to inaccuracies in the discretization. The talk will show examples from linear elasticity and compressible (Navier)Stokes flows at low Mach numbers and in stratified flows. Thus, connections to wellbalanced schemes and WENO schemes will be drawn.
Structure Preserving Reduced Order Methods for Fluidstructure Interaction Parametric Problems: State of the Art and Perspectives
Gianluigi Rozza (SISSA, Trieste, Italy)
We describe the state of the art and perspectives in the developments of efficient structure preserving reduced order methods for parametric nonlinear fluidstructure interaction problems by monolithic and segregated approaches, as well as the use of numerical techniques to enhance the reduction of the Kolmogorov nwidth in order to improve computational performances. Special attention is dedicated to the approximation stability of the reduced order model by supremisers, to the structure preserving property, as well as to the imposition of accurate and efficient coupling conditions to guarantee continuity of quantities at the fluidstructure interface.
Joint work with Monica Nonino, Francesco Ballarin (SISSA) and Yvon Maday (Sorbonne Paris LJLL).
A Fast Solver for Highdegree IGA
Giancarlo Sangalli (University of Pavia, Italy)
The concept of krefinement was proposed as one of the key features of isogeometric analysis, "a new, more efficient, higherorder concept", in the seminal work [1]. The idea of using highdegree and continuity splines/NURBS as a basis for a new highorder method appeared very promising from the beginning, and received confirmations from the next developments. The krefinement leads to several advantages: higher accuracy per degreeoffreedom, improved spectral accuracy, the possibility of structurepreserving smooth discretizations are the most interesting features that have been studied actively in the community. At the same time, the krefinement brings significant challenges at the computational level: using standard finite element routines, its computational cost grows with respect to the degree, making degree raising computationally expensive. This presentation gives an overview of some recent results that extend what we did in [2].
[1] T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement, Comput. Methods Appl. Mech. Engrg., Vol. 194, pp. 41354195 (2005).This research activity is developed with Monica Montardini, Mattia Tani, and other collaborators.