# American Institute of Mathematical Sciences

October  2018, 14(4): 1579-1594. doi: 10.3934/jimo.2018022

## Adjoint-based parameter and state estimation in 1-D magnetohydrodynamic (MHD) flow system

 1 School of Automation, Guangdong University of Technology, Guangzhou, Guangdong, China 2 School of Mathematical Sciences, Zhejiang University, Hangzhou, Zhejiang, China 3 School of Control Science and Engineering, Shandong University, Jinan, Shandong, China

* Corresponding author: Zongze Wu

Received  April 2017 Revised  October 2017 Published  January 2018

Fund Project: This work is supported by the National Natural Science Foundation of China grants (61703114,61673126,61473253) and the Open Research Project of the State Key Laboratory of Industrial Control Technology, Zhejiang University, China (ICT170288, ICT170301).

In this paper, an adjoint-based optimization method is employed to estimate the unknown coefficients and states arising in an one-dimensional (1-D) magnetohydrodynamic (MHD) flow, whose dynamics can be modeled by a coupled partial differential equations (PDEs) under some suitable assumptions. In this model, the coefficients of the Reynolds number and initial conditions as well as state variables are supposed to be unknown and need to be estimated. We first employ the Lagrange multiplier method to connect the dynamics of the 1-D MHD system and the cost functional defined as the least square errors between the measurements in the experiment and the numerical simulation values. Then, we use the adjoint-based method to the augmented Lagrangian cost functional to get an adjoint coupled PDEs system, and the exact gradients of the defined cost functional with respect to these unknown parameters and initial states are further derived. The existed gradient-based optimization technique such as sequential quadratic programming (SQP) is employed for minimizing the cost functional in the optimization process. Finally, we illustrate the numerical examples to verify the effectiveness of our adjoint-based estimation approach.

Citation: Zhigang Ren, Shan Guo, Zhipeng Li, Zongze Wu. Adjoint-based parameter and state estimation in 1-D magnetohydrodynamic (MHD) flow system. Journal of Industrial & Management Optimization, 2018, 14 (4) : 1579-1594. doi: 10.3934/jimo.2018022
##### References:

show all references

