Células somáticas y extracción con CO[2] supercrítico perfeccionamiento para la producción de queso de cabra bajo en grasa

Programa de doctorado: Sanidad Animal. Premio Extraordinario de Doctorado en la rama de Ciencias de la Salud 2011-2012

Modelización teórico-experimental de nuevo molde electroconformado para rotomoldeo con calentamiento/enfriamiento directo mediante fluido térmico

Programa de doctorado: Ingeniería Ambiental y Desalinización

Radiactividad natural y Medio Ambiente en las Islas Canarias

XX Jornadas Blascabrerianas, organizadas por la Academia de Ciencias e Ingenierías de Lanzarote y el Departamento de Física de la Universidad de Las Palmas de Gran Canaria

Estudio termofísico de líquidos iónicos basados en piridínio: posibilidades de aplicación específica en agricultura

[ES] Los líquidos iónicos (LI) son materiales prometedores para la industria química. Sin embargo, para una implementación exitosa, es imprescindible tener datos fiables sobre sus propiedades físicas y conocer su comportamiento en solución con otras sustancias. Por ello, en el laboratorio de Termodinámica y Fisicoquímica de Fluidos de la ULPGC hemos analizado algunas propiedades termofísicas de sustancias puras, como densidades y viscosidades de tres isómeros de un líquido iónico a base de piridinio del que se carece de datos básicos, así como otras propiedades en solución con solventes. Respecto a las últimas, se han determinado los volúmenes de mezcla y las entalpías de mezcla con agua y con una serie de alcoholes. Los resultados para los tres isómeros fueron comparados para deducir información sobre el efecto de la posición de un sustituyente. El conocimiento acumulado permitirá una mejor y más rápida selección del correcto LI para una tarea específica. Estas tareas pueden incluir disolución de azúcares u otros hidratos de carbono para su análisis, purificación o separación.

A 2D BEM-FEM approach for time harmonic fluid-structure interaction analysis of thin elastic bodies.

[EN]This paper deals with two-dimensional time harmonic fluid-structure interaction problems when the fluid is at rest, and the elastic bodies have small thicknesses. A BEM-FEM numerical approach is used, where the BEM is applied to the fluid, and the structural FEM is applied to the thin elastic bodies.

Repertorio compartido y cosificación en una comunidad de prácticas universitaria de enseñanza de física

El presente trabajo forma parte de una investigación más amplia, orientada a describir y caracterizar al Taller de Enseñanza de Física (TEF), una experiencia innovadora de enseñanza de Física universitaria, como comunidad de prácticas (CoP) y como espacio de formación docente. En esta presentación, profundizamos en una de las dinámicas que constituyen la negociación de significados en toda CoP: la cosificación. Para el análisis de este proceso en el TEF y la caracterización de su repertorio compartido, focalizamos en la programación de la enseñanza de las unidades pedagógicas que se vinculaban con el tema Energía (vale decir, "Impulso, trabajo y energía", "Fluidos" y "Termodinámica"), desarrollada en dos períodos históricos diferentes (2002 y 2009). Dichas propuestas de enseñanza fueron una manifestación del carácter dinámico y creativo, no rutinario, del trabajo de esta CoP. Este análisis nos ha permitido identificar algunos elementos sustanciales de su repertorio compartido: a) la valoración de las ideas previas de los estudiantes; b) la humanización de la Física, a través de diversos procedimientos y recursos; y c) la enseñanza de "herramientas metodológicas" para la realización de experiencias, contenidos procedimentales fundamentales para el TEF en el análisis de situaciones físicas y de otros campos de conocimiento

Repertorio compartido y cosificación en una comunidad de prácticas universitaria de enseñanza de física

