Desarrollo, implementación y análisis de un modelo cinemático local de velocidades para Costa Rica basado en la integración y homogenización de datos de estaciones GNSS continuas

El primer procesamiento estricto realizado con el software científico Bernese y contemplando las más estrictas normas de cálculo recomendadas internacionalmente, permitió obtener un campo puntual de alta exactitud, basado en la integración y estandarización de los datos de una red GPS ubicada en Costa Rica. Este procesamiento contempló un total de 119 semanas de datos diarios, es decir unos 2,3 años, desde enero del año 2009 hasta abril del año 2011, para un total de 30 estaciones GPS, de las cuales 22 están ubicadas en el territorio nacional de Costa Rica y 8 internaciones pertenecientes a la red del Sistema Geocéntrico para las Américas (SIRGAS). Las denominadas soluciones semilibres generaron, semana a semana, una red GPS con una alta exactitud interna definida por medio de los vectores entre las estaciones y las coordenadas finales de la constelación satelital. La evaluación semanal dada por la repetibilidad de las soluciones brindó en promedio errores de 1,7 mm, 1,4 mm y 5,1 mm en las componentes [n e u], confirmando una alta consistencia en estas soluciones. Aunque las soluciones semilibres poseen una alta exactitud interna, las mismas no son utilizables para fines de análisis cinemático, pues carecen de un marco de referencia. En Latinoamérica, la densificación del Marco Internacional Terrestre de Referencia (ITRF), está representado por la red de estaciones de operación continua GNSS de SIRGAS, denominada como SIRGAS-CON. Por medio de las denominadas coordenadas semanales finales de las 8 estaciones consideradas como vínculo, se refirió cada una de las 119 soluciones al marco SIRGAS. La introducción del marco de referencia SIRGAS a las soluciones semilibres produce deformaciones en estas soluciones. Las deformaciones de las soluciones semilibres son producto de las cinemática de cada una de las placas en las que se ubican las estaciones de vínculo. Luego de efectuado el amarre semanal a las coordenadas SIRGAS, se hizo una estimación de los vectores de velocidad de cada una de las estaciones, incluyendo las de amarre, cuyos valores de velocidad se conocen con una alta exactitud. Para la determinación de las velocidades de las estaciones costarricenses, se programó una rutina en ambiente MatLab, basada en una ajuste por mínimos cuadrados. Los valores obtenidos en el marco de este proyecto en comparación con los valores oficiales, brindaron diferencias promedio del orden de los 0,06 cm/a, -0,08 cm/a y -0,10 cm/a respectivamente para las coordenadas [X Y Z]. De esta manera se logró determinar las coordenadas geocéntricas [X Y Z]T y sus variaciones temporales [vX vY vZ]T para el conjunto de 22 estaciones GPS de Costa Rica, dentro del datum IGS05, época de referencia 2010,5. Aunque se logró una alta exactitud en los vectores de coordenadas geocéntricas de las 22 estaciones, para algunas de las estaciones el cálculo de las velocidades no fue representativo debido al relativo corto tiempo (menos de un año) de archivos de datos. Bajo esta premisa, se excluyeron las ocho estaciones ubicadas al sur de país. Esto implicó hacer una estimación del campo local de velocidades con solamente veinte estaciones nacionales más tres estaciones en Panamá y una en Nicaragua. El algoritmo usado fue el denominado Colocación por Mínimos Cuadrados, el cual permite la estimación o interpolación de datos a partir de datos efectivamente conocidos, el cual fue programado mediante una rutina en ambiente MatLab. El campo resultante se estimó con una resolución de 30' X 30' y es altamente constante, con una velocidad resultante promedio de 2,58 cm/a en una dirección de 40,8° en dirección noreste. Este campo fue validado con base en los datos del modelo VEMOS2009, recomendado por SIRGAS. Las diferencias de velocidad promedio para las estaciones usadas como insumo para el cálculo del campo fueron del orden los +0,63 cm/a y +0,22 cm/a para los valores de velocidad en latitud y longitud, lo que supone una buena determinación de los valores de velocidad y de la estimación de la función de covarianza empírica, necesaria para la aplicación del método de colocación. Además, la grilla usada como base para la interpolación brindó diferencias del orden de -0,62 cm/a y -0,12 cm/a para latitud y longitud. Adicionalmente los resultados de este trabajo fueron usados como insumo para hacer una aproximación en la definición del límite del llamado Bloque de Panamá dentro del territorio nacional de Costa Rica. El cálculo de las componentes del Polo de Euler por medio de una rutina programa en ambiente MatLab y aplicado a diferentes combinaciones de puntos no brindó mayores aportes a la definición física de este límite. La estrategia lo que confirmó fue simplemente la diferencia en la dirección de todos los vectores velocidad y no permitió reveló revelar con mayor detalle una ubicación de esta zona dentro del territorio nacional de Costa Rica. ABSTRACT The first strict processing performed with the Bernese scientific software and contemplating the highest standards internationally recommended calculation, yielded a precise field of high accuracy, based on the integration and standardization of data from a GPS network located in Costa Rica. This processing watched a total of 119 weeks of daily data, is about 2.3 years from January 2009 to April 2011, for a total of 30 GPS stations, of which 22 are located in the country of Costa Rica and 8 hospitalizations within the network of Geocentric System for the Americas (SIRGAS). The semi-free solutions generated, every week a GPS network with high internal accuracy defined by vectors between stations and the final coordinates of the satellite constellation. The weekly evaluation given by repeatability of the solutions provided in average errors of 1.7 mm 1.4 mm and 5.1 mm in the components [n e u], confirming a high consistency in these solutions. Although semi-free solutions have a high internal accuracy, they are not used for purposes of kinematic analysis, because they lack a reference frame. In Latin America, the densification of the International Terrestrial Reference Frame (ITRF), is represented by a network of continuously operating GNSS stations SIRGAS, known as SIRGAS-CON. Through weekly final coordinates of the 8 stations considered as a link, described each of the solutions to the frame 119 SIRGAS. The introduction of the frame SIRGAS to semi-free solutions generates deformations. The deformations of the semi-free solutions are products of the kinematics of each of the plates in which link stations are located. After SIRGAS weekly link to SIRGAS frame, an estimate of the velocity vectors of each of the stations was done. The velocity vectors for each SIRGAS stations are known with high accuracy. For this calculation routine in MatLab environment, based on a least squares fit was scheduled. The values obtained compared to the official values, gave average differences of the order of 0.06 cm/yr, -0.08 cm/yr and -0.10 cm/yr respectively for the coordinates [XYZ]. Thus was possible to determine the geocentric coordinates [XYZ]T and its temporal variations [vX vY vZ]T for the set of 22 GPS stations of Costa Rica, within IGS05 datum, reference epoch 2010.5. The high accuracy vector for geocentric coordinates was obtained, however for some stations the velocity vectors was not representative because of the relatively short time (less than one year) of data files. Under this premise, the eight stations located in the south of the country were excluded. This involved an estimate of the local velocity field with only twenty national stations plus three stations in Panama and Nicaragua. The algorithm used was Least Squares Collocation, which allows the estimation and interpolation of data from known data effectively. The algorithm was programmed with MatLab. The resulting field was estimated with a resolution of 30' X 30' and is highly consistent with a resulting average speed of 2.58 cm/y in a direction of 40.8° to the northeast. This field was validated based on the model data VEMOS2009 recommended by SIRGAS. The differences in average velocity for the stations used as input for the calculation of the field were of the order of +0.63 cm/yr, +0.22 cm/yr for the velocity values in latitude and longitude, which is a good determination velocity values and estimating the empirical covariance function necessary for implementing the method of application. Furthermore, the grid used as the basis for interpolation provided differences of about -0.62 cm/yr, -0.12 cm/yr to latitude and longitude. Additionally, the results of this investigation were used as input to an approach in defining the boundary of Panama called block within the country of Costa Rica. The calculation of the components of the Euler pole through a routine program in MatLab and applied to different combinations of points gave no further contributions to the physical definition of this limit. The strategy was simply confirming the difference in the direction of all the velocity vectors and not allowed to reveal more detail revealed a location of this area within the country of Costa Rica.

