The meccano method for isogeometric analysis of planar domains

[EN]The authors have recently introduced the meccano method for tetrahedral mesh generation and volume parameterization of solids. In this paper, we present advances of the method for T-spline modelling and analysis of complex geometries. We consider a planar domain composed by several irregular sub-domains. These sub-regions are defined by their boundaries and can represent different materials. The bivariate T-spline representation of the whole physical domain is constructed from a square. In this procedure, a T-mesh optimization method is crucial. We show results of an elliptic problem by using a quadtree local T-mesh refinement technique…

Uso de métodos ensemble para la predicción de la radiación solar

[EN]Se propone un modelo de radiación solar adaptativo como una nueva herramienta para la generación de mapas de radiación solar. Este introduce mejoras a los modelos existentes como la adaptación de la malla a la orografía y al albedo. Esta estrategia adaptativa nos permite generar un código eficiente que reduce el coste computacional para una precisión dada. La radiación global es obtenida como suma de sus tres componentes, la directa, la difusa y la reflejada, sobre una región de estudio bajo condiciones de cielo limpio. En este sentido, las superficies inclinadas tendrán un tratamiento diferente de las horizontales y se tendrá en cuenta el efecto de las zonas en sombra…

Desarrollo de un método ensemble para la predicción del viento a escala local usando elementos finitos

[EN]Los autores han desarrollado un modelo de viento de masa consistente especialmente diseñado para su aplicación en la escala local y en zonas de orografía compleja. Se ha dotado a este modelo de carácter predictivo usando como entrada resultados de HARMONIE. El HARMONIE es un modelo meteorológico predictivo de escala regional usado en la AEMET. Por otra parte, en los últimos años los métodos ensemble se han consolidado en la predicción meteorológica a escala regional…

Application of genetic algorithms for the calibration of an air quality model and its validation using pollutant measures from the surroundings of an electric power plant

[EN]This work presents the calibration and validation of an air quality finite element model applied to emissions from a thermal power plant located in Gran Canaria. The calibration is performed using genetic algorithms. To calibrate and validate the model, the authors use empirical measures of pollutants concentrations from 4 stations located nearby the power plant; an hourly record per station during 3 days is available. Measures from 3 stations will be used to calibrate, while validation will use measures from the remaining station…

Calibration and validation of an air quality finite element model around an electric power plant in Gran Canaria island

[EN]This work presents the calibration and validation of an air quality finite element model applied to the surroundings of Jinamar electric power plant in Gran Canaria island (Spain). The model involves the generation of an adaptive tetrahedral mesh, the computation of an ambient wind field, the inclusion of the plume rise effect in the wind field, and the simulation of transport and reaction of pollutants. The main advantage of the model is the treatment of complex terrains that introduces an alternative to the standard implementation of current models. In addition, it improves the computational cost through the use of unstructured meshes...

T-spline parameterization of 2D geometries based on the meccano method with a new T-mesh optimization algorithm

[EN]We present a new method, based on the idea of the meccano method and a novel T-mesh optimization procedure, to construct a T-spline parameterization of 2D geometries for the application of isogeometric analysis. The proposed method only demands a boundary representation of the geometry as input data. The algorithm obtains, as a result, high quality parametric transformation between 2D objects and the parametric domain, the unit square. First, we define a parametric mapping between the input boundary of the object and the boundary of the parametric domain. Then, we build a T-mesh adapted to the geometric singularities of the domain in order to preserve the features of the object boundary with a desired tolerance…

A new predictive solar radiation numerical model

[EN]A predictive solar radiation numerical model is presented. Starting from the works of, a solar radiation numerical model is developed considering the terrain surface through 2-D adaptive meshes of triangles which are constructed using a refinement/derefinement procedure in accordance with the variations of terrain surface and albedo. The effect of shadows is considered in each time step. Solar radiation is first computed for clear-sky (CS) conditions and then, real-sky values are computed daily in terms of the CS index computed using all the observational data which are available for each day at several points of the studied zone…

Tetrahedral mesh optimization combining boundary and inner node relocation and adaptive local refinement

[EN]This work introduces a new technique for tetrahedral mesh optimization. The procedure relocates boundary and inner nodes without changing the mesh topology. In order to maintain the boundary approximation while boundary nodes are moved, a local refinement of tetrahedra with faces on the solid boundary is necessary in some cases. New nodes are projected on the boundary by using a surface parameterization. In this work, the proposed method is applied to tetrahedral meshes of genus-zero solids that are generated by the meccano method. In this case, the solid boundary is automatically decomposed into six surface patches which are parameterized into the six faces of a cube with the Floater parameterization...

