# American Institute of Mathematical Sciences

April  2018, 14(2): 817-831. doi: 10.3934/jimo.2017077

## Parameter identification techniques applied to an environmental pollution model

 1 School of Mathematics and Statistics, Nanjing University of Information Science and Technology (NUIST), Nanjing, China, 210044 2 Department of Scientific Computing, Florida State University, Tallahassee, FL 32306, USA

* Corresponding author: I. Michael Navon

* Corresponding author: Yuepeng Wang

Received  January 2016 Revised  August 2017 Published  September 2017

Fund Project: The first author is supported by NSFC (41375115,61572015) and (ICT1600262)

The retrieval of parameters related to an environmental model is explored. We address computational challenges occurring due to a significant numerical difference of up to two orders of magnitude between the two model parameters we aim to retrieve. First, the corresponding optimization problem is poorly scaled, causing minimization algorithms to perform poorly (see Gill et al., practical optimization, AP, 1981,401pp). This issue is addressed by proper rescaling. Difficulties also arise from the presence of strong nonlinearity and ill-posedness which means that the parameters do not converge to a single deterministic set of values, but rather there exists a range of parameter combinations that produce the same model behavior. We address these computational issues by the addition of a regularization term in the cost function. All these computational approaches are addressed in the framework of variational adjoint data assimilation. The used observational data are derived from numerical simulation results located at only two spatial points. The effect of different initial guess values of parameters on retrieval results is also considered. As indicated by results of numerical experiments, the method presented in this paper achieves a near perfect parameter identification, and overcomes the indefiniteness that may occur in inversion process even in the case of noisy input data.

Citation: Yuepeng Wang, Yue Cheng, I. Michael Navon, Yuanhong Guan. Parameter identification techniques applied to an environmental pollution model. Journal of Industrial & Management Optimization, 2018, 14 (2) : 817-831. doi: 10.3934/jimo.2017077
##### References:

show all references

##### References:
Flow chart of parameter estimation with NCG method
Evolution of the value of parameters with iterations for the initial guess values (1, 1)(solid line -), (10, 1)(dashed line --), (30, 10)(dotted line :), (20, 3)(dotted-dashed line -.) and (30, 1)(star dotted-dashed line *-.), respectively
The evolution of the cost functional as a function of the number of minimization iterations for the case of initial guess value $(1, 1)$
Retrieval process of $K_1$ and $(\mu_{\ast})_0$ with iteration in different cases: $(a).~ \sigma=0, \gamma=0$; $(b).~\sigma=0, \gamma\neq{0}$; $(c).~\sigma\neq{0}, \gamma=0$; $(d).~\sigma\neq{0}, \gamma\neq{0}$
Change of $(\mu_{\ast})_0$ magnified partially: $(a)$. in Fig.5 (c) without regularization; $(b)$. in Fig.5 (d) with regularization
Comparison between the methods of NCG, ncg, lbfgs and tn with different initial guess. The true parameters are (0.38, 50)
 NCG ncg lbfgs tn Iter. Result Iter. Result Iter. Result Iter. Result [0.3;40] 2 [0.3751; 50.7179] 2 [0.3759; 50.7398] 5 [0.3752; 50.7179] 6 [0.3753; 50.7524] [0.15;18] 3 [0.3794; 50.0722] 2 [0.3962; 52.4113] 11 [0.3791; 50.1205] 4 [0.3797; 50.0491] [0.01;1.2] 3 [0.3771;50.4490] 2 [0.6661; 89.1588] 10 [0.3773; 50.4361] 4 [0.3534; 54.4259]
 NCG ncg lbfgs tn Iter. Result Iter. Result Iter. Result Iter. Result [0.3;40] 2 [0.3751; 50.7179] 2 [0.3759; 50.7398] 5 [0.3752; 50.7179] 6 [0.3753; 50.7524] [0.15;18] 3 [0.3794; 50.0722] 2 [0.3962; 52.4113] 11 [0.3791; 50.1205] 4 [0.3797; 50.0491] [0.01;1.2] 3 [0.3771;50.4490] 2 [0.6661; 89.1588] 10 [0.3773; 50.4361] 4 [0.3534; 54.4259]
