Finite element simulation of a local scale air quality model over complex terrain

Oliver, Albert; Montero, G.; Montenegro, R.; Rodríguez, Eduardo; Escobar, J. M.; Foguet, Agustí Pérez

doi:10.5194/asr-8-105-2012

Cited by 4 publications

(4 citation statements)

References 40 publications

(45 reference statements)

Supporting

Mentioning

Contrasting

Order By: Relevance

“…For example, in certain pollution models a system of diffusion-advection-reaction equations have to be solved for a given velocity field (the velocity field is determined by analytic modeling, measurements, or computer simulation). In this case the stream function corresponding to the velocity field takes on the role of ψ [17]. A popular choice for tokamak magnetic fields are so-called flux or magnetic coordinates [18], in which the magnetic flux ψ or a function of it is the first coordinate.…”

Section: Introductionmentioning

confidence: 99%

Streamline integration as a method for two-dimensional elliptic grid generation

Wiesenberger

Held

Einkemmer

2017

Journal of Computational Physics

View full text Add to dashboard Cite

We propose a new numerical algorithm to construct a structured numerical elliptic grid of a doubly connected domain. Our method is applicable to domains with boundaries defined by two contour lines of a two-dimensional function. The resulting grids are orthogonal to the boundary. Grid points as well as the elements of the Jacobian matrix can be computed efficiently and up to machine precision. In the simplest case we construct conformal grids, yet with the help of weight functions and monitor metrics we can control the distribution of cells across the domain. Our algorithm is parallelizable and easy to implement with elementary numerical methods. We assess the quality of grids by considering both the distribution of cell sizes and the accuracy of the solution to elliptic problems. Among the tested grids these key properties are best fulfilled by the grid constructed with the monitor metric approach.

show abstract

Section: Introductionmentioning

confidence: 99%

Streamline integration as a method for two-dimensional elliptic grid generation

Wiesenberger

Held

Einkemmer

2017

Journal of Computational Physics

View full text Add to dashboard Cite

show abstract

“…Using the resulting wind field, we simulate the transport and reaction of pollutants. A finite element model designed to simulate the microscale() is used in this example. This model uses a finite element method stabilized with least squares to prevent the numerical oscillations that arise with the classical Galerkin formulation.…”

Section: Examplesmentioning

confidence: 99%

Insertion of triangulated surfaces into a meccano tetrahedral discretization by means of mesh refinement and optimization procedures

Ruiz‐Gironés

Oliver

Socorro-Marrero

et al. 2017

Numerical Meth Engineering

View full text Add to dashboard Cite

This is the peer reviewed version of the following article: Ruiz Gironès , E., Oliver , A., Socorro, G., Cascón, J., Escobar, J.M., Montenegro, R., Sarrate, J. Insertion of triangulated surfaces into a meccano tetrahedral discretization by means of mesh refinement and optimization procedures. "International journal for numerical methods in engineering", 2018, which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5706/pdf. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.In this paper, we present a new method for inserting several triangulated surfaces into an existing tetrahedral mesh generated by the meccano method. The result is a conformal mesh where each inserted surface is approximated by a set of faces of the final tetrahedral mesh. First, the tetrahedral mesh is refined around the inserted surfaces to capture their geometric features. Second, each immersed surface is approximated by a set of faces from the tetrahedral mesh. Third, following a novel approach, the nodes of the approximated surfaces are mapped to the corresponding immersed surface. Fourth, we untangle and smooth the mesh by optimizing a regularized shape distortion measure for tetrahedral elements in which we move all the nodes of the mesh, restricting the movement of the edge and surface nodes along the corresponding entity they belong to. The refining process allows approximating the immersed surface for any initial meccano tetrahedral mesh. Moreover, the proposed projection method avoids computational expensive geometric projections. Finally, the applied simultaneous untangling and smoothing process delivers a high-quality mesh and ensures that the immersed surfaces are interpolated. Several examples are presented to assess the properties of the proposed method.Peer ReviewedPostprint (author's final draft

show abstract

“…The RIVAD/ARM3 scheme is a simplified model that predicts the sulfate and nitrate production rates assuming a steady state concentration of the hydroxil radical. This model considers four species and production rates, which are defined as :

r_{{SO}_{2}} = - r_{{SO}_{4}^{2 -}} = α_{1} (u) u_{{SO}_{2}} = \frac{- γ_{1}}{u_{{SO}_{2}} + δ_{1} u_{{NO}_{x}}} u_{{SO}_{2}}

r_{{NO}_{x}} = - r_{{NO}_{3}} = α_{2} (u) u_{{NO}_{x}} = \frac{- γ_{2}}{δ_{2} u_{{SO}_{2}} + u_{{NO}_{x}}} u_{{NO}_{x}}

where δ 1 , δ 2 , γ 1 , and γ 2 are the four constants. The initial and the inflow boundary conditions are zero for all species, and the emission rate of the first primary pollutant is bigger than the second.…”

Section: Test Examplesmentioning

confidence: 99%

“…The RIVAD/ARM3 scheme [34] is a simplified model that predicts the sulfate and nitrate production rates assuming a steady state concentration of the hydroxil radical. This model considers four species and production rates, which are defined as [16,35]:…”

Section: Rivad/arm3 Testsmentioning

confidence: 99%

A multimesh adaptive scheme for air quality modeling with the finite element method

Monforte

Foguet

2013

Numerical Methods in Fluids

View full text Add to dashboard Cite

Abstract. A multi-mesh adaptive scheme for convection-diffusion-reaction problems is presented. The proposal is applied to air quality modeling, especifically to the simulation of a pollutant punctual emissions. The performance of the proposal is analyzed with different nonlinear reaction models, including the photochemical model CB05 implmented within the Comunity Multiscale Air Quality model, which involves sixty-two species and very different characteristic reaction times. The problem is solved with splitting of transport and reaction processes. This allows to discretize the species in distinct computational meshes, adapted to the distribution of the error indicator of each case. A common reference mesh is used for all species and during all problem evolution. A remeshing technique based on imposing the volume of new elements is used to define and update the computational meshes. An error indicator well suited for problems involving large variation of the unknowns is used. A single-mesh strategy, with remeshing adapted to the most demanding specie in each part of the domain, is used for comparison. The results of the examples presented show that the accuracy of single and multi-mesh strategies are similar. Instead, computational cost of multi-mesh is lower than single-mesh in most cases. Reduction increases with the number of species and the number of plumes. An example of a punctual emissor in a three-dimensional domain, with realistic values of CB05 components, is presented.

show abstract

Finite element simulation of a local scale air quality model over complex terrain

Cited by 4 publications

References 40 publications

Streamline integration as a method for two-dimensional elliptic grid generation

Streamline integration as a method for two-dimensional elliptic grid generation

Insertion of triangulated surfaces into a meccano tetrahedral discretization by means of mesh refinement and optimization procedures

A multimesh adaptive scheme for air quality modeling with the finite element method

Contact Info

Product

Resources

About