Abstract
We consider elliptic distributed optimal control problems with energy regularization. Here the standard L 2 -norm regularization is replaced by the H - 1 -norm leading to more focused controls. In this case, the optimality system can be reduced to a single singularly perturbed diffusion-reaction equation known as differential filter in turbulence theory. We investigate the error between the finite element approximation u Qh to the state u and the desired state u ¯ in terms of the mesh-size h and the regularization parameter Q. The choice Q = h 2 ensures optimal convergence the rate of which only depends on the regularity of the target function u ¯. The resulting symmetric and positive definite system of finite element equations is solved by the conjugate gradient (CG) method preconditioned by algebraic multigrid (AMG) or balancing domain decomposition by constraints (BDDC). We numerically study robustness and efficiency of the AMG preconditioner with respect to h, Q, and the number of subdomains (cores) p. Furthermore, we investigate the parallel performance of the BDDC preconditioned CG solver.
Original language | English |
---|---|
Pages (from-to) | 97-111 |
Number of pages | 15 |
Journal | Computational Methods in Applied Mathematics |
Volume | 22 |
Issue number | 1 |
DOIs | |
Publication status | Published - 1 Jan 2022 |
Keywords
- A Priori Error Estimates
- Elliptic Optimal Control Problems
- Energy Regularization
- Fast Solvers
- Finite Element Discretization
- Parallelization
ASJC Scopus subject areas
- Numerical Analysis
- Computational Mathematics
- Applied Mathematics