Direct and adjoint global stability analysis of turbulent transonic flows over a NACA0012 profile

In this work, various turbulent solutions of the two-dimensional (2D) and three-dimensional compressible Reynolds averaged Navier?Stokes equations are analyzed using global stability theory. This analysis is motivated by the onset of flow unsteadiness (Hopf bifurcation) for transonic buffet conditions where moderately high Reynolds numbers and compressible effects must be considered. The buffet phenomenon involves a complex interaction between the separated flow and a shock wave. The efficient numerical methodology presented in this paper predicts the critical parameters, namely, the angle of attack and Mach and Reynolds numbers beyond which the onset of flow unsteadiness appears. The geometry, a NACA0012 profile, and flow parameters selected reproduce situations of practical interest for aeronautical applications. The numerical computation is performed in three steps. First, a steady baseflow solution is obtained; second, the Jacobian matrix for the RANS equations based on a finite volume discretization is computed; and finally, the generalized eigenvalue problem is derived when the baseflow is linearly perturbed. The methodology is validated predicting the 2D Hopf bifurcation for a circular cylinder under laminar flow condition. This benchmark shows good agreement with the previous published computations and experimental data. In the transonic buffet case, the baseflow is computed using the Spalart?Allmaras turbulence model and represents a mean flow where the high frequency content and length scales of the order of the shear-layer thickness have been averaged. The lower frequency content is assumed to be decoupled from the high frequencies, thus allowing a stability analysis to be performed on the low frequency range. In addition, results of the corresponding adjoint problem and the sensitivity map are provided for the first time for the buffet problem. Finally, an extruded three-dimensional geometry of the NACA0012 airfoil, where all velocity components are considered, was also analyzed as a Triglobal stability case, and the outcoming results were compared to the previous 2D limited model, confirming that the buffet onset is well detected.

Analysis of measured and simulated supraglottal acoustic waves

To date, although much attention has been paid to the estimation and modeling of the voice source (ie, the glottal airflow volume velocity), the measurement and characterization of the supraglottal pressure wave have been much less studied. Some previous results have unveiled that the supraglottal pressure wave has some spectral resonances similar to those of the voice pressure wave. This makes the supraglottal wave partially intelligible. Although the explanation for such effect seems to be clearly related to the reflected pressure wave traveling upstream along the vocal tract, the influence that nonlinear source-filter interaction has on it is not as clear. This article provides an insight into this issue by comparing the acoustic analyses of measured and simulated supraglottal and voice waves. Simulations have been performed using a high-dimensional discrete vocal fold model. Results of such comparative analysis indicate that spectral resonances in the supraglottal wave are mainly caused by the regressive pressure wave that travels upstream along the vocal tract and not by source-tract interaction. On the contrary and according to simulation results, source-tract interaction has a role in the loss of intelligibility that happens in the supraglottal wave with respect to the voice wave. This loss of intelligibility mainly corresponds to spectral differences for frequencies above 1500 Hz.

Análisis comparativo y nuevas estimaciones de las tasas de rebase en obras marítimas en grandes profundidades

