# American Institute of Mathematical Sciences

July  2021, 14(7): 2557-2570. doi: 10.3934/dcdss.2020400

## Self-adaptive algorithm based on a posteriori analysis of the error applied to air quality forecasting using the finite volume method

 1 MAE2D laboratory, University Abdelmalek Essaadi, Polydisciplinary Faculty of Larache, Road of Rabat, Larache, Morocco 2 National Meteorological Direction of Morocco, Airport Casa-Anfa, Casa Oasis - Casablanca, Morocco

* Corresponding author: tayeq.hatim@gmail.com (Hatim Tayeq)

Received  July 2019 Revised  October 2019 Published  July 2021 Early access  June 2020

Fund Project: This work was partially supported by the project PICS Mairoc with a corporation with the Jean Leray Laboratory of Nantes university (France) and the National Meteorological Direction (DMN) of Morocco

In this work, we present a self-adaptive algorithm based on the techniques of the a posteriori estimates for the transport equation modeling the dispersion of pollutants in the atmospheric boundary layer at the local scale (industrial accident, urban air quality). The goal is to provide a powerful model for forecasting pollutants concentrations with better manipulation of available computing resources.

This analysis is based on a vertex-centered Finite Volume Method in space and an implicit Euler scheme in time. We apply and validate our model, using a self-adaptive algorithm, with real atmospheric data of the Grand Casablanca area (Morocco).

Citation: Hatim Tayeq, Amal Bergam, Anouar El Harrak, Kenza Khomsi. Self-adaptive algorithm based on a posteriori analysis of the error applied to air quality forecasting using the finite volume method. Discrete & Continuous Dynamical Systems - S, 2021, 14 (7) : 2557-2570. doi: 10.3934/dcdss.2020400
##### References:

show all references

##### References:
Vertex-centred cell in 2D
The measurement stations on the Grand Casablanca area
Numerical simulations obtained using a uniform mesh with 443 elements at two times $t = 50$, (A), and $t = 1400$, (B)
Numerical simulations obtained using a uniform mesh with 22887 elements at two times $t = 50$, (A), and $t = 1400$, (B)
(A) and (B) are the approximate solution using a adaptive mesh at $t = 50$ and $t = 1400$ respectively in level 2
(A), (B) and (C) are the adaptively refined meshes at $t = 50$ in level 1, 2 and 3 respectively
(A), (B) and (C) are the adaptively refined meshes at $t = 1400$ in level 1, 2 and 3 respectively
Forecasted concentration of ozone before and after the implementation of our algorithm in two positions: Jahid station, (A) and Khansae station, (B)
The calculated parameters for the numerical model
 Parameter Signification Value $D$ Diffusion coefficient $95.3\ m^2.s^{-1}$ $V$ Velocity vector $[3.9\ m.s^{-1}, 1.58\ m.s^{-1}]$ $r$ Reaction balance $3.9$ $f$ Second member $2.95\ \mu g.m^{-3}$ $c_0$ Initial concentration $64\ \mu g.m^{-3}$
 Parameter Signification Value $D$ Diffusion coefficient $95.3\ m^2.s^{-1}$ $V$ Velocity vector $[3.9\ m.s^{-1}, 1.58\ m.s^{-1}]$ $r$ Reaction balance $3.9$ $f$ Second member $2.95\ \mu g.m^{-3}$ $c_0$ Initial concentration $64\ \mu g.m^{-3}$
The self-adaptive algorithm at an iteration
 Step 1: Mesh generating and data input. Step 2: Solve the discrete problem. Step 3: Calculate the local error indicators. Step 4: Mark mesh cells and adapt mesh. Step 5: Solve the discrete problem in the new adapted mesh. Step 6: if the stopping test is satisfied, go to the step 7 else, go to the step 3. Step 7: Interpolate solution and visualization.
 Step 1: Mesh generating and data input. Step 2: Solve the discrete problem. Step 3: Calculate the local error indicators. Step 4: Mark mesh cells and adapt mesh. Step 5: Solve the discrete problem in the new adapted mesh. Step 6: if the stopping test is satisfied, go to the step 7 else, go to the step 3. Step 7: Interpolate solution and visualization.
