Catégorie : Nouvelles

Koondanibha Mitra Thursday 13th January at 11:00   ABSTRACT: Nonlinear advection-diffusion-reaction equations are used to model various complex flow processes in porous media, and in biological systems. They also exhibit parabolic-hyperbolic and parabolic-elliptic kinds of degeneracies resulting in the loss of regularity of the solutions. The nonlinear degenerate nature of the equations makes it challenging to provide sharp error bounds to any numerical solutions of the problem. When discretized in time, such equations result in a sequence of nonlinear degenerate elliptic problems which requires linear iterative schemes to solve. The linear iterates can be used to provide upper/lower bounds to the error, and to separate the error contributions due to linearization and discretization. However, the nonlinearity, as before, impedes the derivation of sharp error bounds in the standard error norm. In the first part of this study, we provide reliable, fully computable, and locally space-time efficient a posteriori error bounds for numerical approximations of such nonlinear degenerate parabolic problems. For showing global reliability, a nonlocal-in-time error estimate is derived individually for the time-integrated $H^1(H^{-1})$, $L^2(L^2)$, and the $L^2(H^1)$ errors. A maximum principle and a degeneracy estimator are employed for the last one. Global and local space-time efficiency error bounds are then obtained in a standard $H^1(H^{-1})\cap L^2(H^1)$ norm. The reliability and efficiency norms employed coincide when there is no nonlinearity. Moreover, error contributors such as flux nonconformity, time discretization, quadrature, and data oscillation are identified and separated. Numerical tests are conducted for nondegenerate and degenerate cases having exact solutions, as well as for realistic cases. It is shown that the estimators correctly identify the errors up to a factor of the order of unity. In the second part, using linear iterative schemes, we derive reliable, fully computable, and efficient error bounds for the finite element solution of the elliptic problem which…

Lire la suite

Gregor Gantner Friday 10th December at 11:00   ABSTRACT: Currently, there is a growing interest in simultaneous space-time methods for solving parabolic evolution equations. The main reasons are that, compared to classical time-marching methods, space-time methods are much better suited for massively parallel implementation, are guaranteed to give quasi-optimal approximations from the employed trial space, have the potential to drive optimally converging simultaneously space-time adaptive refinement routines, and they provide enhanced possibilities for reduced order modelling of parameter-dependent problems. On the other hand, space-time methods require more storage. This disadvantage however vanishes for problems of optimal control, for which the solution is needed simultaneously over the whole time interval anyway. While the common space-time variational formulation of a parabolic equation results in a bilinear form that is non-coercive, [1] recently proved the well-posedness of a space-time first-order system least-squares formulation of the heat equation. Least-squares formulations always correspond to a symmetric and coercive bilinear form. In particular, the Galerkin approximation from any conforming trial space exists and is a quasi-best approximation. Additionally, the least-squares functional automatically provides a reliable and efficient error estimator. In [2], we have generalized the least-squares method of [1] to general second-order parabolic PDEs with possibly inhomogeneous Dirichlet or Neumann boundary conditions. For homogeneous Dirichlet conditions, we present in this talk convergence of a standard adaptive finite element method driven by the least-squares estimator, which has also been demonstrated in [2]. The convergence analysis is applicable to a wide range of least-squares formulations for other PDEs, answering a long-standing open question in the literature. Moreover, we employ the space-time least-squares method for parameter-dependent problems as well as optimal control problems. In both cases, the coercivity of the corresponding bilinear form plays a crucial role. [1] T. Führer and M. Karkulik. Space–time least-squares finite elements for parabolic equations.…

Lire la suite

Oliver Sutton: Thursday 3 June at 11:00 am   ABSTRACT: The linear Boltzmann problem is a widely used model for the transport of particles through a scattering medium, such as neutrons in a nuclear reactor, or photons during radiotherapy or in the atmosphere of a star. The key challenge in simulating such phenomena using this model lies in the fact that the problem is an integro-partial-differential equation in 6 independent variables: three position variables and three momentum variables (7 if time is also included). Despite this, there is a long history of this model being successfully applied in practice. A well-studied class of numerical schemes for simulating phenomena governed by this model couples a discontinuous Galerkin spatial discretisation with a multigroup discrete-ordinates discretisation in the momentum variables. Standard multigroup discrete ordinates discretisations may be viewed as employing a piecewise constant approximation space, and possess the particularly attractive characteristic of decoupling the fully-coupled problem into a sequence of three-dimensional linear transport problems which may be solved independently and in parallel. In this talk, we will discuss a new generalisation of these multigroup discrete ordinates schemes. These new schemes employ arbitrarily high order polynomials in the discretisation of the momentum variables, providing high order convergence properties, and offer a familiar Galerkin framework for their analysis. Crucially, moreover, they retain the simple algorithmic structure of their classical counterparts.

Lire la suite