El rebase se define como el transporte de una cantidad importante de agua sobre la coronación de una estructura. Por tanto, es el fenómeno que, en general, determina la cota de coronación del dique dependiendo de la cantidad aceptable del mismo, a la vista de condicionantes funcionales y estructurales del dique. En general, la cantidad de rebase que puede tolerar un dique de abrigo desde el punto de vista de su integridad estructural es muy superior a la cantidad permisible desde el punto de vista de su funcionalidad. Por otro lado, el diseño de un dique con una probabilidad de rebase demasiado baja o nula conduciría a diseños incompatibles con consideraciones de otro tipo, como son las estéticas o las económicas. Existen distintas formas de estudiar el rebase producido por el oleaje sobre los espaldones de las obras marítimas. Las más habituales son los ensayos en modelo físico y las formulaciones empíricas o semi-empíricas. Las menos habituales son la instrumentación en prototipo, las redes neuronales y los modelos numéricos. Los ensayos en modelo físico son la herramienta más precisa y fiable para el estudio específico de cada caso, debido a la complejidad del proceso de rebase, con multitud de fenómenos físicos y parámetros involucrados. Los modelos físicos permiten conocer el comportamiento hidráulico y estructural del dique, identificando posibles fallos en el proyecto antes de su ejecución, evaluando diversas alternativas y todo esto con el consiguiente ahorro en costes de construcción mediante la aportación de mejoras al diseño inicial de la estructura. Sin embargo, presentan algunos inconvenientes derivados de los márgenes de error asociados a los ”efectos de escala y de modelo”. Las formulaciones empíricas o semi-empíricas presentan el inconveniente de que su uso está limitado por la aplicabilidad de las fórmulas, ya que éstas sólo son válidas para una casuística de condiciones ambientales y tipologías estructurales limitadas al rango de lo reproducido en los ensayos. El objetivo de la presente Tesis Doctoral es el contrate de las formulaciones desarrolladas por diferentes autores en materia de rebase en distintas tipologías de diques de abrigo. Para ello, se ha realizado en primer lugar la recopilación y el análisis de las formulaciones existentes para estimar la tasa de rebase sobre diques en talud y verticales. Posteriormente, se llevó a cabo el contraste de dichas formulaciones con los resultados obtenidos en una serie de ensayos realizados en el Centro de Estudios de Puertos y Costas. Para finalizar, se aplicó a los ensayos de diques en talud seleccionados la herramienta neuronal NN-OVERTOPPING2, desarrollada en el proyecto europeo de rebases CLASH (“Crest Level Assessment of Coastal Structures by Full Scale Monitoring, Neural Network Prediction and Hazard Analysis on Permissible Wave Overtopping”), contrastando de este modo la tasa de rebase obtenida en los ensayos con este otro método basado en la teoría de las redes neuronales. Posteriormente, se analizó la influencia del viento en el rebase. Para ello se han realizado una serie de ensayos en modelo físico a escala reducida, generando oleaje con y sin viento, sobre la sección vertical del Dique de Levante de Málaga. Finalmente, se presenta el análisis crítico del contraste de cada una de las formulaciones aplicadas a los ensayos seleccionados, que conduce a las conclusiones obtenidas en la presente Tesis Doctoral. Overtopping is defined as the volume of water surpassing the crest of a breakwater and reaching the sheltered area. This phenomenon determines the breakwater’s crest level, depending on the volume of water admissible at the rear because of the sheltered area’s functional and structural conditioning factors. The ways to assess overtopping processes range from those deemed to be most traditional, such as semi-empirical or empirical type equations and physical, reduced scale model tests, to others less usual such as the instrumentation of actual breakwaters (prototypes), artificial neural networks and numerical models. Determining overtopping in reduced scale physical model tests is simple but the values obtained are affected to a greater or lesser degree by the effects of a scale model-prototype such that it can only be considered as an approximation to what actually happens. Nevertheless, physical models are considered to be highly useful for estimating damage that may occur in the area sheltered by the breakwater. Therefore, although physical models present certain problems fundamentally deriving from scale effects, they are still the most accurate, reliable tool for the specific study of each case, especially when large sized models are adopted and wind is generated Empirical expressions obtained from laboratory tests have been developed for calculating the overtopping rate and, therefore, the formulas obtained obviously depend not only on environmental conditions – wave height, wave period and water level – but also on the model’s characteristics and are only applicable in a range of validity of the tests performed in each case. The purpose of this Thesis is to make a comparative analysis of methods for calculating overtopping rates developed by different authors for harbour breakwater overtopping. First, existing equations were compiled and analysed in order to estimate the overtopping rate on sloping and vertical breakwaters. These equations were then compared with the results obtained in a number of tests performed in the Centre for Port and Coastal Studies of the CEDEX. In addition, a neural network model developed in the European CLASH Project (“Crest Level Assessment of Coastal Structures by Full Scale Monitoring, Neural Network Prediction and Hazard Analysis on Permissible Wave Overtopping“) was also tested. Finally, the wind effects on overtopping are evaluated using tests performed with and without wind in the physical model of the Levante Breakwater (Málaga).

Global Stability Analysis of Turbulent Transonic Flows on Airfoil geometries

