Abstract:In this paper, a fast solver is studied for saddle point system arising from a second-order Crank-Nicolson discretization of an initial-valued parabolic PDE constrained optimal control problem, which is indefinite and ill-conditioned. Different from the saddle point system arising from the first-order Euler discretization, the saddle point system arising from Crank-Nicolson discretization has a dense and non-symmetric Schur complement, which brings challenges to fast solver designing. To remedy this, a novel symmetrization technique is applied to the saddle point system so that the new Schur complement is symmetric definite and the well-known matching-Schur-complement (MSC) preconditioner is applicable to the new Schur complement. Nevertheless, the new Schur complement is still a dense matrix and the inversion of the corresponding MSC preconditioner is not parallel-in-time (PinT) and thus time consuming. For this concern, a modified MSC preconditioner for the new Schur complement system. Our new preconditioner can be implemented in a fast and PinT way via a temporal diagonalization technique. Theoretically, the eigenvalues of the preconditioned matrix by our new preconditioner are proven to be lower and upper bounded by positive constants independent of matrix size and the regularization parameter. With such spectrum, the preconditioned conjugate gradient (PCG) solver for the Schur complement system is proven to have a convergence rate independent of matrix size and regularization parameter. To the best of my knowledge, it is the first time to have an iterative solver with problem-independent convergence rate for the saddle point system arising from Crank-Nicolson discretization of the optimal control problem. Numerical results are reported to show that the performance of the proposed preconditioner.

A Parallel-in-time Method Based on Preconditioner for Biot's Model

A Parallel Simulator for Massive Reservoir Models Utilizing Distributed-Memory Parallel Systems

Parameter-robust discretization and preconditioning of Biot's consolidation model

Some Optimally Convergent Algorithms for Decoupling the Computation of Biot’s Model

Uniformly well-posed hybridized discontinuous Galerkin/hybrid mixed discretizations for Biot's consolidation model

Biot model with generalized eigenvalue problems for scalability and robustness to parameters

Monolithic two-level Schwarz preconditioner for Biot's consolidation model in two space dimensions

Efficient parallel inversion of ParaOpt preconditioners

Comparisons of Some Iterative Algorithms for Biot Equations.

An H(div)-Conforming Finite Element Method for the Biot Consolidation Model

PiTSBiCG: Parallel in Time Stable Bi-Conjugate Gradient Algorithm

Crank-Nicolson-type iterative decoupled algorithms for Biot's consolidation model using total pressure

AN ITERATIVE DECOUPLED ALGORITHM WITH UNCONDITIONAL STABILITY FOR BIOT MODEL

A lattice Boltzmann method for Biot's consolidation model of linear poroelasticity

A priori error estimates of two fully discrete coupled schemes for Biot's consolidation model

Circulant Preconditioned Iterative Methods for Peridynamic Model Simulation

A parallel-in-time preconditioner for Crank-Nicolson discretization of a parabolic optimal control problem

A parallel preconditioner for the all-at-once linear system from evolutionary PDEs with Crank-Nicolson discretization

Oscillation-free numerical schemes for Biot's model and their iterative coupling solution

A Mixed Virtual Element Method for Biot's Consolidation Model

An optimal parallel-in-time preconditioner for parabolic optimal control problems