The conservative characteristic FD methods for atmospheric aerosol transport problems

In the paper, we develop the new conservative characteristic finite difference methods (C-CFD) for the atmospheric aerosol transport problems. We propose the time second-order and spatial high-order conservative characteristic finite difference methods for the aerosol vertical advection-diffusion process and the two-dimensional conservative characteristic finite difference methods for aerosol horizontal transport process in the second-order splitting algorithm. Based on the characteristic form of advection-diffusion equations tracking back along the characteristic curve, we treat the integrals over the tracking cells at the previous time level by the conservative interpolations and propose to treat the diffusion terms by the average along the characteristics, where the high-order discrete fluxes are obtained by approximating the cumulative mass function and are continuous at the tracking points. The important feature is that the proposed C-CFD schemes preserve mass and have second-order accuracy in time and high-order accuracy in space. Numerical tests are taken to show the accuracy in time and space and mass conservation of our C-CFD schemes, compared with the standard CFD method. A real case of air quality modelling during the 2008 Beijing Olympics and a severe haze in North China are further simulated and analyzed by using our C-CFD algorithm. Simulated results are in good agreement with observations. The developed C-CFD algorithm can be used for efficiently solving large scale atmospheric aerosol transport problems.

[1]  Meigen Zhang,et al.  Emission controls versus meteorological conditions in determining aerosol concentrations in Beijing during the 2008 Olympic Games , 2011 .

[2]  Rodolfo Bermejo,et al.  A Galerkin-characteristic algorithm for transport-diffusion equations , 1995 .

[3]  Dong Liang,et al.  A new splitting wavelet method for solving the general aerosol dynamics equation , 2008 .

[4]  Greg Yarwood,et al.  Modeling Europe with CAMx for the Air Quality Model Evaluation International Initiative (AQMEII) , 2012 .

[5]  Edouard Debry,et al.  Solving aerosol coagulation with size-binning methods , 2007 .

[6]  S. Friedlander.  Smoke, Dust, and Haze: Fundamentals of Aerosol Dynamics , 2000 .

[7]  Xiaoliang Wang,et al.  Instruction Manual for the Aerodynamic Lens Calculator , 2006 .

[8]  Georg A. Grell,et al.  Fully coupled “online” chemistry within the WRF model , 2005 .

[9]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[10]  J. Seinfeld,et al.  Development and application of the Model of Aerosol Dynamics, Reaction, Ionization, and Dissolution (MADRID) , 2004 .

[11]  Peter H. McMurry,et al.  Modal Aerosol Dynamics Modeling , 1997 .

[12]  F. Binkowski,et al.  Models-3 community multiscale air quality (cmaq) model aerosol component , 2003 .

[13]  Zong-Liang Yang,et al.  Effects of vegetation canopy processes on snow surface energy and mass balances , 2004 .

[14]  S. Madronich,et al.  Characterizations of chemical oxidants in Mexico City: A regional chemical dynamical model (WRF-Chem) study , 2007 .

[15]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[16]  Jiming Hao,et al.  Quantifying the air pollutants emission reduction during the 2008 Olympic games in Beijing. , 2010, Environmental science & technology.

[17]  Greet Janssens-Maenhout,et al.  Emissions of air pollutants and greenhouse gases over Asian regions during 2000–2008: Regional Emission inventory in ASia (REAS) version 2 , 2013 .

[18]  Qiang Zhang,et al.  The 2013 severe haze over southern Hebei, China: model evaluation, source apportionment, and policy implications , 2013 .

[19]  Dong Liang,et al.  Multi-component atmospheric aerosols prediction by a multi-functional MC-HDMR approach , 2012 .

[20]  M. Memmesheimer,et al.  Modal aerosol dynamics model for Europe: development and first applications , 1998 .

[21]  Christodoulos Pilinis.  Derivation and numerical solution of the species mass distribution equations for multicomponent particulate systems , 1990 .

[22]  Jeffrey S. Scroggs,et al.  A forward-trajectory global semi-Lagrangian transport scheme , 2003 .

[23]  Adrian Sandu,et al.  Piecewise Polynomial Solutions of Aerosol Dynamic Equation , 2006 .

[24]  A. Nenes,et al.  ISORROPIA: A New Thermodynamic Equilibrium Model for Multiphase Multicomponent Inorganic Aerosols , 1998 .

[25]  J. Verwer,et al.  Analysis of operator splitting for advection-diffusion-reaction problems from air pollution modelling , 1999 .

[26]  Bonyoung Koo,et al.  Development and application of a three-dimensional aerosol chemical transport model, PMCAMx , 2007 .

[27]  Peter Hansbo,et al.  The characteristic streamline diffusion method for convection-diffusion problems , 1992 .

[28]  Dong Liang,et al.  An efficient approach of aerosol thermodynamic equilibrium predictions by the HDMR method , 2010 .

[29]  S. Pratsinis.  Simultaneous nucleation, condensation, and coagulation in aerosol reactors , 1988 .

[30]  Bruno Sportisse.  Fundamentals in Air Pollution: From Processes to Modelling , 2009 .

[31]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[32]  Dong Liang,et al.  An efficient second‐order characteristic finite element method for non‐linear aerosol dynamic equations , 2009, International Journal for Numerical Methods in Engineering.

[33]  Nigel Wood,et al.  SLICE: A Semi‐Lagrangian Inherently Conserving and Efficient scheme for transport problems , 2002 .

[34]  Andrew Staniforth,et al.  Application of the parabolic spline method (PSM) to a multi-dimensional conservative semi-Lagrangian transport scheme (SLICE) , 2007, J. Comput. Phys..

[35]  S. Freitas,et al.  The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) – Part 1: Model description and evaluation , 2007 .

[36]  Ralf Wolke,et al.  Implicit-explicit Runge-Kutta methods applied to atmospheric chemistry-transport modelling , 2000, Environ. Model. Softw..

[37]  John H. Seinfeld,et al.  Simulation of multicomponent aerosol condensation by the moving sectional method , 1990 .

[38]  J. Seinfeld,et al.  Atmospheric Chemistry and Physics: From Air Pollution to Climate Change , 1997 .

[39]  Mary F. Wheeler,et al.  Some improved error estimates for the modified method of characteristics , 1989 .