Tenacidad de fractura dinámica de iniciación a diferentes velocidades de solicitación: técnicas experimentales

En esta investigación se presentan diversos dispositivos experimentales diseñados para la realización de ensayos de tenacidad de fractura a diferentes velocidades de solicitación. Se realizan ensayos de flexión en tres puntos de probetas normalizadas (ASTM E399) empleando una máquina servohidráulica, una torre de caída, una barra Hopkinson modificada y un dispositivo especialmente diseñado para trabajar con cargas explosivas. Para poder validar dichos dispositivos experimentales, se realizan estos ensayos sobre un material independiente de la velocidad de deformación; dicho material es la aleación de aluminio 7017-T73. Todos los ensayos realizados sobre dicha aleación dieron el mismo valor de tenacidad de fractura dinámica de iniciación; por lo tanto, los dispositivos experimentales diseñados se consideraron válidos para este fin. Una vez validada la técnica experimental, se llevan a cabo ensayos de flexión en tres puntos en un material dependiente de la velocidad de deformación. Por lo tanto, los dispositivos experimentales diseñados fueron utilizados para determinar la tenacidad de fractura dinámica de iniciación del acero Mars 240. Los resultados obtenidos muestran que, en el caso de este acero, la tenacidad de fractura dinámica de iniciación aumenta con la velocidad de solicitación.

Estudio de la estructura de autofunciones de sistemas caóticos en una base de funciones de scar

