6,758 research outputs found
A stabilized finite element method for inverse problems subject to the convection-diffusion equation. I: diffusion-dominated regime
The numerical approximation of an inverse problem subject to the
convection--diffusion equation when diffusion dominates is studied. We derive
Carleman estimates that are on a form suitable for use in numerical analysis
and with explicit dependence on the P\'eclet number. A stabilized finite
element method is then proposed and analysed. An upper bound on the condition
number is first derived. Combining the stability estimates on the continuous
problem with the numerical stability of the method, we then obtain error
estimates in local - or -norms that are optimal with respect to the
approximation order, the problem's stability and perturbations in data. The
convergence order is the same for both norms, but the -estimate requires
an additional divergence assumption for the convective field. The theory is
illustrated in some computational examples.Comment: 21 pages, 6 figures; in v2 we added two remarks and an appendix on
psiDOs, and made some minor correction
Dual weighted residual based error control for nonstationary convection-dominated equations: potential or ballast?
Even though substantial progress has been made in the numerical approximation
of convection-dominated problems, its major challenges remain in the scope of
current research. In particular, parameter robust a posteriori error estimates
for quantities of physical interest and adaptive mesh refinement strategies
with proved convergence are still missing. Here, we study numerically the
potential of the Dual Weighted Residual (DWR) approach applied to stabilized
finite element methods to further enhance the quality of approximations. The
impact of a strict application of the DWR methodology is particularly focused
rather than the reduction of computational costs for solving the dual problem
by interpolation or localization.Comment: arXiv admin note: text overlap with arXiv:1803.1064
Robust error estimates in weak norms for advection dominated transport problems with rough data
We consider mixing problems in the form of transient convection--diffusion
equations with a velocity vector field with multiscale character and rough
data. We assume that the velocity field has two scales, a coarse scale with
slow spatial variation, which is responsible for advective transport and a fine
scale with small amplitude that contributes to the mixing. For this problem we
consider the estimation of filtered error quantities for solutions computed
using a finite element method with symmetric stabilization. A posteriori error
estimates and a priori error estimates are derived using the multiscale
decomposition of the advective velocity to improve stability. All estimates are
independent both of the P\'eclet number and of the regularity of the exact
solution
Primal dual mixed finite element methods for indefinite advection--diffusion equations
We consider primal-dual mixed finite element methods for the
advection--diffusion equation. For the primal variable we use standard
continuous finite element space and for the flux we use the Raviart-Thomas
space. We prove optimal a priori error estimates in the energy- and the
-norms for the primal variable in the low Peclet regime. In the high
Peclet regime we also prove optimal error estimates for the primal variable in
the norm for smooth solutions. Numerically we observe that the method
eliminates the spurious oscillations close to interior layers that pollute the
solution of the standard Galerkin method when the local Peclet number is high.
This method, however, does produce spurious solutions when outflow boundary
layer presents. In the last section we propose two simple strategies to remove
such numerical artefacts caused by the outflow boundary layer and validate them
numerically.Comment: 25 pages, 6 figures, 5 table
Fast iterative solvers for convection-diffusion control problems
In this manuscript, we describe effective solvers for the optimal control of stabilized convection-diffusion problems. We employ the local projection stabilization, which we show to give the same matrix system whether the discretize-then-optimize or optimize-then-discretize approach for this problem is used. We then derive two effective preconditioners for this problem, the �first to be used with MINRES and the second to be used with the Bramble-Pasciak Conjugate Gradient method. The key components of both preconditioners are an accurate mass matrix approximation, a good approximation of the Schur complement, and an appropriate multigrid process to enact this latter approximation. We present numerical results to demonstrate that these preconditioners result in convergence in a small number of iterations, which is robust with respect to the mesh size h, and the regularization parameter β, for a range of problems
- …