Results obtained with several sets of experiments
 Guessed value Error level $\sigma$ Iters $J$ Estimated Value Regularization.Para.$\gamma$ (1, 1) 0.00 12 2.339E-7 (49.4799, 0.3974) 0 (1, 1) 0.00 13 1.225E-6 (49.5557, 0.3734) 3.5E-10 (1, 1) 0.02 15 3.020E-6 (50.4796, 0.2683) 0 (1, 1) 0.02 14 4.900E-6 (50.2095, 0.3864) 3.5E-10 (20, 20) 0.02 5 4.1237E-6 (50.4308, 0.3902) 3.0E-10 (20, 20) 0.02 5 3.1178E-6 (50.0751, 0.4118) 0 (20, 20) 0.02 50 2.9000E-6 (49.1653, 0.4669) 0
 Guessed value Error level $\sigma$ Iters $J$ Estimated Value Regularization.Para.$\gamma$ (1, 1) 0.00 12 2.339E-7 (49.4799, 0.3974) 0 (1, 1) 0.00 13 1.225E-6 (49.5557, 0.3734) 3.5E-10 (1, 1) 0.02 15 3.020E-6 (50.4796, 0.2683) 0 (1, 1) 0.02 14 4.900E-6 (50.2095, 0.3864) 3.5E-10 (20, 20) 0.02 5 4.1237E-6 (50.4308, 0.3902) 3.0E-10 (20, 20) 0.02 5 3.1178E-6 (50.0751, 0.4118) 0 (20, 20) 0.02 50 2.9000E-6 (49.1653, 0.4669) 0
 [1] C.Y. Wang, M.X. Li. Convergence property of the Fletcher-Reeves conjugate gradient method with errors. Journal of Industrial & Management Optimization, 2005, 1 (2) : 193-200. doi: 10.3934/jimo.2005.1.193 [2] Guanghui Zhou, Qin Ni, Meilan Zeng. A scaled conjugate gradient method with moving asymptotes for unconstrained optimization problems. Journal of Industrial & Management Optimization, 2017, 13 (2) : 595-608. doi: 10.3934/jimo.2016034 [3] El-Sayed M.E. Mostafa. A nonlinear conjugate gradient method for a special class of matrix optimization problems. Journal of Industrial & Management Optimization, 2014, 10 (3) : 883-903. doi: 10.3934/jimo.2014.10.883 [4] Nam-Yong Lee, Bradley J. Lucier. Preconditioned conjugate gradient method for boundary artifact-free image deblurring. Inverse Problems & Imaging, 2016, 10 (1) : 195-225. doi: 10.3934/ipi.2016.10.195 [5] Xing Li, Chungen Shen, Lei-Hong Zhang. A projected preconditioned conjugate gradient method for the linear response eigenvalue problem. Numerical Algebra, Control & Optimization, 2018, 8 (4) : 389-412. doi: 10.3934/naco.2018025 [6] Qinghua Ma, Zuoliang Xu, Liping Wang. Recovery of the local volatility function using regularization and a gradient projection method. Journal of Industrial & Management Optimization, 2015, 11 (2) : 421-437. doi: 10.3934/jimo.2015.11.421 [7] Juan Carlos De los Reyes, Estefanía Loayza-Romero. Total generalized variation regularization in data assimilation for Burgers' equation. Inverse Problems & Imaging, 2019, 13 (4) : 755-786. doi: 10.3934/ipi.2019035 [8] Saman Babaie–Kafaki, Reza Ghanbari. A class of descent four–term extension of the Dai–Liao conjugate gradient method based on the scaled memoryless BFGS update. Journal of Industrial & Management Optimization, 2017, 13 (2) : 649-658. doi: 10.3934/jimo.2016038 [9] Min Li. A three term Polak-Ribière-Polyak conjugate gradient method close to the memoryless BFGS quasi-Newton method. Journal of Industrial & Management Optimization, 2020, 16 (1) : 245-260. doi: 10.3934/jimo.2018149 [10] Yanfei Wang, Dmitry Lukyanenko, Anatoly Yagola. Magnetic parameters inversion method with full tensor gradient data. Inverse Problems & Imaging, 2019, 13 (4) : 745-754. doi: 10.3934/ipi.2019034 [11] Alexandre J. Chorin, Fei Lu, Robert N. Miller, Matthias Morzfeld, Xuemin Tu. Sampling, feasibility, and priors in data assimilation. Discrete & Continuous Dynamical Systems - A, 2016, 36 (8) : 4227-4246. doi: 10.3934/dcds.2016.36.4227 [12] Bernd Hofmann, Barbara Kaltenbacher, Elena Resmerita. Lavrentiev's regularization method in Hilbert spaces revisited. Inverse Problems & Imaging, 2016, 10 (3) : 741-764. doi: 10.3934/ipi.2016019 [13] Hongxia Yin. An iterative method for general variational inequalities. Journal of Industrial & Management Optimization, 2005, 1 (2) : 201-209. doi: 10.3934/jimo.2005.1.201 [14] Jianjun Zhang, Yunyi Hu, James G. Nagy. A scaled gradient method for digital tomographic image reconstruction. Inverse Problems & Imaging, 2018, 12 (1) : 239-259. doi: 10.3934/ipi.2018010 [15] Stefan Kindermann. Convergence of the gradient method for ill-posed problems. Inverse Problems & Imaging, 2017, 11 (4) : 703-720. doi: 10.3934/ipi.2017033 [16] Daniela Saxenhuber, Ronny Ramlau. A gradient-based method for atmospheric tomography. Inverse Problems & Imaging, 2016, 10 (3) : 781-805. doi: 10.3934/ipi.2016021 [17] Wanyou Cheng, Zixin Chen, Donghui Li. Nomonotone spectral gradient method for sparse recovery. Inverse Problems & Imaging, 2015, 9 (3) : 815-833. doi: 10.3934/ipi.2015.9.815 [18] Yigui Ou, Yuanwen Liu. A memory gradient method based on the nonmonotone technique. Journal of Industrial & Management Optimization, 2017, 13 (2) : 857-872. doi: 10.3934/jimo.2016050 [19] José Antonio Carrillo, Yanghong Huang, Francesco Saverio Patacchini, Gershon Wolansky. Numerical study of a particle method for gradient flows. Kinetic & Related Models, 2017, 10 (3) : 613-641. doi: 10.3934/krm.2017025 [20] Débora A. F. Albanez, Maicon J. Benvenutti. Continuous data assimilation algorithm for simplified Bardina model. Evolution Equations & Control Theory, 2018, 7 (1) : 33-52. doi: 10.3934/eect.2018002

2018 Impact Factor: 1.025