El problema central que aborda esta tesis doctoral es el estudio de la correspondencia entre la mecanica clasica y la mecanica cuantica en sistemas hamiltonianos clasicamente caoticos, tema que se enmarca dentro del llamado caos cuantico. En concreto, en este trabajo proponemos un nuevo y efectivo metodo para calcular las autofunciones de sistemas caoticos usando una base de funciones de scar. Dichas funciones de scar juegan un papel fundamental en el estudio de las manifestaciones cuanticas del caos, ya que se trata de funciones de onda semiclasicas con una dispersion muy peque~na y localizadas a lo largo de las variedades invariantes de las orbitas periodicas inestables del sistema que conforman la estructura organizativa del caos clasico. El metodo de calculo desarrollado se ha denominado Metodo de Gram- Schmidt Selectivo (MGSS), dado que construye la base haciendo uso del metodo de Gram{Schmidt convencional pero teniendo en cuenta, ademas, la dispersi on de las funciones de scar y la longitud de la orbita periodica a lo largo de la cual se localizan. El MGSS nos ha permitido calcular con gran precision las 2400 autofunciones con menor energa de un oscilador cuartico altamente caotico con dos grados de libertad acoplados, as como autofunciones muy excitadas en una ventana de energa, utilizando en ambos casos una base mucho mas eciente que las descritas en la literatura. Ademas, hemos empleado el MGSS para calcular las autofunciones del sistema molecular LiNC/LiCN, que presenta un espacio de fases con zonas de regularidad y con regiones en las que el movimiento es altamente caotico; en este ultimo sistema, hemos calculado de forma muy eciente las autofunciones asociadas a las primeras 66 energas. Finalmente, hemos propuesto un metodo perturbativo para calcular velocidades de reaccion en sistemas abiertos descritos por potenciales anarmonicos. Con este metodo, hemos calculado la velocidad de reaccion de distintos potenciales de uno y dos grados de libertad, as como la velocidad de isomerizacion del sistema molecular LiNC/LiCN, tanto sometido a un ruido blanco (sin correlaciones) como en presencia de un ba~no de atomos de argon, lo que constituye un entorno con correlaciones. El metodo desarrollado es independiente de la supercie divisoria y nos ha permitido obtener correcciones analticas a la famosa formula de Kramers, lo que posibilita el calculo exacto de velocidades de reaccion en potenciales anarmonicos que interaccionan con el entorno.

Metodología para la evaluación de la consistencia del trazado de carreteras interurbanas de dos carriles

La consistencia del trazado puede interpretarse como la relación entre las características geométricas de una carretera y lo que espera encontrar el conductor que circula por ella. Si hay una correspondencia entre estos dos aspectos, la conducción puede hacerse de modo continuo, sin sobresaltos, lo que incide favorablemente sobre la seguridad en la circulación. Si bien hay una serie de recomendaciones desde el punto de vista geométrico para obtener trazados consistentes, esto no siempre se logra, y sólo en los últimos años se ha iniciado el estudio de metodologías para evaluar ésto, tanto en vías existentes como en vías proyectadas. La mayor parte de estas metodologías sólo considera el trazado en planta, olvidándose del trazado en alzado y de la coordinación entre los mismos. En esta Tesis doctoral se ha desarrollado una metodología para evaluar la consistencia del trazado en carreteras interurbanas de dos carriles que considera dichos aspectos. Para ello, se hizo un análisis exhaustivo de los índices de trazado, los cuales evalúan las características geométricas en planta y en alzado. Los índices se correlacionaron con la accidentalidad, para determinar cuál de ellos tiene mayor incidencia, encontrándose que es el cambio de curvatura vertical (VCCR); a este índice se le estableció un rango de calificación. Como elemento de evaluación complementario de análisis se seleccionó el perfil de velocidades de operación, procedimiento que ha sido probado en diferentes investigaciones, y del cual se desarrolló un modelo aplicado a Colombia. Para la coordinación de trazados en planta y alzado se evaluaron diferentes combinaciones geométricas, algunas de las cuales generaron reapariciones del trazado. Se ha definido un nuevo índice (Irt) que permite determinar numéricamente la posibilidad de que se presente esta situación, indeseable desde el punto de vista de la seguridad vial. La combinación de estos tres elementos permite una evaluación integral de los diferentes aspectos que inciden sobre la consistencia del trazado de una carretera. La metodología desarrollada se aplicó en el estudio de consistencia del trazado en algunas carreteras españolas y colombianas, ubicadas en distintos tipos de terreno. ABSTRACT Geometric Design Consistency can be defined as the relationship between the geometric characteristics of a road and what the driver expects to find when driving. If there is a correspondence between these two aspects, driving is smoother and unexpected events are minimized, which increases traffic safety conditions. Although from the geometric point of view there are several recommendations to ensure consistent designs, this is not always successfully applied. The study of methods to evaluate design consistency in existing and future routes has only begun in recent years. Most existing methods only consider the horizontal alignment of the road and overlook both the vertical alignment and the coordination that must exist between the vertical and the horizontal. The present Doctoral Thesis proposes a method to evaluate the geometric design consistency of a two-lane rural highway which considers all three of these aspects: the horizontal alignment, the vertical alignment and the coordination that must exist between them. In order to achieve this, several different alignment indices, that evaluate horizontal and vertical geometric characteristics, were thoroughly analyzed to determine their correlation with traffic accidents. The Vertical Curvature Change Rate (VCCR) index showed the highest correlation, and rating thresholds for this index have been established. To complement the evaluation, the operating speed profile, was chosen. This procedure has been extensively tested by several researchers. An operating speed prediction model adapted to Colombia was developed. To study the coordination between the horizontal and the vertical alignments of the road, several geometric combinations of the two were used. Some of these combinations generate undesirable losses of visibility. For this reason, a new index (Irt) was defined to numerically detect those cases, which are undesirable from the point of view of traffic safety. The combination of these three factors allows a comprehensive evaluation of the different aspects that affect the geometric design consistency of a highway. The methodology was applied to some Spanish and Colombian roads located in different types of terrain.

Evaluación de la madera estructural para la rehabilitación de las Naves Cros (Valencia, España) por medio de diferentes técnicas no destructivas

El proyecto de rehabilitación de una de las naves del complejo fabril de la industria química ?CROS? en Valencia se llevó a cabo con el criterio de mantener, en la medida de lo posible, los elementos estructurales presentes en la nave. Con este objetivo se realizaron una serie de ensayos no destructivos (END) in situ. Estos ensayos permitieron evaluar la calidad de la madera, determinar qué elementos estructurales debían ser sustituidos y comprobar la aptitud de los que iban a ser reutilizados. Los END empleados en este estudio fueron los siguientes: (1) Identificación de la especie por técnicas anatómicas, (2) Clasificación resistente por método visual, (3) Estimación de humedad por la técnica de resistencia eléctrica; (4) Obtención de velocidades de propagación ultrasónicas (5) Resistógrafía y (6) Alteración de la propagación de ondas electromagnéticas por medio de Georradar. Para la calibración de estos END se tomó una muestra de piezas y se hicieron ensayos destructivos bajo condiciones controladas en laboratorio. En el trabajo que aquí se presenta se muestra la metodología empleada durante el proceso de toma de datos, de análisis de resultados y de cruce de la información obtenida a partir de cada uno de los ensayos hasta llegar a un diagnóstico para los elementos analizados. The assessment of structural timber was requested in the rehabilitation project of the Naves of the chemical industry "CROS". The criterion was to maintain as much as possible timber of the structure and to make only partial replacements. In order not to damage the existing structure and to assess the quality of the existing timber, a series of non-destructive testing (NDT) in the entire structure were performed: (1) Identification of the species by anatomical techniques, (2) Strength grading by visual method, (3) Estimation of moisture content by the technique of electrical resistance, (4) Acquisition of ultrasonic propagation velocities (5) Resistography and (6) Record of the propagation of electromagnetic waves by means of Ground-penetrating radar. Following, a sample group was chose to carry out destructive testing in the lab and compare the NDT results with those obtained with the standard UNE-EN408 (modules of strength, stiffness and density). In the present work, the results provided by each of the NDT techniques are detailed and above all, what is more important, the validity of these after they have been contrasted with the destructive standard tests.

Comportamiento frente a la corrosión de barras de acero de armar embebidas en morteros con incorporación de escorias blancas de horno cuchara

El objetivo del presente trabajo es estudiar el comportamiento frente a la corrosión de barras de acero de armar embebidas en probetas de mortero, en las que se ha sustituido parcialmente el árido y el cemento por escorias blancas de horno cuchara (LFS) y se ha añadido en el momento del amasado distintos porcentajes de ión cloruro, en peso de cemento. Se han fabricado probetas prismáticas de 6x8x2 cm3, en las que se han embebido 3 barras de acero de diámetro 6 mm, con el fin de comparar el potencial de corrosión y la velocidad de corrosión, en las probetas patrón y en las probetas con escorias. Los resultados alcanzados permiten afirmar que el comportamiento frente a la corrosión inducida por iones cloruro en barras de acero embebidas en probetas patrón y en probetas fabricadas parcialmente con escorias blancas de horno cuchara son semejantes, tanto en las velocidades de corrosión alcanzadas como en el umbral de cloruros que produce la despasivación del acero.

Estudio dinámico de puentes de ferrocarril sometidos a nuevos modelos de carga en el marco europeo

Más allá de la alta velocidad, y debido en parte a su éxito, en los últimos años han aparecido trenes de pasajeros capaces de circular a velocidades más elevadas y ser más eficientes (como, por ejemplo, los de tipo automotor). Además, algunos de los antiguos modelos recogidos en las normativas parecen no cubrir determinados casos, lo que ha llevado a la reciente modificación de éstas. A todo esto se le añade la necesidad de crear una normativa común europea que establezca las bases para lograr la interoperabilidad. Con este proyecto se pretende estudiar los efectos dinámicos en puentes de estas nuevas tipologías de trenes, haciendo especial hincapié en los trenes automotores. El trabajo se centra principalmente en puentes isostáticos de luces cortas, por ser éstos especialmente sensibles a este tipo de acciones. El parámetro crítico, en estos casos, son las aceleraciones en centro de vano. Se realizara el cálculo dinámico por diferentes métodos: elementos finitos, análisis modal, así como programas de integración numérica en Matlab de los modos de vibración, analizando y comparando la efectividad de estos. Las diferencias entre la consideración o no de la distribución de las cargas transmitidas al tablero de la estructura de la vía a serán discutidas en el proyecto, así como la interacción dinámica entre estructura y vehículo. Finalmente se aplicarán los conocimientos obtenidos a varios casos reales de puentes incluidos en la red ferroviaria española, como por ejemplo, el viaducto de alta velocidad del Jalón. Los resultados de estos estudios tienen una repercusión en diversas normativas europeas actualmente en revisión y en las cuales está involucrado el Grupo de Mecánica Computacional: EN15528 European Committee for Standarization (2008), EN1991-2 European Committee for Standarization (2003) y las Especificaciones Técnicas de Interoperabilidad Ferroviaria (TSI).

Estudios de confort de vehículos ferroviarios mediante técnicas combinadas de simulación FEM y MBS

La elasticidad de la caja de viajeros de los vehículos ferroviarios tiene una gran influencia sobre el confort. Por esta razón, cuando se desea simular el comportamiento dinámico del vehículo para estudios de confort, resulta conveniente construir un modelo elástico de caja, a fin de obtener resultados más precisos. La construcción de este tipo de modelos pasa por el desarrollo de dos etapas fundamentales, que comprenden la generación de un modelo de caja mediante la técnica de los elementos finitos (FEM) para su caracterización dinámica como cuerpo elástico y la definición de un modelo de sistema multicuerpo (MBS) que englobe los restantes componentes del vehículo. En este artículo se presentan los resultados obtenidos en un estudio comparativo llevado a cabo con modelos de caja rígida y elástica, en los que se ha valorado el nivel de confort obtenido con ambas configuraciones. Para ello, se ha simulado el comportamiento del vehículo a dos velocidades distintas, de 70 y 110km/h, y con dos niveles de irregularidades. Se han analizado las aceleraciones de la caja, que se han procesado de acuerdo a las especificaciones de la norma EN12299, a fin de obtener el índice de comodidad. Este parámetro se ha utilizado para comparar el nivel de confort obtenido con ambos modelos, habiéndose encontrado una gran diferencia en los índices calculados con caja rígida y con caja elástica, lo que confirma la gran influencia de la elasticidad de la caja en los estudios de confort llevados a cabo mediante técnicas de simulación dinámica.

An Intercomparison of Modelling Approaches for Simulating the Atmospheric Dispersion of Ammonia Emitted by Agricultural Sources

RESUMEN La dispersión del amoniaco (NH3) emitido por fuentes agrícolas en medias distancias, y su posterior deposición en el suelo y la vegetación, pueden llevar a la degradación de ecosistemas vulnerables y a la acidificación de los suelos. La deposición de NH3 suele ser mayor junto a la fuente emisora, por lo que los impactos negativos de dichas emisiones son generalmente mayores en esas zonas. Bajo la legislación comunitaria, varios estados miembros emplean modelos de dispersión inversa para estimar los impactos de las emisiones en las proximidades de las zonas naturales de especial conservación. Una revisión reciente de métodos para evaluar impactos de NH3 en distancias medias recomendaba la comparación de diferentes modelos para identificar diferencias importantes entre los métodos empleados por los distintos países de la UE. En base a esta recomendación, esta tesis doctoral compara y evalúa las predicciones de las concentraciones atmosféricas de NH3 de varios modelos bajo condiciones, tanto reales como hipotéticas, que plantean un potencial impacto sobre ecosistemas (incluidos aquellos bajo condiciones de clima Mediterráneo). En este sentido, se procedió además a la comparación y evaluación de varias técnicas de modelización inversa para inferir emisiones de NH3. Finalmente, se ha desarrollado un modelo matemático simple para calcular las concentraciones de NH3 y la velocidad de deposición de NH3 en ecosistemas vulnerables cercanos a una fuente emisora. La comparativa de modelos supuso la evaluación de cuatro modelos de dispersión (ADMS 4.1; AERMOD v07026; OPS-st v3.0.3 y LADD v2010) en un amplio rango de casos hipotéticos (dispersión de NH3 procedente de distintos tipos de fuentes agrícolas de emisión). La menor diferencia entre las concentraciones medias estimadas por los distintos modelos se obtuvo para escenarios simples. La convergencia entre las predicciones de los modelos fue mínima para el escenario relativo a la dispersión de NH3 procedente de un establo ventilado mecánicamente. En este caso, el modelo ADMS predijo concentraciones significativamente menores que los otros modelos. Una explicación de estas diferencias podríamos encontrarla en la interacción de diferentes “penachos” y “capas límite” durante el proceso de parametrización. Los cuatro modelos de dispersión fueron empleados para dos casos reales de dispersión de NH3: una granja de cerdos en Falster (Dinamarca) y otra en Carolina del Norte (EEUU). Las concentraciones medias anuales estimadas por los modelos fueron similares para el caso americano (emisión de granjas ventiladas de forma natural y balsa de purines). La comparación de las predicciones de los modelos con concentraciones medias anuales medidas in situ, así como la aplicación de los criterios establecidos para la aceptación estadística de los modelos, permitió concluir que los cuatro modelos se comportaron aceptablemente para este escenario. No ocurrió lo mismo en el caso danés (nave ventilada mecánicamente), en donde el modelo LADD no dio buenos resultados debido a la ausencia de procesos de “sobreelevacion de penacho” (plume-rise). Los modelos de dispersión dan a menudo pobres resultados en condiciones de baja velocidad de viento debido a que la teoría de dispersión en la que se basan no es aplicable en estas condiciones. En situaciones de frecuente descenso en la velocidad del viento, la actual guía de modelización propone usar un modelo que sea eficaz bajo dichas condiciones, máxime cuando se realice una valoración que tenga como objeto establecer una política de regularización. Esto puede no ser siempre posible debido a datos meteorológicos insuficientes, en cuyo caso la única opción sería utilizar un modelo más común, como la versión avanzada de los modelos Gausianos ADMS o AERMOD. Con el objetivo de evaluar la idoneidad de estos modelos para condiciones de bajas velocidades de viento, ambos modelos fueron utilizados en un caso con condiciones Mediterráneas. Lo que supone sucesivos periodos de baja velocidad del viento. El estudio se centró en la dispersión de NH3 procedente de una granja de cerdos en Segovia (España central). Para ello la concentración de NH3 media mensual fue medida en 21 localizaciones en torno a la granja. Se realizaron también medidas de concentración de alta resolución en una única localización durante una campaña de una semana. En este caso, se evaluaron dos estrategias para mejorar la respuesta del modelo ante bajas velocidades del viento. La primera se basó en “no zero wind” (NZW), que sustituyó periodos de calma con el mínimo límite de velocidad del viento y “accumulated calm emissions” (ACE), que forzaban al modelo a calcular las emisiones totales en un periodo de calma y la siguiente hora de no-calma. Debido a las importantes incertidumbres en los datos de entrada del modelo (inputs) (tasa de emisión de NH3, velocidad de salida de la fuente, parámetros de la capa límite, etc.), se utilizó el mismo caso para evaluar la incertidumbre en la predicción del modelo y valorar como dicha incertidumbre puede ser considerada en evaluaciones del modelo. Un modelo dinámico de emisión, modificado para el caso de clima Mediterráneo, fue empleado para estimar la variabilidad temporal en las emisiones de NH3. Así mismo, se realizó una comparativa utilizando las emisiones dinámicas y la tasa constante de emisión. La incertidumbre predicha asociada a la incertidumbre de los inputs fue de 67-98% del valor medio para el modelo ADMS y entre 53-83% del valor medio para AERMOD. La mayoría de esta incertidumbre se debió a la incertidumbre del ratio de emisión en la fuente (50%), seguida por la de las condiciones meteorológicas (10-20%) y aquella asociada a las velocidades de salida (5-10%). El modelo AERMOD predijo mayores concentraciones que ADMS y existieron más simulaciones que alcanzaron los criterios de aceptabilidad cuando se compararon las predicciones con las concentraciones medias anuales medidas. Sin embargo, las predicciones del modelo ADMS se correlacionaron espacialmente mejor con las mediciones. El uso de valores dinámicos de emisión estimados mejoró el comportamiento de ADMS, haciendo empeorar el de AERMOD. La aplicación de estrategias destinadas a mejorar el comportamiento de este último tuvo efectos contradictorios similares. Con el objeto de comparar distintas técnicas de modelización inversa, varios modelos (ADMS, LADD y WindTrax) fueron empleados para un caso no agrícola, una colonia de pingüinos en la Antártida. Este caso fue empleado para el estudio debido a que suponía la oportunidad de obtener el primer factor de emisión experimental para una colonia de pingüinos antárticos. Además las condiciones eran propicias desde el punto de vista de la casi total ausencia de concentraciones ambiente (background). Tras el trabajo de modelización existió una concordancia suficiente entre las estimaciones obtenidas por los tres modelos. De este modo se pudo definir un factor de emisión de para la colonia de 1.23 g NH3 por pareja criadora por día (con un rango de incertidumbre de 0.8-2.54 g NH3 por pareja criadora por día). Posteriores aplicaciones de técnicas de modelización inversa para casos agrícolas mostraron también un buen compromiso estadístico entre las emisiones estimadas por los distintos modelos. Con todo ello, es posible concluir que la modelización inversa es una técnica robusta para estimar tasas de emisión de NH3. Modelos de selección (screening) permiten obtener una rápida y aproximada estimación de los impactos medioambientales, siendo una herramienta útil para evaluaciones de impactos en tanto que permite eliminar casos que presentan un riesgo potencial de daño bajo. De esta forma, lo recursos del modelo pueden Resumen (Castellano) destinarse a casos en donde la posibilidad de daño es mayor. El modelo de Cálculo Simple de los Límites de Impacto de Amoniaco (SCAIL) se desarrolló para obtener una estimación de la concentración media de NH3 y de la tasa de deposición seca asociadas a una fuente agrícola. Está técnica de selección, basada en el modelo LADD, fue evaluada y calibrada con diferentes bases de datos y, finalmente, validada utilizando medidas independientes de concentraciones realizadas cerca de las fuentes. En general SCAIL dio buenos resultados de acuerdo a los criterios estadísticos establecidos. Este trabajo ha permitido definir situaciones en las que las concentraciones predichas por modelos de dispersión son similares, frente a otras en las que las predicciones difieren notablemente entre modelos. Algunos modelos nos están diseñados para simular determinados escenarios en tanto que no incluyen procesos relevantes o están más allá de los límites de su aplicabilidad. Un ejemplo es el modelo LADD que no es aplicable en fuentes con velocidad de salida significativa debido a que no incluye una parametrización de sobreelevacion del penacho. La evaluación de un esquema simple combinando la sobreelevacion del penacho y una turbulencia aumentada en la fuente mejoró el comportamiento del modelo. Sin embargo más pruebas son necesarias para avanzar en este sentido. Incluso modelos que son aplicables y que incluyen los procesos relevantes no siempre dan similares predicciones. Siendo las razones de esto aún desconocidas. Por ejemplo, AERMOD predice mayores concentraciones que ADMS para dispersión de NH3 procedente de naves de ganado ventiladas mecánicamente. Existe evidencia que sugiere que el modelo ADMS infraestima concentraciones en estas situaciones debido a un elevado límite de velocidad de viento. Por el contrario, existen evidencias de que AERMOD sobreestima concentraciones debido a sobreestimaciones a bajas Resumen (Castellano) velocidades de viento. Sin embrago, una modificación simple del pre-procesador meteorológico parece mejorar notablemente el comportamiento del modelo. Es de gran importancia que estas diferencias entre las predicciones de los modelos sean consideradas en los procesos de evaluación regulada por los organismos competentes. Esto puede ser realizado mediante la aplicación del modelo más útil para cada caso o, mejor aún, mediante modelos múltiples o híbridos. ABSTRACT Short-range atmospheric dispersion of ammonia (NH3) emitted by agricultural sources and its subsequent deposition to soil and vegetation can lead to the degradation of sensitive ecosystems and acidification of the soil. Atmospheric concentrations and dry deposition rates of NH3 are generally highest near the emission source and so environmental impacts to sensitive ecosystems are often largest at these locations. Under European legislation, several member states use short-range atmospheric dispersion models to estimate the impact of ammonia emissions on nearby designated nature conservation sites. A recent review of assessment methods for short-range impacts of NH3 recommended an intercomparison of the different models to identify whether there are notable differences to the assessment approaches used in different European countries. Based on this recommendation, this thesis compares and evaluates the atmospheric concentration predictions of several models used in these impact assessments for various real and hypothetical scenarios, including Mediterranean meteorological conditions. In addition, various inverse dispersion modelling techniques for the estimation of NH3 emissions rates are also compared and evaluated and a simple screening model to calculate the NH3 concentration and dry deposition rate at a sensitive ecosystem located close to an NH3 source was developed. The model intercomparison evaluated four atmospheric dispersion models (ADMS 4.1; AERMOD v07026; OPS-st v3.0.3 and LADD v2010) for a range of hypothetical case studies representing the atmospheric dispersion from several agricultural NH3 source types. The best agreement between the mean annual concentration predictions of the models was found for simple scenarios with area and volume sources. The agreement between the predictions of the models was worst for the scenario representing the dispersion from a mechanically ventilated livestock house, for which ADMS predicted significantly smaller concentrations than the other models. The reason for these differences appears to be due to the interaction of different plume-rise and boundary layer parameterisations. All four dispersion models were applied to two real case studies of dispersion of NH3 from pig farms in Falster (Denmark) and North Carolina (USA). The mean annual concentration predictions of the models were similar for the USA case study (emissions from naturally ventilated pig houses and a slurry lagoon). The comparison of model predictions with mean annual measured concentrations and the application of established statistical model acceptability criteria concluded that all four models performed acceptably for this case study. This was not the case for the Danish case study (mechanically ventilated pig house) for which the LADD model did not perform acceptably due to the lack of plume-rise processes in the model. Regulatory dispersion models often perform poorly in low wind speed conditions due to the model dispersion theory being inapplicable at low wind speeds. For situations with frequent low wind speed periods, current modelling guidance for regulatory assessments is to use a model that can handle these conditions in an acceptable way. This may not always be possible due to insufficient meteorological data and so the only option may be to carry out the assessment using a more common regulatory model, such as the advanced Gaussian models ADMS or AERMOD. In order to assess the suitability of these models for low wind conditions, they were applied to a Mediterranean case study that included many periods of low wind speed. The case study was the dispersion of NH3 emitted by a pig farm in Segovia, Central Spain, for which mean monthly atmospheric NH3 concentration measurements were made at 21 locations surrounding the farm as well as high-temporal-resolution concentration measurements at one location during a one-week campaign. Two strategies to improve the model performance for low wind speed conditions were tested. These were ‘no zero wind’ (NZW), which replaced calm periods with the minimum threshold wind speed of the model and ‘accumulated calm emissions’ (ACE), which forced the model to emit the total emissions during a calm period during the first subsequent non-calm hour. Due to large uncertainties in the model input data (NH3 emission rates, source exit velocities, boundary layer parameters), the case study was also used to assess model prediction uncertainty and assess how this uncertainty can be taken into account in model evaluations. A dynamic emission model modified for the Mediterranean climate was used to estimate the temporal variability in NH3 emission rates and a comparison was made between the simulations using the dynamic emissions and a constant emission rate. Prediction uncertainty due to model input uncertainty was 67-98% of the mean value for ADMS and between 53-83% of the mean value for AERMOD. Most of this uncertainty was due to source emission rate uncertainty (~50%), followed by uncertainty in the meteorological conditions (~10-20%) and uncertainty in exit velocities (~5-10%). AERMOD predicted higher concentrations than ADMS and more of the simulations met the model acceptability criteria when compared with the annual mean measured concentrations. However, the ADMS predictions were better correlated spatially with the measurements. The use of dynamic emission estimates improved the performance of ADMS but worsened the performance of AERMOD and the application of strategies to improved model performance had similar contradictory effects. In order to compare different inverse modelling techniques, several models (ADMS, LADD and WindTrax) were applied to a non-agricultural case study of a penguin colony in Antarctica. This case study was used since it gave the opportunity to provide the first experimentally-derived emission factor for an Antarctic penguin colony and also had the advantage of negligible background concentrations. There was sufficient agreement between the emission estimates obtained from the three models to define an emission factor for the penguin colony (1.23 g NH3 per breeding pair per day with an uncertainty range of 0.8-2.54 g NH3 per breeding pair per day). This emission estimate compared favourably to the value obtained using a simple micrometeorological technique (aerodynamic gradient) of 0.98 g ammonia per breeding pair per day (95% confidence interval: 0.2-2.4 g ammonia per breeding pair per day). Further application of the inverse modelling techniques for a range of agricultural case studies also demonstrated good agreement between the emission estimates. It is concluded, therefore, that inverse dispersion modelling is a robust technique for estimating NH3 emission rates. Screening models that can provide a quick and approximate estimate of environmental impacts are a useful tool for impact assessments because they can be used to filter out cases that potentially have a minimal environmental impact allowing resources to be focussed on more potentially damaging cases. The Simple Calculation of Ammonia Impact Limits (SCAIL) model was developed as a screening model to provide an estimate of the mean NH3 concentration and dry deposition rate downwind of an agricultural source. This screening tool, based on the LADD model, was evaluated and calibrated with several experimental datasets and then validated using independent concentration measurements made near sources. Overall SCAIL performed acceptably according to established statistical criteria. This work has identified situations where the concentration predictions of dispersion models are similar and other situations where the predictions are significantly different. Some models are simply not designed to simulate certain scenarios since they do not include the relevant processes or are beyond the limits of their applicability. An example is the LADD model that is not applicable to sources with significant exit velocity since the model does not include a plume-rise parameterisation. The testing of a simple scheme combining a momentum-driven plume rise and increased turbulence at the source improved model performance, but more testing is required. Even models that are applicable and include the relevant process do not always give similar predictions and the reasons for this need to be investigated. AERMOD for example predicts higher concentrations than ADMS for dispersion from mechanically ventilated livestock housing. There is evidence to suggest that ADMS underestimates concentrations in these situations due to a high wind speed threshold. Conversely, there is also evidence that AERMOD overestimates concentrations in these situations due to overestimation at low wind speeds. However, a simple modification to the meteorological pre-processor appears to improve the performance of the model. It is important that these differences between the predictions of these models are taken into account in regulatory assessments. This can be done by applying the most suitable model for the assessment in question or, better still, using multiple or hybrid models.

Simulación dinámica y térmica de una habitación de acuerdo con el Código Técnico de la Edificación

El Código Técnico de la edificación (CTE), es el marco normativo por el que se regulan las exigencias básicas de calidad que deben de cumplir los edificios. Dichas exigencias establecen los requisitos imprescindibles en el ámbito de seguridad estructural, seguridad de utilización ahorro de energía, aislamiento térmico, higiene, salud, protección del medio ambiente…. y se rigen, cada una de ellas, por su norma específica (Norma UNE). La norma ISO 7730:2006 es una norma internacional que aborda la evaluación del medio ambiente, en el interior de un edificio, tal y como se recoge en el CTE. Con este propósito la norma establece la predicción de índices de calidad del aire IAQ (Indoor Air cuality) así como de la predicción del campo de velocidades, de temperatura y, en general, de todas las magnitudes físicas, características, tanto del aire como de los elementos constructivos utilizados en la edificación. Basándonos en estos criterios el objetivo de este trabajo es presentar los primeros resultados obtenidos en la simulación de parámetros de confort en una habitación, de una vivienda construida con criterios bioclimáticos, mediante el programa de simulación STAR_CCM+ de Computational Fluid Dynamics (CFD). Utilizando STAR_CCM+ diseñamos la habitación en estudio y simulamos el campo de velocidades y de temperatura, en el interior de la misma, para distintas configuraciones de aberturas utilizando los criterios de velocidad y de temperatura del aire dado por la ISO 7730:2006. En esta línea se está investigando actualmente y el trabajo que se presenta muestra que el programa de simulación utilizado permite, por un lado, diseñar la disposición más adecuada de las aberturas en una habitación en particular, y de una vivienda en general, y, por otro lado, hacer un estudio cualitativo y cuantitativo de los parámetros de confort en diferentes situaciones meteorológicas, de acuerdo con el CTE, incluso antes de que la vivienda esté construida

Una solución exacta en magnetogasdinámica transónica

Se plantean las ecuaciones de la Magnetogasdinámica desde un punto de vista puramente macroscópico,bajo la hipótesis fundamental de que las ecuaciones de MAXWELL no son modificadas por el movimiento gasdinámico y de que el gas es un medio continuo que se comporta como un fluido neutro y no como una mezcla de partículas individuales positivas y negativas. A continuación se consideran únicamente movimientos estacionarios y no disipativos. Entre este tipo particular de corrientes se eligen aquellas en las que el campo magnético es, en cada punto, paralelo al campo de velocidades; debiendo investigar a posteriori las condiciones de contorno que garanticen la existencia de tal tipo de soluciones.Después de comprobar que, en nuestro caso particular, el problema magnetogasdinámico se puede plantear en una forma muy similar a la de ciertos problemas puramente gasdinámicos.Se estudia el problema transónico de RlNGLEB para un gas conductor.La conclusión fundamental deducida de la solución obtenida, es la de que la presencia del campo magnético permite, bajo ciertas condiciones, obtener corrientes continuas isentrópicas a números de Mach mucho mayores que los que se obtienen en en el problema gasdinámico clásico.

Efectos de desequilibrio químico en aerodinámica hipersónica

Se considera la influencia del desequilibrio químico de las corrientes hipersónicas, en los dos casos siguientes: a) Recombinación y disociación en la capa límite. — Se pone de manifiesto que, puesto que en el dominio de vuelo de máximo calentamiento convectivo la corriente anterior no viscosa está en equilibrio, el alejamiento del equilibrio tiene lugar únicamente cerca de la pared. Esto facilita la obtención de una solución analítica sencilla en la que se muestran claramente los efectos de las velocidades de reacción, forma del obstáculo y regímenes de vuelo. b) Reacciones químicas entre los constituyentes del aire y una substancia inyectada o vaporizada desde la pared.—En este caso, se supone que la reacción química es tan rápida que tiene lugar en una zona muy estrecha. Esto simplifica la obtención de una solución analítica con velocidad de reacción finita, la cual también proporciona un criterio para la extinción de la llama.

