A posteriori finite element output bounds with adaptive mesh refinement: application to a heat transfer problem in a three-dimensional rectangular duct

Numerical simulations based on an a posteriori finite element bound method with adaptive mesh refinement are presented for the three-dimensional convection–diffusion equation. The bound method provides relevant, quantitative, inexpensive, and rigorous lower and upper bounds to the output on a very fine discretization (“truth” discretization) at a cost close to the coarse mesh calculation (“working” discretization). To achieve a desired bound gap (i.e., difference between upper and lower bounds) at the lowest cost, an adaptive mesh refinement technique is used to refine the mesh only where needed. An optimal stabilization parameter is also applied to improve the sharpness of the bound gap. In this paper, the output of a heat transfer problem in a rectangular duct with a given velocity field is investigated. The average temperature at one section of the duct is bounded for a given inlet temperature and heat flux. For this problem, the adaptive mesh refinement strategy provides the same bound gap with only half the number of elements required by an uniform mesh refinement strategy. The hybrid flux calculation on the coarse mesh introduced for the domain decomposition approach is compared with the hybrid flux calculation on the fine mesh to analyze the contribution of the hybrid flux to the bound gap.

[1]  Anthony T. Patera,et al.  Asymptotic a Posteriori Finite Element Bounds for the Outputs of Noncoercive Problems: the Helmholtz , 1999 .

[2]  J. Peraire,et al.  A posteriori finite element bounds for linear-functional outputs of elliptic partial differential equations , 1997 .

[3]  Anthony T. Patera,et al.  A general formulation for a posteriori bounds for output functionals of partial differential equations; application to the eigenvalue problem* , 1999 .

[4]  Rolf Rannacher,et al.  Adaptive Galerkin finite element methods for partial differential equations , 2001 .

[5]  Anthony T. Patera,et al.  A posteriori finite-element output bounds for the incompressible Navier-Stokes equations: application to a natural convection problem , 2001 .

[6]  Michael B. Giles,et al.  Adjoint Recovery of Superconvergent Functionals from PDE Approximations , 2000, SIAM Rev..

[7]  Charbel Farhat,et al.  Implicit parallel processing in structural mechanics , 1994 .

[8]  Anthony T. Patera,et al.  Bounds for Linear–Functional Outputs of Coercive Partial Differential Equations : Local Indicators and Adaptive Refinement , 1998 .

[9]  Marius Paraschivoiu,et al.  A posteriori finite element bounds for linear-functional outputs of coercive partial differential equations and of the Stokes problem , 1997 .

[10]  Pierre Ladevèze,et al.  Error Estimate Procedure in the Finite Element Method and Applications , 1983 .

[11]  C. Farhat,et al.  A method of finite element tearing and interconnecting and its parallel solution algorithm , 1991 .

[12]  Marius Paraschivoiu.  A posteriori finite element output bounds in three space dimensions using the FETI method , 2001 .

[13]  J. Oden,et al.  A unified approach to a posteriori error estimation using element residual methods , 1993 .

[14]  Anthony T. Patera,et al.  A hierarchical duality approach to bounds for the outputs of partial differential equations , 1998 .

[15]  Eugene M. Cliff,et al.  Computational Methods for Optimal Design and Control , 1998 .

[16]  Anthony T. Patera,et al.  A Posteriori Bounds for Linear-Functional Outputs of Crouzeix-Raviart Finite Element Discretizations of the Incompressible Stokes Problem , 2000 .