La aparición de inestabilidades en un flujo es un problema importante que puede afectar a algunas aplicaciones aerodinámicas. De hecho existen diferentes tipos de fenómenos no-estacionarios que actualmente son tema de investigación; casos como la separación a altos ángulos de ataque o el buffet transónico son dos ejemplos de cierta relevancia. El análisis de estabilidad global permite identificar la aparición de dichas condiciones inestables, proporcionando información importante sobre la región donde la inestabilidad es dominante y sobre la frecuencia del fenómeno inestable. La metodología empleada es capaz de calcular un flujo base promediado mediante una discretización con volúmenes finitos y posteriormente la solución de un problema de autovalores asociado a la linealización que aparece al perturbar el flujo base. El cálculo numérico se puede dividir en tres pasos: primero se calcula una solución estacionaria para las ecuaciones RANS, luego se extrae la matriz del Jacobiano que representa el problema linealizado y finalmente se deriva y se resuelve el problema de autovalores generalizado mediante el método iterativo de Arnoldi. Como primer caso de validación, la técnica descrita ha sido aplicada a un cilindro circular en condiciones laminares para detectar el principio de las oscilaciones de los vórtices de von Karman, y se han comparado los resultados con experimentos y cálculos anteriores. La parte más importante del estudio se centra en el análisis de flujos compresibles en régimen turbulento. La predicción de la aparición y la progresión de flujo separado a altos ángulos de ataque se han estudiado en el perfil NACA0012 en condiciones tanto subsónicas como supersónicas y en una sección del ala del A310 en condiciones de despegue. Para todas las geometrías analizadas, se ha podido observar que la separación gradual genera la aparición de un modo inestable específico para altos ángulos de ataque siempre mayores que el ángulo asociado al máximo coeficiente de sustentación. Además, se ha estudiado el problema adjunto para obtener información sobre la zona donde una fuerza externa provoca el máximo cambio en el campo fluido. El estudio se ha completado calculando el mapa de sensibilidad estructural y localizando el centro de la inestabilidad. En el presente trabajo de tesis se ha analizado otro importante fenómeno: el buffet transónico. En condiciones transónicas, la interacción entre la onda de choque y la capa límite genera una oscilación de la posición de la onda de choque y, por consiguiente, de las fuerzas aerodinámicas. El conocimiento de las condiciones críticas y su origen puede ayudar a evitar la oscilación causada por estas fuerzas. Las condiciones para las cuales comienza la inestabilidad han sido calculadas y comparadas con trabajos anteriores. Por otra parte, los resultados del correspondiente problema adjunto y el mapa de sensibilidad se han obtenido por primera vez para el buffet, indicando la región del dominio que sera necesario modificar para crear el mayor cambio en las propiedades del campo fluido. Dado el gran consumo de memoria requerido para los casos 3D, se ha realizado un estudio sobre la reducción del domino con la finalidad de reducirlo a la región donde está localizada la inestabilidad. La eficacia de dicha reducción de dominio ha sido evaluada investigando el cambio en la dimensión de la matriz del Jacobiano, no resultando muy eficiente en términos del consumo de memoria. Dado que el buffet es un problema en general tridimensional, el análisis TriGlobal de una geometría 3D podría considerarse el auténtico reto futuro. Como aproximación al problema, un primer estudio se ha realizado empleando una geometría tridimensional extruida del NACA00f2. El cálculo del flujo 3D y, por primera vez en casos tridimensionales compresibles y turbulentos, el análisis de estabilidad TriGlobal, se han llevado a cabo. La comparación de los resultados obtenidos con los resultados del anterior modelo 2D, ha permitido, primero, verificar la exactitud del cálculo 2D realizado anteriormente y también ha proporcionado una estimación del consumo de memoria requerido para el caso 3D. ABSTRACT Flow unsteadiness is an important problem in aerodynamic applications. In fact, there are several types of unsteady phenomena that are still at the cutting edge of research in the field; separation at high angles of attack and transonic buffet are two important examples. Global Stability Analysis can identify the unstable onset conditions, providing important information about the instability location in the domain and the frequency of the unstable phenomenon. The methodology computes a base flow averaged state based on a finite volume discretization and a solution for a generalized eigenvalue problem corresponding to the perturbed linearized equations. The numerical computation is then performed in three steps: first, a steady solution for the RANS equation is computed; second, the Jacobian matrix that represents the linearized problem is obtained; and finally, the generalized eigenvalue problem is derived and solved with an Arnoldi iterative method. As a first validation test, the technique has been applied on a laminar circular cylinder in order to detect the von Karman vortex shedding onset, comparing the results with experiments and with previous calculations. The main part of the study focusses on turbulent and compressible cases. The prediction of the origin and progression of separated flows at high angles of attack has been studied on the NACA0012 airfoil at subsonic and transonic conditions and for the A310 airfoil in take-off configuration. For all the analyzed geometries, it has been found that gradual separation generates the appearance of one specific unstable mode for angles of attack always greater than the ones related to the maximum lift coefficient. In addition, the adjoint problem has been studied to suggest the location of an external force that results in the largest change to the flow field. From the direct and the adjoint analysis the structural sensitivity map has been computed and the core of the instability has been located. The other important phenomenon analyzed in this work is the transonic buffet. In transonic conditions, the interaction between the shock wave and the boundary layer leads to an oscillation of the shock location and, consequently, of the aerodynamic forces. Knowing the critical operational conditions and its origin can be helpful in preventing such fluctuating forces. The instability onset has then been computed and compared with the literature. Moreover, results of the corresponding adjoint problem and a sensitivity map have been provided for the first time for the buffet problem, indicating the region that must be modified to create the biggest change in flow field properties. Because of the large memory consumption required when a 3D case is approached, a domain reduction study has been carried out with the aim of limiting the domain size to the region where the instability is located. The effectiveness of the domain reduction has been evaluated by investigating the change in the Jacobian matrix size, not being very efficient in terms of memory consumption. Since buffet is a three-dimensional problem, TriGlobal stability analysis can be seen as a future challenge. To approximate the problem, a first study has been carried out on an extruded three-dimensional geometry of the NACA0012 airfoil. The 3D flow computation and the TriGlobal stability analysis have been performed for the first time on a compressible and turbulent 3D case. The results have been compared with a 2D model, confirming that the buffet onset evaluated in the 2D case is well detected. Moreover, the computation has given an indication about the memory consumption for a 3D case.

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.

Desarrollo de un Procedimiento Basado en Algoritmos de Optimización para el Dimensionado de Absorbedores Puntuales Aplicados a la Conversión de Energía Undimotriz

