1 / 19

Local scale air quality modelling based on CMAQ forecast data

Local scale air quality modelling based on CMAQ forecast data. A. Oliver, A. Perez-Foguet,. Laboratori de Càlcul Numèric (LaCàN) Departament de Matemàtica Aplicada III Universitat Politècnica de Catalunya - Barcelonatech. Motivation.

kevinyeager
Download Presentation

Local scale air quality modelling based on CMAQ forecast data

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Local scale air quality modelling based on CMAQ forecast data A. Oliver, A. Perez-Foguet, Laboratori de Càlcul Numèric (LaCàN) Departament de Matemàtica Aplicada III Universitat Politècnica de Catalunya - Barcelonatech

  2. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 2 Motivation Use CMAQ model with prediction capabilities to predict the air quality in a local scale using a Finite Element model Local scale domain: 16x16km2

  3. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 3 Algorithm • Construction of a mesh • Mesh adapted to the terrain • Adapted to the CMAQ layers • CMAQ data interpolation • Horizontal and vertical interpolation from data • Wind field adjustment • Mass consistent model • Pollutant dispersion modeling • Wind field plume rise perturbation • Transport and reaction pollutant simulation

  4. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 4 Mesh construction • Two dimensional triangular mesh of terrain • Stack discretization • Vertical spacing following layers • Three dimensional tetrahedral mesh creation • Smoothing and untangling in order to guarantee a minimum quality

  5. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 5 Mesh construction Layer construction • Layer by later two-dimensional mesh • Element size depends on curvature • Layers are flattened in order to get coarser meshes in the upper layer

  6. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 6 Mesh construction Stack inclusion

  7. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 7 Mesh construction Generated mesh

  8. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 8 CMAQ data interpolation • Bilinear interpolation in all the layers • Height correction • Vertical linear interpolation ____ CMAQ terrain ······· Real terrain

  9. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 9 CMAQ data interpolation Interpolated SO2 concenctration Interpolated SO4 concenctration

  10. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 10 Wind field modeling • Mass-consistent model • Lagrange multiplier

  11. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 11 Wind field modeling

  12. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 12 Plume rise modeling • Briggs formula • Buoyant (wc < 4Vo)‏ • Driving-force: gas temperature difference • Curved trajectory • Momentum (wc > 4Vo)‏ • Driving-force: Gas velocity • Vertical straight trajectory

  13. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 13 Plume rise modeling

  14. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 14 Air quality modeling Stack outflow Inlet wind boundaries Outlet wind boundaries Initial condition

  15. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 15 Air quality modeling RIVAD reactive model (4 species)

  16. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 16 Air quality modeling • Splitting • Strang Splitting • Rosembrock 2 J = Jacobian s(c)

  17. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 17 Air quality modeling • Temporal discretization: Cranck-Nicolson • Spatial discretization: Least Squares FEM • System solver: Conjugate gradient preconditioned with an Incomplete Cholesky Factorization • Matrix storage: sparse MCS

  18. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 18 Air quality modeling Plume evolution Concentration slices

  19. 12th EMS Annual Meeting & 9th European Conference on Applied Climatology | 10-14 September 2012 | 19 Conclusions • CMAQ forecast data can be used as the initial and boundary conditions in a local scale air quality model • At least a mass consistent model is essential • Splitting and FEM resolution is an efficient technique to solve the problem • A more realistic reactive model can be used

More Related