Numerical study of the coefficient identification algorithm based on ensembles of adjoint problem solutions for a production-destruction model

A numerical algorithm for the solution of an inverse coefficient problem for nonstationary, nonlinear production-destruction type model is proposed and tested on an example of the Lorenz'63 system. With an ensemble of adjoint problem solutions, the inverse problem is transformed into a quasi-linear matrix problem and solved with Newton-type algorithm. Two different ways of the adjoint ensemble construction are compared. In the first one, a trigonometric basis is used. In the second one in situ measurements are taken into account. Local convergence properties of the algorithm are studied numerically to find out when the use of more data can lead to the degradation of the reconstruction results.