Application of regularization technique to variational adjoint method: A case for nonlinear convection-diffusion problem

Abstract A discrete assimilation system for a one-dimensional variable coefficient convection–diffusion equation is constructed. The variational adjoint method combined with the regularization technique is employed to retrieve the initial condition and diffusion coefficient with the aid of a set of simulated observations. Several numerical experiments are performed: (a) retrieving both the initial condition and diffusion coefficient jointly (Experiment JR), (b) retrieving either of them separately (Experiment SR), (c) retrieving only the diffusion coefficient with the iteration count increased to 800 (Experiment NoR-SR), and (d) retrieving only the diffusion coefficient with the consideration of a regularization term based on the Experiment NoR-SR (Experiment AdR-SR). The results indicate that within the limit of 100 iterations, the retrieval quality of the Experiment SR is better than those from the Experiment JR. Compared with the initial condition, the diffusion coefficient is a little difficult to retrieve, whereas we still achieve the desired result by increasing the iterations or integrating the regularization term into the cost functional for the improvement with respect to the diffusion coefficient. Further comparisons between the Experiment NoR-SR and AdR-SR show that the regularization term can really help not only improve the precision of retrieval to a large extent, but also speed up the convergence of solution, even if some perturbations are imposed on those observations.

[1]  William H. Press,et al.  Numerical Recipes 3rd Edition: The Art of Scientific Computing , 2007 .

[2]  Milija Zupanski,et al.  Regional Four-Dimensional Variational Data Assimilation in a Quasi-Operational Forecasting Environment , 1993 .

[3]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[4]  Dusanka Zupanski,et al.  Model Error Estimation Employing an Ensemble Data Assimilation Approach , 2006 .

[5]  Fathalla A. Rihan,et al.  Four-dimensional variational data assimilation for Doppler radar wind data , 2005 .

[6]  F. L. Dimet,et al.  Variational algorithms for analysis and assimilation of meteorological observations: theoretical aspects , 1986 .

[7]  D. Daescu.  On the Sensitivity Equations of Four-Dimensional Variational (4D-Var) Data Assimilation , 2008 .

[8]  Jimy Dudhia,et al.  Four-Dimensional Variational Data Assimilation for WRF : Formulation and Preliminary Results , 2009 .

[9]  D. Zupanski,et al.  Soil model parameter estimation with ensemble data assimilation , 2009 .

[10]  P. Courtier,et al.  Variational Assimilation of Meteorological Observations With the Adjoint Vorticity Equation. I: Theory , 2007 .

[11]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[12]  Jun Zou,et al.  Numerical identifications of parameters in parabolic systems , 1998 .

[13]  S. Cohn,et al.  An Introduction to Estimation Theory , 1997 .

[14]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[15]  L. Jiang,et al.  An analysis of 4D variational data assimilation and its application , 2008, Computing.

[16]  W. Thacker,et al.  An Optimal-Control/Adjoint-Equations Approach to Studying the Oceanic General Circulation , 1989 .

[17]  Clifford H. Thurber,et al.  Parameter estimation and inverse problems , 2005 .

[18]  M. Zupanski.  Maximum Likelihood Ensemble Filter: Theoretical Aspects , 2005 .

[19]  P. Courtier,et al.  The ECMWF implementation of three‐dimensional variational assimilation (3D‐Var). I: Formulation , 1998 .

[20]  R. Todling,et al.  Adjoint Estimation of the Variation in Model Functional Output due to the Assimilation of Data , 2009 .

[21]  P. Houtekamer,et al.  A Sequential Ensemble Kalman Filter for Atmospheric Data Assimilation , 2001 .

[22]  A. Lorenc.  A Global Three-Dimensional Multivariate Statistical Interpolation Scheme , 1981 .

[23]  J. Lions.  Optimal Control of Systems Governed by Partial Differential Equations , 1971 .

[24]  Stephen E. Cohn,et al.  An Introduction to Estimation Theory (gtSpecial IssueltData Assimilation in Meteology and Oceanography: Theory and Practice) , 1997 .

[25]  Ionel Michael Navon,et al.  Practical and theoretical aspects of adjoint parameter estimation and identifiability in meteorology and oceanography , 1998 .

[26]  Xuecheng Tai,et al.  Parameter Estimation with the Augmented Lagrangian Method for a Parabolic Equation , 2005 .

[27]  P. Courtier,et al.  Variational Assimilation of Meteorological Observations With the Adjoint Vorticity Equation. Ii: Numerical Results , 2007 .

[28]  Wei Han,et al.  Theoretical analyses and numerical experiments of variational assimilation for one-dimensional ocean temperature model with techniques in inverse problems , 2004 .

[29]  Dusanka Zupanski,et al.  Four-Dimensional Variational Data Assimilation for the Blizzard of 2000 , 2002 .

[30]  Jun Zou,et al.  A new approach to convergence rate analysis of Tikhonov regularization for parameter identification in heat conduction , 2000 .

[31]  P. Cochat,et al.  Et al , 2008, Archives de pediatrie : organe officiel de la Societe francaise de pediatrie.

[32]  I. Rutherford.  Data Assimilation by Statistical Interpolation of Forecast Error Fields , 1972 .

[33]  Thomas Kaminski,et al.  Recipes for adjoint code construction , 1998, TOMS.

[34]  Gilbert Strang,et al.  Computational Science and Engineering , 2007 .

[35]  Y. Sasaki,et al.  An ObJective Analysis Based on the Variational Method , 1958 .