Mención Internacional en el título de doctorNature offers a vast amount of examples of efficient locomotion. Millions of years of evolution have allowed animals —such as fish, insects and birds—, and even plants —such as winged-seeds or dandelions— to achieve outstanding locomotive skills. Therefore, it is not a surprise that scientists and engineers have tried to replicate the flight and swimming capabilities of the former examples in order to develop efficient aerial and nautical robots. In fact, these efforts have led to the design and development of several successful bioinspired robots. However, their performance is still far below their living counterparts. One of the main reasons is that the understanding of the physics underlying biological locomotion is still limited. This is due to the complexity of the problem under consideration: the locomotion of a body through a fluid medium. This can be considered fluid structure interaction (FSI) problem where the dynamics of the specimens is the result from the hydrodynamic interaction with the surrounding fluid, which in turn is modified by the motion of the specimens. Consequently, the resulting problem is highly nonlinear and complex from a mathematical standpoint. This dissertation attempts to contribute to further understand the fluid structure interactions in bioinspired locomotion problems. To that end, direct numerical simulations of several examples of bioinspired FSI problems are performed. These examples include the auto-rotation of a winged-seed, the flow interactions between the wings of
a dragonfly, and the schooling patterns that emerge between two fish. In the first part of this dissertation, the algorithm which has been developed to perform part of the aforementioned studies is presented. The proposed algorithm allows the study of the FSI of systems of connected rigid bodies —which serve as a model for the actual specimens— immersed in an incompressible fluid. It is built based on a preexisting flow solver, coupled with a robotic algorithm for the computation of the dynamic equations of the bodies. The use of robotic algorithms endows the proposed methodology with a great fiexibility, allowing to simulate a large variety of problems with different geometries and configurations. The second part of the thesis is devoted to the analysis of the aforementioned examples. In this regard, we first consider the flight of a winged-seed. This is a very interesting, yet complex, problem of fluid-dynamic interaction; in which the auto-rotative motion is the result of a subtle equilibrium between the aerodynamic forces and the inertia properties of the winged-seed. In our study, the dynamics and the flow surrounding the auto-rotating seed are characterized in a range of Reynolds numbers, Re. Specifically, we focus on the study of the leading edge vortex (LEV) that is developed on the upper surface of the seed's wing as it auto-rotates. Our findings suggest that, in the explored range Re = [80 — 240], LEV's stability is not driven by vorticity transport along the spanwise direction nor viscous effects, as reported in the literature of rotating wings. Instead, fictitious accelerations (i.e., Coriolis and centrifugal accelerations) are the most suitable candidates to stabilize the LEV over the seed's wing. In the second example, we study the effect of the three-dimensional (3D) interactions in the performance of two tandem wings, resembling those of a dragonfly. To that end, the wings undergo a two-dimensional (2D) optimum kinematics which is a combination of heaving and pitching. We first analyze the effect of wings' aspect ratio, AR, by comparing the 3D and 2D simulations. The results show that 3D vertical interactions are detrimental for the thrust production of the hindwing, but they do not significantly affect the propulsive efficiency of the tandem arrangement. Next, a more realistic flapping kinematics of the 3D is considered and compared to the previous heaving kinematics. We find a decrease in the propulsive efficiency of the flapping wings compared to their heaving counterparts, which has been linked to a non-desired shedding of vorticity on the inboard region of the wings. The last bioinspired example corresponds to the collective motion of two self-propelled three-dimensional bodies. These bodies are idealized as rectangular, flat plates with flexibility along their chordwise direction, and that self-propels thanks to a prescribed vertical motion of their leading edges. We observe that tandem configurations emerge where both plates swim at a constant mean horizontal velocity and with a mean equilibrium horizontal distance. These configurations can be classified, attending to the resulting flow interactions, into compact and regular configurations.
In the former, the performance of the upstream flapper is modified due to the close interaction with the downstream flapper. However, in the regular configurations, the performance of the upstream flapper is similar to that of an isolated flapper. Conversely, the performance of the downstream flapper is affected in both configurations by the interaction with the wake of the upstream flapper. We are able to link the changes in the downstream flapper's performance to its interaction with the vertical jet induced by vortex rings of the upstream flapper's wake. Finally, we propose a model to qualitatively predict the performance of a hypothetical downstream flapper based on the flow field of and isolated flapper, showing good agreement with the actual simulations.La naturaleza ofrece una gran cantidad de ejemplos de locomoción eficiente. Millones de años de evolución han permitido a animales —tales como peces, insectos o pájaros— e incluso plantas —como sainaras o dientes de león— lograr unas habilidades de lomoción excepcionales. Por lo tanto, no es una sorpresa que científicos e ingenieros hayan intentado replicar la capacidades de vuelo y nado de los anteriores ejemplos, con el objetivo de desarrollar robots aéreos y nadadores más eficientes. De hecho, estos esfuerzos han dado lugar al diseño y desarrollo exitoso de varios robots bioinsipirados. Sin embargo, el rendimiento de éstos es todavía muy inferior al de sus referentes biológicos. Una de las principales razones es que la comprensión de la física subyacente de la lomococión de sistemas biológicos es aún limitada. Esto es debido a la complejidad del problema, a saber, el movimiento de un cuerpo a través de un medio fluido. Este se puede considerar como un problema de interacción fluido estructura (FSI) donde la dinámica del espécimen es el resultado de la interacción fluidodinámica con el fluido de alrededor, el cual es a su vez modificado por el movimiento del cuerpo. Consecuentemente, el problema resultante es altamente no lineal y complejo desde un punto de vista matemático. Con esta disertación se pretende contribuir a una mayor comprensión de la interacción fluido estructura en problemas de locomoción bioninspirados. Con tal propósito, se han realizado simulaciones numéricas directas de varios ejemplos bioinspirados de interacción fluido estructura. Estos ejemplos incluyen la autorrotación de una sámara, las interaccionés fluidas entre las alas de una libélula y los patrones de nado que surgen
entre dos peces. Durante la primera parte de esta disertación, se describe el algoritmo que ha sido desarrollado con el propósito de simular alguno de los problemas anteriormente citados. El algoritmo propuesto permite el estudio de la interacción fluido estructura de sistemas de cuerpos rígidos conectados —los cuales sirven como modelo de los especímenes reales— que están sumergidos en un fluido incompresible. Está construido sobre un solver fluido pre-existente, acoplado a un algoritmo robótico que se encarga de calcular las ecuaciones dinámicas de los cuerpos. El uso de algoritmos robóticos proporciona a la metodología propuesta una gran flexibilidad, permitiendo simular una gran variedad de problemas con diversas geometrías y configuraciones. La segunda parte de esta tesis está dedicada al análisis de los ejemplos mencionados anteriormente. En este respecto, consideramos primero el vuelo de una sámara, el cual es un problema muy interesante, aunque complejo, de interacción fluido dinámica
en el cual el movimiento autorrotativo es el resultado de un sutil equilibrio entre las fuerzas aerodinámicas y las propiedades inerciales de la semilla. En nuestro estudio, caracterizamos la dinámica y el flujo alrededor de la semilla autorrotante en un rango de números de Reynolds, Re. En concreto, nos centramos en el estudio del vórtice del borde de ataque (LEV) que se forma en la parte superior del ala de la sámara cuando ésta autorrota. Nuestros hallazgos sugieren que, en el rango explorado de Re = [80 — 240], la estabilidad del LEV no se debe a un transporte de vorticidad a lo largo de la dirección de la envergadura del ala, ni a efectos viscosos —como se ha mencionado en la literatura de alas rotativas—, sino que las aceleraciones ficticias (es decir, las aceleraciones centrífugas y de Coriolis), son las candidatas más probables responsables de la estabilización del LEV. En el segundo ejemplo, se estudia el efecto de las interacciones tridimensionales (3D) en el rendimiento de dos alas en configuración tándem, basadas en las de una libélula. Para ello, se prescribe que el movimiento de las alas sea una combinación de cabeceo y oscilación vertical que es óptimo en 2 dimensiones (2D). Primero analizamos el efecto de la relación de aspecto de las alas, A% comparando los resultados de las simulaciones en 3D y en 2D. Los resultados revelan que las interacciones vorticales en 3D son perjudiciales para la generación de empuje del ala trasera, pero estas interacciones no afectan de forma significativa a la eficiencia propulsiva del conjunto. Posteriormente, se considera un movimiento de batimiento más realista de las alas, y se compara su eficiencia con la obtenida previamente para las alas en movimiento oscilatorio vertical. Se observa una menor eficiencia de las alas en batimiento en comparación con las mismas alas en movimiento oscilatorio vertical. Este deterioro es asociado a un desprendimiento de estructuras vorticales cerca de los bordes marginales de las alas en batimiento. El último ejemplo bioinspirado es el del movimiento colectivo de dos cuerpos tridimiensionales que se auto propulsan. Estos cuerpos se idealizan como placas planas rectangulares, siendo flexibles a lo largo de su cuerda, y que se auto propulsan gracias
a un movimiento vertical impuesto de sus bordes de ataque. Los resultados muestran la aparición de configuraciones tándem donde sendas placas nadan con una velocidad inedia constante y separadas a una distancia de equilibrio que es también constante. Estas configuraciones son clasificadas —atendiendo a las interacciones fluidas— entre compactas y regulares. En las primeras, el rendimiento de la placa que nada aguas arriba (a la que llamaremos líder) se ve afectado por las interacciones cercanas con
el cuerpo que nada aguas abajo (al que denominaremos seguidor). En cambio, en las configuraciones regulares el redimiento del líder es el mismo que el de una placa similar nadando de forma aislada. Por el contrario, el rendimiento del seguidor se ve afectado en ambas configuraciones debido a las interacciones con la estela del líder. Se ha podido relacionar estos cambios en la eficiencia del seguidor con la interacción con el chorro inducido por los anillos vorticales de la estela del líder. Finalmente, hemos
propuesto un modelo que permite predecir, de forma cualitativa, el rendimiento de un seguidor hipotético basándonos en el campo fluido de una placa aislada. El modelo muestra una buena correlación con los datos obtenidos de las simulaciones numéricas.This thesis has been carried out in the Bioengineering and Aerospace Engineering Department at Universidad Carlos III de Madrid. The financial support has been provided by the Spanish Ministry of Economy and Competitiveness through grant DPI2016-76151-C2-2-R (AEI/FEDER, UE).Programa de Doctorado en Mecánica de Fluidos por la Universidad Carlos III de Madrid; la Universidad de Jaén; la Universidad de Zaragoza; la Universidad Nacional de Educación a Distancia; la Universidad Politécnica de Madrid y la Universidad Rovira i VirgiliPresidente: Francisco Javier Huera-Huarte.- Secretario: Javier Rodríguez Rodríguez.- Vocal: Ignazio María Viol