La energía transportada por el oleaje a través de los océanos (energía undimotriz) se enmarca dentro de las denominadas energías oceánicas. Su aprovechamiento para generar energía eléctrica (o ser aprovechada de alguna otra forma) es una idea reflejada ya hace más de dos siglos en una patente (1799). Desde entonces, y con especial intensidad desde los años 70, ha venido despertando el interés de instituciones ligadas al I+D+i y empresas del sector energético y tecnológico, debido principalmente a la magnitud del recurso disponible. Actualmente se puede considerar al sector en un estado precomercial, con un amplio rango de dispositivos y tecnologías en diferente grado de desarrollo en los que ninguno destaca sobre los otros (ni ha demostrado su viabilidad económica), y sin que se aprecie una tendencia a converger un único dispositivo (o un número reducido de ellos). El recurso energético que se está tratando de aprovechar, pese a compartir la característica de no-controlabilidad con otras fuentes de energía renovable como la eólica o la solar, presenta una variabilidad adicional. De esta manera, diferentes localizaciones, pese a poder presentar recursos de contenido energético similar, presentan oleajes de características muy diferentes en términos de alturas y periodos de oleaje, y en la dispersión estadística de estos valores. Esta variabilidad en el oleaje hace que cobre especial relevancia la adecuación de los dispositivos de aprovechamiento de energía undimotriz (WEC: Wave Energy Converter) a su localización, de cara a mejorar su viabilidad económica. Parece razonable suponer que, en un futuro, el proceso de diseño de un parque de generación undimotriz implique un rediseño (en base a una tecnología conocida) para cada proyecto de implantación en una nueva localización. El objetivo de esta tesis es plantear un procedimiento de dimensionado de una tecnología de aprovechamiento de la energía undimotriz concreta: los absorbedores puntuales. Dicha metodología de diseño se plantea como un problema de optimización matemático, el cual se resuelve utilizando un algoritmo de optimización bioinspirado: evolución diferencial. Este planteamiento permite automatizar la fase previa de dimensionado implementando la metodología en un código de programación. El proceso de diseño de un WEC es un problema de ingería complejo, por lo que no considera factible el planteamiento de un diseño completo mediante un único procedimiento de optimización matemático. En vez de eso, se platea el proceso de diseño en diferentes etapas, de manera que la metodología desarrollada en esta tesis se utilice para obtener las dimensiones básicas de una solución de referencia de WEC, la cual será utilizada como punto de partida para continuar con las etapas posteriores del proceso de diseño. La metodología de dimensionado previo presentada en esta tesis parte de unas condiciones de contorno de diseño definidas previamente, tales como: localización, características del sistema de generación de energía eléctrica (PTO: Power Take-Off), estrategia de extracción de energía eléctrica y concepto concreto de WEC). Utilizando un algoritmo de evolución diferencial multi-objetivo se obtiene un conjunto de soluciones factibles (de acuerdo con una ciertas restricciones técnicas y dimensionales) y óptimas (de acuerdo con una serie de funciones objetivo de pseudo-coste y pseudo-beneficio). Dicho conjunto de soluciones o dimensiones de WEC es utilizado como caso de referencia en las posteriores etapas de diseño. En el documento de la tesis se presentan dos versiones de dicha metodología con dos modelos diferentes de evaluación de las soluciones candidatas. Por un lado, se presenta un modelo en el dominio de la frecuencia que presenta importantes simplificaciones en cuanto al tratamiento del recurso del oleaje. Este procedimiento presenta una menor carga computacional pero una mayor incertidumbre en los resultados, la cual puede traducirse en trabajo adicional en las etapas posteriores del proceso de diseño. Sin embargo, el uso de esta metodología resulta conveniente para realizar análisis paramétricos previos de las condiciones de contorno, tales como la localización seleccionada. Por otro lado, la segunda metodología propuesta utiliza modelos en el domino estocástico, lo que aumenta la carga computacional, pero permite obtener resultados con menos incertidumbre e información estadística muy útil para el proceso de diseño. Por este motivo, esta metodología es más adecuada para su uso en un proceso de dimensionado completo de un WEC. La metodología desarrollada durante la tesis ha sido utilizada en un proyecto industrial de evaluación energética preliminar de una planta de energía undimotriz. En dicho proceso de evaluación, el método de dimensionado previo fue utilizado en una primera etapa, de cara a obtener un conjunto de soluciones factibles de acuerdo con una serie de restricciones técnicas básicas. La selección y refinamiento de la geometría de la solución geométrica de WEC propuesta fue realizada a posteriori (por otros participantes del proyecto) utilizando un modelo detallado en el dominio del tiempo y un modelo de evaluación económica del dispositivo. El uso de esta metodología puede ayudar a reducir las iteraciones manuales y a mejorar los resultados obtenidos en estas últimas etapas del proyecto. ABSTRACT The energy transported by ocean waves (wave energy) is framed within the so-called oceanic energies. Its use to generate electric energy (or desalinate ocean water, etc.) is an idea expressed first time in a patent two centuries ago (1799). Ever since, but specially since the 1970’s, this energy has become interesting for R&D institutions and companies related with the technological and energetic sectors mainly because of the magnitude of available energy. Nowadays the development of this technology can be considered to be in a pre-commercial stage, with a wide range of devices and technologies developed to different degrees but with none standing out nor economically viable. Nor do these technologies seem ready to converge to a single device (or a reduce number of devices). The energy resource to be exploited shares its non-controllability with other renewable energy sources such as wind and solar. However, wave energy presents an additional short-term variability due to its oscillatory nature. Thus, different locations may show waves with similar energy content but different characteristics such as wave height or wave period. This variability in ocean waves makes it very important that the devices for harnessing wave energy (WEC: Wave Energy Converter) fit closely to the characteristics of their location in order to improve their economic viability. It seems reasonable to assume that, in the future, the process of designing a wave power plant will involve a re-design (based on a well-known technology) for each implementation project in any new location. The objective of this PhD thesis is to propose a dimensioning method for a specific wave-energy-harnessing technology: point absorbers. This design methodology is presented as a mathematical optimization problem solved by using an optimization bio-inspired algorithm: differential evolution. This approach allows automating the preliminary dimensioning stage by implementing the methodology in programmed code. The design process of a WEC is a complex engineering problem, so the complete design is not feasible using a single mathematical optimization procedure. Instead, the design process is proposed in different stages, so the methodology developed in this thesis is used for the basic dimensions of a reference solution of the WEC, which would be used as a starting point for the later stages of the design process. The preliminary dimensioning methodology presented in this thesis starts from some previously defined boundary conditions such as: location, power take-off (PTO) characteristic, strategy of energy extraction and specific WEC technology. Using a differential multi-objective evolutionary algorithm produces a set of feasible solutions (according to certain technical and dimensional constraints) and optimal solutions (according to a set of pseudo-cost and pseudo-benefit objective functions). This set of solutions or WEC dimensions are used as a reference case in subsequent stages of design. In the document of this thesis, two versions of this methodology with two different models of evaluation of candidate solutions are presented. On the one hand, a model in the frequency domain that has significant simplifications in the treatment of the wave resource is presented. This method implies a lower computational load but increased uncertainty in the results, which may lead to additional work in the later stages of the design process. However, use of this methodology is useful in order to perform previous parametric analysis of boundary conditions such as the selected location. On the other hand, the second method uses stochastic models, increasing the computational load, but providing results with smaller uncertainty and very useful statistical information for the design process. Therefore, this method is more suitable to be used in a detail design process for full dimensioning of the WEC. The methodology developed throughout the thesis has been used in an industrial project for preliminary energetic assessment of a wave energy power plant. In this assessment process, the method of previous dimensioning was used in the first stage, in order to obtain a set of feasible solutions according to a set of basic technical constraints. The geometry of the WEC was refined and selected subsequently (by other project participants) using a detailed model in the time domain and a model of economic evaluation of the device. Using this methodology can help to reduce the number of design iterations and to improve the results obtained in the last stages of the project.

