Adaptive finite element simulation of stack pollutant emissions over complex terrains
TL;DR: In this article, a three-dimensional finite element model for the pollutant dispersion is presented, where the first stage consists on the construction of an adaptive tetrahedral mesh of a rectangular region bounded in its lower part by the terrain and in its upper part by a horizontal plane.
About: This article is published in Energy.The article was published on 2013-01-01 and is currently open access. It has received 26 citations till now. The article focuses on the topics: Finite element method & Discretization.
Citations
More filters
TL;DR: In this paper, the authors presented a numerical model to study pollutant dispersion in the atmospheric boundary layer (ABL) by accounting for the mechanisms of advection by the mean wind in the horizontal direction, turbulent diffusion in the vertical direction to ground surface, dry deposition, and radioactive decay.
20 citations
TL;DR: A new method for wind field forecasting over complex terrain using the predictions of the HARMONIE meso-scale model as the input data for an adaptive finite element mass-consistent wind model with a minimal user intervention.
Abstract: In this paper, we introduce a new method for wind field forecasting over complex terrain. The main idea is to use the predictions of the HARMONIE meso-scale model as the input data for an adaptive finite element mass-consistent wind model. The HARMONIE results (obtained with a maximum resolution of about 1 km) are refined in a local scale (about a few metres). An interface between both models is implemented in such a way that the initial wind field is obtained by a suitable interpolation of the HARMONIE results. Genetic algorithms are used to calibrate some parameters of the local wind field model in accordance to the HARMONIE data. In addition, measured data are considered to improve the reliability of the simulations. An automatic tetrahedral mesh generator, based on the meccano method, is applied to adapt the discretization to complex terrains. The main characteristic of the framework is a minimal user intervention. The final goal is to validate our model in several realistic applications on Gran Canaria island, Spain, with some experimental data obtained by the AEMET in their meteorological stations. The source code of the mass-consistent wind model is available online at http://www.dca.iusiani.ulpgc.es/Wind3D/
.
18 citations
TL;DR: The method mixes the benefits resulting from isogeometric analysis, implicit dynamics, residual minimization, and alternating direction solver with direction splitting, and shows that the resulting system of linear equations has a Kronecker product structure, which results in a linear computational cost of the direct solver.
Abstract: In this paper, we propose a novel computational implicit method, which we call Isogeometric Residual Minimization (iGRM) with direction splitting. The method mixes the benefits resulting from isogeometric analysis, implicit dynamics, residual minimization, and alternating direction solver. We utilize tensor product B-spline basis functions in space, implicit second order time integration schemes, residual minimization in every time step, and we exploit Kronecker product structure of the matrix to employ linear computational cost alternating direction solver. We implement an implicit time integration scheme and apply, for each space-direction, a stabilized mixed method based on residual minimization. We show that the resulting system of linear equations has a Kronecker product structure, which results in a linear computational cost of the direct solver, even using implicit time integration schemes together with the stabilized mixed formulation. We test our method on three advection–diffusion computational examples, including model “membrane” problem, the circular wind problem, and the simulations modeling pollution propagating from a chimney.
16 citations
Cites methods from "Adaptive finite element simulation ..."
...We use the pollution model based on [39], with the following simplifying assumptions: we consider only one component of the pollution vector field, we neglect the chemical interactions between different components, and we assume cube shape of the computational domain....
[...]
TL;DR: In this paper, a multimesh adaptive scheme for convection-diffusion-reaction problems for a large number of components is presented, where the evaluation of the nonreactive part for each component and the reaction at each node constitute independent tasks.
Abstract: SUMMARY
A multimesh adaptive scheme for convection–diffusion–reaction problems for a large number of components is presented. The problem is solved by splitting transport and reaction processes. This way, the evaluation of the nonreactive part for each component and the reaction at each node constitute independent tasks. This allows to discretize each component of the solution on a distinct computational mesh, adapted on the basis of its error indicator. The standard single-mesh strategy is used for comparison. Simulations of a point emission in a 3D domain are presented. Low remeshing periods of the adaptive scheme are found to be optimal, in terms of computational cost and accuracy, for the nonreactive problem. Examples with several reaction terms, with an increase of the complexity, are then presented. Results show that the accuracy of single-mesh and multimesh strategies are similar. Instead, the computational cost of the multimesh strategy is lower than the single-mesh in the majority of the examples; this process is controlled by the stiff behavior of the reactive term. The problem size of the multimesh scheme is much lower, and therefore, larger spatial discretizations can be simulated for a given available memory. The efficiency of the multimesh strategy increases with the number of species and the number of species that develop a plume. Finally, an example of a punctual emission considering realistic values of the initial concentrations and using the Community Multiscale Air Quality-CBO5 reaction model, which involves 62 components, is presented; the small-scale structure of the different nitrogen components near the emitter is captured. Copyright © 2013 John Wiley & Sons, Ltd.
11 citations
Cites background from "Adaptive finite element simulation ..."
...The extension to complex geometries (topography and build elements) [16], real wind fields [17], and more realistic stack emissions models [16,18] is left for future developments....
[...]
...This model considers four species and production rates, which are defined as [16, 35]:...
[...]
References
More filters
TL;DR: The literature is now becoming overpopulated with plume rise formulas, and, although many more data are now available, comparisons of formulas and data have not kept up as discussed by the authors.
Abstract: While reliable plume rise formulas have long been needed for dispersion calculation, up to ten years ago there were very few data to test the half dozen or so published formulas. The literature is now becoming overpopulated with plume rise formulas, and, although many more data are now available, comparisons of formulas and data have not kept up. The author undertook a more comprehensive set of comparisons as part of a state-of-the-art report on plume rise for the Nuclear Safety Information Center of the Oak Ridge National Laboratory. Some of the results are summarized here.
33 citations
"Adaptive finite element simulation ..." refers background in this paper
...The trajectory of the plume rise has been widely studied in the past [52, 53, 54]....
[...]
TL;DR: The problem of obtaining an incompressible wind field that adjusts to an experimental one is solved by two different adaptive finite element methods, using linear interpolation of a potential function and a mixed finite element method.
33 citations
21 Apr 2002
TL;DR: A code for generating a tetrahedral mesh from an "optimal" node distribution in the domain based on Delaunay triangulation and a procedure to optimise the resulting mesh is developed.
Abstract: In the finite element simulation of environmental processes that occur in a three-dimensional domain defined over an irregular terrain, a mesh generator capable of adapting itself to the topographic characteristics is essen tial. The present study develops a code for generating a tetrahedral mesh from an "optimal" node distribution in the domain. The main ideas for the construction of the initial mesh combine the use of a refinement/derefinement algorithm for two-dimensional domains and a tetrahedral mesh generator algorithm based on Delaunay triangulation. Moreover, we propose a procedure to optimise the resulting mesh. A function to define the vertical distance betw een nodes distributed in the domain is also analysed. Finally, these tec hniques are applied to the construction of meshes adapted to the topography of the southern section of La Palma (Canary Islands).
32 citations
"Adaptive finite element simulation ..." refers background in this paper
...A detailed description of the mesh generation procedure can be seen in [44, 45]....
[...]
01 Nov 2003
TL;DR: This work optimize the sequential program of the STEM-II (Sulphur Transport Eulerian Model 2) program, a large-scale pollution modeling application, with the aim of increasing data locality, and parallelized the program using OpenMP directives for shared memory systems, and the MPI library for distributed memory machines.
Abstract: The aim of this work is to provide a high performance air quality simulation using the STEM-II (Sulphur Transport Eulerian Model 2) program, a large-scale pollution modeling application First, we optimize the sequential program with the aim of increasing data locality Then, we parallelized the program using OpenMP directives for shared memory systems, and the MPI library for distributed memory machines Performance results are presented for a SGI O2000 multiproccessor, a Fujitsu AP3000 multicomputer and a Cluster of PCs Experimental results show that the parallel versions of the code achieve important reductions in the CPU time needed by each simulation This will allow us to obtain results with adequate speed and reliability for the industrial environment where it is intended to be applied
30 citations
"Adaptive finite element simulation ..." refers methods in this paper
...The convection, diffusion and reaction problem is usually solved using splitting schemes [31, 32] and specific numerical solvers for time integration of photochemical reaction terms [33, 34, 23]....
[...]
TL;DR: In this article, the comprehensive Regional Atmospheric Modeling Systems (RAMS) model, coupled online with a chemical module including 29 species, is used to follow the chemistry of a polluted zone over Southern France.
30 citations
"Adaptive finite element simulation ..." refers background or methods in this paper
...However, because of the great coupling among key components as ozone, nitrogen compounds, and Volatile Organic Components [15, 16], and due to the awareness about the socio-economic impacts of their immission [17, 18, 19], some references about the need of coupling local emissions using regional planning with Plume in Grid photochemical models can be found [20, 9, 10]....
[...]
...The modelling domains usually vary from about few thousands to tens of kilometres, and the grid spatial discretization varies from about tens of kilometres to one kilometre using nested sub-grids [9, 10, 11]....
[...]