Global Stability Analysis of a Compressible Turbulent Flow around a High-Lift Configuration

The purpose of this work is to analyze a complex high lift configuration for which significant regions of separated flow are present. Current state of the art methods have some diffculty to predict the origin and the progression of this separated flow when increasing the angle of attack. The mechanisms responsible for the maximum lift limit on multi-element wing con?gurations are not clear; this stability analysis could help to understand the physics behind the phenomenon and to find a relation between the flow separation and the instability onset. The methodology presented herein consists in the computation of a steady base flow solution based on a finite volume discretization and a proposal of the solution for a generalized eigenvalue problem corresponding to the perturbed and linearized problem. The eigenvalue problem has been solved with the Arnoldi iterative method, one of the Krylov subspace projection methods. The described methodology was applied to the NACA0012 test case in subsonic and in transonic conditions and, finally, for the first time to the authors knowledge, on an industrial multi-component geometry, such as the A310 airfoil, in order to identify low frequency instabilities related to the separation. One important conclusion is that for all the analyzed geometries, one unstable mode related to flow separation appears for an angle of attack greater than the one correspondent to the maximum lift coe?cient condition. Finally, an adjoint study was carried out in order to evaluate the receptivity and the structural sensitivity of the geometries, giving an indication of the domain region that could be modified resulting in the biggest change of the flowfield.

Dynamic interactions between solids and viscous liquids with free surface = Interacciones dinámicas entre sólidos y líquidos viscosos con superficie libre

En esta tesis se investiga la interacción entre un fluido viscoso y un cuerpo sólido en presencia de una superficie libre. El problema se expresa teóricamente poniendo especial atención a los aspectos de conservación de energía y de la interacción del fluido con el cuerpo. El problema se considera 2D y monofásico, y un desarrollo matemático permite una descomposición de los términos disipativos en términos relacionados con la superficie libre y términos relacionados con la enstrofía. El modelo numérico utilizado en la tesis se basa en el método sin malla Smoothed Particle Hydrodynamics (SPH). De manera análoga a lo que se hace a nivel continuo, las propiedades de conservación se estudian en la tesis con el sistema discreto de partículas. Se tratan también las condiciones de contorno de un cuerpo que se mueve en un flujo viscoso, implementadas con el método ghost-fluid. Se ha desarrollado un algoritmo explícito de interacción fluido / cuerpo. Se han documentado algunos casos de modo detallado con el objetivo de comprobar la capacidad del modelo para reproducir correctamente la disipación de energía y el movimiento del cuerpo. En particular se ha investigado la atenuación de una onda estacionaria, comparando la simulación numérica con predicciones teóricas. Se han realizado otras pruebas para monitorizar la disipación de energía para flujos más violentos que implican la fragmentación de la superficie libre. La cantidad de energía disipada con los diferentes términos se ha evaluado en los casos estudiados con el modelo numérico. Se han realizado otras pruebas numéricas para verificar la técnica de modelización de la interacción fluido / cuerpo, concretamente las fuerzas ejercidas por las olas en cuerpos con formas simples, y el equilibrio de un cuerpo flotante con una forma compleja. Una vez que el modelo numérico ha sido validado, se han realizado simulaciones numéricas para obtener una comprensión más completa de la física implicada en casos (casi) realistas sobre los había aspectos que no se conocían suficientemente. En primer lugar se ha estudiado el el flujo alrededor de un cilindro bajo la superficie libre. El estudio se ha realizado con un número de Reynolds moderado, para un rango de inmersiones del cilindro y números de Froude. La solución numérica permite una investigación de los patrones complejos que se producen. La estela del cilindro interactúa con la superficie libre. Se han identificado algunos inestabilidades características. El segundo estudio se ha realizado sobre el problema de sloshing, tanto experimentalmente como numéricamente. El análisis se restringe a aguas poco profundas y con oscilación horizontal, pero se ha estudiado un gran número de condiciones, lo que lleva a una comprensión bastante completa de los sistemas de onda involucradas. La última parte de la tesis trata también sobre un problema de sloshing pero esta vez el tanque está oscilando con rotación y hay acoplamiento con un sistema mecánico. El sistema se llama pendulum-TLD (Tuned Liquid Damper - con líquido amortiguador). Este tipo de sistema se utiliza normalmente para la amortiguación de las estructuras civiles. El análisis se ha realizado analíticamente, numéricamente y experimentalmente utilizando líquidos con viscosidades diferentes, centrándose en características no lineales y mecanismos de disipación. ABSTRA C T The subject of the present thesis is the interaction between a viscous fluid and a solid body in the presence of a free surface. The problem is expressed first theoretically with a particular focus on the energy conservation and the fluid-body interaction. The problem is considered 2D and monophasic, and some mathematical development allows for a decomposition of the energy dissipation into terms related to the Free Surface and others related to the enstrophy. The numerical model used on the thesis is based on Smoothed Particle Hydrodynamics (SPH): a computational method that works by dividing the fluid into particles. Analogously to what is done at continuum level, the conservation properties are studied on the discrete system of particles. Additionally the boundary conditions for a moving body in a viscous flow are treated and discussed using the ghost-fluid method. An explicit algorithm for handling fluid-body coupling is also developed. Following these theoretical developments on the numerical model, some test cases are devised in order to test the ability of the model to correctly reproduce the energy dissipation and the motion of the body. The attenuation of a standing wave is used to compare what is numerically simulated to what is theoretically predicted. Further tests are done in order to monitor the energy dissipation in case of more violent flows involving the fragmentation of the free-surface. The amount of energy dissipated with the different terms is assessed with the numerical model. Other numerical tests are performed in order to test the fluid/body interaction method: forces exerted by waves on simple shapes, and equilibrium of a floating body with a complex shape. Once the numerical model has been validated, numerical tests are performed in order to get a more complete understanding of the physics involved in (almost) realistic cases. First a study is performed on the flow passing a cylinder under the free surface. The study is performed at moderate Reynolds numbers, for various cylinder submergences, and various Froude numbers. The capacity of the numerical solver allows for an investigation of the complex patterns which occur. The wake from the cylinder interacts with the free surface, and some characteristical flow mechanisms are identified. The second study is done on the sloshing problem, both experimentally and numerically. The analysis is restrained to shallow water and horizontal excitation, but a large number of conditions are studied, leading to quite a complete understanding of the wave systems involved. The last part of the thesis still involves a sloshing problem but this time the tank is rolling and there is coupling with a mechanical system. The system is named pendulum-TLD (Tuned Liquid Damper). This kind of system is normally used for damping of civil structures. The analysis is then performed analytically, numerically and experimentally for using liquids with different viscosities, focusing on non-linear features and dissipation mechanisms.