Análisis y reducción del riesgo de explosión por incendio en los almacenes de pirotecnia

Un incendio en el interior de un almacén de pirotecnia puede provocar la iniciación más o menos simultánea de los artículos existentes en su interior, con su consecuente onda de choque potenciada por la presión de los propios gases generados en su interior. La finalidad de esta tesis doctoral consiste en analizar los fenómenos que se desarrollan como consecuencia de un incendio interior en los almacenes de pirotecnia y establecer posibles medidas eficaces de prevención y protección con objeto de reducir el riesgo de explosión.. Inicialmente se llevó a cabo un análisis del riesgo de explosión, evaluando las consecuencias en caso de producirse un fuego interior. Los resultados fueron catastróficos. Finalmente, y con objeto de minimizar la posibilidad de producirse un incendio en el interior de un almacén, se evaluaron diferentes sistemas de detección y extinción automática de incendios para almacenes pirotécnicos, realizándose pruebas a escala real. En particular, la presente tesis consta de los siguientes trabajos: 1) Evaluación de la sobrepresión generada por la explosión de artificios pirotécnicos en el interior de almacenes, mediante pruebas de explosión de artículos pirotécnicos, con el fin de evaluar las posibles consecuencias en caso de materializarse una explosión, y disponer de información que permita conocer el grado de resistencia de las estructuras de los almacenes, así como el diseño adecuado de las mismas para soportar los efectos de tales explosiones. 2) Evaluación del efecto de un incendio en el interior de un almacén de pirotecnia mediante pruebas a escala real, con el fin de evaluar la posibilidad de la iniciación simultánea de los artículos pirotécnicos almacenados en su interior, y su consecuente explosión, así como las consecuencias en caso de materializarse la explosión. 3) Análisis de sistemas de detección y extinción automática de incendios para almacenes de pirotecnia, mediante la realización de pruebas a escala real para verificar la eficiencia y eficacia de tales sistemas y agentes para sofocar un incendio pirotécnico de forma tal que se evite la posible explosión asociada. 4) Evaluación de nuevos agentes extintores utilizados en los sistemas de extinción automática de incendios, que actualmente se están utilizando en el sector pirotécnico nacional, mediante la realización de pruebas a escala real. Los resultados obtenidos han permitido al Ministerio de Industria, Energía y Turismo actualizar la reglamentación vigente aplicable (Instrucción Técnica Complementaria número 17 del Reglamento de artículos pirotécnicos y cartuchería, aprobado por el Real Decreto 563/2010, de 7 de mayo). ABSTRACT A fire inside a pyrotechnics magazine can cause the more or less simultaneous initiation of the stored articles, with the wave of shock that is produced promoted by the pressure of the gases generated in its interior. The purpose of this PhD thesis is to analyze the phenomena that develop as a result of a fire inside the pyrotechnics stores and establish possible effective prevention and protection measures to reduce the risk of explosionAn analysis of the risk of explosion and fire consequences in a pyrotechnic magazine, testing to real scale, was initially carried out. The results were catastrophic. Finally, and in order to minimize the possibility of a fire inside a pyrotechnic magazine, different systems of detection and automatic extinction of fires in pyrotechnic magazines were evaluated, testing to real scale. Particularly, this thesis has the following parts: 1) Evaluation of the pressure generated by the pyrotechnic explosion in a magazine, with explosión tests of pyrotechnic articles, in order to evaluate the possible consequences in case of materialize an explosion, and provide information to know the resistance grade of the magazine structures, as well as the right design of those to withstand the effects of such explosions. 2) Evaluation of the effect of a fire inside a pyrotechnic magazine by full-scale tests, in order to evaluate the possibility of simultaneous initiation of the pyrotechnic articles stored inside, and the resulting explosion, as well as the consequences in case of materialize an explosion. 3) Analysis and evaluation of detection and automatic fire extinguishing in pyrotechnic magazines, by full-scale tests in order to verify the efficiency and effectiveness of such systems and agents to quell a pyrotechnic fire and avoid the posible mass explosion. 4) Evaluation of new extinguishing agent used in automatic fire extinguishers, which are currently being used in the national fireworks industry, by full-scale tests. The results obtained have allowed the Ministry of Industry, Energy and Tourism, update the current regulation (Complementary Technical Instruction 17 of the Regulation of Pyrotecnic Articles and Ammunition, approved by the Royal Decree 563/2010, of May 7).

