Simulación numérica de un problema de contaminación atmosférica

La ecuación en derivadas parciales de advección difusión con reacción química es la base de los modelos de dispersión de contaminantes en la atmósfera, y los diferentes métodos numéricos empleados para su resolución han sido objeto de amplios estudios a lo largo de su desarrollo. En esta Tesis se presenta la implementación de un nuevo método conservativo para la resolución de la parte advectiva de la ecuación en derivadas parciales que modela la dispersión de contaminantes dentro del modelo mesoescalar de transporte químico CHIMERE. Este método está basado en una técnica de volúmenes finitos junto con una interpolación racional. La ventaja de este método es la conservación exacta de la masa transportada debido al empleo de la ley de conservación de masas. Para ello emplea una formulación de flujo basado en el cálculo de la integral ponderada dentro de cada celda definida para la discretización del espacio en el método de volúmenes finitos. Los resultados numéricos obtenidos en las simulaciones realizadas (implementando el modelo conservativo para la advección en el modelo CHIMERE) se han comparado con los datos observados de concentración de contaminantes registrados en la red de estaciones de seguimiento y medición distribuidas por la Península Ibérica. Los datos estadísticos de medición del error, la media normalizada y la media absoluta normalizada del error, presentan valores que están dentro de los rangos propuestos por la EPA para considerar el modelo preciso. Además, se introduce un nuevo método para resolver la parte advectivadifusiva de la ecuación en derivadas parciales que modeliza la dispersión de contaminantes en la atmósfera. Se ha empleado un método de diferencias finitas de alto orden para resolver la parte difusiva de la ecuación de transporte de contaminantes junto con el método racional conservativo para la parte advectiva en una y dos dimensiones. Los resultados obtenidos de la aplicación del método a diferentes situaciones incluyendo casos académicos y reales han sido comparados con la solución analítica de la ecuación de advección-difusión, demostrando que el nuevo método proporciona un resultado preciso para aproximar la solución. Por último, se ha desarrollado un modelo completo que contempla los fenómenos advectivo y difusivo con reacción química, usando los métodos anteriores junto con una técnica de diferenciación regresiva (BDF2). Esta técnica consiste en un método implícito multipaso de diferenciación regresiva de segundo orden, que nos permite resolver los problemas rígidos típicos de la química atmosférica, modelizados a través de sistemas de ecuaciones diferenciales ordinarias. Este método hace uso de la técnica iterativa Gauss- Seidel para obtener la solución de la parte implícita de la fórmula BDF2. El empleo de la técnica de Gauss-Seidel en lugar de otras técnicas comúnmente empleadas, como la iteración por el método de Newton, nos proporciona rapidez de cálculo y bajo consumo de memoria, ideal para obtener modelos operativos para la resolución de la cinética química atmosférica. ABSTRACT Extensive research has been performed to solve the atmospheric chemicaladvection- diffusion equation and different numerical methods have been proposed. This Thesis presents the implementation of an exactly conservative method for the advection equation in the European scale Eulerian chemistry transport model CHIMERE based on a rational interpolation and a finite volume algorithm. The advantage of the method is that the cell-integrated average is predicted via a flux formulation, thus the mass is exactly conserved. Numerical results are compared with a set of observation registered at some monitoring sites in Spain. The mean normalized bias and the mean normalized absolute error present values that are inside the range to consider an accurate model performance. In addition, it has been introduced a new method to solve the advectiondiffusion equation. It is based on a high-order accurate finite difference method to solve de diffusion equation together with a rational interpolation and a finite volume to solve the advection equation in one dimension and two dimensions. Numerical results obtained from solving several problems include academic and real atmospheric problems have been compared with the analytical solution of the advection-diffusion equation, showing that the new method give an efficient algorithm for solving such problems. Finally, a complete model has been developed to solve the atmospheric chemical-advection-diffusion equation, adding the conservative method for the advection equation, the high-order finite difference method for the diffusion equation and a second-order backward differentiation formula (BDF2) to solve the atmospheric chemical kinetics. The BDF2 is an implicit, second order multistep backward differentiation formula used to solve the stiff systems of ordinary differential equations (ODEs) from atmospheric chemistry. The Gauss-Seidel iteration is used for approximately solving the implicitly defined BDF solution, giving a faster tool than the more commonly used iterative modified Newton technique. This method implies low start-up costs and a low memory demand due to the use of Gauss-Seidel iteration.

Speed-independent gait identification for mobile devices

