Mathematical Modeling of Atmospheric Flow and Computation of Convex Envelopes

A. Caboussat

Mathematical Modelling of Natural Phenomena (2011)

  • Volume: 6, Issue: 5, page 44-66
  • ISSN: 0973-5348

Abstract

top
Atmospheric flow equations govern the time evolution of chemical concentrations in the atmosphere. When considering gas and particle phases, the underlying partial differential equations involve advection and diffusion operators, coagulation effects, and evaporation and condensation phenomena between the aerosol particles and the gas phase. Operator splitting techniques are generally used in global air quality models. When considering organic aerosol particles, the modeling of the thermodynamic equilibrium of each particle leads to the determination of the convex envelope of the energy function. Two strategies are proposed to address the computation of convex envelopes. The first one is based on a primal-dual interior-point method, while the second one relies on a variational formulation, an appropriate augmented Lagrangian, an Uzawa iterative algorithm, and finite element techniques. Numerical experiments are presented for chemical systems of atmospheric interest, in order to compute convex envelopes in various space dimensions.

How to cite

top

Caboussat, A.. "Mathematical Modeling of Atmospheric Flow and Computation of Convex Envelopes." Mathematical Modelling of Natural Phenomena 6.5 (2011): 44-66. <http://eudml.org/doc/222410>.