Mies van der Rohe's Illinois Institute of Technology: Analysis and History of a Compositive Development = El IIT de Mies van der Rohe: análisis e historia de un proceso compositivo

The Illinois Institute of Technology (iit) campus, Chicago, by architect Ludwig Mies van der Rohe, is often considered as a transitional work, usually acknowledged as significant for the reorientation of his professional career after he emigrated to the United States. Moreover, its favorable recognition today is somehow indicative of its relevance as a model for urban intervention in the contemporary American city and for contemporary city planning in general, not to mention the profound impact that it had on the cityscape of Chicago. However, today we know it was rather the result of a close collaboration between he and Ludwig Hilberseimer —later on, to be completed with Alfred Caldwell— who merged their personal ideas and expertise in the design for the first time. In addition to this, when one tries to locate the design within its own historical context and evaluate the sources of its approach to it, some contradictions arise. The major impact of the images produced by Mies to promote its realization —widely disseminated in most contemporary architectural periodicals— probably outshined the particular circumstances in which the design was conceived. In fact, it would never be materialized as originally presented, but it was, instead, continuously reworked according to land availability in the site —a circumstance often ignored by subsequent architectural critic, that enthusiastically praised the design even before it was fully completed. One of the main consequences of looking at iit from such a standpoint is that, when historically contextualized, one can appreciate that, due to the urban scale of its implementation process, the design had to face a complex reality very different to that initially planned by the architect, often far from his actual possibilities of intervention. Such approach is in contradiction with the common description of the design as a ‘tabula rasa’ that allegedly would have been formulated on the basis of a full denial of its context. On the contrary, the ever-changing circumstances of the design motivated a necessary re-interpretation of the relation between its executed fragments, in order to keep the original identity of the whole in an ever-changing context. This situation implied a continuous transformation of the design by means of a steady re-composition of its elements: as the number of completed buildings increased in its successive stages, their relation to their site-specific context changed, in a very particular process that these lines try to delineate. Requiring decades to be erected, neither of its authors would ever see the design finished as planned, partially because of the difficulties in acquiring the extension of land that it required. Considering the study of this process as able to provide a valuable gateway to understand the urban discourse that the architects entailed, the aim of these lines is to analyze the problems that the iit campus design had to face. As a starting point, a relationship between practice and theory in the activity of the authors implied in iit campus design has been assumed. Far from being interrupted during World War ii, strong historical evidence can be found to infer that both were developed in parallel. Consequently, the historical sequence of the preserved testimonies has been put into context, as well as their transformation while Mies remained in charge for the campus Master Plan. Notably, when seen from this perspective, some ideas already expressed during his previous European practice were still present during the design process. Particularly, Mies's particular understanding of certain architectural concepts — such as those of ‘order’ and ‘structure’—can be traced paralleling the theories about urban planning from his collaborators, a fact that possibly facilitated the campus successful development. The study of the way these ideas were actually redeveloped and modified in the American urban context, added to the specific process of the implementation of iit campus design, sheds a new light for a critical interpretation of the reasons that made it possible, and of the actual responsibility of Mies's collaborators in its overall development and final completion. RESUMEN El campus del Illinois Institute of Technology (iit) de Chicago, obra del arquitecto Ludwig Mies van der Rohe, es a menudo considerado como una obra de transición que, por lo general, ha venido siendo reconocida como relevante para la reorientación de su carrera profesional posterior a su exilio en los Estados Unidos. El reconocimiento del que goza el proyecto es indicativo, de algún modo, de su importancia como modelo para la intervención urbana en la ciudad norteamericana contemporánea y el planeamiento de la ciudad contemporánea en general, sin olvidar el profundo impacto que ha tenido sobre el paisaje urbano de Chicago. Sin embargo, hoy sabemos que el resultado se benefició de su estrecha colaboración con Ludwig Hilberseimer y se completaría más tarde con la de Alfred Caldwell, quienes unieron sus ideas y experiencia profesional en el proyecto por primera vez. Asimismo, cuando se intenta ubicar el proyecto dentro de su propio contexto histórico y evaluar los criterios de su manera de abordarlo, surgen algunas contradicciones. El considerable impacto de las imágenes producidas por Mies para impulsar su ejecución —ampliamente difundidas en la mayoría de publicaciones de arquitectura de la época— probablemente eclipsó las particulares circunstancias en las que el proyecto fue concebido. De hecho, nunca llegó a materializarse tal y como fue inicialmente presentado. Por contra, fue reelaborado de manera continua, de acuerdo a la disponibilidad de suelo en el emplazamiento; una circunstancia a menudo ignorada por la crítica posterior, que elogió con entusiasmo el proyecto antes siquiera de que fuese terminado. Una de las principales consecuencias de contemplar el iit desde semejante punto de vista es que, una vez contextualizada históricamente su puesta en obra, se puede apreciar que el arquitecto tuvo que enfrentarse a una compleja realidad urbana muy diferente a la inicialmente prevista —probablemente debido a la escala del proyecto— a menudo lejos de sus posibilidades reales de intervención. Este enfoque contradice la descripción habitual del proyecto como una ‘tabula rasa’, que supuestamente se habría formulado sobre la base de una negación completa de su contexto. Por el contrario, las circunstancias cambiantes del proyecto obligaron una necesaria reinterpretación de la relación entre sus frag mentos ejecutados, con el fin de mantener la identidad original del conjunto en un contexto en constante cambio. Esta situación implicó una continua transformación del proyecto por medio de una permanente re-composición de sus elementos: según se incrementaba el número de edificios construidos en las etapas sucesivas de desarrollo del conjunto, variaba su relación con el contexto específico en que se emplazaban, en un proceso muy particular que estas líneas tratan de perfilar. Al necesitar décadas para ser levantado, ninguno de sus autores vería el conjunto terminado según lo planificado, en parte debido a las dificultades para la adquisición de la extensión de suelo que demandaba. Asumiendo que el estudio de este proceso es capaz de proporcionar una valiosa puerta de entrada para elucidar el discurso urbano asumido por los Mies, el objetivo de estas líneas es analizar los problemas a los que el proyecto del campus del iit tuvo que enfrentarse. Como punto de partida, se ha supuesto una relación entre la práctica y la teoría en la actividad de los autores implicados en el proyecto del campus del iit. Lejos de interrumpirse durante la Segunda Guerra Mundial, existen evidencias históricas sólidas para deducir que ambas vertientes se desarrollaron en paralelo. En consecuencia, se ha contextualizado la secuencia histórica de los testimonios conservados, así como su transformación durante el periodo en que Mies estuvo a cargo del Plan General del campus. Significativamente, al ser contempladas bajo esta perspectiva, algunas ideas ya expresadas durante su práctica europea anterior resultan aún presentes durante la redacción del proyecto. En concreto, se puede trazar un paralelismo entre la comprensión particular de Mies de ciertos conceptos arquitectónicos —como los de ‘orden’ y ‘estructura’— y las teorías sobre el urbanismo de sus colaboradores, hecho que posiblemente facilitó el exitoso desarrollo del proyecto. El estudio de la manera en que estas ideas fueron reelaboradas y modificadas en el contexto urbano estadounidense, sumado al proceso específico de su aplicación en el proyecto del campus del iit, arroja una nueva luz para una interpretación crítica tanto de las razones que lo hicieron posible, como del papel real que los colaboradores de Mies tuvieron en su desarrollo y ejecución final.

