Numerical method for a dynamic optimization problem arising in the modeling of a population of aerosol particles

A model coupling differential equations and a sequence of constrained optimization problems is proposed for the simulation of the evolution of a population of particles at equilibrium interacting through a common medium.The first order optimality conditions of the optimization problems relaxed with barrier functions are coupled with the differential equations into a system of differential-algebraic equations that is discretized in time with an implicit first order scheme. The resulting system of nonlinear algebraic equations is solved at each time step with an interior-point/Newton method. The Newton system is block-structured and solved with Schur complement techniques, in order to take advantage of its sparsity. Application to the dynamics of a population of organic atmospheric aerosol particles is given to illustrate the evolution of particles of different sizes. To cite this article: A. Caboussat, A. Leonard, C. R. Acad. Sci. Paris, Ser. I 346 (2008).

RésuméNous proposons un modèle pour le calcul de l'évolution d'une population de particules. Ce modèle couple un système d'équations différentielles et une séquence de problèmes d'optimisation sous contraintes.Les conditions de premier ordre de chacun des problèmes d'optimisation et une discrétisation implicite des équations différentielles forment un système d'équations non linéaires qui est traité avec une méthode de point intérieur, couplé à une itération de Newton. Le système linéaire correspondant a une structure par blocs. Une méthode de résolution directe basée sur le complément de Schur prend en compte la structure creuse de la matrice et permet de découpler les différentes particules du système. Des résultats numériques pour une population de particules organiques montrent l'évolution de particules de différentes tailles. Pour citer cet article : A. Caboussat, A. Leonard, C. R. Acad. Sci. Paris, Ser. I 346 (2008).