Comportamiento aerodinámico de barreras cortavientos 1

El trabajo presenta resultados de la exploración del campo de velocidades en estelas de distintos modelos de barreras porosas, con vistas a suposible utilización para proteger plantaciones o zonas habitadas de los efectos perniciosos del viento. Los ensayos se han hecho en túnel aerodinámico y, para las medidas, se ha utilizado un anemómetro de hilo caliente. Se ha medido: la velocidad media temporal, el nivel de turbulencia y, en algunos casos, el transporte turbulento y la inclinación del vector velocidad, en cuatro secciones distintas corriente abajo de cada barrera. El trabajo está dividido en ocho capítulos. En el primero se discute la posible utilidad de las barreras en agricultura. En el segundo se revisan ciertas características de las estelas turbulentas de diferentes obstáculos,con el fin de contrastar los resultados obtenidos con otros análogos. En el capítulo tercero se describen los métodos de medida. El cuarto presenta el grueso de los resultados, y los cuatro siguientes: la influencia del número de Reynolds, rugosidad del suelo, perfil del viento incidente y presencia de las paredes y techo del túnel, respectivamente. Los ensayos muestran que los perfiles de velocidades correspondientes a secciones situadas a cierta distancia de la barrera (mayor de seis a ocho veces la altura de ésta) se asemejan a los de un semichorro de baja velocidad que descarga paralelamente a una corriente más rápida. Los perfiles próximos y las condiciones iniciales del semichorro dependen de la configuración de la barrera. El modelo del semichorro equivalente permite calcular las características de la estela a distancias mayores que las que es posible reproducir en los experimentos. Por otra parte, sugiere ciertas modificaciones de la forma de la barrera para aumentar la longitud de la zona protegida del viento. La rugosidad del suelo y la existencia de perfiles de viento distintos del uniforme y más ajustados a la realidad, contribuyen a disminuir la longitud de la zona protegida. Esta observación está de acuerdo con los resultados de otros autores. Se observa que el nivel de turbulencia es muy sensible a las características geométricas de la barrera, lo que sugiere la posibilidad de controlar la capa límite sobre el terreno y, por tanto, el transporte de calor y masa en provecho de la productividad de cultivo. El trabajo que se presenta es parte de un programa más amplio, que tiene por objeto transmitir tecnología avanzada a ciertos dominios de interés para la agricultura y la industria.