Design Parameter Analysis of Point Absorber WEC via an Evolutionary-Algorithm-Based Dimensioning Tool

Wave energy conversion has an essential difference from other renewable energies since the dependence between the devices design and the energy resource is stronger. Dimensioning is therefore considered a key stage when a design project of Wave Energy Converters (WEC) is undertaken. Location, WEC concept, Power Take-Off (PTO) type, control strategy and hydrodynamic resonance considerations are some of the critical aspects to take into account to achieve a good performance. The paper proposes an automatic dimensioning methodology to be accomplished at the initial design project stages and the following elements are described to carry out the study: an optimization design algorithm, its objective functions and restrictions, a PTO model, as well as a procedure to evaluate the WEC energy production. After that, a parametric analysis is included considering different combinations of the key parameters previously introduced. A variety of study cases are analysed from the point of view of energy production for different design-parameters and all of them are compared with a reference case. Finally, a discussion is presented based on the results obtained, and some recommendations to face the WEC design stage are given.

Analysis of the behavior of multijunction solar cells under high irradiance Gaussian light profiles showing chromatic aberration with emphasis on tunnel junction performance

In this work, we explain the behavior of multijunction solar cells under non-uniform (spatially and in spectral content) light profiles in general and in particular when Gaussian light profiles cause a photo-generated current density, which exceeds locally the peak current density of the tunnel junction. We have analyzed the implications on the tunnel junction's limitation, that is, in the loss of efficiency due to the appearance of a dip in the I–V curve. For that, we have carried out simulations with our three-dimensional distributed model for multijunction solar cells, which contemplates a full description of the tunnel junction and also takes into account the lateral resistances in the tunnel junction. The main findings are that the current density photo-generated spreads out through the lateral resistances of the device, mainly through the tunnel junction layers and the back contact. Therefore, under non-uniform light profiles these resistances are determinant not only to avoid the tunnel junction's limitation but also for mitigating losses in the fill factor. Therefore, taking into account these lateral resistances could be the key for jointly optimizing the concentrator photovoltaic system (concentrator optics, front grid layout and semiconductor structure)