El presente trabajo forma parte de una investigación más amplia, orientada a describir y caracterizar al Taller de Enseñanza de Física (TEF), una experiencia innovadora de enseñanza de Física universitaria, como comunidad de prácticas (CoP) y como espacio de formación docente. En esta presentación, profundizamos en una de las dinámicas que constituyen la negociación de significados en toda CoP: la cosificación. Para el análisis de este proceso en el TEF y la caracterización de su repertorio compartido, focalizamos en la programación de la enseñanza de las unidades pedagógicas que se vinculaban con el tema Energía (vale decir, "Impulso, trabajo y energía", "Fluidos" y "Termodinámica"), desarrollada en dos períodos históricos diferentes (2002 y 2009). Dichas propuestas de enseñanza fueron una manifestación del carácter dinámico y creativo, no rutinario, del trabajo de esta CoP. Este análisis nos ha permitido identificar algunos elementos sustanciales de su repertorio compartido: a) la valoración de las ideas previas de los estudiantes; b) la humanización de la Física, a través de diversos procedimientos y recursos; y c) la enseñanza de "herramientas metodológicas" para la realización de experiencias, contenidos procedimentales fundamentales para el TEF en el análisis de situaciones físicas y de otros campos de conocimiento

Repertorio compartido y cosificación en una comunidad de prácticas universitaria de enseñanza de física

El presente trabajo forma parte de una investigación más amplia, orientada a describir y caracterizar al Taller de Enseñanza de Física (TEF), una experiencia innovadora de enseñanza de Física universitaria, como comunidad de prácticas (CoP) y como espacio de formación docente. En esta presentación, profundizamos en una de las dinámicas que constituyen la negociación de significados en toda CoP: la cosificación. Para el análisis de este proceso en el TEF y la caracterización de su repertorio compartido, focalizamos en la programación de la enseñanza de las unidades pedagógicas que se vinculaban con el tema Energía (vale decir, "Impulso, trabajo y energía", "Fluidos" y "Termodinámica"), desarrollada en dos períodos históricos diferentes (2002 y 2009). Dichas propuestas de enseñanza fueron una manifestación del carácter dinámico y creativo, no rutinario, del trabajo de esta CoP. Este análisis nos ha permitido identificar algunos elementos sustanciales de su repertorio compartido: a) la valoración de las ideas previas de los estudiantes; b) la humanización de la Física, a través de diversos procedimientos y recursos; y c) la enseñanza de "herramientas metodológicas" para la realización de experiencias, contenidos procedimentales fundamentales para el TEF en el análisis de situaciones físicas y de otros campos de conocimiento

Reduced order models for industrial thermo-fluid problems

