Comptes Rendus
Numerical Analysis
Numerical method for a dynamic optimization problem arising in the modeling of a population of aerosol particles
[Une méthode numérique pour un problème d'optimisation dynamique lié à une population de particules]
Comptes Rendus. Mathématique, Volume 346 (2008) no. 11-12, pp. 677-680.

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.

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.

Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crma.2008.04.016
Alexandre Caboussat 1 ; Allison Leonard 1

1 Department of Mathematics, University of Houston, 4800 Calhoun Rd, Houston, TX 77204-3008, USA
@article{CRMATH_2008__346_11-12_677_0,
     author = {Alexandre Caboussat and Allison Leonard},
     title = {Numerical method for a dynamic optimization problem arising in the modeling of a population of aerosol particles},
     journal = {Comptes Rendus. Math\'ematique},
     pages = {677--680},
     publisher = {Elsevier},
     volume = {346},
     number = {11-12},
     year = {2008},
     doi = {10.1016/j.crma.2008.04.016},
     language = {en},
}
TY  - JOUR
AU  - Alexandre Caboussat
AU  - Allison Leonard
TI  - Numerical method for a dynamic optimization problem arising in the modeling of a population of aerosol particles
JO  - Comptes Rendus. Mathématique
PY  - 2008
SP  - 677
EP  - 680
VL  - 346
IS  - 11-12
PB  - Elsevier
DO  - 10.1016/j.crma.2008.04.016
LA  - en
ID  - CRMATH_2008__346_11-12_677_0
ER  - 
%0 Journal Article
%A Alexandre Caboussat
%A Allison Leonard
%T Numerical method for a dynamic optimization problem arising in the modeling of a population of aerosol particles
%J Comptes Rendus. Mathématique
%D 2008
%P 677-680
%V 346
%N 11-12
%I Elsevier
%R 10.1016/j.crma.2008.04.016
%G en
%F CRMATH_2008__346_11-12_677_0
Alexandre Caboussat; Allison Leonard. Numerical method for a dynamic optimization problem arising in the modeling of a population of aerosol particles. Comptes Rendus. Mathématique, Volume 346 (2008) no. 11-12, pp. 677-680. doi : 10.1016/j.crma.2008.04.016. https://comptes-rendus.academie-sciences.fr/mathematique/articles/10.1016/j.crma.2008.04.016/

[1] N.R. Amundson; A. Caboussat; J.W. He; J.H. Seinfeld Primal-dual interior-point algorithm for chemical equilibrium problems related to modeling of atmospheric organic aerosols, J. Optim. Theory Appl., Volume 130 (2006) no. 3, pp. 375-407

[2] N.R. Amundson; A. Caboussat; J.W. He; C. Landry; J.H. Seinfeld A dynamic optimization problem related to organic aerosols, C. R. Acad. Sci. Paris, Ser. I, Volume 344 (2007) no. 8, pp. 519-522

[3] A. Caboussat, C. Landry, A second order scheme for solving optimization-constrained differential equations with discontinuities, in: Numerical Mathematics and Advanced Applications: Proceedings of Enumath 2007 the 7th European Conference on Numerical Mathematics and Advanced Applications, Graz, Austria, September 2007, submitted for publication

[4] A.V. Fiacco; G.P. McCormick Nonlinear Programming: Sequential Unconstrained Minimization Techniques, John Wiley and Sons, Inc., 1968

[5] J. Gondzio; A. Grothey Direct solution of linear systems of size 109 arising in optimization with interior point methods, Lecture Notes in Computer Science, Springer-Verlag, Berlin, 2006, pp. 513-525

[6] J. Gondzio; A. Grothey Solving nonlinear portfolio optimization problems with the primal-dual interior point method, Eur. J. Operational Res., Volume 181 (2007) no. 3, pp. 1019-1029

[7] E. Hairer; G. Wanner Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems, Springer Series in Computational Mathematics, vol. 14, Springer-Verlag, 1996

[8] J.H. Seinfeld; S.N. Pandis Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, Wiley, New York, 1998

[9] R.A. Zaveri, R.C. Easter, J.D. Fast, L.K. Peters, Model for simulating aerosol interactions and chemistry (MOSAIC), J. Geophys. Res. (2008), , in press | DOI

Cité par Sources :

Commentaires - Politique


Ces articles pourraient vous intéresser

A dynamic optimization problem related to organic aerosols

Neal R. Amundson; Alexandre Caboussat; Jiwen He; ...

C. R. Math (2007)


An optimization problem related to the modeling of atmospheric inorganic aerosols

Neal R. Amundson; Alexandre Caboussat; Jiwen He; ...

C. R. Math (2005)


An optimization problem related to the modeling of atmospheric organic aerosols

Neal R. Amundson; Alexandre Caboussat; Jiwen He; ...

C. R. Math (2005)