Adaptive finite elements for a linear parabolic problem

Marco Picasso
1998
Article

A posteriori error estimates for the heat equation in two space dimensions are presented. A classical discretization is used, Euler backward in time, and continuous, piecewise linear triangular finite elements in space. The error is bounded above and below by an explicit error estimator based on the residual. Numerical results are presented for uniform triangulations and constant time steps. The quality of our error estimator is discussed. An adaptive algorithm is then proposed. Successive Delaunay triangulations are generated, so that the estimated relative error is close to a preset tolerance. Again, numerical results demonstrate the efficiency of our approach. (C) 1998 Elsevier Science S.A. All rights reserved.

À propos de ce résultat

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés

Chargement

Publications associées

Chargement

Marco Picasso
1998
Article

Résumé

Official source

https://infoscience.epfl.ch/record/89263?ln=fr

À propos de ce résultat

Concepts associés (19)

Concepts associés

Chargement

Publications associées (64)

Chargement

Afficher plus

Publications associées

Chargement

Méthode des éléments finis

En analyse numérique, la méthode des éléments finis (MEF, ou FEM pour finite element method en anglais) est utilisée pour résoudre numériquement des équations aux dérivées partielles. Celles-ci pe

Analyse numérique

L’analyse numérique est une discipline à l'interface des mathématiques et de l'informatique. Elle s’intéresse tant aux fondements qu’à la mise en pratique des méthodes permettant de résoudre, par des

Maximum a posteriori

L'estimateur du maximum a posteriori (MAP), tout comme la méthode du maximum de vraisemblance, est une méthode pouvant être utilisée afin d'estimer un certain nombre de paramètres inconnus, comme le

In this paper we derive two a posteriori upper bounds for the heat equation. A continuous, piecewise linear finite element discretization in space and the Crank-Nicolson method for the time discretization are used. The error due to the space discretization is derived using anisotropic interpolation estimates and a postprocessing procedure. The error due to the time discretization is obtained using two different continuous, piecewise quadratic time reconstructions. The first reconstruction is developed following G. Akrivis, C. Makridakis, and R. H. Nochetto [Math. Comp., 75 (2006), pp. 511-531], while the second one is new. Moreover, in the case of isotropic meshes only, upper and lower bounds are provided as in [R. Verfurth, Calcolo, 40 (2003), pp. 195 212]. An adaptive algorithm is developed. Numerical studies are reported for several test cases and show that the second error estimator is more efficient than the first one. In particular, the second error indicator is of optimal order with respect to both the mesh size and the time step when using our adaptive algorithm.

2009

Chargement

Space-time adaptive algorithms for parabolic problems

Virabouth Prachittham

We developed a space and time adaptive method to simulate electroosmosis and mass transport of a sample concentration within a network of microchannels. The space adaptive criteria is based on an error estimator derived using anisotropic interpolation estimates and a post-processing procedure. In order to improve the accuracy of the numerical solution and to reduce even further the computational cost of the numerical simulation, a time adaptive procedure is combined with the one in space. To do so, a time error estimator is derived for a first model problem, the linear heat equation discretized in time with the Crank-Nicolson scheme. The main difficulty is then to obtain an optimal second order error estimator. Applying standard energy techniques with a continuous, piecewise linear approximation in time fail in recovering the optimal order. To restore the appropriate rate of convergence, a continuous piecewise quadratic polynomial function in time is needed. For this purpose, two different quadratic functions are introduced and two different time error estimators are then derived. It turns out that the second error estimator is more efficient than the first one when considering our adaptive algorithm. Thus, using the second quadratic polynomial, an upper bound for the error is derived for a second model problem, the time-dependent convection-diffusion problem discretized in time with the Crank-Nicolson scheme. The corresponding space and time error estimators are finally used for the numerical simulation of mass transport of a sample concentration within a complex network of microchannels driven by an electroosmotic flow and/or by a pressure-driven flow. Numerical results presented show the efficiency and the robustness of this approach.

EPFL2009

Chargement

An adaptive algorithm for the Stokes problem using continuous, piecewise linear stabilized finite elements and meshes with high aspect ratio

Marco Picasso

An adaptive algorithm is proposed for solving the Stokes problem with continuous, piecewise linear stabilized finite elements and meshes with high aspect ratio. Anisotropic, a posteriori error estimates in the H-1 x L-2 norm are derived, the constant being independent from the mesh aspect ratio. The estimates involve the first derivatives of the velocity error and the velocity of a dual problem having the pressure error in the right side. An explicit error indicator is then proposed. The first derivatives of the velocity error are approximated using Zienkiewicz-Zhu error estimator (post-processing). The pressure error in the right side of the dual problem is approximated by smoothing. Numerical results on meshes with large aspect ratio show that the error indicator is robust whenever the pressure error in the dual problem is approximated with sufficient accuracy. Finally, an adaptive algorithm is proposed, with goal to build an anisotropic triangulation such that the relative estimated error is close to a preset tolerance. (c) 2004 IMACS. Published by Elsevier B.V. All rights reserved.

2005

Chargement

Afficher plus