Field-winding fault detection in synchronous machines with static excitation through frequency response analysis.

Frequency Response Analysis is a well-known technique for the diagnosis of power transformers. Currently, this technique is under research for its application in rotary electrical machines. This paper presents significant results on the application of Frequency Response Analysis to fault detection in field winding of synchronous machines with static excitation. First, the influence of the rotor position on the frequency response is evaluated. Secondly, some relevant test results are shown regarding ground fault and inter-turn fault detection in field windings at standstill condition. The influence of the fault resistance value is also taken into account. This paper also studies the applicability of Frequency Response Analysis in fault detection in field windings while rotating. This represents an important feature because some defects only appear with the machine rated speed. Several laboratory test results show the applicability of this fault detection technique in field windings at full speed with no excitation current.

Transcripts from a single full-length cDNA clone of hepatitis C virus are infectious when directly transfected into the liver of a chimpanzee

We have succeeded in constructing a stable full-length cDNA clone of strain H77 (genotype 1a) of hepatitis C virus (HCV). We devised a cassette vector with fixed 5′ and 3′ termini and constructed multiple full-length cDNA clones of H77 in a single step by cloning of the entire ORF, which was amplified by long reverse transcriptase–PCR, directly into this vector. The infectivity of two complete full-length cDNA clones was tested by the direct intrahepatic injection of a chimpanzee with RNA transcripts. However, we found no evidence for HCV replication. Sequence analysis of these and 16 additional full-length clones revealed that seven clones were defective for polyprotein synthesis, and the remaining nine clones had 6–28 amino acid mutations in the predicted polyprotein compared with the consensus sequence of H77. Next, we constructed a consensus chimera from four of the full-length cDNA clones with just two ligation steps. Injection of RNA transcripts from this consensus clone into the liver of a chimpanzee resulted in viral replication. The sequence of the virus recovered from the chimpanzee was identical to that of the injected RNA transcripts. This stable infectious molecular clone should be an important tool for developing a better understanding of the molecular biology and pathogenesis of HCV.

Invariance between subjects of brain wave representations of language

In three experiments, electric brain waves of 19 subjects were recorded under several different experimental conditions for two purposes. One was to test how well we could recognize which sentence, from a set of 24 or 48 sentences, was being processed in the cortex. The other was to study the invariance of brain waves between subjects. As in our earlier work, the analysis consisted of averaging over trials to create prototypes and test samples, to both of which Fourier transforms were applied, followed by filtering and an inverse transformation to the time domain. A least-squares criterion of fit between prototypes and test samples was used for classification. In all three experiments, averaging over subjects improved the recognition rates. The most significant finding was the following. When brain waves were averaged separately for two nonoverlapping groups of subjects, one for prototypes and the other for test samples, we were able to recognize correctly 90% of the brain waves generated by 48 different sentences about European geography.

Analysis of the Drosophila host defense in domino mutant larvae, which are devoid of hemocytes

We have analyzed the Drosophila immune response in domino mutant larvae, which are devoid of blood cells. The domino mutants have a good larval viability, but they die as prepupae. We show that, on immune challenge, induction of the genes encoding antimicrobial peptides in the fat body is not affected significantly in the mutant larvae, indicating that hemocytes are not essential in this process. The hemocoele of domino larvae contains numerous live microorganisms, the presence of which induces a weak antimicrobial response in the fat body. A full response is observed only after septic injury. We propose that the fat body cells are activated both by the presence of microorganisms and by injury and that injury potentiates the effect of microorganisms. Survival experiments after an immune challenge showed that domino mutants devoid of blood cells maintain a wild-type resistance to septic injury. This resistance was also observed in mutant larvae in which the synthesis of antibacterial peptides is impaired (immune deficiency larvae) and in mutants that are deficient for humoral melanization (Black cells larvae). However, if domino was combined with either the immune deficiency or the Black cell mutation, the resistance to septic injury was reduced severely. These results establish the relevance of the three immune reactions: phagocytosis, synthesis of antibacterial peptides, and melanization. By working in synergy, they provide Drosophila a highly effective defense against injury and/or infection.

Invariance of brain-wave representations of simple visual images and their names

In two experiments, electric brain waves of 14 subjects were recorded under several different conditions to study the invariance of brain-wave representations of simple patches of colors and simple visual shapes and their names, the words blue, circle, etc. As in our earlier work, the analysis consisted of averaging over trials to create prototypes and test samples, to both of which Fourier transforms were applied, followed by filtering and an inverse transformation to the time domain. A least-squares criterion of fit between prototypes and test samples was used for classification. The most significant results were these. By averaging over different subjects, as well as trials, we created prototypes from brain waves evoked by simple visual images and test samples from brain waves evoked by auditory or visual words naming the visual images. We correctly recognized from 60% to 75% of the test-sample brain waves. The general conclusion is that simple shapes such as circles and single-color displays generate brain waves surprisingly similar to those generated by their verbal names. These results, taken together with extensive psychological studies of auditory and visual memory, strongly support the solution proposed for visual shapes, by Bishop Berkeley and David Hume in the 18th century, to the long-standing problem of how the mind represents simple abstract ideas.