A new method is presented to generate reduced order models (ROMs) in Fluid Dynamics problems of industrial interest. The method is based on the expansion of the flow variables in a Proper Orthogonal Decomposition (POD) basis, calculated from a limited number of snapshots, which are obtained via Computational Fluid Dynamics (CFD). Then, the POD-mode amplitudes are calculated as minimizers of a properly defined overall residual of the equations and boundary conditions. The method includes various ingredients that are new in this field. The residual can be calculated using only a limited number of points in the flow field, which can be scattered either all over the whole computational domain or over a smaller projection window. The resulting ROM is both computationally efficient(reconstructed flow fields require, in cases that do not present shock waves, less than 1 % of the time needed to compute a full CFD solution) and flexible(the projection window can avoid regions of large localized CFD errors).Also, for problems related with aerodynamics, POD modes are obtained from a set of snapshots calculated by a CFD method based on the compressible Navier Stokes equations and a turbulence model (which further more includes some unphysical stabilizing terms that are included for purely numerical reasons), but projection onto the POD manifold is made using the inviscid Euler equations, which makes the method independent of the CFD scheme. In addition, shock waves are treated specifically in the POD description, to avoid the need of using a too large number of snapshots. Various definitions of the residual are also discussed, along with the number and distribution of snapshots, the number of retained modes, and the effect of CFD errors. The method is checked and discussed on several test problems that describe (i) heat transfer in the recirculation region downstream of a backwards facing step, (ii) the flow past a two-dimensional airfoil in both the subsonic and transonic regimes, and (iii) the flow past a three-dimensional horizontal tail plane. The method is both efficient and numerically robust in the sense that the computational effort is quite small compared to CFD and results are both reasonably accurate and largely insensitive to the definition of the residual, to CFD errors, and to the CFD method itself, which may contain artificial stabilizing terms. Thus, the method is amenable for practical engineering applications. Resumen Se presenta un nuevo método para generar modelos de orden reducido (ROMs) aplicado a problemas fluidodinámicos de interés industrial. El nuevo método se basa en la expansión de las variables fluidas en una base POD, calculada a partir de un cierto número de snapshots, los cuales se han obtenido gracias a simulaciones numéricas (CFD). A continuación, las amplitudes de los modos POD se calculan minimizando un residual global adecuadamente definido que combina las ecuaciones y las condiciones de contorno. El método incluye varios ingredientes que son nuevos en este campo de estudio. El residual puede calcularse utilizando únicamente un número limitado de puntos del campo fluido. Estos puntos puede encontrarse dispersos a lo largo del dominio computacional completo o sobre una ventana de proyección. El modelo ROM obtenido es tanto computacionalmente eficiente (en aquellos casos que no presentan ondas de choque reconstruir los campos fluidos requiere menos del 1% del tiempo necesario para calcular una solución CFD) como flexible (la ventana de proyección puede escogerse de forma que evite contener regiones con errores en la solución CFD localizados y grandes). Además, en problemas aerodinámicos, los modos POD se obtienen de un conjunto de snapshots calculados utilizando un código CFD basado en la versión compresible de las ecuaciones de Navier Stokes y un modelo de turbulencia (el cual puede incluir algunos términos estabilizadores sin sentido físico que se añaden por razones puramente numéricas), aunque la proyección en la variedad POD se hace utilizando las ecuaciones de Euler, lo que hace al método independiente del esquema utilizado en el código CFD. Además, las ondas de choque se tratan específicamente en la descripción POD para evitar la necesidad de utilizar un número demasiado grande de snapshots. Varias definiciones del residual se discuten, así como el número y distribución de los snapshots,el número de modos retenidos y el efecto de los errores debidos al CFD. El método se comprueba y discute para varios problemas de evaluación que describen (i) la transferencia de calor en la región de recirculación aguas abajo de un escalón, (ii) el flujo alrededor de un perfil bidimensional en regímenes subsónico y transónico y (iii) el flujo alrededor de un estabilizador horizontal tridimensional. El método es tanto eficiente como numéricamente robusto en el sentido de que el esfuerzo computacional es muy pequeño comparado con el requerido por el CFD y los resultados son razonablemente precisos y muy insensibles a la definición del residual, los errores debidos al CFD y al método CFD en sí mismo, el cual puede contener términos estabilizadores artificiales. Por lo tanto, el método puede utilizarse en aplicaciones prácticas de ingeniería.

Aplicación de las inclusions fluidas en la búsqueda de yacimientos de cobre. Ejemplo: el pórfido de cobre Lara. Ayacucho – Perú

El pórfido de cobre Lara se encuentra en el flanco pacífico de la Cordillera Occidental de los Andes del Perú, dentro de la faja metalogenética Cu-Mo del Cretáceo superior. Está alineado con los pórfidos de cobre del Paleoceno-Eoceno más grandes del sur del Perú: Toquepala, Cuajone y Quellaveco. Geográficamente, esta localizado al sureste de la costa del Perú, a 400 km de Lima. Se ubica en la provincia de Laramate del departamento de Ayacucho, a 40 km del poblado de Palpa que se encuentra en la Panamericana Sur. Su elevación es de 1850 msnm. Los objetivos de este estudio fueron: caracterizar las inclusiones de fluidos de los eventos hidrotermales relacionados con la mineralización económica, estudiar la naturaleza de los fluidos hidrotermales que formaron la alteración y mineralización, determinar la temperatura de formación de los fluidos mineralizantes y establecer una relación entre las inclusiones de fluidos con el tonelaje y la edad de mineralización de los pórfidos de cobre.