Due to the intensive use of mobile phones for diferent purposes, these devices usually contain condential information which must not be accessed by another person apart from the owner of the device. Furthermore, the new generation phones commonly incorporate an accelerometer which may be used to capture the acceleration signals produced as a result of owner s gait. Nowadays, gait identication in basis of acceleration signals is being considered as a new biometric technique which allows blocking the device when another person is carrying it. Although distance based approaches as Euclidean distance or dynamic time warping have been applied to solve this identication problem, they show di±culties when dealing with gaits at diferent speeds. For this reason, in this paper, a method to extract an average template from instances of the gait at diferent velocities is presented. This method has been tested with the gait signals of 34 subjects while walking at diferent motion speeds (slow, normal and fast) and it has shown to improve the performance of Euclidean distance and classical dynamic time warping.

On the Use of Matrix-Free Shift-Invert Strategies For Global Flow Instability Analysis

A novel time-stepping shift-invert algorithm for linear stability analysis of laminar flows in complex geometries is presented. This method, based on a Krylov subspace iteration, enables the solution of complex non-symmetric eigenvalue problems in a matrix-free framework. Validations and comparisons to the classical exponential method have been performed in three different cases: (i) stenotic flow, (ii) backward-facing step and (iii) lid-driven swirling flow. Results show that this new approach speeds up the required Krylov subspace iterations and has the capability of converging to specific parts of the global spectrum. It is shown that, although the exponential method remains the method of choice if leading eigenvalues are sought, the performance of the present method could be dramatically improved with the use of a preconditioner. In addition, as opposed to other methods, this strategy can be directly applied to any time-stepper, regardless of the temporal or spatial discretization of the latter.

Comparison of micrometeorological techinques for estimating ammonia emission from covered slurry storage and land spreading of cattle slurry

The micrometeorological mass-balance integrated horizontal flux (IHF) technique has been commonly employed for measuring ammonia (NH3) emissions inon-field experiments. However, the inverse-dispersion modeling technique, such as the backward Lagrangian stochastic (bLS) modeling approach, is currently highlighted as offering flexibility in plot design and requiring a minimum number of samplers (Ro et al., 2013). The objective of this study was to make a comparison between the bLS technique with the IHF technique for estimating NH3 emission from flexible bag storage and following landspreading of dairy cattle slurry. Moreover, considering that NH3 emission in storage could have been non uniform, the effect on bLS estimates of a single point and multiple downwind concentration measurements was tested, as proposed by Sanz et al. (2010).

Methodology of quantifying curvature of Fresnel lenses and its effect on CPV module performance

Fresnel lenses used as primary optics in concentrating photovoltaic modules may show warping produced by lens manufacturing or module assembly (e.g., stress during molding or weight load) or due to stress during operation (e.g., mismatch of thermal expansion between different materials). To quantify this problem, a simple method called “checkerboard method” is presented. The proposed method identifies shape errors on the front surface of primary lenses by analyzing the Fresnel reflections. This paper also deals with the quantification of the effects these curvatures have on their optical performance and on the electrical performance of concentrating modules incorporating them. This method can be used to perform quality control of Fresnel lenses in scenarios of high volume production.

Parallel Computer Vision Algorithms for Graphics Processing Units

La evolución de los teléfonos móviles inteligentes, dotados de cámaras digitales, está provocando una creciente demanda de aplicaciones cada vez más complejas que necesitan algoritmos de visión artificial en tiempo real; puesto que el tamaño de las señales de vídeo no hace sino aumentar y en cambio el rendimiento de los procesadores de un solo núcleo se ha estancado, los nuevos algoritmos que se diseñen para visión artificial han de ser paralelos para poder ejecutarse en múltiples procesadores y ser computacionalmente escalables. Una de las clases de procesadores más interesantes en la actualidad se encuentra en las tarjetas gráficas (GPU), que son dispositivos que ofrecen un alto grado de paralelismo, un excelente rendimiento numérico y una creciente versatilidad, lo que los hace interesantes para llevar a cabo computación científica. En esta tesis se exploran dos aplicaciones de visión artificial que revisten una gran complejidad computacional y no pueden ser ejecutadas en tiempo real empleando procesadores tradicionales. En cambio, como se demuestra en esta tesis, la paralelización de las distintas subtareas y su implementación sobre una GPU arrojan los resultados deseados de ejecución con tasas de refresco interactivas. Asimismo, se propone una técnica para la evaluación rápida de funciones de complejidad arbitraria especialmente indicada para su uso en una GPU. En primer lugar se estudia la aplicación de técnicas de síntesis de imágenes virtuales a partir de únicamente dos cámaras lejanas y no paralelas—en contraste con la configuración habitual en TV 3D de cámaras cercanas y paralelas—con información de color y profundidad. Empleando filtros de mediana modificados para la elaboración de un mapa de profundidad virtual y proyecciones inversas, se comprueba que estas técnicas son adecuadas para una libre elección del punto de vista. Además, se demuestra que la codificación de la información de profundidad con respecto a un sistema de referencia global es sumamente perjudicial y debería ser evitada. Por otro lado se propone un sistema de detección de objetos móviles basado en técnicas de estimación de densidad con funciones locales. Este tipo de técnicas es muy adecuada para el modelado de escenas complejas con fondos multimodales, pero ha recibido poco uso debido a su gran complejidad computacional. El sistema propuesto, implementado en tiempo real sobre una GPU, incluye propuestas para la estimación dinámica de los anchos de banda de las funciones locales, actualización selectiva del modelo de fondo, actualización de la posición de las muestras de referencia del modelo de primer plano empleando un filtro de partículas multirregión y selección automática de regiones de interés para reducir el coste computacional. Los resultados, evaluados sobre diversas bases de datos y comparados con otros algoritmos del estado del arte, demuestran la gran versatilidad y calidad de la propuesta. Finalmente se propone un método para la aproximación de funciones arbitrarias empleando funciones continuas lineales a tramos, especialmente indicada para su implementación en una GPU mediante el uso de las unidades de filtraje de texturas, normalmente no utilizadas para cómputo numérico. La propuesta incluye un riguroso análisis matemático del error cometido en la aproximación en función del número de muestras empleadas, así como un método para la obtención de una partición cuasióptima del dominio de la función para minimizar el error. ABSTRACT The evolution of smartphones, all equipped with digital cameras, is driving a growing demand for ever more complex applications that need to rely on real-time computer vision algorithms. However, video signals are only increasing in size, whereas the performance of single-core processors has somewhat stagnated in the past few years. Consequently, new computer vision algorithms will need to be parallel to run on multiple processors and be computationally scalable. One of the most promising classes of processors nowadays can be found in graphics processing units (GPU). These are devices offering a high parallelism degree, excellent numerical performance and increasing versatility, which makes them interesting to run scientific computations. In this thesis, we explore two computer vision applications with a high computational complexity that precludes them from running in real time on traditional uniprocessors. However, we show that by parallelizing subtasks and implementing them on a GPU, both applications attain their goals of running at interactive frame rates. In addition, we propose a technique for fast evaluation of arbitrarily complex functions, specially designed for GPU implementation. First, we explore the application of depth-image–based rendering techniques to the unusual configuration of two convergent, wide baseline cameras, in contrast to the usual configuration used in 3D TV, which are narrow baseline, parallel cameras. By using a backward mapping approach with a depth inpainting scheme based on median filters, we show that these techniques are adequate for free viewpoint video applications. In addition, we show that referring depth information to a global reference system is ill-advised and should be avoided. Then, we propose a background subtraction system based on kernel density estimation techniques. These techniques are very adequate for modelling complex scenes featuring multimodal backgrounds, but have not been so popular due to their huge computational and memory complexity. The proposed system, implemented in real time on a GPU, features novel proposals for dynamic kernel bandwidth estimation for the background model, selective update of the background model, update of the position of reference samples of the foreground model using a multi-region particle filter, and automatic selection of regions of interest to reduce computational cost. The results, evaluated on several databases and compared to other state-of-the-art algorithms, demonstrate the high quality and versatility of our proposal. Finally, we propose a general method for the approximation of arbitrarily complex functions using continuous piecewise linear functions, specially formulated for GPU implementation by leveraging their texture filtering units, normally unused for numerical computation. Our proposal features a rigorous mathematical analysis of the approximation error in function of the number of samples, as well as a method to obtain a suboptimal partition of the domain of the function to minimize approximation error.

Estudo cefalométrico-radiográfico das alterações dento-esqueléticas em pacientes com má oclusão de Classe III submetidos a tratamento ortodôntico-cirúrgico

Este estudo teve como finalidade avaliar cefalometricamente, por meio de telerradiografias em norma lateral, as alterações dento-esqueléticas em pacientes Classe III submetidos a tratamento ortodôntico-cirúrgico. A amostra experimental constituiu-se de 16 pacientes Brasileiros, dos sexos masculino e feminino, na faixa etária pré-cirúrgica média de 21 anos e 11 meses, apresentando má oclusão de Classe III com indicação de tratamento cirúrgico representado por recuo mandibular isolado. Para cada paciente foram realizadas telerradiografias nas fases inicial, pré-cirúrgica e pós-cirúrgica, sendo comparadas a um grupo controle, constituído de telerradiografias de indivíduos com oclusão normal. Segundo a metodologia empregada e pela análise dos resultados obtidos, avaliados estatisticamente, constatou-se que os pacientes Classe III com indicação de recuo mandibular foram caracterizados por um mau relacionamento entre as bases esqueléticas representado por um bom posicionamento da maxila associado a prognatismo mandibular, aumento da altura facial ântero-inferior, incisivos inferiores e sínfise mandibular lingualizados e incisivos superiores vestibularizados. A partir do preparo ortodôntico pré-cirúrgico, observou-se uma rotação mandibular no sentido horário, descompensação dentária representada por lingualização e extrusão dos incisivos superiores e vestibularização dos incisivos inferiores, acompanhada por uma remodelação da cortical óssea vestibular da sínfise mandibular. Esta descompensação ortodôntica definiu características dento-alveolares semelhantes às dos indivíduos com oclusão normal. O comportamento das variáveis dento-esqueléticas após a cirurgia ortognática, a partir do deslocamento póstero-superior das estruturas dento-esqueléticas da mandíbula, proporcionou um equilíbrio destas estruturas, em relação à oclusão normal.

Estudo cefalométrico-radiográfico das alterações dento-esqueléticas em pacientes com má oclusão de Classe III submetidos a tratamento ortodôntico-cirúrgico

Este estudo teve como finalidade avaliar cefalometricamente, por meio de telerradiografias em norma lateral, as alterações dento-esqueléticas em pacientes Classe III submetidos a tratamento ortodôntico-cirúrgico. A amostra experimental constituiu-se de 16 pacientes Brasileiros, dos sexos masculino e feminino, na faixa etária pré-cirúrgica média de 21 anos e 11 meses, apresentando má oclusão de Classe III com indicação de tratamento cirúrgico representado por recuo mandibular isolado. Para cada paciente foram realizadas telerradiografias nas fases inicial, pré-cirúrgica e pós-cirúrgica, sendo comparadas a um grupo controle, constituído de telerradiografias de indivíduos com oclusão normal. Segundo a metodologia empregada e pela análise dos resultados obtidos, avaliados estatisticamente, constatou-se que os pacientes Classe III com indicação de recuo mandibular foram caracterizados por um mau relacionamento entre as bases esqueléticas representado por um bom posicionamento da maxila associado a prognatismo mandibular, aumento da altura facial ântero-inferior, incisivos inferiores e sínfise mandibular lingualizados e incisivos superiores vestibularizados. A partir do preparo ortodôntico pré-cirúrgico, observou-se uma rotação mandibular no sentido horário, descompensação dentária representada por lingualização e extrusão dos incisivos superiores e vestibularização dos incisivos inferiores, acompanhada por uma remodelação da cortical óssea vestibular da sínfise mandibular. Esta descompensação ortodôntica definiu características dento-alveolares semelhantes às dos indivíduos com oclusão normal. O comportamento das variáveis dento-esqueléticas após a cirurgia ortognática, a partir do deslocamento póstero-superior das estruturas dento-esqueléticas da mandíbula, proporcionou um equilíbrio destas estruturas, em relação à oclusão normal.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

Game theory and reciprocity in some extensive form experimental games

We examine decision making in two-person extensive form game trees using nine treatments that vary matching protocol, payoffs, and payoff information. Our objective is to establish replicable principles of cooperative versus noncooperative behavior that involve the use of signaling, reciprocity, and backward induction strategies, depending on the availability of dominated direct punishing strategies and the probability of repeated interaction with the same partner. Contrary to the predictions of game theory, we find substantial support for cooperation under complete information even in various single-play treatments.

Borrelia burgdorferi periplasmic flagella have both skeletal and motility functions

Bacterial shape usually is dictated by the peptidoglycan layer of the cell wall. In this paper, we show that the morphology of the Lyme disease spirochete Borrelia burgdorferi is the result of a complex interaction between the cell cylinder and the internal periplasmic flagella. B. burgdorferi has a bundle of 7–11 helically shaped periplasmic flagella attached at each end of the cell cylinder and has a flat-wave cell morphology. Backward moving, propagating waves enable these bacteria to swim in both low viscosity media and highly viscous gel-like media. Using targeted mutagenesis, we inactivated the gene encoding the major periplasmic flagellar filament protein FlaB. The resulting flaB mutants not only were nonmotile, but were rod-shaped. Western blot analysis indicated that FlaB was no longer synthesized, and electron microscopy revealed that the mutants were completely deficient in periplasmic flagella. Wild-type cells poisoned with the protonophore carbonyl cyanide-m-chlorophenylhydrazone retained their flat-wave morphology, indicating that the periplasmic flagella do not need to be energized for the cell to maintain this shape. Our results indicate that the periplasmic flagella of B. burgdorferi have a skeletal function. These organelles dynamically interact with the rod-shaped cell cylinder to enable the cell to swim, and to confer in part its flat-wave morphology.

A FYVE-finger-containing protein, Rabip4, is a Rab4 effector involved in early endosomal traffic

The small GTPase Rab4 is implicated in endocytosis in all cell types, but also plays a specific role in some regulated processes. To better understand the role of Rab4 in regulation of vesicular trafficking, we searched for an effector(s) that specifically recognizes its GTP-bound form. We cloned a ubiquitous 69-kDa protein, Rabip4, that behaves as a Rab4 effector in the yeast two-hybrid system and in the mammalian cell. Rabip4 contains two coiled-coil domains and a FYVE-finger domain. When expressed in CHO cells, Rabip4 is present in early endosomes, because it is colocated with endogenous Early Endosome Antigen 1, although it is absent from Rab11-positive recycling endosomes and Rab-7 positive late endosomes. The coexpression of Rabip4 with active Rab4, but not with inactive Rab4, leads to an enlargement of early endosomes. It strongly increases the degree of colocalization of markers of sorting (Rab5) and recycling (Rab11) endosomes with Rab4. Furthermore, the expression of Rabip4 leads to the intracellular retention of a recycling molecule, the glucose transporter Glut 1. We propose that Rabip4, an effector of Rab4, controls early endosomal traffic possibly by activating a backward transport step from recycling to sorting endosomes.

Functional and structural mapping of human cerebral cortex: Solutions are in the surfaces

The human cerebral cortex is notorious for the depth and irregularity of its convolutions and for its variability from one individual to the next. These complexities of cortical geography have been a chronic impediment to studies of functional specialization in the cortex. In this report, we discuss ways to compensate for the convolutions by using a combination of strategies whose common denominator involves explicit reconstructions of the cortical surface. Surface-based visualization involves reconstructing cortical surfaces and displaying them, along with associated experimental data, in various complementary formats (including three-dimensional native configurations, two-dimensional slices, extensively smoothed surfaces, ellipsoidal representations, and cortical flat maps). Generating these representations for the cortex of the Visible Man leads to a surface-based atlas that has important advantages over conventional stereotaxic atlases as a substrate for displaying and analyzing large amounts of experimental data. We illustrate this by showing the relationship between functionally specialized regions and topographically organized areas in human visual cortex. Surface-based warping allows data to be mapped from individual hemispheres to a surface-based atlas while respecting surface topology, improving registration of identifiable landmarks, and minimizing unwanted distortions. Surface-based warping also can aid in comparisons between species, which we illustrate by warping a macaque flat map to match the shape of a human flat map. Collectively, these approaches will allow more refined analyses of commonalities as well as individual differences in the functional organization of primate cerebral cortex.

A new view of language acquisition

At the forefront of debates on language are new data demonstrating infants' early acquisition of information about their native language. The data show that infants perceptually “map” critical aspects of ambient language in the first year of life before they can speak. Statistical properties of speech are picked up through exposure to ambient language. Moreover, linguistic experience alters infants' perception of speech, warping perception in the service of language. Infants' strategies are unexpected and unpredicted by historical views. A new theoretical position has emerged, and six postulates of this position are described.

Modeling back propagating action potential in weakly excitable dendrites of neocortical pyramidal cells.

Simultaneous recordings from the soma and apical dendrite of layer V neocortical pyramidal cells of young rats show that, for any location of current input, an evoked action potential (AP) always starts at the axon and then propagates actively, but decrementally, backward into the dendrites. This back-propagating AP is supported by a low density (-gNa = approximately 4 mS/cm2) of rapidly inactivating voltage-dependent Na+ channels in the soma and the apical dendrite. Investigation of detailed, biophysically constrained, models of reconstructed pyramidal cells shows the following. (i) The initiation of the AP first in the axon cannot be explained solely by morphological considerations; the axon must be more excitable than the soma and dendrites. (ii) The minimal Na+ channel density in the axon that fully accounts for the experimental results is about 20-times that of the soma. If -gNa in the axon hillock and initial segment is the same as in the soma [as recently suggested by Colbert and Johnston [Colbert, C. M. & Johnston, D. (1995) Soc. Neurosci. Abstr. 21, 684.2]], then -gNa in the more distal axonal regions is required to be about 40-times that of the soma. (iii) A backward propagating AP in weakly excitable dendrites can be modulated in a graded manner by background synaptic activity. The functional role of weakly excitable dendrites and a more excitable axon for forward synaptic integration and for backward, global, communication between the axon and the dendrites is discussed.