The results obtained before and after the application of the self-adaptive algorithm of the mesh at $t = 50$
 Level Number of jump residual CPU elements indicator indicator time primal mesh $443$ $5.0132e-03$ $1.6088e-04$ $1.4053\ s$ $t=50$ $1$ $196$ $9.8885e-05$ $9.4178e-04$ $6.0905\ s$ $2$ $278$ $3.1335e-05$ $2.6348e-04$ $7.2734\ s$ $3$ $629$ $7.8928e-06$ $3.1045e-06$ $8.2617\ s$ uniform mesh $22887$ $1.3548e-07$ $1.9964e-06$ $53.4408\ s$
 Level Number of jump residual CPU elements indicator indicator time primal mesh $443$ $5.0132e-03$ $1.6088e-04$ $1.4053\ s$ $t=50$ $1$ $196$ $9.8885e-05$ $9.4178e-04$ $6.0905\ s$ $2$ $278$ $3.1335e-05$ $2.6348e-04$ $7.2734\ s$ $3$ $629$ $7.8928e-06$ $3.1045e-06$ $8.2617\ s$ uniform mesh $22887$ $1.3548e-07$ $1.9964e-06$ $53.4408\ s$
The results obtained before and after the application of the self-adaptive algorithm of the mesh at $t = 1400$
 Level Number of Jump Residual CPU elements indicator indicator time Primal mesh $443$ $1.6483e-05$ $1.9601e-05$ $1.4053\ s$ Adaptive mesh $1$ $828$ $1.3212e-07$ $1.9516e-06$ $18.1240\ s$ $2$ $1374$ $7.2050e-08$ $1.5499e-07$ $25.0560\ s$ uniform mesh $22887$ $1.1405e-08$ $1.2477e-07$ $95.6868\ s$
 Level Number of Jump Residual CPU elements indicator indicator time Primal mesh $443$ $1.6483e-05$ $1.9601e-05$ $1.4053\ s$ Adaptive mesh $1$ $828$ $1.3212e-07$ $1.9516e-06$ $18.1240\ s$ $2$ $1374$ $7.2050e-08$ $1.5499e-07$ $25.0560\ s$ uniform mesh $22887$ $1.1405e-08$ $1.2477e-07$ $95.6868\ s$
 [1] Zhao-Han Sheng, Tingwen Huang, Jian-Guo Du, Qiang Mei, Hui Huang. Study on self-adaptive proportional control method for a class of output models. Discrete & Continuous Dynamical Systems - B, 2009, 11 (2) : 459-477. doi: 10.3934/dcdsb.2009.11.459 [2] Ya-Zheng Dang, Zhong-Hui Xue, Yan Gao, Jun-Xiang Li. Fast self-adaptive regularization iterative algorithm for solving split feasibility problem. Journal of Industrial & Management Optimization, 2020, 16 (4) : 1555-1569. doi: 10.3934/jimo.2019017 [3] Gang Qian, Deren Han, Lingling Xu, Hai Yang. Solving nonadditive traffic assignment problems: A self-adaptive projection-auxiliary problem method for variational inequalities. Journal of Industrial & Management Optimization, 2013, 9 (1) : 255-274. doi: 10.3934/jimo.2013.9.255 [4] Anouar El Harrak, Hatim Tayeq, Amal Bergam. A posteriori error estimates for a finite volume scheme applied to a nonlinear reaction-diffusion equation in population dynamics. Discrete & Continuous Dynamical Systems - S, 2021, 14 (7) : 2183-2197. doi: 10.3934/dcdss.2021062 [5] Stefan Berres, Ricardo Ruiz-Baier, Hartmut Schwandt, Elmer M. Tory. An adaptive finite-volume method for a model of two-phase pedestrian flow. Networks & Heterogeneous Media, 2011, 6 (3) : 401-423. doi: 10.3934/nhm.2011.6.401 [6] Alberto Bressan, Wen Shen. A posteriori error estimates for self-similar solutions to the Euler equations. Discrete & Continuous Dynamical Systems, 2021, 41 (1) : 113-130. doi: 10.3934/dcds.2020168 [7] Ruiqi Li, Yifan Chen, Xiang Zhao, Yanli Hu, Weidong Xiao. Time series based urban air quality predication. Big Data & Information Analytics, 2016, 1 (2&3) : 171-183. doi: 10.3934/bdia.2016003 [8] Hsueh-Chen Lee, Hyesuk Lee. An a posteriori error estimator based on least-squares finite element solutions for viscoelastic fluid flows. Electronic Research Archive, 2021, 29 (4) : 2755-2770. doi: 10.3934/era.2021012 [9] John D. Towers. An explicit finite volume algorithm for vanishing viscosity solutions on a network. Networks & Heterogeneous Media, 2021  doi: 10.3934/nhm.2021021 [10] Lili Ju, Wensong Wu, Weidong Zhao. Adaptive finite volume methods for steady convection-diffusion equations with mesh optimization. Discrete & Continuous Dynamical Systems - B, 2009, 11 (3) : 669-690. doi: 10.3934/dcdsb.2009.11.669 [11] Hao Wang, Wei Yang, Yunqing Huang. An adaptive edge finite element method for the Maxwell's equations in metamaterials. Electronic Research Archive, 2020, 28 (2) : 961-976. doi: 10.3934/era.2020051 [12] Timilehin Opeyemi Alakoya, Lateef Olakunle Jolaoso, Oluwatosin Temitope Mewomo. A self adaptive inertial algorithm for solving split variational inclusion and fixed point problems with applications. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020152 [13] Benedict Geihe, Martin Rumpf. A posteriori error estimates for sequential laminates in shape optimization. Discrete & Continuous Dynamical Systems - S, 2016, 9 (5) : 1377-1392. doi: 10.3934/dcdss.2016055 [14] Walter Allegretto, Yanping Lin, Ningning Yan. A posteriori error analysis for FEM of American options. Discrete & Continuous Dynamical Systems - B, 2006, 6 (5) : 957-978. doi: 10.3934/dcdsb.2006.6.957 [15] Caterina Calgaro, Meriem Ezzoug, Ezzeddine Zahrouni. Stability and convergence of an hybrid finite volume-finite element method for a multiphasic incompressible fluid model. Communications on Pure & Applied Analysis, 2018, 17 (2) : 429-448. doi: 10.3934/cpaa.2018024 [16] Christos V. Nikolopoulos, Georgios E. Zouraris. Numerical solution of a non-local elliptic problem modeling a thermistor with a finite element and a finite volume method. Conference Publications, 2007, 2007 (Special) : 768-778. doi: 10.3934/proc.2007.2007.768 [17] Nan Li, Song Wang, Shuhua Zhang. Pricing options on investment project contraction and ownership transfer using a finite volume scheme and an interior penalty method. Journal of Industrial & Management Optimization, 2020, 16 (3) : 1349-1368. doi: 10.3934/jimo.2019006 [18] Tan Bui-Thanh, Omar Ghattas. A scalable algorithm for MAP estimators in Bayesian inverse problems with Besov priors. Inverse Problems & Imaging, 2015, 9 (1) : 27-53. doi: 10.3934/ipi.2015.9.27 [19] Gang Bao, Mingming Zhang, Bin Hu, Peijun Li. An adaptive finite element DtN method for the three-dimensional acoustic scattering problem. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 61-79. doi: 10.3934/dcdsb.2020351 [20] Michael Hintermüller, Monserrat Rincon-Camacho. An adaptive finite element method in $L^2$-TV-based image denoising. Inverse Problems & Imaging, 2014, 8 (3) : 685-711. doi: 10.3934/ipi.2014.8.685

2020 Impact Factor: 2.425