Schittkowski, K. Parameter identification in one-dimensional partial differential algebraic equations. (English) Zbl 1158.65335 GAMM-Mitt. 30, No. 2, 352-375 (2007). Summary: We discuss a couple of situations, where algebraic equations are to be attached to a system of one-dimensional partial differential equations. Besides of models leading directly to algebraic equations because of the underlying practical background, for example in case of stationary equations, there are many others where the specific mathematical structure requires a certain reformulation leading to time-independent equations. To be able to apply our approach to a large class of real-life problems, we have to take into account flux formulations, constraints, switching points, different integration areas with transition conditions, and coupled ordinary differential algebraic equations (DAEs), for example. The system of partial differential algebraic equations (PDAEs) is discretized by the method of lines leading to a large system of differential algebraic equations which can be solved by any available implicit integration method. Standard difference formulas are applied to discretize first and second partial derivatives, and upwind formulae are used for transport equations. Proceeding from given experimental data, i.e., observation times and measurements, the minimum least squares distance of measured data from a fitting criterion is computed, which depends on the solution of the system of PDAEs. Parameters to be identified can be part of the differential equations, initial, transition, or boundary conditions, coupled DAEs, constraints, fitting criterion, etc. Also the switching points can become optimization variables. The resulting least squares problem is solved by an adapted sequential quadratic programming (SQP) algorithm which retains typical features of a classical Gauss-Newton method by retaining robustness and fast convergence speed of SQP methods. The mathematical structure of the identification problems is outlined in detail, and we present a number of case studies to illustrate the different model classes which can be treated by our approach. Cited in 3 Documents MSC: 65M32 Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs 35R10 Partial functional-differential equations 35R30 Inverse problems for PDEs 65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs 65M20 Method of lines for initial value and initial-boundary value problems involving PDEs 65L80 Numerical methods for differential-algebraic equations 90C55 Methods of successive quadratic programming type 49M15 Newton-type methods 65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs Keywords:inverse problem; finite difference method; parameter estimation; data fitting; least squares optimization; partial differential algebraic equations; method of lines; differential algebraic equations; transport equations; sequential quadratic programming algorithm; Gauss-Newton method; convergence Software:DASSL; Diffpack; EASY - FIT; PDEFIT; PCOMP; Schittkowski; NL2SOL PDFBibTeX XMLCite \textit{K. Schittkowski}, GAMM-Mitt. 30, No. 2, 352--375 (2007; Zbl 1158.65335) Full Text: DOI References: [1] Argentine, Physical Reviews E 56 pp 2359– (1997) [2] (1998): Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations, SIAM, Philadelphia · Zbl 0908.65055 [3] Bär, Physical Reviews E 59 pp 337– (1999) [4] Bazeze, Numerical Methods for Partial Differential Equations 15 pp 558– (1999) [5] Birk, Journal of Process Control 9 pp 325– · Zbl 1239.05187 [6] (1990): Least Squares Methods, Elsevier · Zbl 0875.65055 [7] Blom, ACM Transactions on Mathematical Software 20 pp 194– (1994) [8] Caracotsios, Computers and Chemical Engineering 19 pp 1019– (1995) [9] Chakravarthy, SIAM Journal on Numerical Analysis 21 pp 955– (1984) [10] Chakravarthy, Lectures in Applied Mathematics 22 pp 57– (1985) [11] , , , , , (2003): Index analysis of a nonlinear PDAE system describing a molten carbonate fuel cell, to appear: PAMM · Zbl 1354.35169 [12] Chudej, PAMM-Proceedings of Applied Mathematics and Mechanics 3 pp 521– (2003) [13] (1977): Nonlinear least squares, in: The State of the Art in Numerical Analysis, D. Jacobs ed., Academic Press, New York, London [14] Dennis, ACM Transactions on Mathematical Software 7 pp 348– (1981) [15] Dennis, ACM Transactions on Mathematical Software 7 pp 369– (1981) [16] Dobmann, ACM Transactions on Mathematical Software 21 pp 233– (1995) [17] (1999): Dynamische Gebietszerlegung für Optimalsteuerungsprobleme auf vernetzten Gebieten unter Verwendung vonMehrgitterverfahren, Diploma Thesis, Dept. ofMathematics, University of Bayreuth, Germany [18] Gill, SIAM Journal on Numerical Analysis 15 pp 977– (1978) [19] , (1981): Practical Optimization, Academic Press · Zbl 0503.90062 [20] (1978): Identification of parameters in distributed systems, in: Distributed Parameter Systems, W.H. Ray, D.G. Lainiotis eds., Marcel Dekker, New York, Basel, 47-134 [21] , , (2001): Modelling, stabilization and control of flow in networks of open channels, in: Online Optimization of Large Scale Systems, M. Groetschel, S.O. Krumke, J. Rambau eds., Springer, 251-270 [22] (1991): Solving Ordinary Differential Equations II. Stiff and Differential- Algebraic Problems, Springer Series Computational Mathematics, Vol. 14, Springer · Zbl 0729.65051 [23] , (2001): Numerical solutions of the equal width wave equation using an adaptive method of lines, in: Adaptive Methods of Lines, A. Vande Wouwer, Ph. Saucec Ph., W. Schiesser eds., Chapman and Hall/CRC, Boca Raton [24] Heidebrecht, Chemical Engineering Sciences 58 pp 1029– [25] Johnson, Journal of Fluid Dynamics 42 pp 49– (1970) [26] (1999): Computational Partial Differential Equations, Lecture Notes in Computational Science and Engineering, Vol. 2, Springer, Berlin, Heidelberg [27] Levenberg, Quarterly Applied Mathematics 2 pp 164– (1944) [28] Lewis, Finite Elements in Design 34 pp 271– (2000) [29] (1982): A stabilized Gauß-Newton algorithm for unconstrained least squares problems, Report UMINF-102.82, Institute of Information Processing, University of Umea, Umea, Sweden [30] (2001): Transport Modeling in Hydrochemical Systems, Interdisciplinary Applied Mathematics, Springer, New York [31] (1996): Models of quasi-linear PDAEs with convection, Report, Dept. of Mathematics and Computer Science, University of Halle, Germany [32] Marquardt, SIAM Journal of Applied Mathematics 11 pp 431– (1963) [33] (1990): Computer aided modelling of distributed parameter process, Technical Report No. 193, School of Electrical and Computer Engineering, Chalmers University of Technology, Göteborg, Sweden [34] Pesch, PAMM-Proceedings of Applied Mathematics and Mechanics 3 pp 521– (2003) [35] (1982): A description of DASSL: A differential/algebraic system solver, in: Proceedings of the 10th IMACS World Congress, Montreal, Canada [36] Pfeiffer, Mathematics and Computers in Simulation 42 pp 617– (1996) [37] (1991): Biplicit numerical integration of partial differential equations with the transversal method of lines, Report No. 279, DFG SPP Anwendungsbezogene Optimierung und Steuerung, Technical University, Dept. of Mathematics, Munich [38] Gu, The Astrophysical Journal 534 pp 380– (2000) [39] Ramsin, Nordisk Tidstr. Informationsbehandlung (BIT) 17 pp 72– (1977) [40] (1982): The Method of Discretization in Time and Partial Differential Equations, Reidel, Dordrecht [41] (1993): An Introduction to Partial Differential Equations, Texts in Applied Mathematics, Vol. 13, Springer, Berlin [42] (1991): The Numerical Method of Lines, Academic Press, San Diego [43] Schiesser, Computers in Mathematics and Applications 28 pp 147– (1994) · Zbl 0857.65083 [44] Schittkowski, Annals of Operations Research 5 pp 485– (1985/86) [45] (1988): Solving nonlinear least squares problems by a general purpose SQPmethod, in: Trends in Mathematical Optimization, K.-H. Hoffmann, J.-B. Hiriart-Urruty, C. Lemarechal, J. Zowe eds., International Series of Numerical Mathematics, Vol. 84, Birkhäuser, pp. 295-309. [46] Schittkowski, Optimization Methods and Software 7 pp 165– (1997) [47] (2000): Parameter estimation in a mathematical model for substrate diffusion in a metabolically active cutaneous tissue, Progress in Optimization, X. Yang et al. eds., Kluwer Academic Publishers, 329 - 342 · Zbl 0953.92018 [48] Schittkowski, Optimization Methods and Software 10 pp 539– (1999) [49] Schittkowski, Structural and Multidisciplinary Optimization 23 pp 153– (2001) [50] (2002): Numerical Data Fitting in Dynamical Systems - A Practical Introduction with Applications and Software, Kluwer Academic Publishers, Dordrecht, Boston, London · Zbl 1018.65077 [51] (2004): PCOMP: A modeling language for nonlinear programs with automatic differentiation, in: Modeling Languages in Mathematical Optimization, J. Kallrath ed., Kluwer Academic Publishers, 349-367 [52] (2004): NLPQLP20: A Fortran implementation of a sequential quadratic programming algorithm with distributed and non-monotone line search-User’s guide, Report, Department of Mathematics, University of Bayreuth [53] (1985): Numerical Solution of Partial Differential Equations: Finite Difference Methods, Clarendon Press, Oxford Applied Mathematics and Computing Science Series [54] (1995): Numerical Partial Differential Equations, Texts in Applied Mathematics, Vol. 22, Springer, Berlin [55] (1970): Theory of Elasticity, McGraw Hill, New York [56] (1999): Some remarks on second order sufficient optimality conditions for nonlinear elliptic and parabolic control problems, in: Proceedings of the Workshop ’Stabilität und Sensitivit ät von Optimierungs- und Steuerungsproblemen’, Burg (Spreewald), Germany, 21.-23.4.99 [57] Segal A. (1995): Numerik partieller Differentialgleichungen für Ingenieure, Teubner [58] Wansbrough, Chemical Engineering 5 pp 95– (1985) [59] Yang, Journal of Computational Physics 102 pp 139– (1992) [60] (1985): Construction of a class of symmetric TVD schemes, Lectures in Applied Mathematics, Vol. 22, 381-395, Springer, Berlin [61] (1972): Differential equations for the haertbeat and nerve impulse, in: Towards a Theoretical Biology, C.H.Waddington ed., Edinburgh University Press, Vol. 4, 8-67 This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.