Pulsatile blood flow interaction with arterial walls of aorta : autoregulation and impedance pressure boundary condition and its biomedical applications

Para las decisiones urgentes sobre intervenciones quirúrgicas en el sistema cardiovascular se necesitan simulaciones computacionales con resultados fiables y que consuman un tiempo de cálculo razonable. Durante años los investigadores han trabajado en diversos métodos numéricos de cálculo que resulten atractivos para los cirujanos. Estos métodos, precisos pero costosos desde el punto de vista del coste computacional, crean un desajuste entre la oferta de los ingenieros que realizan las simulaciones y los médicos que operan en el quirófano. Por otra parte, los métodos de cálculo más simplificados reducen el tiempo de cálculo pero pueden proporcionar resultados no realistas. El objetivo de esta tesis es combinar los conceptos de autorregulación e impedancia del sistema circulatorio, la interacción flujo sanguíneo-pared arterial y modelos geométricos idealizados tridimensionales de las arterias pero sin pérdida de realismo, con objeto de proponer una metodología de simulación que proporcione resultados correctos y completos, con tiempos de cálculo moderados. En las simulaciones numéricas, las condiciones de contorno basadas en historias de presión presentan inconvenientes por ser difícil conocerlas con detalle, y porque los resultados son muy sensibles ante pequeñas variaciones de dichas historias. La metodología propuesta se basa en los conceptos de autorregulación, para imponer la demanda de flujo aguas abajo del modelo en el ciclo cardiaco, y la impedancia, para representar el efecto que ejerce el flujo en el resto del sistema circulatorio sobre las arterias modeladas. De este modo las historias de presión en el contorno son resultados del cálculo, que se obtienen de manera iterativa. El método propuesto se aplica en una geometría idealizada del arco aórtico sin patologías y en otra geometría correspondiente a una disección Stanford de tipo A, considerando la interacción del flujo pulsátil con las paredes arteriales. El efecto de los tejidos circundantes también se incorpora en los modelos. También se hacen aplicaciones considerando la interacción en una geometría especifica de un paciente anciano que proviene de una tomografía computarizada. Finalmente se analiza una disección Stanford tipo B con tres modelos que incluyen la fenestración del saco. Clinicians demand fast and reliable numerical results of cardiovascular biomechanic simulations for their urgent pre-surgery decissions. Researchers during many years have work on different numerical methods in order to attract the clinicians' confidence to their colorful contours. Though precise but expensive and time-consuming methodologies create a gap between numerical biomechanics and hospital personnel. On the other hand, simulation simplifications with the aim of reduction in computational time may cause in production of unrealistic outcomes. The main objective of the current investigation is to combine ideas such as autoregulation, impedance, fluid-solid interaction and idealized geometries in order to propose a computationally cheap methodology without excessive or unrealistic simplifications. The pressure boundary conditions are critical and polemic in numerical simulations of cardiovascular system, in which a specific arterial site is of interest and the rest of the netwrok is neglected but represented by a boundary condition. The proposed methodology is a pressure boundary condition which takes advantage of numerical simplicity of application of an imposed pressure boundary condition on outlets, while it includes more sophisticated concepts such as autoregulation and impedance to gain more realistic results. Incorporation of autoregulation and impedance converts the pressure boundary conditions to an active and dynamic boundary conditions, receiving feedback from the results during the numerical calculations and comparing them with the physiological requirements. On the other hand, the impedance boundary condition defines the shapes of the pressure history curves applied at outlets. The applications of the proposed method are seen on idealized geometry of the healthy arotic arch as well as idealized Stanford type A dissection, considering the interaction of the arterial walls with the pulsatile blood flow. The effect of surrounding tissues is incorporated and studied in the models. The simulations continue with FSI analysis of a patient-specific CT scanned geometry of an old individual. Finally, inspiring of the statistic results of mortality rates in Stanford type B dissection, three models of fenestrated dissection sac is studied and discussed. Applying the developed boundary condition, an alternative hypothesis is proposed by the author with respect to the decrease in mortality rates in patients with fenestrations.

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.

