CONVERGENCE RESULTS FOR THE FLUX IDENTIFICATION IN A SCALAR CONSERVATION LAW

Here we study an inverse problem for a quasilinear hyperbolic equation. We start by proving the existence of solutions to the problem which is posed as the minimization of a suitable cost function. Then we use a Lagrangian formulation in order to formally compute the gradient of the cost function introducing an adjoint equation. Despite the fact that the Lagrangian formulation is formal and that the cost function is not necessarily differentiable, a viscous perturbation and a numerical approximation of the problem allow us to justify this computation. When the adjoint problem for the quasi-linear equation admits a smooth solution, then the perturbed adjoint states can be proved to converge to that very solution. The sequences of gradients for both perturbed problems are also proved to converge to the same element of the subdifferential of the cost function. We evidence these results for a large class of numerical schemes and particular cost functions which can be applied to the identification of isotherms for conservation laws modeling distillation or chromatography. They are illustrated by numerical examples.

[1]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[2]  O. Oleinik Discontinuous solutions of non-linear differential equations , 1963 .

[3]  E. Conway Generalized solutions of linear differential equations with discontinuous coefficients and the uniqueness question for multidimensional quasilinear conservation laws , 1967 .

[4]  R. Tyrrell Rockafellar,et al.  Convex Analysis , 1970, Princeton Landmarks in Mathematics and Physics.

[5]  G. Chavent Identification of Distributed Parameter Systems: About the Output Least Square Method, its Implementation, and Identifiability , 1979 .

[6]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[7]  J. Smoller Shock Waves and Reaction-Diffusion Equations , 1983 .

[8]  D. Hoff The sharp form of Oleĭnik’s entropy condition in several space variables , 1983 .

[9]  B. Lucier A moving mesh numerical method for hyperbolic conservation laws , 1986 .

[10]  Pierre Rouchon,et al.  Numerical Simulation of Band Propagation in Nonlinear Chromatography , 1987, Preparative-Scale Chromatography.

[11]  E. Tadmor Local error estimates for discontinuous solutions of nonlinear hyperbolic equations , 1991 .

[12]  Mauricio Sepulveda Cortes Identification de paramètres pour un système hyperbolique : application à l'estimation des isothermes en chromatographie , 1993 .

[13]  Mauricio Sepúlveda,et al.  Parameter identification for a model of chromatographic column , 1994 .

[14]  François Bouchut,et al.  Equations de transport unidimensionnelles à coefficients discontinus , 1995 .

[15]  Huijiang Zhao FIRST ORDER QUASILINEAR EQUATIONS IN SEVERAL INDEPENDENT VARIABLES WITH SINGULAR INITIAL DATA Lp(p , 1996 .

[16]  Statistical Thermodynamics Models for Multicomponent Isothermal Diphasic Equilibria , 1997 .

[17]  F. James,et al.  One-dimensional transport equations with discontinuous coefficients , 1998 .

[18]  A. Bressan Hyperbolic Systems of Conservation Laws , 1999 .

[19]  Mauricio Sepúlveda,et al.  Determination of binary competitive equilibrium isotherms from the individual chromatographic band profiles , 1999 .