Feedback stabilization of parabolic coupled system and its numerical study

Wasim Akram,Debanjana Mitra,Neela Nataraj,Mythily Ramaswamy
2023-03-20
Abstract:In the first part of this article, we study feedback stabilization of a parabolic coupled system by using localized interior controls. The system is feedback stabilizable with exponential decay $-\omega<0$ for any $\omega>0$. A stabilizing control is found in feedback form by solving a suitable algebraic Riccati equation. In the second part, a conforming finite element method is employed to approximate the continuous system by a finite dimensional discrete system. The approximated system is also feedback stabilizable (uniformly) with exponential decay $-\omega+\epsilon$, for any $\epsilon>0$ and the feedback control is obtained by solving a discrete algebraic Riccati equation. The error estimate of stabilized solutions as well as stabilizing feedback controls are obtained. We validate the theoretical results by numerical implementations.
Analysis of PDEs,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the feedback stabilization of parabolic coupled systems and their numerical study. Specifically, the authors study how to use local internal control to achieve the feedback stabilization of parabolic coupled systems and explore the exponential decay characteristics of the systems. In addition, the paper also discusses in detail the application of numerical methods in this process, especially using the finite - element method to approximate the continuous system and analyzes the error estimates of the approximate solution and feedback control. ### Main problems 1. **Feedback stabilization**: - Study how to achieve the feedback stabilization of parabolic coupled systems through local internal control. - Prove that the system can achieve the stabilization of any exponential decay \(-\omega < 0\). 2. **Numerical analysis**: - Use the finite - element method to discretize the continuous system to obtain a finite - dimensional discrete system. - Analyze the feedback stabilization properties of the discrete system, including uniform exponential decay. - Provide error estimates of the stabilized solution and feedback control. ### Formulas and methods - **Continuous system**: - The form of the system is: \[ Y'(t) = AY(t)+Bu(t)\quad\text{for all }t > 0,\quad Y(0)=Y_0 \] - Where \(A\) and \(B\) are linear operators and control operators, defined as follows: \[ A=\begin{pmatrix} \eta_0\Delta-\nu_0I&-\eta_1I\\ I&\beta_0\Delta - (\kappa+\nu_0)I \end{pmatrix},\quad D(A)=(H^2(\Omega)\cap H^1_0(\Omega))^2 \] \[ Bf=\begin{pmatrix} f\chi_O\\ 0 \end{pmatrix}\quad\text{for all }f\in L^2(\Omega) \] - **Discrete system**: - Use the finite - element method to obtain discrete operators \(A_h\) and \(B_h\), defined as follows: \[ A_{\omega h}=A_h+\omega I_h,\quad B_h = \Pi_hB \] - The form of the discrete system is: \[ Y'_h(t)=A_hY_h(t)+B_h u_h(t)\quad\text{for all }t > 0,\quad Y_h(0)=Y_{0h} \] ### Main results - **Stabilization of the continuous system**: - By solving the algebraic Riccati equation, the feedback control \(u^{\sharp}(t)=-B^{*}PY^{\sharp}(t)\) is obtained, which stabilizes the system. - It is proved that the generated semigroup \(\{e^{tA_{\omega, P}}\}_{t\geq0}\) is analytic and exponentially stable. - **Stabilization of the discrete system**: - For each \(h\), there exists a unique non - negative self - adjoint Riccati operator \(P_h\) that satisfies the discrete Riccati equation. - By solving the discrete Riccati equation, the feedback control \(u^{\sharp}_h(t)=-B_h^{*}P_hY^{\sharp}_h(t)\) of the discrete system is obtained. - It is proved that the stabilization of the discrete system is uniform, that is, the generated semigroup \(\{e^{tA_{\omega h, P_h}}\}_{t\geq0}\) is analytic and exponentially stable. - **Error estimates**: