Wind tunnel analysis of the flow topology over complex topographies

El principal objetivo de este trabajo es aportar conocimiento para contestar la pregunta: ¿hasta que punto los ensayos en túnel aerodinámico pueden contribuir a determinar las características que afectan la respuesta dinámica de los aerogeneradores operando en terreno complejo?. Esta pregunta no es nueva, de hecho, el debate en la comunidad científica comenzó en el primer tercio del siglo pasado y aún está intensamente vivo. El método generalmente aceptado para enfrentar el mencionado problema consiste en analizar un caso de estudio determinado en el cual se aplican tanto ensayos a escala real como análisis computacionales y ensayos en túnel aerodinámico. Esto no es ni fácil ni barato. Esta es la razón por la cual desde el experimento de Askervein en 1988, los modelizadores del flujo atmosférico tuvieron que esperar hasta 2007 a que el experimento de Bolund fuese puesto en marcha con un despliegue de medios técnicos equivalentes (teniendo en cuenta la evolución de las tecnologías de sensores y computación). El problema contempla tantos aspectos que ambas experiencias fueron restringidas a condiciones de atmósfera neutra con efectos de Coriolis despreciables con objeto de reducir la complejidad. Este es el contexto en el que se ha desarrollado la presente tesis doctoral. La topología del flujo sobre la isla de Bolund ha sido estudiada mediante la reproducción del experimento de Bolund en los túneles aerodinámicos A9 y ACLA16 del IDR. Dos modelos de la isla de Bolund fueron fabricados a dos escalas, 1:230 y 1:115. El flujo de entrada en el túnel aerodinámico simulando la capa límite sin perturbar correspondía a régimen de transición (transitionally rough regime) y fue usado como situación de referencia. El modelo a escala 1:230 fue ensayado en el túnel A9 para determinar la presión sobre su superficie. La distribución del coeficiente de presión sobre la isla proporcionó una visualización y estimación de una región de desprendimiento sobre el pequeño acantilado situado al frente de la misma. Las medidas de presión instantánea con suficiente grado de resolución temporal pusieron de manifiesto la no estacionariedad en la región de desprendimiento. El modelo a escala 1:115 fue ensayado utilizando hilo caliente de tres componentes y un sistema de velocimetría por imágenes de partículas de dos componentes. El flujo fue caracterizado por el ratio de aceleración, el incremento normalizado de energía cinética turbulenta y los ángulos de inclinación y desviación horizontal. Los resultados a lo largo de la dirección 270°y alturas de 2 m y 5 m presentaron una gran similitud con los resultados a escala real del experimento de Bolund. Los perfiles verticales en las localizaciones de las torres meteorológicas mostraron un acuerdo significativo con los resultados a escala real. El análisis de los esfuerzos de Reynolds y el análisis espectral en las localizaciones de los mástiles meteorológicos presentaron niveles de acuerdo variados en ciertas posiciones, mientras que en otras presentaron claras diferencias. El mapeo horizontal del flujo, para una dirección de viento de 270°, permitió caracterizar el comportamiento de la burbuja intermitente de recirculación sobre el pequeño acantilado existente al frente de la isla así como de la región de relajación y de la capa de cortadura en la región corriente abajo de Bolund. Se realizaron medidas de velocidad con alta resolución espacial en planos perpendiculares a la dirección del flujo sin perturbar. Estas medidas permitieron detectar y caracterizar una estructura de flujo similar a un torbellino longitudinal con regiones con altos gradientes de velocidad y alta intensidad de turbulencia. Esta estructura de flujo es, sin duda, un reto para los modelos computacionales y puede considerarse un factor de riesgo para la operación de los aerogeneradores. Se obtuvieron y analizaron distribuciones espaciales de los esfuerzos de Reynolds mediante 3CHW y PIV. Este tipo de parámetros no constituyen parte de los resultados habituales en los ensayos en túnel sobre topografías y son muy útiles para los modelizadores que utilizan simulación de grades torbellinos (LES). Se proporciona una interpretación de los resultados obtenidos en el túnel aerodinámico en términos de utilidad para los diseñadores de parques eólicos. La evolución y variación de los parámetros del flujo a lo largo de líneas, planos y superficies han permitido identificar como estas propiedades del flujo podrían afectar la localización de los aerogeneradores y a la clasificación de emplazamientos. Los resultados presentados sugieren, bajo ciertas condiciones, la robustez de los ensayos en túnel para estudiar la topología sobre terreno complejo y su comparabilidad con otras técnicas de simulación, especialmente considerando el nivel de acuerdo del conjunto de resultados presentados con los resultados a escala real. De forma adicional, algunos de los parámetros del flujo obtenidos de las medidas en túnel son difícilmente determinables en ensayos a escala real o por medios computacionales, considerado el estado del arte. Este trabajo fue realizado como parte de las actividades subvencionadas por la Comisión Europea como dentro del proyecto FP7-PEOPLE-ITN-2008WAUDIT (Wind Resource Assessment Audit and Standardization) dentro de la FP7 Marie-Curie Initial Training Network y por el Ministerio Español de Economía y Competitividad dentro del proyecto ENE2012-36473, TURCO (Determinación en túnel aerodinámico de la distribución espacial de parámetros estadísticos de la turbulencia atmosférica sobre topografías complejas) del Plan Nacional de Investigación (Subprograma de investigación fundamental no orientada 2012). El informe se ha organizado en siete capítulos y un conjunto de anexos. En el primer capítulo se introduce el problema. En el capítulo dos se describen los medios experimentales utilizados. Seguidamente, en el capítulo tres, se analizan en detalle las condiciones de referencia del principal túnel aerodinámico utilizado en esta investigación. En el capítulo tres se presentan resultados de ensayos de presión superficial sobre un modelo de la isla. Los principales resultados del experimento de Bolund se reproducen en el capítulo cinco. En el capítulo seis se identifican diferentes estructuras del flujo sobre la isla y, finalmente, en el capitulo siete, se recogen las conclusiones y una propuesta de lineas de trabajo futuras. ABSTRACT The main objective of this work is to contribute to answer the question: to which extend can the wind tunnel testing contribute to determine the flow characteristics that affect the dynamic response of wind turbines operating in highly complex terrains?. This question is not new, indeed, the debate in the scientific community was opened in the first third of the past century and it is still intensely alive. The accepted approach to face this problem consists in analysing a given case study where full-scale tests, computational modelling and wind tunnel testing are applied to the same topography. This is neither easy nor cheap. This is is the reason why since the Askervein experience in 1988, the atmospheric flow modellers community had to wait till 2007 when the Bolund experiment was setup with a deployment of technical means equivalent (considering the evolution of the sensor and computing techniques). The problem is so manifold that both experiences were restricted to neutral conditions without Coriolis effects in order to reduce the complexity. This is the framework in which this PhD has been carried out. The flow topology over the Bolund Island has been studied by replicating the Bolund experiment in the IDR A9 and ACLA16 wind tunnels. Two mock-ups of the Bolund island were manufactured at two scales of 1:230 and 1:115. The in-flow in the empty wind tunnel simulating the incoming atmospheric boundary layer was in the transitionally rough regime and used as a reference case. The 1:230 model was tested in the A9 wind tunnel to measure surface pressure. The mapping of the pressure coefficient across the island gave a visualisation and estimation of a detachment region on the top of the escarpment in front of the island. Time resolved instantaneous pressure measurements illustrated the non-steadiness in the detachment region. The 1:115 model was tested using 3C hot-wires(HW) and 2C Particle Image Velocimetry(PIV). Measurements at met masts M3, M6, M7 and M8 and along Line 270°were taken to replicate the result of the Bolund experiment. The flow was characterised by the speed-up ratio, normalised increment of the turbulent kinetic energy, inclination angle and turning angle. Results along line 270°at heights of 2 m and 5 m compared very well with the full-scale results of the Bolund experiment. Vertical profiles at the met masts showed a significant agreement with the full-scale results. The analysis of the Reynolds stresses and the spectral analysis at the met mast locations gave a varied level of agreement at some locations while clear mismatch at others. The horizontal mapping of the flow field, for a 270°wind direction, allowed to characterise the behaviour of the intermittent recirculation bubble on top of the front escarpment followed by a relaxation region and the presence of a shear layer in the lee side of the island. Further detailed velocity measurements were taken at cross-flow planes over the island to study the flow structures on the island. A longitudinal vortex-like structure with high mean velocity gradients and high turbulent kinetic energy was characterised on the escarpment and evolving downstream. This flow structure is a challenge to the numerical models while posing a threat to wind farm designers when siting wind turbines. Spatial distribution of Reynold stresses were presented from 3C HW and PIV measurements. These values are not common results from usual wind tunnel measurements and very useful for modellers using large eddy simulation (LES). An interpretation of the wind tunnel results in terms of usefulness to wind farm designers is given. Evolution and variation of the flow parameters along measurement lines, planes and surfaces indicated how the flow field could affect wind turbine siting. Different flow properties were presented so compare the level of agreement to full-scale results and how this affected when characterising the site wind classes. The results presented suggest, under certain conditions, the robustness of the wind tunnel testing for studying flow topology over complex terrain and its capability to compare to other modelling techniques especially from the level of agreement between the different data sets presented. Additionally, some flow parameters obtained from wind tunnel measurements would have been quite difficult to be measured at full-scale or by computational means considering the state of the art. This work was carried out as a part of the activities supported by the EC as part of the FP7- PEOPLE-ITN-2008 WAUDIT project (Wind Resource Assessment Audit and Standardization) within the FP7 Marie-Curie Initial Training Network and by the Spanish Ministerio de Economía y Competitividad, within the framework of the ENE2012-36473, TURCO project (Determination of the Spatial Distribution of Statistic Parameters of Flow Turbulence over Complex Topographies in Wind Tunnel) belonging to the Spanish National Program of Research (Subprograma de investigación fundamental no orientada 2012). The report is organised in seven chapters and a collection of annexes. In chapter one, the problem is introduced. In chapter two the experimental setup is described. Following, in chapter three, the inflow conditions of the main wind tunnel used in this piece of research are analysed in detail. In chapter three, preliminary pressure tests results on a model of the island are presented. The main results from the Bolund experiment are replicated in chapter five. In chapter six, an identification of specific flow strutures over the island is presented and, finally, in chapter seven, conclusions and lines for future works related to the presented one are included.

Determination of Maximum Mechanical Energy Efficiency in Energy Galloping Systems

Transverse galloping is a type of aeroelastic instability characterized by large amplitude, low frequency, normal to wind oscillations that appear in some elastic two-dimensional bluff bodies when subjected to a fluid flow, provided that the flow velocity exceeds a threshold critical value. Such an oscillatory motion is explained because of the energy transfer from the flow to the two-dimensional bluff body. The 7 amount of energy that can be extracted depends on the cross section of the galloping prism. Assuming that the Glauert-Den Hartog quasistatic criterion for galloping instability is satisfied in a first approximation, the suitability of a given cross section for energy harvesting is evaluated by analyzing the lateral aerodynamic force coefficient, fitting a function with a power series in tan a (a being the angle of attack) to 10 available experimental data. In this paper, a fairly large number of simple prisms (triangle, ellipse, biconvex, and rhombus cross sections, as well 11 as D-shaped bodies) is analyzed for suitability as energy harvesters. The influence of the fitting process in the energy harvesting efficiency evaluation is also demonstrated. The analysis shows that the more promising bodies are those with isosceles or approximate isosceles cross sections.

Global instability analysis and control of three-dimensional long open cavity flow

The three-dimensional wall-bounded open cavity may be considered as a simplified geometry found in industrial applications such as leading gear or slotted flats on the airplane. Understanding the three-dimensional complex flow structure that surrounds this particular geometry is therefore of major industrial interest. At the light of the remarkable former investigations in this kind of flows, enough evidences suggest that the lateral walls have a great influence on the flow features and hence on their instability modes. Nevertheless, even though there is a large body of literature on cavity flows, most of them are based on the assumption that the flow is two-dimensional and spanwise-periodic. The flow over realistic open cavity should be considered. This thesis presents an investigation of three-dimensional wall-bounded open cavity with geometric ratio 6:2:1. To this aim, three-dimensional Direct Numerical Simulation (DNS) and global linear instability have been performed. Linear instability analysis reveals that the onset of the first instability in this open cavity is around Recr 1080. The three-dimensional shear layer mode with a complex structure is shown to be the most unstable mode. I t is noteworthy that the flow pattern of this high-frequency shear layer mode is similar to the observed unstable oscillations in supercritical unstable case. DNS of the cavity flow carried out at different Reynolds number from steady state until a nonlinear saturated state is obtained. The comparison of time histories of kinetic energy presents a clearly dominant energetic mode which shifts between low-frequency and highfrequency oscillation. A complete flow patterns from subcritical cases to supercritical case has been put in evidence. The flow structure at the supercritical case Re=1100 resembles typical wake-shedding instability oscillations with a lateral motion existed in the subcritical cases. Also, This flow pattern is similar to the observations in experiments. In order to validate the linear instability analysis results, the topology of the composite flow fields reconstructed by linear superposition of a three-dimensional base flow and its leading three-dimensional global eigenmodes has been studied. The instantaneous wall streamlines of those composited flows display distinguish influence region of each eigenmode. Attention has been focused on the leading high-frequency shear layer mode; the composite flow fields have been fully recognized with respect to the downstream wave shedding. The three-dimensional shear layer mode is shown to give rise to a typical wake-shedding instability with a lateral motions occurring downstream which is in good agreement with the experiment results. Moreover, the spanwise-periodic, open cavity with the same length to depth ratio has been also studied. The most unstable linear mode is different from the real three-dimensional cavity flow, because of the existence of the side walls. Structure sensitivity of the unstable global mode is analyzed in the flow control context. The adjoint-based sensitivity analysis has been employed to localized the receptivity region, where the flow is more sensible to momentum forcing and mass injection. Because of the non-normality of the linearized Navier-Stokes equations, the direct and adjoint field has a large spatial separation. The strongest sensitivity region is locate in the upstream lip of the three-dimensional cavity. This numerical finding is in agreement with experimental observations. Finally, a prototype of passive flow control strategy is applied.

Motivación, flow e intención de práctica físico-deportiva en estudiantes de educación física de España y Latinoamérica

A pesar de los innegables beneficios físicos y psicosociales que la práctica de actividad física tiene para la población en las últimas décadas se observa un importante descenso en el nivel de práctica de actividad física durante la adolescencia. La Educación Física (EF) es una asignatura presente en las primeras etapas educativas de los currículums de todos los países y supone para la población infantil y adolescente, durante el periodo de escolarización obligatoria, un entorno de contacto cotidiano con la actividad física y deportiva lo que explica el creciente interés que existe por investigar en este contexto y extraer conclusiones que puedan ayudar a revertir el proceso de progresivo abandono de la práctica en estas edades. El presente trabajo se centra en un periodo, la adolescencia, especialmente sensible al abandono de la práctica de actividad física. Apoyándonos en tres de los marcos teóricos más importantes de los últimos años - Teoría de Metas de Logro, Teoría de la Autodeterminación y Teoría del Flow - hemos realizado un diseño con tres estudios complementarios que nos permitan avanzar en el conocimiento de las variables que afectan a esta realidad. En el primero de los estudios, en el que participaron 3990 estudiantes de EF de España, Argentina, Colombia y Ecuador con edades comprendidas entre 12 y 18 años, se fijaron como objetivos el establecimiento de perfiles motivacionales en estudiantes de EF considerando las orientaciones motivacionales, la motivación intrínseca y el flow disposicional; la comprobación de la posible estabilidad de los perfiles establecidos en cuatro países diferentes; y el análisis de la relación entre los perfiles motivacionales establecidos con la práctica de actividad física realizada en la actualidad y la intención de ser activo en el futuro. Los resultados revelaron la existencia de tres perfiles motivacionales en los estudiantes: perfil de motivación baja-moderada, perfil de motivación de alta motivación con bajo ego y perfil de alta motivación. Se encontraron similitudes en los perfiles motivacionales obtenidos en los cuatro países analizados, sugiriéndose la posibilidad de que la lengua materna compartida pueda ser una variable de influencia en la existencia de patrones conductuales similares. La orientación al ego resultó ser un elemento diferenciador en la relación de los dos perfiles que mostraron niveles más altos de motivación intrínseca, flow disposicional y orientación a la tarea con las variables de práctica. Mayores niveles de orientación al ego determinaron mayores valores de práctica e intención de práctica. En el segundo estudio, en el que participaron 365 estudiantes de EF de España, Argentina y Colombia con edades comprendidas entre 12 y 16 años, se testó un modelo explicativo de la intención de ser físicamente activo en el futuro en estudiantes de EF poco motivados introduciendo el flow disposicional como variable explicativa junto a las orientaciones motivacionales y la motivación intrínseca. El modelo de ecuaciones estructurales testado muestra que tanto la orientación al ego como la orientación a la tarea y la motivación intrínseca predijeron el flow disposicional. No se encontró un efecto significativo de la orientación a la tarea sobre la motivación intrínseca, ni de esta sobre la intención de práctica de AFD en el futuro. El flow disposicional resultó ser un mediador entre ambas orientaciones motivacionales y dicha intención, adquiriendo un papel relevante en la predicción de la intención de realizar actividad física en el futuro en los sujetos con poca motivación hacia la clase de EF. En el último estudio, con una muestra de 53 estudiantes de EF pertenecientes a dos clases de segundo curso de ESO con edades comprendidas entre 13 y 15 años, se analizó el efecto de una intervención de apoyo a las NPB sobre la satisfacción de las mismas, la motivación intrínseca, la disposición a experimentar flow y la intención de ser físicamente activo en el futuro; estudiándose también el efecto de la misma sobre ciertas variables presentes en el proceso de enseñanza-aprendizaje tales como la diversión, el aprendizaje y el gusto hacia el contenido y la metodología empleada. La intervención resultó efectiva en la mejora de todas las variables motivacionales analizadas a excepción de la satisfacción de la necesidad de relaciones sociales, variable que podría ser más difícilmente modificable en tanto que los grupos de clase suelen permanecer estables a lo largo de esta etapa educativa. Los alumnos de la condición experimental presentaron niveles más elevados de diversión, aprendizaje y gusto hacia el contenido y la metodología empleada. ABSTRACT Despite the undeniable physical and psychological benefits that physical activity has on health, a significant drop in the level of physical activity during adolescence has been observed over the last few decades. Physical education (PE) is present in the early stages of educational curricula of all countries, creating an environment in which children and adolescents can maintain daily contact with sport and physical activity during the period of compulsory education. There is thus a growing interest in research within this context, helping to draw conclusions aimed at reversing the progressive decline of physical activity in this age bracket. The work presented here focuses on adolescence – a period particularly sensitive when faced with declining rates of physical activity. Relying on three of the most important recent theoretical frameworks – Achievement Goal Theory, Self-Determination Theory and Flow Theory – we have carried out a design using three complementary studies that allow us to advance the understanding of variables affecting this reality. The first two studies employ a correlational methodology while the third follows a quasi-experimental design. The objectives of the first study, with the participation of 3990 PE students aged 12 to 18 from Spain, Argentina, Colombia and Ecuador, were: a) to establish motivational profiles in physical education students taking into consideration motivational orientations, intrinsic motivation and dispositional flow; b) to check the possible stability of the established profiles in the four countries; and c) to analyze the relationship between motivational profiles and both current levels of physical activity and intention to be physically active in the future. Results revealed the existence of three motivational profiles: low-moderate motivation profile, high motivation and low ego profile and high motivation profile. Similarities were found across the four countries, which presents the possibility that sharing the same mother tongue might influence the existence of similar behavioral patterns. Amongst the profiles that showed the highest levels of intrinsic motivation, dispositional flow and task orientation with practice variables, ego orientation became apparent as a differentiating element. Higher levels of ego orientation correlated with higher levels both in current level of physical activity and intention to be physically active in the future. In the second study, with the participation of 365 PE students aged 12 to 16 from Spain, Argentina and Colombia, an explanatory model of the intention to be physically active within the low-moderate motivation profile was tested, including dispositional flow as an influence variable together with motivational orientations and intrinsic motivation. The structural model equation in the low-moderate motivation profile shows some differences with respect to findings of similar studies using heterogeneous samples in motivational profiles. It was found that ego orientation, task orientation and intrinsic motivation predicted dispositional flow. No significant effects of task orientation on intrinsic motivation, or of intrinsic motivation on intention to be physically active in the future were found. Dispositional flow emerged as a mediator between both motivational orientations and intention, acquiring an important role in the prediction of future intentions to be physically active in adolescents who show low levels of motivation towards physical education. In the last of the three studies, with the participation of 53 PE students drawn from two classes in the second year of secondary education aged 13 to 15, the effect of an intervention to support basic psychological needs was tested against the satisfaction of these needs, intrinsic motivation, disposition to experience flow and intention to be physically active in the future. Alongside this, the aforementioned intervention was tested on certain variables present in the teaching-learning process such as fun, learning and taste for the activity and methodology. Supportive intervention was effective in improving all the motivational variables, except for the satisfaction of the relatedness variable, which may be more difficult to modify as class groups tend to remain stable throughout this educational state. Students in the experimental condition showed higher levels of fun, learning and taste for the content and methodology.

Automated wordlength optimization framework for multi-source statistical interval-based analysis of nonlinear systems with control-flow structures

El uso de aritmética de punto fijo es una opción de diseño muy extendida en sistemas con fuertes restricciones de área, consumo o rendimiento. Para producir implementaciones donde los costes se minimicen sin impactar negativamente en la precisión de los resultados debemos llevar a cabo una asignación cuidadosa de anchuras de palabra. Encontrar la combinación óptima de anchuras de palabra en coma fija para un sistema dado es un problema combinatorio NP-hard al que los diseñadores dedican entre el 25 y el 50 % del ciclo de diseño. Las plataformas hardware reconfigurables, como son las FPGAs, también se benefician de las ventajas que ofrece la aritmética de coma fija, ya que éstas compensan las frecuencias de reloj más bajas y el uso más ineficiente del hardware que hacen estas plataformas respecto a los ASICs. A medida que las FPGAs se popularizan para su uso en computación científica los diseños aumentan de tamaño y complejidad hasta llegar al punto en que no pueden ser manejados eficientemente por las técnicas actuales de modelado de señal y ruido de cuantificación y de optimización de anchura de palabra. En esta Tesis Doctoral exploramos distintos aspectos del problema de la cuantificación y presentamos nuevas metodologías para cada uno de ellos: Las técnicas basadas en extensiones de intervalos han permitido obtener modelos de propagación de señal y ruido de cuantificación muy precisos en sistemas con operaciones no lineales. Nosotros llevamos esta aproximación un paso más allá introduciendo elementos de Multi-Element Generalized Polynomial Chaos (ME-gPC) y combinándolos con una técnica moderna basada en Modified Affine Arithmetic (MAA) estadístico para así modelar sistemas que contienen estructuras de control de flujo. Nuestra metodología genera los distintos caminos de ejecución automáticamente, determina las regiones del dominio de entrada que ejercitarán cada uno de ellos y extrae los momentos estadísticos del sistema a partir de dichas soluciones parciales. Utilizamos esta técnica para estimar tanto el rango dinámico como el ruido de redondeo en sistemas con las ya mencionadas estructuras de control de flujo y mostramos la precisión de nuestra aproximación, que en determinados casos de uso con operadores no lineales llega a tener tan solo una desviación del 0.04% con respecto a los valores de referencia obtenidos mediante simulación. Un inconveniente conocido de las técnicas basadas en extensiones de intervalos es la explosión combinacional de términos a medida que el tamaño de los sistemas a estudiar crece, lo cual conlleva problemas de escalabilidad. Para afrontar este problema presen tamos una técnica de inyección de ruidos agrupados que hace grupos con las señales del sistema, introduce las fuentes de ruido para cada uno de los grupos por separado y finalmente combina los resultados de cada uno de ellos. De esta forma, el número de fuentes de ruido queda controlado en cada momento y, debido a ello, la explosión combinatoria se minimiza. También presentamos un algoritmo de particionado multi-vía destinado a minimizar la desviación de los resultados a causa de la pérdida de correlación entre términos de ruido con el objetivo de mantener los resultados tan precisos como sea posible. La presente Tesis Doctoral también aborda el desarrollo de metodologías de optimización de anchura de palabra basadas en simulaciones de Monte-Cario que se ejecuten en tiempos razonables. Para ello presentamos dos nuevas técnicas que exploran la reducción del tiempo de ejecución desde distintos ángulos: En primer lugar, el método interpolativo aplica un interpolador sencillo pero preciso para estimar la sensibilidad de cada señal, y que es usado después durante la etapa de optimización. En segundo lugar, el método incremental gira en torno al hecho de que, aunque es estrictamente necesario mantener un intervalo de confianza dado para los resultados finales de nuestra búsqueda, podemos emplear niveles de confianza más relajados, lo cual deriva en un menor número de pruebas por simulación, en las etapas iniciales de la búsqueda, cuando todavía estamos lejos de las soluciones optimizadas. Mediante estas dos aproximaciones demostramos que podemos acelerar el tiempo de ejecución de los algoritmos clásicos de búsqueda voraz en factores de hasta x240 para problemas de tamaño pequeño/mediano. Finalmente, este libro presenta HOPLITE, una infraestructura de cuantificación automatizada, flexible y modular que incluye la implementación de las técnicas anteriores y se proporciona de forma pública. Su objetivo es ofrecer a desabolladores e investigadores un entorno común para prototipar y verificar nuevas metodologías de cuantificación de forma sencilla. Describimos el flujo de trabajo, justificamos las decisiones de diseño tomadas, explicamos su API pública y hacemos una demostración paso a paso de su funcionamiento. Además mostramos, a través de un ejemplo sencillo, la forma en que conectar nuevas extensiones a la herramienta con las interfaces ya existentes para poder así expandir y mejorar las capacidades de HOPLITE. ABSTRACT Using fixed-point arithmetic is one of the most common design choices for systems where area, power or throughput are heavily constrained. In order to produce implementations where the cost is minimized without negatively impacting the accuracy of the results, a careful assignment of word-lengths is required. The problem of finding the optimal combination of fixed-point word-lengths for a given system is a combinatorial NP-hard problem to which developers devote between 25 and 50% of the design-cycle time. Reconfigurable hardware platforms such as FPGAs also benefit of the advantages of fixed-point arithmetic, as it compensates for the slower clock frequencies and less efficient area utilization of the hardware platform with respect to ASICs. As FPGAs become commonly used for scientific computation, designs constantly grow larger and more complex, up to the point where they cannot be handled efficiently by current signal and quantization noise modelling and word-length optimization methodologies. In this Ph.D. Thesis we explore different aspects of the quantization problem and we present new methodologies for each of them: The techniques based on extensions of intervals have allowed to obtain accurate models of the signal and quantization noise propagation in systems with non-linear operations. We take this approach a step further by introducing elements of MultiElement Generalized Polynomial Chaos (ME-gPC) and combining them with an stateof- the-art Statistical Modified Affine Arithmetic (MAA) based methodology in order to model systems that contain control-flow structures. Our methodology produces the different execution paths automatically, determines the regions of the input domain that will exercise them, and extracts the system statistical moments from the partial results. We use this technique to estimate both the dynamic range and the round-off noise in systems with the aforementioned control-flow structures. We show the good accuracy of our approach, which in some case studies with non-linear operators shows a 0.04 % deviation respect to the simulation-based reference values. A known drawback of the techniques based on extensions of intervals is the combinatorial explosion of terms as the size of the targeted systems grows, which leads to scalability problems. To address this issue we present a clustered noise injection technique that groups the signals in the system, introduces the noise terms in each group independently and then combines the results at the end. In this way, the number of noise sources in the system at a given time is controlled and, because of this, the combinato rial explosion is minimized. We also present a multi-way partitioning algorithm aimed at minimizing the deviation of the results due to the loss of correlation between noise terms, in order to keep the results as accurate as possible. This Ph.D. Thesis also covers the development of methodologies for word-length optimization based on Monte-Carlo simulations in reasonable times. We do so by presenting two novel techniques that explore the reduction of the execution times approaching the problem in two different ways: First, the interpolative method applies a simple but precise interpolator to estimate the sensitivity of each signal, which is later used to guide the optimization effort. Second, the incremental method revolves on the fact that, although we strictly need to guarantee a certain confidence level in the simulations for the final results of the optimization process, we can do it with more relaxed levels, which in turn implies using a considerably smaller amount of samples, in the initial stages of the process, when we are still far from the optimized solution. Through these two approaches we demonstrate that the execution time of classical greedy techniques can be accelerated by factors of up to ×240 for small/medium sized problems. Finally, this book introduces HOPLITE, an automated, flexible and modular framework for quantization that includes the implementation of the previous techniques and is provided for public access. The aim is to offer a common ground for developers and researches for prototyping and verifying new techniques for system modelling and word-length optimization easily. We describe its work flow, justifying the taken design decisions, explain its public API and we do a step-by-step demonstration of its execution. We also show, through an example, the way new extensions to the flow should be connected to the existing interfaces in order to expand and improve the capabilities of HOPLITE.

Determination of the lifetime and force dependence of interactions of single bonds between surface-attached CD2 and CD48 adhesion molecules

We studied single molecular interactions between surface-attached rat CD2, a T-lymphocyte adhesion receptor, and CD48, a CD2 ligand found on antigen-presenting cells. Spherical particles were coated with decreasing densities of CD48–CD4 chimeric molecules then driven along CD2-derivatized glass surfaces under a low hydrodynamic shear rate. Particles exhibited multiple arrests of varying duration. By analyzing the dependence of arrest frequency and duration on the surface density of CD48 sites, it was concluded that (i) arrests were generated by single molecular bonds and (ii) the initial bond dissociation rate was about 7.8 s−1. The force exerted on bonds was increased from about 11 to 22 pN; the detachment rate exhibited a twofold increase. These results agree with and extend studies on the CD2–CD48 interaction by surface plasmon resonance technology, which yielded an affinity constant of ≈104 M−1 and a dissociation rate of ≥6 s−1. It is concluded that the flow chamber technology can be an useful complement to atomic force microscopy for studying interactions between isolated biomolecules, with a resolution of about 20 ms and sensitivity of a few piconewtons. Further, this technology might be extended to actual cells.

Desenvolvimento e avaliação de métodos de extração e separação cromatográfica em colunas monolíticas e superficialmente porosas para determinação de herbicidas triazínicos em solos e águas

O uso de pesticidas levou ao aumento da produtividade e qualidade dos produtos agrícolas, porém o seu uso acarreta na intoxicação dos seres vivos pela ingestão gradativa de seus resíduos que contaminam o solo, a água e os alimentos. Dessa forma, há a necessidade do monitoramento constante de suas concentrações nos compartimentos ambientais. Para isto, busca-se o desenvolvimento de métodos de extração e enriquecimento de forma rápida, com baixo custo, gerando um baixo volume de resíduos, contribuindo com a química verde. Dentre estes métodos destacam-se a extração por banho de ultrassom e a extração por ponto nuvem. Após o procedimento de extração, o extrato obtido pode ser analisado por técnicas de Cromatografia a Líquido de Alta Eficiência (HPLC) e a Cromatografia por Injeção Sequencial (SIC), empregando fases estacionárias modernas, tais como as monolíticas e as partículas superficialmente porosas. O emprego de SIC com coluna monolítica (C18, 50 x 4,6 mm) e empacotada com partículas superficialmente porosas (C18, 30 x 4,6 mm, tamanho de partícula 2,7 µm) foi estudado para separação de simazina (SIM) e atrazina (ATR), e seus metabólitos, desetilatrazina (DEA), desisopropilatrazina (DIA) e hidroxiatrazina (HAT). A separação foi obtida por eluição passo-a-passo, com fases móveis compostas de acetonitrila (ACN) e tampão Acetato de Amônio/Ácido acético (NH4Ac/HAc) 2,5 mM pH 4,2. A separação na coluna monolítica foi realizada com duas fases móveis: MP1= 15:85 (v v-1) ACN:NH4Ac/HAc e MP2= 35:65 (v v-1) ACN:NH4Ac/HAc a uma vazão de 35 µL s-1. A separação na coluna com partículas superficialmente porosas foi efetivada com as fases móveis MP1= 13:87 (v v-1) ACN: NH4Ac/HAc e MP2= 35:65 (v v-1) ACN:NH4Ac/HAc à vazão de 8 µL s-1. A extração por banho de ultrassom em solo fortificado com os herbicidas (100 e 1000 µg kg-1) resultou em recuperações entre 42 e 160%. A separação de DEA, DIA, HAT, SIM e ATR empregando HPLC foi obtida por um gradiente linear de 13 a 35% para a coluna monolítica e de 10 a 35% ACN na coluna com partículas superficialmente porosas, sendo a fase aquosa constituída por tampão NH4Ac/HAc 2,5 mM pH 4,2. Em ambas as colunas a vazão foi de 1,5 mL min-1 e o tempo de análise 15 min. A extração por banho de ultrassom das amostras de solo com presença de ATR, fortificadas com concentrações de 250 a 1000 µg kg-1, proporcionou recuperações entre 40 e 86%. A presença de ATR foi confirmada por espectrometria de massas. Foram realizados estudos de fortificação com ATR e SIM em amostras de água empregando a extração por ponto nuvem com o surfactante Triton-X114. A separação empregando HPLC foi obtida por um gradiente linear de 13 a 90% de ACN para a coluna monolítica e de 10 a 90% de ACN para a coluna empacotada, sempre em tampão NH4Ac/HAc 2,5 mM pH 4,2. Em ambas as colunas a vazão foi de 1,5 mL min-1 e o tempo de análise 16 min. Fortificações entre 1 e 50 µg L-1 resultaram em recuperações entre 65 e 132%.

Coagent Modified Polypropylene Prepared by Reactive Extrusion: A New Look Into the Structure-Property Relations of Injection Molded Parts

Peroxide-mediated reactive extrusion of linear isotactic polypropylene (L-PP) was conducted in the presence of trimethylolpropane trimethacrylate (TMPTMA) and triallyl trimesate (TAM) coagents, using a twin screw extruder. The resulting coagent-modified polypropylenes (CM-PP) had higher viscosities and elasticities, as well as increased crystallization temperature compared to PP reacted only with peroxide (DCP-PP). Additionally, deviations from terminal flow, and strain hardening were observed in PP modified with TAM, signifying the presence of long chain branching (LCB). The CM-PP formulations retained the modulus and tensile strength of the parent L-PP, in spite of their lower molar mass and viscosities, whereas their elongation at break and the impact strength were better. This was attributed to the finer spherulitic structure of these materials, and to the disappearance of the skin-core layer in the injection molded specimens.

(Table 1) Age determination of sediment core MD02-2589

Sortable silt mean grain sizes together with oxygen and carbon isotopic data produced on the benthic foraminiferal species Fontbotia wuellerstorfi are used to construct high-resolution records of near-bottom flow vigour and deep water ventilation at a core site MD02-2589 located at 2660 m water depth on the southern Agulhas Plateau. The results suggest that during glacial periods (marine oxygen isotope stages 2 and 6, MIS 2 and MIS 6, respectively), there was a persistent contribution of a well-ventilated water mass within the Atlantic to Indian oceanic gateway with a d13C signature similar to present-day Northern Component Water (NCW), e.g., North Atlantic Deep Water (NADW). The records of chemical ventilation and near-bottom flow vigor reflect changes in the advection of northern source waters and meridional variability in the location of the Antarctic Circumpolar Current and its associated fronts. We suggest that during Termination II (TII), changes in chemical ventilation are largely decoupled from near-bottom physical flow speeds. A mid-TII climate optimum is associated with a low-flow speed plateau concurrent with a period of increased ventilation shown in the benthic d13C of other Southern Ocean records but not in our benthic d13C of MD02-2589. The climate optimum is followed by a period of southern cooling around 128 ka coincident with a stronger influence of NCW to interglacial levels at around 124 ka. All proxy records show a near synchronous and rapid shift during the transition from MIS 5a-4 (73 ka). This large event is attributed to a rapid decrease in NADW influence and replacement over the Agulhas Plateau by southern source waters.

(Table 2) Age determination of sediment cores from the Panama Basin

Application of the 230Th normalization method to estimate sediment burial fluxes in six cores from the eastern equatorial Pacific (EEP) reveals that bulk sediment and organic carbon fluxes display a coherent regional pattern during the Holocene that is consistent with modern oceanographic conditions, in contrast with estimates of bulk mass accumulation rates (MARs) derived from core chronologies. Two nearby sites (less than 10 km apart), which have different MARs, show nearly identical 230Th-normalized bulk fluxes. Focusing factors derived from the 230Th data at the foot of the Carnegie Ridge in the Panama Basin are >2 in the Holocene, implying that lateral sediment addition is significant in this part of the basin. New geochemical data and existing literature provide evidence for a hydrothermal source of sediment in the southern part of the Panama Basin and for downslope transport from the top of the Carnegie Ridge. The compilation of core records suggests that sediment focusing is spatially and temporally variable in the EEP. During oxygen isotope stage 2 (OIS 2, from 13-27 ka BP), focusing appears even higher compared to the Holocene at most sites, similar to earlier findings in the eastern and central equatorial Pacific. The magnitude of the glacial increase in focusing factors, however, is strongly dependent on the accuracy of age models. We offer two possible explanations for the increase in glacial focusing compared to the Holocene. The first one is that the apparent increase in lateral sediment redistribution is partly or even largely an artifact of insufficient age control in the EEP, while the second explanation, which assumes that the observed increase is real, involves enhanced deep sea tidal current flow during periods of low sea level stand.

(Table 1) Summary of basalt unit thickness, vesicularity, number of flow units, and average porosity at ODP Site 135-839

An inflatable drill-string packer was used at Site 839 to measure the bulk in-situ permeability within basalts cored in Hole 839B. The packer was inflated at two depths, 398.2 and 326.9 mbsf; all on-board information indicated that the packer mechanically closed off the borehole, although apparently the packer hydraulically sealed the borehole only at 398.2 mbsf. Two pulse tests were run at each depth, two constant-rate injection tests were run at the first set, and four were run at the second. Of these, only the constant-rate injection tests at the first set yielded a permeability, calculated as ranging from 1 to 5 * 10**-12 m**2. Pulse tests and constant-rate injection tests for the second set did not yield valid data. The measured permeability is an upper limit; if the packer leaked during the experiments, the basalt would be less permeable. In comparison, permeabilities measured at other Deep Sea Drilling Project and Ocean Drilling Program sites in pillow basalts and flows similar to those measured in Hole 839B are mainly about 10**-13 to 10**-14 m**2. Thus, if our results are valid, the basalts at Site 839 are more permeable than ocean-floor basalts investigated elsewhere. Based on other supporting evidence, we consider these results to be a valid measure of the permeability of the basalts. Temperature data and the geochemical and geotechnical properties of the drilled sediments all indicate that the site is strongly affected by fluid flow. The heat flow is very much less than expected in young oceanic basalts, probably a result of rapid fluid circulation through the crust. The geochemistry of pore fluids is similar to that of seawater, indicating seawater flow through the sediments, and sediments are uniformly underconsolidated for their burial depth, again indicating probable fluid flow. The basalts are highly vesicular. However, the vesicularity can only account for part of the average porosity measured on the neutron porosity well log; the remainder of the measured porosity is likely present as voids and fractures within and between thin-bedded basalts. Core samples, together with porosity, density, and resistivity well-log data show locations where the basalt section is thin bedded and probably has from 15% to 35% void and fracture porosity. Thus, the measured permeability seems reasonable with respect to the high measured porosity. Much of the fluid flow at Site 839 could be directed through highly porous and permeable zones within and between the basalt flows and in the sediment layer just above the basalt. Thus, the permeability measurements give an indication of where and how fluid flow may occur within the oceanic crust of the Lau Basin.

(Table 1) Age determination of sediment cores of the Labrador Sea

A close look at the sedimentology of Heinrich event 4 from the northwest Labrador Sea indicates that an extended ice margin, perhaps greater than before Heinrich events 1 or 2 (H-1 and H-2), existed in the Hudson Strait region pre-Heinrich event 4 (H-4) and, that on the basis of characteristics of the sediment unit, Heinrich event-4 was different than Heinrich events 1 or 2 (i.e., larger ice sheet collapse(?), longer duration(?), "dirtier" icebergs(?)). Other data from across the southern and eastern margin of the Laurentide Ice Sheet, as well as Greenland and the North Atlantic, support this interpretation, possibly indicating a relative mid-Wisconsin glacial maximum pre-Heinrich event 4. Many of these data also indicate that Heinrich event 4 (35 ka) resulted in serious climatic and oceanographic reorganizations. We suggest that Heinrich event 4 gutted the Hudson Strait, leaving it devoid of ice for Heinrich event 3. We further hypothesize that Heinrich event 3 did not originate from axial ice transport along the Hudson Strait; thus Heinrich event 3 may be more analogous to the proposed northward advancing ice from Ungava Bay during Heinrich event 0 than to the more typical down-the-strait flow during H-1, H-2, and H-4. Consequently, the climatic and oceanographic impacts resulting from Heinrich events are highly susceptible to the type, origin, and magnitude of ice sheet collapse, something which varied per Heinrich event during the last glacial period.

Simplified techniques for determination of theoretical and effective self-purification coefficients of polluted streams /

Mode of access: Internet.

Electrochemical deposition of cupric ions for the determination of mass transfer rates to a stepped rotating cylinder electrode

The stepped rotating cylinder electrode (SRCE) geometry has been developed as a simple aid to the practical study of the flow-enhanced corrosion and applied electrochemistry problems commonly observed under conditions of disturbed, turbulent flow. The electrodeposition of cupric ions from an acid sulphate plating bath has been used to characterise differential rates of mass transfer to the SRCE. The variation in thickness of electrodeposited copperfilms has allowed the mapping of local rates of mass transfer over the active surface of this geometry. Both optical and scanning electron microscopy were used for the examination of metallographic sections to provide a high resolution evaluation of the distribution of mass transfer coefficient. Results are also discussed using the convective-diffusion model in combination with the existing direct numerical flow simulation (DNS) data for this geometry.

Determination of lifter design, speed and filling effects in AG mills by 3D DEM

The power required to operate large gyratory mills often exceeds 10 MW. Hence, optimisation of the power consumption will have a significant impact on the overall economic performance and environmental impact of the mineral processing plant. In most of the published models of tumbling mills (e.g. [Morrell, S., 1996. Power draw of wet tumbling mills and its relationship to charge dynamics, Part 2: An empirical approach to modelling of mill power draw. Trans. Inst. Mining Metall. (Section C: Mineral Processing Ext. Metall.) 105, C54-C62. Austin, L.G., 1990. A mill power equation for SAG mills. Miner. Metall. Process. 57-62]), the effect of lifter design and its interaction with mill speed and filling are not incorporated. Recent experience suggests that there is an opportunity for improving grinding efficiency by choosing the appropriate combination of these variables. However, it is difficult to experimentally determine the interactions of these variables in a full scale mill. Although some work has recently been published using DEM simulations, it was basically. limited to 2D. The discrete element code, Particle Flow Code 3D (PFC3D), has been used in this work to model the effects of lifter height (525 cm) and mill speed (50-90% of critical) on the power draw and frequency distribution of specific energy (J/kg) of normal impacts in a 5 m diameter autogenous (AG) mill. It was found that the distribution of the impact energy is affected by the number of lifters, lifter height, mill speed and mill filling. Interactions of lifter design, mill speed and mill filling are demonstrated through three dimensional distinct element methods (3D DEM) modelling. The intensity of the induced stresses (shear and normal) on lifters, and hence the lifter wear, is also simulated. (C) 2004 Elsevier Ltd. All rights reserved.