Tesis:
Instability Analysis of Vortex Systems in Lifting-Body Wakes
- Autor: TENDERO VENTANAS, Juan Ángel
- Título: Instability Analysis of Vortex Systems in Lifting-Body Wakes
- Fecha: 2016
- Materia: Sin materia definida
- Escuela: E.T.S. DE INGENIEROS AERONAUTICOS
- Departamentos: MECANICA DE FLUIDOS Y PROPULSIÓN AEROESPACIAL
- Acceso electrónico: http://oa.upm.es/43277/
- Director/a 1º: THEOFILIS, Vassilios
- Resumen: The present work deals with vortical flow instabilities, in particular with those that appear in the wake left by aircraft, wind turbines or any other body surrounded by fluid, generating a force perpendicular to the flow, namely lift. The case of the aircraft wake destruction caused by the growth of vortex instabilities is of great importance, as the selection of safe distances between aircraft is in most of the cases very conservative due to the lack of knowledge of the phenomena involved in this process. Therefore, any improvement on their knowledge could help to reduce these distances without concerning the safety of the flight, which is of great importance when large aircraft, like Airbus 380 or Boeing 747, are flying followed by small aeroplanes. Other relevant case where the analysis of vortical wakes is crucial, is that of the wind farms, where the interaction of the wake left from a turbine upwind could degrade the performance of the wind turbine downstream by a great amount or even concern its structural safety. A wide variety of methodologies are discussed and employed in the present thesis to perform the stability analysis of such flows. Two groups can be distinguished among these methods. These belonging to one of them, designed vortex methods, are, in most of the cases, incompressible and inviscid, in contrast with standard stability analysis methods, that are usually viscous. However, the main difference between these two approaches is the fact that vortex methods are grid-less (Lagrangian description), in contrast to standard methods, which need a grid (Eulerian description). Both, vortex and grid methods, could be used to solve the fluid motion as an evolution problem (Initial Value Problems, IVP). This approach is usually denominated DNS (Direct Numerical Simulation), although here this denomination will be used only for the grid IVP. The study of the flow behaviour to instabilities could be done using IVP of the two types, where imposed perturbations or random noise will activate the instabilities, whose growth or decay will be observed. However, in many cases, it is more interesting to perform the stability analysis by means of linearisation of the equations and analysis of the eigenvalues and eigenmodes of the resulting EVP (EigenValue Problem). This second approach could be also used with both of the aforementioned methods, based on the Lagrangian or the Eulerian descriptions of the flow equations, which constitutes the core of the present research work. Some canonical cases are studied first, and the complexity is increased gradually to add more realistic features. However, even the results of the most simple cases have their application to explain phenomena related to the behaviour of vortices in the real world. The case of the counter-rotating vortex pair is a clear example. Despite its simplicity, the results extracted from it are of great relevance, as it is the instability discovered by Crow [35]. In this thesis, the movement of the vortices in a two-dimensional plane is studied: from the simplification of point vortices by a vortex simulation, to the scenario where the vortices are close and there is viscosity in the flow by means of a two-dimensional DNS. Homogeneous instabilities in the axial direction are also studied, first for two point vortices, following the idea of Crow [35], and later for the viscous vortex dipole, using the BiGlobal methodology. Comparisons of these two approaches could be found in Tendero et al. [129, 130, 131]. In addition, axial diffusion could be added and the resulting problem studied employing the three dimensional Parabolised Stability Equations (PSE-3D) analysis concept. All the methodologies considered have their range of validity. Nevertheless, for a true counter-rotating vortex pair with a great amount of axial diffusion, only the PSE-3D analysis proves that can give trustworthy results, as the axial diffusion leads to saturation of the modal growth, in contrast to parallel vortex results. This thesis collects results from several publications: Paredes et al. [113], Tendero et al. [129, 130, 131], He et al. [63], and Suryanarayanan et al. [128], as well as some original work. RESUMEN El presente trabajo trata inestabilidades en flujos con vórtices, en particular, aquellas que aparecen en la estela dejada por aeronaves, aerogeneradores o cualquier otro cuerpo inmerso en un fluido que genere una fuerza perpendicular al mismo, conocida como sustentación. El caso de la destrucción de las estelas de aviones causadas por el crecimiento de inestabilidades es de gran importancia, ya que la selección de distancias de seguridad entre aeronaves es en la mayoría de los casos muy conservadora debido a la falta de conocimiento de los fenómenos involucrados en este proceso. Por lo tanto, cualquier mejora en su conocimiento podría ayudar a reducir esas distancias sin afectar la seguridad del vuelo, lo que es de gran importancia cuando grandes aviones, como el Airbus 380 o el Boeing 747, vuelan seguidos de aviones pequeños. Otro caso relevante en el que el análisis de estelas de vórtices es crucial, es el de los campos de aerogeneradores, donde la interacción de la estela dejada por un aerogenerador aguas arriba podría degradar de una forma importante las actuaciones de una turbina aguas abajo, incluso llegando a afectar a su integridad estructural. Una gran cantidad de metodologías se discuten y emplean en la presente tesis para proceder con el análisis de estabilidad de este tipo de flujos. Pueden diferenciarse dos grupos entre estos métodos. Los pertenecientes a uno de ellos, denominados métodos de vórtices, son, en la mayoría de los casos, incompresibles y no viscosos, al contrario que los métodos típicos para el análisis de estabilidad, que normalmente son viscosos. Sin embargo, la principal distinción entre estas dos aproximaciones es el hecho de que los métodos de vórtices son sin malla (descripción lagrangiana), al contrario que los métodos tradicionales, que sí la necesitan (descripción euleriana). Ambos, métodos de vórtices y con malla, pueden usarse para resolver el movimiento del fluido como un problema de evolución (Problema de Valores Iniciales, IVP en inglés). Este enfoque es normalmente denominado DNS (Simulación Numérica Directa por sus siglas en inglés), aunque aquí se usará esa denominación sólo para los IVP con malla. El estudio del comportamiento del fluido en cuanto a inestabilidades puede hacerse usando IVP de los dos tipos analizados, donde la imposición de perturbaciones o el ruido aleatorio activarían las inestabilidades, cuyo crecimiento o amortiguamiento sería observado. Sin embargo, en muchos casos es más interesante ejecutar el análisis de estabilidad por medio de la linealización de las ecuaciones y el análisis de autovalores y autovectores del resultante EVP (Problema de AutoValores, por sus siglas en inglés). Esta segunda estrategia puede utilizarse con los dos métodos mencionados previamente, tanto con los basados en la descripción lagrangiana como en la euleriana de las ecuaciones del fluido, lo que constituye el núcleo del presente trabajo de investigación. Primero se estudian algunos casos canónicos y luego se va incrementando la complejidad gradualmente para ir añadiendo características más realistas. Sin embargo, incluso los resultados de los casos más simples tienen su aplicación para explicar fenómenos relativos al comportamiento de vórtices en casos reales. El caso del par de vórtices contra-rotatorios es un claro ejemplo. A pesar de su simplicidad, los resultados que se extraen de él son de gran relevancia, como lo es la inestabilidad descubierta por Crow [35]. En esta tesis se estudia el movimiento de los vórtices en un plano bidimensional: desde la simplificación de vórtices puntuales mediante una simulación de vórtices, hasta el escenario donde los vórtices son cercanos y hay viscosidad en el flujo mediante una simulación DNS bidimensional. Se estudian inestabilidades con homogeneidad en la dirección axial, primero para dos vórtices puntuales, siguiendo la idea de Crow [35], y luego para el dipolo de vórtices, usando la metodología BiGobal. Se pueden encontrar comparaciones entre estas dos aproximaciones en Tendero et al. [129, 130, 131]. Además, puede añadirse difusión axial y estudiarse el problema resultante empleando el concepto de análisis por medio de las Ecuaciones de Estabilidad Parabolizadas en tres dimensiones (PSE-3D por sus siglas en inglés). Todas las metodologías consideradas tienen su propio rango de validez. No obstante, para un verdadero par de vórtices contra-rotatorios con una gran cantidad de difusión axial, sólo el análisis por medio de PSE-3D prueba que puede proporcionar resultados de confianza, ya que la difusión axial provoca una saturación del crecimiento del modo, en contra de lo que ocurre en el caso de vórtices paralelos. Esta tesis recopila resultados de varias publicaciones: Paredes et al. [113], Tendero et al. [129, 130, 131], He et al. [63], y Suryanarayanan et al. [128], así como trabajos originales.