Red GPS Fiduciaria Nacional de España

El Instituto Geográfico Nacional, por medio del Área de Geodesia, está llevando a cabo el establecimiento de una Red de Estaciones Permanentes GPS que permitan obtener coordenadas muy precisas, así como sus campos de velocidades en un Sistema de Referencia Global (ITRFxx). Dichas estaciones pertenencen a la Red de Estaciones Permanentes de EUREF (EUropean REference Framen) y constituyen el órden cero de la Geodesia Española.

Modelado de materiales piezoeléctricos para aplicaciones de autoalimentación en entornos viales

Este proyecto continua con el estudio iniciado en los proyectos: Alimentación “Energy Harvesting” basada en fuentes piezoeléctricas, Alimentación autónoma: aplicación a fuente piezoeléctrica y Banco de pruebas para sistemas autoalimentados en aplicaciones viales, que se realizaron anteriormente. Este estudio se basa en la búsqueda de soluciones de alimentación Energy Harvesting centrados en fuentes piezoeléctricas. El objetivo de este estudio es conocer las posibilidades que pueden ofrecer los materiales piezoeléctricos para alimentación autónoma cuando son excitados por vehículos en diferentes circunstancias y entornos viales. Para la realización de este estudio se han utilizado señales provenientes del banco de ensayos, así como señales obtenidas de una instalación construida en un vial del Campus. Para poder realizar el estudio de los diversos materiales piezoeléctricos utilizados en este proyecto, se necesitan conocer las características de los modelos eléctricos. Para obtener estas características se ha diseñado un programa basado en LabView, este programa además puede modificar la señal para que simule las diversas circunstancias a las que pueden estar sometidos los diferentes materiales piezoeléctricos. En este estudio se ha tenido en cuenta la distancia de seguridad entre vehículos para diferentes velocidades y la posibilidad de poder amplificar la excitación que pueden recibir los diversos materiales. Utilizando el programa PSpice OrCAD® se realizan los diferentes modelos con el objetivo de observar el comportamiento de los materiales bajo diversas circunstancias. También se han utilizado varias topologías rectificadoras diferentes para poder comprobar en que situación se puede obtener la máxima tensión acumulada posible. Además también se ha evaluado la máxima potencia que puede entregar cada material en cada situación. Por ultimo se han realizado varias pruebas prácticas con el banco de ensayos para comprobar la acumulación energética real que se puede obtener de los materiales instalados en él. También se ha realizado el estudio de una aplicación práctica que sea capaz de alimentar pequeños circuitos o sensores de poco consumo. Abstract This Project goes on with the study started in the projects Feeding “Energy Harvesting” based on piezoelectric sources, independent feeding: application to piezoelectric source and testing bench to self fed systems in roads applications that were previously done. This study is based on the search of feeding solutions Energy Harvesting focus on piezoelectric sources. The aim of this study is know the different piezoelectric materials possibilities that can offer to independent feeding when different circumstances and road surroundings are excited by vehicles in. To develop this study several signals of the testing bench have been used and also signals get from an installation constructed on a Campus road. As far as the study of the different piezoelectric materials used in this project is concerned we should know the electrical models features. In order to obtain these characteristics a program based on LabView was designed, this program can modify the signal so that it simulates the diverse circumstances to which they can be put under the different piezoelectric materials. In this study I have bearded in mind the safe distance between vehicles for different speeds and the possibility of being able to amplify the excitation that can receive the diverse materials. Using the program PSpice OrCAD® the different models with the objective are made to observe the behaviour of the materials under diverse circumstances. Also several different rectifying topologies have been used to be able to verify the maximum overvoltage accumulated. In addition I have assessed the maximum power that can give each material in each situation. Last but not least different practical proofs with the testing bench have been done to check the actual energy accumulation that can be obtained from the materials installed on it. Also the study of a practical application has been made that is able to feed small circuits or low consumption sensor.