@article{Caboussat2011,
abstract = {Atmospheric flow equations govern the time evolution of chemical concentrations in the atmosphere. When considering gas and particle phases, the underlying partial differential equations involve advection and diffusion operators, coagulation effects, and evaporation and condensation phenomena between the aerosol particles and the gas phase. Operator splitting techniques are generally used in global air quality models. When considering organic aerosol particles, the modeling of the thermodynamic equilibrium of each particle leads to the determination of the convex envelope of the energy function. Two strategies are proposed to address the computation of convex envelopes. The first one is based on a primal-dual interior-point method, while the second one relies on a variational formulation, an appropriate augmented Lagrangian, an Uzawa iterative algorithm, and finite element techniques. Numerical experiments are presented for chemical systems of atmospheric interest, in order to compute convex envelopes in various space dimensions.},
author = {Caboussat, A.},
journal = {Mathematical Modelling of Natural Phenomena},
keywords = {atmospheric flow; air quality; convex envelopes; interior-point method; variational problem; augmented Lagrangian},
language = {eng},
month = {8},
number = {5},
pages = {44-66},
publisher = {EDP Sciences},
title = {Mathematical Modeling of Atmospheric Flow and Computation of Convex Envelopes},
url = {http://eudml.org/doc/222410},
volume = {6},
year = {2011},
}

TY - JOUR
AU - Caboussat, A.
TI - Mathematical Modeling of Atmospheric Flow and Computation of Convex Envelopes
JO - Mathematical Modelling of Natural Phenomena
DA - 2011/8//
PB - EDP Sciences
VL - 6
IS - 5
SP - 44
EP - 66
AB - Atmospheric flow equations govern the time evolution of chemical concentrations in the atmosphere. When considering gas and particle phases, the underlying partial differential equations involve advection and diffusion operators, coagulation effects, and evaporation and condensation phenomena between the aerosol particles and the gas phase. Operator splitting techniques are generally used in global air quality models. When considering organic aerosol particles, the modeling of the thermodynamic equilibrium of each particle leads to the determination of the convex envelope of the energy function. Two strategies are proposed to address the computation of convex envelopes. The first one is based on a primal-dual interior-point method, while the second one relies on a variational formulation, an appropriate augmented Lagrangian, an Uzawa iterative algorithm, and finite element techniques. Numerical experiments are presented for chemical systems of atmospheric interest, in order to compute convex envelopes in various space dimensions.
LA - eng
KW - atmospheric flow; air quality; convex envelopes; interior-point method; variational problem; augmented Lagrangian
UR - http://eudml.org/doc/222410
ER -

References

top
  1. 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., 344 (2007), No. 8, 519–522.  
  2. N. R. Amundson, A. Caboussat, J. W. He, C. Landry, C. Tong, J. H. Seinfeld. A new atmospheric aerosol phase equilibrium model (UHAERO): organic systems. Atmos. Chem. Phys., 7 (2007), 4675–4698.  
  3. N. R. Amundson, A. Caboussat, J. W. He, A. V. Martynenko, J. H. Seinfeld, K. Y. Yoo. A new inorganic atmospheric aerosol phase equilibrium model (UHAERO). Atmos. Chem. Phys., 6 (2006), 975–992.  
  4. 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., 130 (2006), No. 3, 375–407.  
  5. H. Y. Benson, D. F. Shanno. Interior-point methods for nonconvex nonlinear programming: regularization and warmstarts. Comput. Optim. Appl., 40 (2008), No. 2, 143–189.  
  6. A. Caboussat. Primal-dual interior-point method for thermodynamic gas-particle partitioning. Computational Optimization and Applications, 48 (2011), No. 3, 717–745.  
  7. A. Caboussat, R. Glowinski. A numerical method for a non-smooth advection-diffusion problem arising in sand mechanics. Com. Pure. Appl. Anal, 8 (2008), No. 1, 161–178.  
  8. A. Caboussat, R. Glowinski, V. Pons. An augmented Lagrangian approach to the numerical solution of a non-smooth eigenvalue problem. J. Numer. Math, 17 (2009), No. 1, 3–26.  
  9. A. Caboussat, C. Landry, J. Rappaz. Optimization problem coupled with differential equations: A numerical algorithm mixing an interior-point method and event detection. J. Optim. Theory Appl., 147 (2010), No. 1, 141–156.  
  10. G. R. Carmichael, L. K. Peters, T. Kitada. A second generation model for the regional-scale transport/chemistry/deposition. Atm. Env., 20 (1986), 173.  
  11. E. J. Dean, R. Glowinski. An augmented Lagrangian approach to the numerical solution of the Dirichlet problem for the elliptic Monge-Ampère equation in two dimensions. Electronic Transactions in Numerical Analysis, 22 (2006), 71–96.  
  12. E. J. Dean, R. Glowinski, G. Guidoboni. On the numerical simulation of Bingham visco-plastic flow: old and new results. Journal of Non Newtonian Fluid Mechanics, 142 (2007), 36–62.  
  13. A. V. Fiacco, G. P. McCormick. Nonlinear programming : sequential unconstrained minimization techniques, Wiley, New York, 1968.  
  14. M. Fortin, R. Glowinski. Augmented Lagrangian Methods: Applications to the Numerical Solution of Boundary-Value Problems, Studies in Mathematics and Its Applications. Elsevier Science Ltd, 1983.  
  15. R. Glowinski. Numerical Methods for Nonlinear Variational Problems, Springer-Verlag, New York, NY, 1984.  
  16. R. Glowinski, P. Le Tallec. Augmented Lagrangians and Operator-Splitting Methods in Nonlinear Mechanics, SIAM, Philadelphia, 1989.  
  17. J. Gondzio, A. Grothey. A new unblocking technique to warmstart interior point methods based on sensitivity analysis. SIAM Journal on Optimization, 19 (2008), No. 3, 1184–1210.  
  18. M. Z. Jacobson. Fundamentals of Atmospheric Modeling, Cambridge, second edition, 2005.  
  19. C. Landry. Numerical Analysis of Optimization-Constrained Differential Equations: Applications to Atmospheric Chemistry. PhD thesis, Ecole Polytechnique Fédérale de Lausanne, 2009. Available at .  URIhttp://library.epfl.ch/theses/?nr=4345
  20. C. Landry, A. Caboussat, E. Hairer. Solving optimization-constrained differential equations with discontinuity points, with application to atmospheric chemistry. SIAM J. Sci. Comp., 31 (2009), No. 5, 3806–3826.  
  21. D. Lanser, J. Verwer. Analysis of operator splitting for advection-diffusion-reaction problems from air pollution modelling. J. Comput. Appl. Math., 111 (1999), 201–216.  
  22. C. M. McDonald, C. A. Floudas. GLOPEQ: A new computational tool for the phase and chemical equilibrium problem. Computers and Chemical Engineering, 21 (1996), No. 1, 1–23.  
  23. G. J. McRae, W. R. Goodin, J. H. Seinfeld. Numerical solution of the atmospheric diffusion equation for chemically reacting flows. J. Comput. Phys., 45 (1982), No. 1, 1–42.  
  24. Z. Meng, D. Dabdub, J. H. Seinfeld. Size-resolved and chemically resolved model of atmospheric aerosol dynamics. J. Geophys. Res., 103 (1998), 3419–3436.  
  25. K. Nguyen, A. Caboussat, D. Dabdub. Mass conservative, positive definite integrator for atmospheric chemical dynamics. Atmos. Env., 43 (2009), No. 40, 6287–6295.  
  26. K. Nguyen, D. Dabdub. Semi-lagrangian flux scheme for the solution of the aerosol condensation/evaporation equation. Aerosol Science & Technology, 36 (2002), 407–418.  
  27. R. T. Rockafellar. Convex analysis, Princeton University Press, Princeton, NJ, 1970.  
  28. J. H. Seinfeld, S. N. Pandis. Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, Wiley, New York, 1998.  
  29. S. Solomon, D. Qin, M. Manning, Z. Chen, M. Marquis, K.B. Averyt, M. Tignor, H.L. Miller, editors. Intergovernmental Panel on Climate Change: Fourth Assessment Report: Climate Change 2007, The Physical Science Basis, Cambridge University Press, 2007.  
  30. B. Sportisse. A review of current issues in air pollution modeling and simulation. Comput. Geosci., 11 (2007), 159–181.  
  31. J.G. Verwer, W. Hundsdorfer, J.G. Blom. Numerical time integration for air pollution models. Surveys Math. Ind., 10 (2002), 107–174.  
  32. R. A. Zaveri, R. C. Easter, J. D. Fast, L. K. Peters. Model for simulating aerosol interactions and chemistry (MOSAIC). J. Geophys. Res. D (Atmospheres), 113 (2008), No. D13, D13204.  
  33. R. A. Zaveri, R. C. Easter, L. K. Peters. A computationally efficient multicomponent equilibrium solver for aerosols (MESA). J. Geophys. Res. D (Atmospheres), 110 (2005), No. D24, D24203.  

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.