Lorenzo Mascotto: Thursday 29 April at 11:00 am ABSTRACT: Solutions to elliptic partial differential equations (PDEs) on smooth domains with smooth data are smooth. However, when solving a PDE on a Lipschitz polygon, its solution is singular at the vertices of the domain. The singular behaviour is known a priori: the solution can be split as a sum of a smooth term, plus series of singular terms that belong to the kernel of the differential operator appearing in the PDE. The virtual element method (VEM) is a generalization of the finite element method (FEM) to polygonal/polyhedral meshes and is based on approximation spaces consisting of solutions to local problems mimicking the target PDE and has been recently generalized to the extended VEM (XVEM). Here, the approximation spaces are enriched with suitable singular functions. A partition of unity (PU) is used to patch local spaces. The approach of the XVEM is close to that of the extended FEM. In this talk, we present a new paradigm for enriching virtual element (VE) spaces. Instead of adding special functions to the global space and eventually patch local spaces with a PU, we modify the definition of the local spaces by tuning the boundary conditions of local problems. By doing this, local VE spaces contain the desired singular functions, but there is no need to patch them with a PU. This results in an effective and slight modification of the already existing implementation of VE codes, as well as in a natural extension for existing theoretical results.

Lire la suite

André Harnist: Thursday 1 April at 11:00 am Abstract: We consider Hybrid High-Order (HHO) approximations of Leray-Lions problems set in W^(1,p) with p in (1,2]. For this class of problems, negative powers of the gradient of the solution can appear in the flux. Depending on the expression of the latter, this can lead to a degeneracy of the problem when the gradient of the solution vanishes or becomes large. The goal of this presentation is to derive novel error estimates depending on the degeneracy of the problem inspired by [1,2,3]. Specifically, we show that, for the globally non-degenerate case, the energy-norm of the error has a convergence rate of (k+1), with k denoting the degree of the HHO approximation. In the globally degenerate case, on the other hand, the energy-norm of the error converges with a rate of (k+1)(p-1), coherently with the estimate originally proved in [4]. We additionally introduce, for each mesh element, a dimensionless number that captures the local degeneracy of the model and identifies the contribution of the element to the global error: from the fully degenerate regime, corresponding to a contribution in (k+1)(p-1), to the non-degenerate regime, corresponding to a contribution in (k+1), through all intermediate regimes. These regime-dependent error estimates are illustrated by a complete panel of numerical experiments. [1] M. Botti, D. Castanon Quiroz, D. A. Di Pietro and A. Harnist, A Hybrid High-Order method for creeping flows of non-Newtonian fluids, Submitted. 2020. URL: https://hal.archives-ouvertes.fr/hal-02519233. [2] D. A. Di Pietro and J. Droniou, The Hybrid High-Order Method for Polytopal Meshes. Modelling, Simulation and Application 19. Springer International Publishing, 2020. ISBN: 978-3-030-37202-6 (Hardcover) 978-3-030-37203-3 (eBook). DOI: 10.1007/978-3-030-37203-3. [3] D. A. Di Pietro and J. Droniou, A Hybrid High-Order method for Leray–Lions elliptic equations on general meshes, Math. Comp., volume 86, 2017, number 307, pages 2159–2191,…

Lire la suite

Omar Duran: Thursday 11 March at 15:00 There are two main approaches to derive a fully discrete method for solving the second-order acoustic wave equation in the time regime. One is to discretize directly the second-order time derivative and the Laplacian operator in space. The other approach is to transform the second-order equation into a first-order hyperbolic system. Firstly, we consider the time second-order form. We devise, analyze the energy-conservation properties, and evaluate numerically a hybrid high-order (HHO) scheme for the space discretization combined with a Newmark-like time-marching scheme. The HHO method uses as discrete unknowns cell- and face-based polynomials of some order 0 ≤ k, yielding for steady problems optimal convergence of order (k + 1) in the energy norm [1]. Secondly, inspired by ideas presented in [2] for hybridizable discontinuous Galerkin (HDG) method and the link between HDG and HHO methods in the steady case [3], first-order explicit or implicit time-marching schemes combined with the HHO method for space discretization are considered. We discuss the selection of the stabilization term and energy conservation and present numerical examples. Extension to the unfitted meshes is contemplated for the acoustic wave equation. We observe that the unfitted approach combined with local cell agglomeration leads to a comparable CFL condition as when using fitted meshes [4]. [1] D.A. Di Pietro, A. Ern, and S. Lemaire. An arbitrary-order and compact-stencil discretization of diffusion on general meshes based on local reconstruction operators. Computational Methods in Applied Mathematics. 14 (2014) 461-472. [2] M. Stanglmeier, N.C. Nguyen, J. Peraire, and B. Cockburn. An explicit hybridizable discontinuous Galerkin method for the acoustic wave equation. Computer Methods in Applied Mechanics and Engineering, 300:748–769, March 2016. [3] B. Cockburn, D. A. Di Pietro, and A. Ern. Bridging the hybrid high-order and hybridizable discontinuous Galerkin methods. ESAIM: Mathematical Modelling and…

Lire la suite