The effects of parametric uncertainties in simulations of a reactive plume using a Lagrangian stochastic model

Abstract A combined Lagrangian stochastic model with micro-mixing and chemical sub-models is used to investigate a reactive plume of nitrogen oxides (NOx) released into a turbulent grid flow doped with ozone (O3). Sensitivities to the model input parameters are explored for different source NOx scenarios. The wind tunnel experiments of Brown and Bilger (1996) provide the simulation conditions for the first case study where photolysis reactions are not included and the main uncertainties occur in parameters defining the turbulence scales, source size and reaction rate of NO with O3. Using nominal values of the parameters from previous studies, the model gives a good representation of the radial profile of the conserved mean scalar Γ ¯ NO x although slightly over predicts peak mean NO2 concentrations Γ ¯ NO 2 compared to the experiments. The high dimensional model representation (HDMR) method is used to investigate the effects of uncertainties in model inputs on the simulation of chemical species concentrations. For this scenario, the Lagrangian velocity structure function coefficient has the largest impact on simulated Γ ¯ NO x profiles. Photolysis reactions are then included in a chemical scheme consisting of eight reactions between species NO, O, O3 and NO2. Independent and interactive effects of 22 input parameters are studied for two source NOx scenarios using HDMR, including turbulence parameters, temperature dependant rate parameters, photolysis rates, temperature, fraction of NO in total NOx at the source and background ozone concentration [O3]. For this reactive case, the variance in the predicted mean plume centre Γ ¯ O 3 is caused by parameters describing both physical (mixing time-scale coefficient) and chemical processes (activation energy for the reaction O3+NO). The variance in predicted plume centre Γ ¯ NO 2 and root mean square NO2 concentration γ NO 2 ′ , is strongly influenced by the fraction of NO in the source NOx, and to a lesser extent the mixing time-scale coefficient. Adjusting the latter gives improved agreement with the Brown and Bilger experiment. Some weak parameter interactions are observed.

[1]  S. Madronich,et al.  Actinometric and radiometric measurement and modeling of the photolysis rate coefficient of ozone to O (^1D) during Mauna Loa Observatory Photochemistry Experiment 2 , 1996 .

[2]  D. Carslaw Evidence of an increasing NO2/NOX emissions ratio from road traffic emissions , 2005 .

[3]  J. Eichhorn Validation of a Microscale Pollution Dispersal Model , 1996 .

[4]  Margaret Bell,et al.  Risks of exceeding the hourly EU limit value for nitrogen dioxide resulting from increased road transport emissions of primary nitrogen dioxide , 2007 .

[5]  R. A. Cox,et al.  Evaluated kinetic and photochemical data for atmospheric chemistry: Volume III - gas phase reactions of inorganic halogens , 2006 .

[6]  A. O'Hagan,et al.  Probabilistic sensitivity analysis of complex models: a Bayesian approach , 2004 .

[7]  R.A.S. Ratchliffe,et al.  Book reviewBoundary-layer meteorology: A new quarterly journal of physical and biological processes in the atmospheric boundary layer. R. E. Munn (Editor). D. Reidel Publishing Company, Dordrecht, 1969, Volume 1(1). , 1971 .

[8]  A. Reynolds Comments on the 'Universality of the Lagrangian Velocity Structure Function Constant (C0) ACROSS DIFFERENT KINDS OF TURBULENCE’ , 1998, Boundary-Layer Meteorology.

[9]  E. D’Asaro,et al.  The Kolmogorov constant for the Lagrangian velocity spectrum and structure function , 2002 .

[10]  D. Thomson Criteria for the selection of stochastic models of particle trajectories in turbulent flows , 1987, Journal of Fluid Mechanics.

[11]  D. Anfossi,et al.  Estimation of the Lagrangian Velocity Structure Function Constant C0 by Large-Eddy Simulation , 2006 .

[12]  Alison S. Tomlin,et al.  Global sensitivity analysis of a 3D street canyon model—Part II: Application and physical insight using sensitivity analysis , 2008 .

[13]  Roger Atkinson,et al.  Evaluated kinetic and photochemical data for atmospheric chemistry: Volume I - gas phase reaxtions of Ox, HOx, NOx and SOx species , 2004 .

[14]  H. Rabitz,et al.  Efficient input-output model representations , 1999 .

[15]  B. L. Sawford,et al.  Micro-Mixing Modelling of Scalar Fluctuations for Plumes in Homogeneous Turbulence , 2004 .

[16]  T. Ziehn,et al.  A global sensitivity study of sulfur chemistry in a premixed methane flame model using HDMR , 2008 .

[17]  Peter C. Young,et al.  State Dependent Parameter metamodelling and sensitivity analysis , 2007, Comput. Phys. Commun..

[18]  David J. Thomson,et al.  Lagrangian modelling of plume chemistry for secondary pollutants in large industrial plumes , 2008 .

[19]  S. Du UNIVERSALITY OF THE LAGRANGIAN VELOCITY STRUCTURE FUNCTION CONSTANT (C0) ACROSS DIFFERENT KINDS OF TURBULENCE , 1997 .

[20]  Alison S. Tomlin,et al.  A Lagrangian stochastic model for predicting concentration fluctuations in urban areas , 2007 .

[21]  Robert W. Bilger,et al.  An experimental study of a reactive plume in grid turbulence , 1996, Journal of Fluid Mechanics.

[22]  Michael Schatzmann,et al.  Wind tunnel measurements of concentration fluctuations in an urban street canyon , 1999 .

[23]  S. Madronich,et al.  Actinometer and Eppley radiometer measurements of the NO2 photolysis rate coefficient during the Mauna Loa Observatory photochemistry experiment , 1992 .

[24]  A. Robins,et al.  Concentration fluctuations and fluxes in plumes from point sources in a turbulent boundary layer , 1982, Journal of Fluid Mechanics.

[25]  Epaminondas Mastorakos,et al.  Micromixing effects in a reacting plume by the Stochastic Fields method , 2006 .

[26]  Alison S. Tomlin,et al.  GUI-HDMR - A software tool for global sensitivity analysis of complex models , 2009, Environ. Model. Softw..

[27]  Herschel Rabitz,et al.  Global uncertainty assessments by high dimensional model representations (HDMR) , 2002 .

[28]  Adrian Sandu,et al.  Benchmarking stiff ode solvers for atmospheric chemistry problems II: Rosenbrock solvers , 1997 .

[29]  Alison S. Tomlin,et al.  Global sensitivity analysis of a 3D street canyon model—Part I: The development of high dimensional model representations , 2008 .

[30]  I. Sobol On the distribution of points in a cube and the approximate evaluation of integrals , 1967 .

[31]  H. Rabitz,et al.  Practical Approaches To Construct RS-HDMR Component Functions , 2002 .

[32]  Pasquale Franzese,et al.  A PDF micromixing model of dispersion for atmospheric flow. Part I: development of the model, application to homogeneous turbulence and to neutral boundary layer , 2005 .