Adaptive Finite Element Method with a local element parametrization nested meses

[EN]This work presents a novel approach to solve a two dimensional problem by using an adaptive finite element approach. The most common strategy to deal with nested adaptivity is to generate a mesh that represents the geometry and the input parameters correctly, and to refine this mesh locally to obtain the most accurate solution. As opposed to this approach, the authors propose a technique using independent meshes : geometry, input data and the unknowns. Each particular mesh is obtained by a local nested refinement of the same coarse mesh at the parametric space…

El método del mecano y análisis isogeométrico

[EN]En este trabajo presentamos avances del método del mecano para la generación de mallas de tetraedros y la parametrización volumétrica de sólidos. El método combina varios procedimientos: una transformación desde la frontera del mecano (espacio paramétrico sencillo) a la superficie del sólido (espacio físico), un algoritmo de refinamiento local de mallas de tetraedros, y un desenredo y suavizado simultáneo de mallas. La base del método consiste en la definición de una transformación volumétrica biyectiva entre el mecano y el sólido. Se mostrarán resultados usando elementos finitos adaptativos en varios problemas medioambientales : campos de viento, contaminación atmósférica, etc…

Construction of polynomial spline spaces over quadtree and octree T-meshes

[EN]We present a new strategy for constructing tensor product spline spaces over quadtree and octree T-meshes. The proposed technique includes some simple rules for inferring local knot vectors to define spline blending functions. These rules allow to obtain for a given T-mesh a set of cubic spline functions that span a space with nice properties: it can reproduce cubic polynomials, the functions are C2-continuous, linearly independent, and spaces spanned by nested T-meshes are also nested. In order to span spaces with these properties applying the proposed rules, the T-mesh should fulfill the only requirement of being a 0-balanced quadtree or octree. ..

Ensemble wind forecasting based on the HARMONIE model and adaptive finite elements in complex orography

[EN]Ensemble forecasting is a methodology to deal with uncertainties in the numerical wind prediction. In this work we propose to apply ensemble methods to the adaptive wind forecasting model presented in. The wind field forecasting is based on a mass-consistent model and a log-linear wind profile using as input data the resulting forecast wind from Harmonie, a Non-Hydrostatic Dynamic model used experimentally at AEMET with promising results. The mass-consistent model parameters are estimated by using genetic algorithms. The mesh is generated using the meccano method and adapted to the geometry…

Local scale air quality model with several pollutant sources

[EN]This paper shows a finite element method for pollutant transport with several pollutant sources. An Eulerian convection–diffusion–reaction model to simulate the pollutant dispersion is used. The discretization of the different sources allows to impose the emissions as boundary conditions. The Eulerian description can deal with the coupling of several plumes. An adaptive stabilized finite element formulation, specifically Least-Squares, with a Crank-Nicolson temporal integration is proposed to solve the problem. An splitting scheme has been used to treat separately the transport and the reaction. A mass-consistent model has been used to compute the wind field of the problem…

3-D wind field simulation over complex terrain

[EN]A new methodology for wind field simulation or forecasting over complex terrain is introduced. The idea is to use wind measurements or predictions of the HARMONIE mesoscale model as the input data for an adaptive finite element mass consistent wind model. The method has been recently implemented in the freely-available Wind3D code. A description of the HARMONIE Non-Hydrostatic Dynamics can be found in. HARMONIE provides wind prediction with a maximum resolution about 1 Km that is refined by the finite element model in a local scale (about a few meters). An interface between both models is implemented such that the initial wind field approximation is obtained by a suitable interpolation of the HARMONIE results…

Surface insertion in a tetrahedral mesh using the meccano method

[EN]This work presents an innovative method to insert an open surface in a tetrahedral mesh. The insertion of a surface in a mesh can be done with 2 different approaches: introduce the surface to the geometry before generating the mesh, or insert the surface once the mesh is generated. This work uses the second approach. Essentially, the surface is first approximated by a set of faces of the existing mesh. This set is refined to obtain a more accurate approximation. Finally, the set is processed to satisfy some topological properties and projected to the actual surface. The strategy is based on a mesh generated by the Meccano Method…