##### References:
Estimation of initial condition $u(x,0)$ (scenario 1).
Estimation of initial condition $B(x,0)$ (scenario 1).
Iteration process of cost function value (scenario 1).
Estimation of state $u(x,t)$ (scenario 1).
Estimation of state $B(x,t)$ (scenario 1).
Estimation of initial condition $u(x,0)$ (scenario 2).
Estimation of initial condition $B(x,0)$ (scenario 2).
Iteration process of cost function value (scenario 2).
Estimation of state $u(x,t)$ (scenario 2).
Estimation of state $B(x,t)$ (scenario 2).
 Algorithm 1. Gradient-based optimization procedure for solving Problem P$_0$. Step 1: Choose the initial guess $\nu$, $\nu_m$ and $u_{0}^{in}(x)$, $B_{0}^{in}(x)$. Step 2: Solve the 1-D MHD model (1) forward in time corresponding to initial guess values from $t=0$ to $t=T$, and solve the adjoint coupled PDE systems (7) backward in time corresponding to initial guess values from $t=T$ to $t=0$ to obtain $u(x, t)$, $B(x, t)$, $\lambda_1(x, t)$, $\lambda_2(x, t)$. Step 3: Compute cost functional and its gradients according to (4) and (5)-(6). Step 4: Use the gradients information obtained in $\mathbf{Step 3}$ to perform an optimality test. If $\nu$, $\nu_m$ and $u_{0}^{in}(x), B_0^{in}(x)$ are optimal, then exit; Otherwise, go to $\mathbf{Step 5}$. Step 5: Use the gradient information obtained in $\mathbf{Step 3}$ to calculate a searching direction. Step 6: Perform a line search to determine the optimal step length. Step 7: Compute the new points $\nu$, $\nu_m$ and $u_{0}^{in}(x)$, $B_{0}^{in}(x)$ and return to $\mathbf{Step 2}$.
 Algorithm 1. Gradient-based optimization procedure for solving Problem P$_0$. Step 1: Choose the initial guess $\nu$, $\nu_m$ and $u_{0}^{in}(x)$, $B_{0}^{in}(x)$. Step 2: Solve the 1-D MHD model (1) forward in time corresponding to initial guess values from $t=0$ to $t=T$, and solve the adjoint coupled PDE systems (7) backward in time corresponding to initial guess values from $t=T$ to $t=0$ to obtain $u(x, t)$, $B(x, t)$, $\lambda_1(x, t)$, $\lambda_2(x, t)$. Step 3: Compute cost functional and its gradients according to (4) and (5)-(6). Step 4: Use the gradients information obtained in $\mathbf{Step 3}$ to perform an optimality test. If $\nu$, $\nu_m$ and $u_{0}^{in}(x), B_0^{in}(x)$ are optimal, then exit; Otherwise, go to $\mathbf{Step 5}$. Step 5: Use the gradient information obtained in $\mathbf{Step 3}$ to calculate a searching direction. Step 6: Perform a line search to determine the optimal step length. Step 7: Compute the new points $\nu$, $\nu_m$ and $u_{0}^{in}(x)$, $B_{0}^{in}(x)$ and return to $\mathbf{Step 2}$.
 [1] Dominique Chapelle, Philippe Moireau, Patrick Le Tallec. Robust filtering for joint state-parameter estimation in distributed mechanical systems. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 65-84. doi: 10.3934/dcds.2009.23.65 [2] Chongyang Liu, Meijia Han, Zhaohua Gong, Kok Lay Teo. Robust parameter estimation for constrained time-delay systems with inexact measurements. Journal of Industrial & Management Optimization, 2021, 17 (1) : 317-337. doi: 10.3934/jimo.2019113 [3] Peter Frolkovič, Viera Kleinová. A new numerical method for level set motion in normal direction used in optical flow estimation. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 851-863. doi: 10.3934/dcdss.2020347 [4] Lekbir Afraites, Chorouk Masnaoui, Mourad Nachaoui. Shape optimization method for an inverse geometric source problem and stability at critical shape. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021006 [5] Weihong Guo, Yifei Lou, Jing Qin, Ming Yan. IPI special issue on "mathematical/statistical approaches in data science" in the Inverse Problem and Imaging. Inverse Problems & Imaging, 2021, 15 (1) : I-I. doi: 10.3934/ipi.2021007 [6] Marion Darbas, Jérémy Heleine, Stephanie Lohrengel. Numerical resolution by the quasi-reversibility method of a data completion problem for Maxwell's equations. Inverse Problems & Imaging, 2020, 14 (6) : 1107-1133. doi: 10.3934/ipi.2020056 [7] Aihua Fan, Jörg Schmeling, Weixiao Shen. $L^\infty$-estimation of generalized Thue-Morse trigonometric polynomials and ergodic maximization. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 297-327. doi: 10.3934/dcds.2020363 [8] Geir Evensen, Javier Amezcua, Marc Bocquet, Alberto Carrassi, Alban Farchi, Alison Fowler, Pieter L. Houtekamer, Christopher K. Jones, Rafael J. de Moraes, Manuel Pulido, Christian Sampson, Femke C. Vossepoel. An international initiative of predicting the SARS-CoV-2 pandemic using ensemble data assimilation. Foundations of Data Science, 2020  doi: 10.3934/fods.2021001 [9] Kien Trung Nguyen, Vo Nguyen Minh Hieu, Van Huy Pham. Inverse group 1-median problem on trees. Journal of Industrial & Management Optimization, 2021, 17 (1) : 221-232. doi: 10.3934/jimo.2019108 [10] Woocheol Choi, Youngwoo Koh. On the splitting method for the nonlinear Schrödinger equation with initial data in $H^1$. Discrete & Continuous Dynamical Systems - A, 2021  doi: 10.3934/dcds.2021019 [11] Shumin Li, Masahiro Yamamoto, Bernadette Miara. A Carleman estimate for the linear shallow shell equation and an inverse source problem. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 367-380. doi: 10.3934/dcds.2009.23.367 [12] Jianli Xiang, Guozheng Yan. The uniqueness of the inverse elastic wave scattering problem based on the mixed reciprocity relation. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2021004 [13] Shahede Omidi, Jafar Fathali. Inverse single facility location problem on a tree with balancing on the distance of server to clients. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2021017 [14] Stanislav Nikolaevich Antontsev, Serik Ersultanovich Aitzhanov, Guzel Rashitkhuzhakyzy Ashurova. An inverse problem for the pseudo-parabolic equation with p-Laplacian. Evolution Equations & Control Theory, 2021  doi: 10.3934/eect.2021005 [15] Ying Liu, Yanping Chen, Yunqing Huang, Yang Wang. Two-grid method for semiconductor device problem by mixed finite element method and characteristics finite element method. Electronic Research Archive, 2021, 29 (1) : 1859-1880. doi: 10.3934/era.2020095 [16] Hui Gao, Jian Lv, Xiaoliang Wang, Liping Pang. An alternating linearization bundle method for a class of nonconvex optimization problem with inexact information. Journal of Industrial & Management Optimization, 2021, 17 (2) : 805-825. doi: 10.3934/jimo.2019135 [17] Xingyue Liang, Jianwei Xia, Guoliang Chen, Huasheng Zhang, Zhen Wang. $\mathcal{H}_{\infty}$ control for fuzzy markovian jump systems based on sampled-data control method. Discrete & Continuous Dynamical Systems - S, 2021, 14 (4) : 1329-1343. doi: 10.3934/dcdss.2020368 [18] José Luiz Boldrini, Jonathan Bravo-Olivares, Eduardo Notte-Cuello, Marko A. Rojas-Medar. Asymptotic behavior of weak and strong solutions of the magnetohydrodynamic equations. Electronic Research Archive, 2021, 29 (1) : 1783-1801. doi: 10.3934/era.2020091 [19] Gang Bao, Mingming Zhang, Bin Hu, Peijun Li. An adaptive finite element DtN method for the three-dimensional acoustic scattering problem. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 61-79. doi: 10.3934/dcdsb.2020351 [20] Xianbo Sun, Zhanbo Chen, Pei Yu. Parameter identification on Abelian integrals to achieve Chebyshev property. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020375

2019 Impact Factor: 1.366