Fenómenos de resonancia asociados al transporte de mercancías líquidas por ferrocarril

Los vagones de mercancías se unen entre sí mediante acoplamientos UIC estándar. Aunque cada plataforma está equipada con un tensor de enganche en ambos extremos, sólo uno de ellos se utiliza, mientras que el otro descansa apoyado en un gancho. Debido a los efectos de oleaje del líquido contenido en las cisternas, este tensor en reposo puede entrar en resonancia con la frecuencia natural de algunos fluidos, llegando a oscilar de forma apreciable. En ocasiones llega a saltar fuera de su alojamiento, provocando ciertas molestias durante la circulación, aunque sin entrañar ningún riesgo para la seguridad de la marcha. Con el propósito de simular la influencia de diversos factores sobre el comportamiento dinámico asociado a estos fenómenos de oleaje, se ha desarrollado un modelo de la interacción fluido-vehículo, empleado técnicas de modelado de sistemas multicuerpo y de elementos finitos.

Métodos para la estimación de la temperatura estática de la formación en yacimientos de petróleo y geotérmicos

Las empresas explotadoras de yacimientos de petróleo y geotérmicos diseñan y construyen sus instalaciones optimizándolas según dos características especialmente representativas: su temperatura y su presión. Por ello, tener información precisa sobre ambos aspectos lo antes posible puede permitir ahorros importantes en las inversiones en infraestructuras y también acortar los plazos de comienzo de explotación; ambos, aspectos esenciales en la rentabilidad del yacimiento. La temperatura estática de una formación es la temperatura existente antes del inicio de las tareas de perforación del yacimiento. Las operaciones de perforación tienen como objetivo perforar, evaluar y terminar un pozo para que produzca de una forma rentable. Durante estas tareas se perturba térmicamente la formación debido a la duración de la perforación, la circulación de fluidos de refrigeración, la diferencia de temperaturas entre la formación y el fluido que perfora, el radio, la difusividad térmica de las formaciones y la tecnología de perforación [7]. La principal herramienta para determinar la temperatura estática de formación son las medidas de temperaturas en los pozos. Estas medidas de temperatura se realizan una vez cesados los trabajos de perforación y circulación de fluidos. El conjunto de medidas de temperatura obtenidas es una serie de valores ascendentes, pero no llega a alcanzar el valor de la temperatura estática de la formación. Para estimar esta temperatura se plantean las ecuaciones diferenciales que describen los fenómenos de transmisión de calor en el yacimiento [14]. Estos métodos se emplean tanto en yacimientos geotérmicos como en yacimientos de petróleo indistintamente [16]. Cada uno de ellos modela el problema de transmisión de calor asociado de una forma distinta, con hipótesis simplificadoras distintas. Se ha comprobado que la aplicación de los distintos métodos a un yacimiento en concreto presenta discrepancias en cuanto a los resultados obtenidos [2,16]. En muchos de los yacimientos no se dispone de información sobre los datos necesarios para aplicar el método de estimación. Esto obliga a adoptar una serie de hipótesis que comprometen la precisión final del método utilizado. Además, puede ocurrir que el método elegido no sea el que mejor reproduce el comportamiento del yacimiento. En el presente trabajo se han analizado los distintos métodos de cálculo. De la gran variedad de métodos existentes [9] se han seleccionado los cuatro más representativos [2,7]: Horner (HM), flujo radial y esférico (SRM), de las dos medidas (TLM) y de fuente de calor cilíndrica (CSM). Estos métodos se han aplicado a una serie de yacimientos de los cuales se conoce la temperatura estática de la formación. De cada yacimiento se disponía de datos tanto de medidas como de las propiedades termo-físicas. Estos datos, en ocasiones, eran insuficientes para la aplicación de todos los métodos, lo cual obligo a adoptar hipótesis sobre los mismos. Esto ha permitido evaluar la precisión de cada método en cada caso. A partir de los resultados obtenidos se han formulado una colección de criterios que permiten seleccionar qué método se adapta mejor para el caso de un yacimiento concreto, de tal manera que el resultado final sea lo más preciso posible. Estos criterios se han fijado en función de las propiedades termo-físicas del yacimiento, que determinan el tipo de yacimiento, y del grado de conocimiento que se dispone del yacimiento objeto de estudio.