MAL, an Integral Element of the Apical Sorting Machinery, Is an Itinerant Protein That Cycles between the Trans-Golgi Network and the Plasma Membrane

The MAL proteolipid is a nonglycosylated integral membrane protein found in glycolipid-enriched membrane microdomains. In polarized epithelial Madin-Darby canine kidney cells, MAL is necessary for normal apical transport and accurate sorting of the influenza virus hemagglutinin. MAL is thus part of the integral machinery for glycolipid-enriched membrane–mediated apical transport. At steady state, MAL is predominantly located in perinuclear vesicles that probably arise from the trans-Golgi network (TGN). To act on membrane traffic and to prevent their accumulation in the target compartment, integral membrane elements of the protein-sorting machinery should be itinerant proteins that cycle between the donor and target compartments. To establish whether MAL is an itinerant protein, we engineered the last extracellular loop of MAL by insertion of sequences containing the FLAG epitope or with sequences containing residues that became O-glycosylated within the cells or that displayed biotinylatable groups. The ectopic expression of these modified MAL proteins allowed us to investigate the surface expression of MAL and its movement through different compartments after internalization with the use of a combination of assays, including surface biotinylation, surface binding of anti-FLAG antibodies, neuraminidase sensitivity, and drug treatments. Immunofluorescence and flow cytometric analyses indicated that, in addition to its Golgi localization, MAL was also expressed on the cell surface, from which it was rapidly internalized. This retrieval implies transport through the endosomal pathway and requires endosomal acidification, because it can be inhibited by drugs such as chloroquine, monensin, and NH4Cl. Resialylation experiments of surface MAL treated with neuraminidase indicated that ∼30% of the internalized MAL molecules were delivered to the TGN, probably to start a new cycle of cargo transport. Together, these observations suggest that, as predicted for integral membrane members of the late protein transport machinery, MAL is an itinerant protein cycling between the TGN and the plasma membrane.

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.

Denaturant-induced movement of the transition state of protein folding revealed by high-pressure stopped-flow measurements

The small all-β protein tendamistat folds and unfolds with two-state kinetics. We determined the volume changes associated with the folding process by performing kinetic and equilibrium measurements at variable pressure between 0.1 and 100 MPa (1 to 1,000 bar). GdmCl-induced equilibrium unfolding transitions reveal that the volume of the native state is increased by 41.4 ± 2.0 cm3/mol relative to the unfolded state. This value is virtually independent of denaturant concentration. The use of a high-pressure stopped-flow instrument enabled us to measure the activation volumes for the refolding (ΔVf0‡) and unfolding reaction (ΔVu0‡) over a broad range of GdmCl concentrations. The volume of the transition state is 60% native-like (ΔVf0‡ = 25.0 ± 1.2 cm3/mol) in the absence of denaturant, indicating partial solvent accessibility of the core residues. The volume of the transition state increases linearly with denaturant concentration and exceeds the volume of the native state above 6 M GdmCl. This result argues for a largely desolvated transition state with packing deficiencies at high denaturant concentrations and shows that the structure of the transition state depends strongly on the experimental conditions.

Anterograde flow of cargo across the Golgi stack potentially mediated via bidirectional “percolating” COPI vesicles

How do secretory proteins and other cargo targeted to post-Golgi locations traverse the Golgi stack? We report immunoelectron microscopy experiments establishing that a Golgi-restricted SNARE, GOS 28, is present in the same population of COPI vesicles as anterograde cargo marked by vesicular stomatitis virus glycoprotein, but is excluded from the COPI vesicles containing retrograde-targeted cargo (marked by KDEL receptor). We also report that GOS 28 and its partnering t-SNARE heavy chain, syntaxin 5, reside together in every cisterna of the stack. Taken together, these data raise the possibility that the anterograde cargo-laden COPI vesicles, retained locally by means of tethers, are inherently capable of fusing with neighboring cisternae on either side. If so, quanta of exported proteins would transit the stack in GOS 28–COPI vesicles via a bidirectional random walk, entering at the cis face and leaving at the trans face and percolating up and down the stack in between. Percolating vesicles carrying both post-Golgi cargo and Golgi residents up and down the stack would reconcile disparate observations on Golgi transport in cells and in cell-free systems.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].