Eliminacion de dioxido de carbono en la propulsion de submarinos

Una planta AIP es cualquier sistema propulsivo capaz de posibilitar la navegación de un vehículo submarino bajo la superficie del mar de forma completamente independiente de la atmósfera terrestre. El uso a bordo de submarinos de plantas AIP basadas en la reacción química entre un hidrocarburo y oxígeno (ambos almacenados en el interior del submarino) da lugar a la producción en grandes cantidades de agua y CO2, residuos que necesitan ser eliminados. En concreto, la producción de CO2 en grandes cantidades (y en estado gaseoso) constituye un auténtico problema en un submarino navegando en inmersión, ya que actualmente no resulta viable almacenarlo a bordo, y su eliminación tiene que llevarse a cabo de forma discreta y con un coste energético reducido. Actualmente, hay varias alternativas para eliminar el CO2 producido en la propulsión de un submarino navegando en inmersión, siendo la más ventajosa la disolución de dicha sustancia en agua de mar y su posterior expulsión al exterior del submarino. Esta alternativa consta básicamente de 3 etapas bien definidas: • Etapa 1.- Introducir agua de mar a bordo del submarino, haciendo bajar su presión desde la existente en el exterior hasta la presión a la que se quiere realizar el proceso de disolución. • Etapa 2.- Llevar a cabo el proceso de disolución a presión constante e independiente de la existente en el exterior del submarino. • Etapa 3.- Expulsar fuera del submarino el agua de mar saturada de CO2 haciendo subir su presión desde la correspondiente al proceso de disolución hasta la existente en el exterior. Para ejecutar las etapas 1 y 3 con un coste energético aceptable, resulta necesaria la instalación de un sistema de recuperación de energía, el cual basa su funcionamiento en aprovechar la energía producida en la caída de presión del flujo de agua entrante para elevar la presión del flujo de agua saliente saturada de CO2. El sistema arriba citado puede implementarse de 3 formas alternativas: • Recuperación de doble salto mediante máquinas hidráulicas de desplazamiento positivo. • Recuperación directa mediante cilindros estacionarios dotados de pistones internos. • Recuperación directa mediante cilindros rotativos sin pistones internos. Por otro lado, para ejecutar la etapa 2 de forma silenciosa, y sin ocupar excesivo volumen, resulta necesaria la instalación de un sistema de disolución de CO2 en agua de mar a baja presión, existiendo actualmente 2 principios funcionales viables: • Dispersión de finas burbujas de gas en el seno de una masa de agua. • Difusión directa de CO2 a través de una inter-fase líquido/gas estable sin procesos de dispersión previos. Una vez dicho todo esto, el objetivo de la tesis consiste en llevar a cabo dos estudios comparativos: uno para analizar las ventajas/inconvenientes que presentan las 3 alternativas de recuperación de energía citadas y otro para analizar las ventajas/inconvenientes que presentan los sistemas de disolución de CO2 en agua de mar basados en los 2 principios funcionales mencionados. En ambos estudios se van a tener en cuenta las singularidades propias de una instalación a bordo de submarinos. Para finalizar este resumen, cabe decir que la ejecución de los estudios arriba citados ha exigido el desarrollo de un código software específico (no disponible en la bibliografía) para llevar a cabo la simulación numérica de los distintos sistemas presentados en la tesis. Este código software se ha desarrollado bajo una serie de restricciones importantes, las cuales se listan a continuación: • Ha sido necesario tener en cuenta fluidos de trabajo multi-componente: agua de mar con CO2 disuelto. • El fluido de trabajo se encuentra normalmente en estado líquido, habiendo sido necesario considerar fenómenos de cambio de fase únicamente en etapas incipientes. • La algoritmia se ha diseñado de la forma más simple posible, al objeto de facilitar el subsiguiente proceso de programación y reducir al máximo el tiempo de ejecución en máquina. • La algoritmia arriba citada se ha diseñado para llevar a cabo análisis de tipo comparativo solamente, y no para obtener resultados extremadamente precisos en términos absolutos.

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos

Numerical and Experimental Studies of Sail Aerodynamics

The purpose of this investigation was the determination of the aerodynamic performance of sails and gain knowledge of the phenomena involved in order to improve the aerody¬namic characteristics. In this research, the airflow around different sails in four scenarios was studied. The method to analyze these scenarios was the combination of numerical simulations and experimental tests by taking advantage of the best of each tool. Two different Com¬putational Fluid Dynamic codes were utilized: the ANSYS-CFX and the CD-Adapco’s STAR-CCM+. The experimental tests were conducted in the Atmospheric Boundary Layer Wind Tunnel at the Universidad de Granada (Spain), the Twisted Flow Wind Tunnel at the University of Auckland (New Zealand) and the A9 Wind Tunnel at the Universidad Polit´ecnica de Madrid (Spain). Through this research, it was found the three-dimensional effect of the mast on the aerodynamic performance of an IMS Class boat. The pressure distribution on a Transpac 52 Class mainsail was also determined. Moreover, the aerodynamic perfor¬mance of the 43ft and 60ft Dhow Classes was obtained. Finally, a feasibility study was conducted to use a structural wing in combination with conventional propulsions systems. The main conclusion was that this research clarified gaps on the knowledge of the aerodynamic performance of sails. Moreover, since commercial codes were not specifically designed to study sails, a procedure was developed. On the other hand, innovative experimental techniques were used and applied to model-scale sails. The achievements of this thesis are promising and some of the results are already in use by the industry on a daily basis. El propósito de este estudio era determinar el comportamiento aerodinámico de unas velas y mejorar el conocimiento de los fenómenos que suceden para optimizar las características aerodinámicas de dichas velas. En esta investigación se estudió el flujo de aire alrededor de diferentes velas en cuatro escenarios. El método para analizar estos escenarios fue la combinación de simulaciones numéricas y ensayos experimentales mediante el aprovechamiento de las ventajas de cada herramienta. Se utilizaron dos códigos de dinámica de fluidos computacional: el ANSYS-CFX y el STAR-CCM+ de la empresa CD-Adapco. Los ensayos experimentales se desarrollaron en el túnel de viento de capa límite de la Universidad de Granada (España), el túnel de viento de la Universidad de Auckland (Nueva Zelanda) y en el túnel A9 de la Universidad Politécnica de Madrid (España). Mediante esta investigación, se determinó el efecto tridimensional del mástil en un velero de la clase IMS. También se describió la distribución de presiones sobre una mayor de un Transpac 52. Además, se obtuvo el comportamiento aerodinámico de las clases 43ft y 60ft de los veleros Dhows. Finalmente, se llevó a cabo un estudio de viabilidad de la utilización de un ala estructural en combinación con sistemas de propulsión convencionales. La conclusión principal de esta investigación fue la capacidad de explicar ciertas lagunas en el conocimiento del comportamiento aerodinámico de las velas en diferentes escenarios. Además, dado que los códigos comerciales no están específicamente diseñados para el estudio de velas, se desarrolló un procedimiento a tal efecto. Por otro lado, se han utilizado innovadoras técnicas experimentales y se han aplicado a modelos de velas a escala. Los logros de esta investigación son prometedores y algunos de los resultados obtenidos ya están siendo utilizados por la industria en su día a día.