Determinación de entalpías de mezcla de sistemas muy viscosos usando un ITC y un microcalorímetro de flujo de construcción propia

Programa de doctorado: Física, Matemáticas, Geología y Clima. La fecha de publicación es la fecha de lectura

Diagramas de flujo para el diseño de un sistema de control de calidad : proceso de elaboración de vino blanco

Para la implementación de un sistema de control de calidad es necesario tener una clara visión de los procesos involucrados, las responsabilidades, organización, registros de la calidad, acciones preventivas y correctivas. En la actualidad se exige un enfoque integrador, donde se estructuran actividades y recursos físicos. Es por ello que los diagramas de flujo tienen una importancia fundamental en las decisiones de dirección. Entre los distintos diagramas de flujo posibles, en este trabajo se propone utilizar los de procesos que permiten ordenar secuencialmente una serie de actividades y recursos físicos. Se presenta un diagrama de flujo en la elaboración de vino blanco en una bodega experimental. A través de él se puede conocer la lógica del proceso de elaboración, realizar un diagnóstico del proceso, detectar los puntos críticos de control que puedan provocar deterioro en la calidad del producto o del ambiente. También es un punto de partida para armar y poner en funcionamiento un plan de inspección y ensayos, para la documentación necesaria en las auditorías y una técnica para detectar puntos de contaminación ambiental.

Uso de DES y RANS no estacionario para la simulación del flujo alrededor de una carena Wigley - Use of Unsteady RANS and DES to simulate flows around a Wigley Hull)

Computer Fluid Dynamics tools have already become a valuable instrument for Naval Architects during the ship design process, thanks to their accuracy and the available computer power. Unfortunately, the development of RANSE codes, generally used when viscous effects play a major role in the flow, has not reached a mature stage, being the accuracy of the turbulence models and the free surface representation the most important sources of uncertainty. Another level of uncertainty is added when the simulations are carried out for unsteady flows, as those generally studied in seakeeping and maneuvering analysis and URANS equations solvers are used. Present work shows the applicability and the benefits derived from the use of new approaches for the turbulence modeling (Detached Eddy Simulation) and the free surface representation (Level Set) on the URANS equations solver CFDSHIP-Iowa. Compared to URANS, DES is expected to predict much broader frequency contents and behave better in flows where boundary layer separation plays a major role. Level Set methods are able to capture very complex free surface geometries, including breaking and overturning waves. The performance of these improvements is tested in set of fairly complex flows, generated by a Wigley hull at pure drift motion, with drift angle ranging from 10 to 60 degrees and at several Froude numbers to study the impact of its variation. Quantitative verification and validation are performed with the obtained results to guarantee their accuracy. The results show the capability of the CFDSHIP-Iowa code to carry out time-accurate simulations of complex flows of extreme unsteady ship maneuvers. The Level Set method is able to capture very complex geometries of the free surface and the use of DES in unsteady simulations highly improves the results obtained. Vortical structures and instabilities as a function of the drift angle and Fr are qualitatively identified. Overall analysis of the flow pattern shows a strong correlation between the vortical structures and free surface wave pattern. Karman-like vortex shedding is identified and the scaled St agrees well with the universal St value. Tip vortices are identified and the associated helical instabilities are analyzed. St using the hull length decreases with the increase of the distance along the vortex core (x), which is similar to results from other simulations. However, St scaled using distance along the vortex cores shows strong oscillations compared to almost constants for those previous simulations. The difference may be caused by the effect of the free-surface, grid resolution, and interaction between the tip vortex and other vortical structures, which needs further investigations. This study is exploratory in the sense that finer grids are desirable and experimental data is lacking for large α, especially for the local flow. More recently, high performance computational capability of CFDSHIP-Iowa V4 has been improved such that large scale computations are possible. DES for DTMB 5415 with bilge keels at α = 20º were conducted using three grids with 10M, 48M and 250M points. DES analysis for flows around KVLCC2 at α = 30º is analyzed using a 13M grid and compared with the results of DES on the 1.6M grid by. Both studies are consistent with what was concluded on grid resolution herein since dominant frequencies for shear-layer, Karman-like, horse-shoe and helical instabilities only show marginal variation on grid refinement. The penalties of using coarse grids are smaller frequency amplitude and less resolved TKE. Therefore finer grids should be used to improve V&V for resolving most of the active turbulent scales for all different Fr and α, which hopefully can be compared with additional EFD data for large α when it becomes available.

Flujo de potencia entre estructuras complejas sometidas a excitaciones de alta frecuencia y aplicaciones a la detección de daño estructural

El tema central de investigación en esta Tesis es el estudio del comportamientodinámico de una estructura mediante modelos que describen la distribución deenergía entre los componentes de la misma y la aplicación de estos modelos parala detección de daños incipientes.Los ensayos dinámicos son un modo de extraer información sobre las propiedadesde una estructura. Si tenemos un modelo de la estructura se podría ajustar éstepara que, con determinado grado de precisión, tenga la misma respuesta que elsistema real ensayado. Después de que se produjese un daño en la estructura,la respuesta al mismo ensayo variará en cierta medida; actualizando el modelo alas nuevas condiciones podemos detectar cambios en la configuración del modeloestructural que nos condujeran a la conclusión de que en la estructura se haproducido un daño.De este modo, la detección de un daño incipiente es posible si somos capacesde distinguir una pequeña variación en los parámetros que definen el modelo. Unrégimen muy apropiado para realizar este tipo de detección es a altas frecuencias,ya que la respuesta es muy dependiente de los pequeños detalles geométricos,dado que el tamaño característico en la estructura asociado a la respuesta esdirectamente proporcional a la velocidad de propagación de las ondas acústicas enel sólido, que para una estructura dada es inalterable, e inversamente proporcionala la frecuencia de la excitación. Al mismo tiempo, esta característica de la respuestaa altas frecuencias hace que un modelo de Elementos Finitos no sea aplicable enla práctica, debido al alto coste computacional.Un modelo ampliamente utilizado en el cálculo de la respuesta de estructurasa altas frecuencias en ingeniería es el SEA (Statistical Energy Analysis). El SEAaplica el balance energético a cada componente estructural, relacionando la energíade vibración de estos con la potencia disipada por cada uno de ellos y la potenciatransmitida entre ellos, cuya suma debe ser igual a la potencia inyectada a cadacomponente estructural. Esta relación es lineal y viene caracterizada por los factoresde pérdidas. Las magnitudes que intervienen en la respuesta se consideranpromediadas en la geometría, la frecuencia y el tiempo.Actualizar el modelo SEA a datos de ensayo es, por lo tanto, calcular losfactores de pérdidas que reproduzcan la respuesta obtenida en éste. Esta actualización,si se hace de manera directa, supone la resolución de un problema inversoque tiene la característica de estar mal condicionado. En la Tesis se propone actualizarel modelo SEA, no en término de los factores de pérdidas, sino en términos deparámetros estructurales que tienen sentido físico cuando se trata de la respuestaa altas frecuencias, como son los factores de disipación de cada componente, susdensidades modales y las rigideces características de los elementos de acoplamiento.Los factores de pérdidas se calculan como función de estos parámetros. Estaformulación es desarrollada de manera original en esta Tesis y principalmente sefunda en la hipótesis de alta densidad modal, es decir, que en la respuesta participanun gran número de modos de cada componente estructural.La teoría general del método SEA, establece que el modelo es válido bajounas hipótesis sobre la naturaleza de las excitaciones externas muy restrictivas,como que éstas deben ser de tipo ruido blanco local. Este tipo de carga es difícil dereproducir en condiciones de ensayo. En la Tesis mostramos con casos prácticos queesta restricción se puede relajar y, en particular, los resultados son suficientementebuenos cuando la estructura se somete a una carga armónica en escalón.Bajo estas aproximaciones se desarrolla un algoritmo de optimización por pasosque permite actualizar un modelo SEA a un ensayo transitorio cuando la carga esde tipo armónica en escalón. Este algoritmo actualiza el modelo no solamente parauna banda de frecuencia en particular sino para diversas bandas de frecuencia demanera simultánea, con el objetivo de plantear un problema mejor condicionado.Por último, se define un índice de daño que mide el cambio en la matriz depérdidas cuando se produce un daño estructural en una localización concreta deun componente. Se simula numéricamente la respuesta de una estructura formadapor vigas donde producimos un daño en la sección de una de ellas; como se tratade un cálculo a altas frecuencias, la simulación se hace mediante el Método delos Elementos Espectrales para lo que ha sido necesario desarrollar dentro de laTesis un elemento espectral de tipo viga dañada en una sección determinada. Losresultados obtenidos permiten localizar el componente estructural en que se haproducido el daño y la sección en que éste se encuentra con determinado grado deconfianza.AbstractThe main subject under research in this Thesis is the study of the dynamic behaviourof a structure using models that describe the energy distribution betweenthe components of the structure and the applicability of these models to incipientdamage detection.Dynamic tests are a way to extract information about the properties of astructure. If we have a model of the structure, it can be updated in order toreproduce the same response as in experimental tests, within a certain degree ofaccuracy. After damage occurs, the response will change to some extent; modelupdating to the new test conditions can help to detect changes in the structuralmodel leading to the conclusión that damage has occurred.In this way incipient damage detection is possible if we are able to detect srnallvariations in the model parameters. It turns out that the high frequency regimeis highly relevant for incipient damage detection, because the response is verysensitive to small structural geometric details. The characteristic length associatedwith the response is proportional to the propagation speed of acoustic waves insidethe solid, but inversely proportional to the excitation frequency. At the same time,this fact makes the application of a Finite Element Method impractical due to thehigh computational cost.A widely used model in engineering when dealing with the high frequencyresponse is SEA (Statistical Energy Analysis). SEA applies the energy balance toeach structural component, relating their vibrational energy with the dissipatedpower and the transmitted power between the different components; their summust be equal to the input power to each of them. This relationship is linear andcharacterized by loss factors. The magnitudes considered in the response shouldbe averaged in geometry, frequency and time.SEA model updating to test data is equivalent to calculating the loss factorsthat provide a better fit to the experimental response. This is formulated as an illconditionedinverse problem. In this Thesis a new updating algorithm is proposedfor the study of the high frequency response regime in terms of parameters withphysical meaning such as the internal dissipation factors, modal densities andcharacteristic coupling stiffness. The loss factors are then calculated from theseparameters. The approach is developed entirely in this Thesis and is mainlybased on a high modal density asumption, that is to say, a large number of modescontributes to the response.General SEA theory establishes the validity of the model under the asumptionof very restrictive external excitations. These should behave as a local white noise.This kind of excitation is difficult to reproduce in an experimental environment.In this Thesis we show that in practical cases this assumption can be relaxed, inparticular, results are good enough when the structure is excited with a harmonicstep function.Under these assumptions an optimization algorithm is developed for SEAmodel updating to a transient test when external loads are harmonic step functions.This algorithm considers the response not only in a single frequency band,but also for several of them simultaneously.A damage index is defined that measures the change in the loss factor matrixwhen a damage has occurred at a certain location in the structure. The structuresconsidered in this study are built with damaged beam elements; as we are dealingwith the high frequency response, the numerical simulation is implemented witha Spectral Element Method. It has therefore been necessary to develop a spectralbeam damaged element as well. The reported results show that damage detectionis possible with this algorithm, moreover, damage location is also possible withina certain degree of accuracy.

Desarrollo de una Interfaz de sustitución sensorial para el estudio de la percepción del flujo vibrotáctil (retina táctil)

En el presente trabajo fin de máster se ha concebido, diseñado e utilizado una interfaz háptica, adecuada para ser utilizada como dispositivo de sustitución sensorial, la cual hemos llamado retina táctil. Por cuanto trata de proporcionar información propia del sentido de la vista a través del sentido del tacto. Durante este trabajo, que fue desarrollado en el grupo de robótica y cibernética CAR UPM-CSIC, se ha trabajado en estrecha colaboración con el departamento de la facultad de psicología de la universidad autónoma de Madrid, los cuales han definido las bases de la información de alto orden, como podrían ser, gradientes de intensidades de vibración, mediante las cuales el individuo llega a tener una mejor comprensión del ambiente. El proyecto maneja teorías psicológicas recientes, como las teorías ecológicas y dinámicas que entienden que la percepción se basa en variables informacionales de alto orden. Ejemplos de tales variables son el flujo óptico, gradientes de movimiento, gradientes de intensidades, cambios en gradientes, etc. Sorprendentemente, nuestra percepción visual es mucho más sensible a variables de alto orden que a variables de bajo orden, lo cual descarta que variables de alto orden se infieran o calculen en base a variables de bajo orden. La hipótesis que maneja la teoría ecológica es que las variables de alto orden se detectan como unidades básicas, sin descomponerlas en variables de bajo orden. Imaginemos el caso de un objeto acercándose, intuitivamente pensaríamos que calculamos la distancia y la velocidad del objeto para determinar el momento en el cual este nos impactaría, ¿pero es este realmente el modo en el que actúa nuestro cerebro?, ¿no seremos capaces en determinar directamente el tiempo de contacto como una variable de alto orden presente en el entorno?, por ejemplo, determinar directamente la relación entre el tamaño del objeto y la tasa de crecimiento. También cabe preguntarse si todas estas suposiciones son válidas para estimulaciónes a través de los receptores táctiles en la piel. El dispositivo desarrollado está conformado por 13 módulos cada uno de los cuales maneja 6 tactores o vibradores, para hacer un total de 78 vibradores (ampliables al agregar módulos adicionales), cada uno de los tactores tiene 8mm de diámetro y proporciona información del flujo óptico asociado al entorno que rodea al usuario a través de información táctil, él mismo puede ser utilizado inalámbricamente a pesar de que el procesamiento de los datos se este realizando en una computadora de mesa, lo cual es muy útil al trabajar con ambientes virtuales. También se presenta la integración de la interfaz con el sistema operativo de robots ROS para usarlo en conjunto con las librerías que han sido desarrolladas para el control de la cámara Microsoft Kinect con la cual se puede obtener una matriz de distancias de puntos en el espacio, permitiendo de esta manera utilizar la interfaz en ambientes reales. Finalmente se realizaron experimentos para comprobar hipótesis sobre la variable de percepción del tiempo de contacto además de verificar el correcto funcionamiento del dispositivo de sustitución sensorial tanto en ambientes reales como en ambientes simulados así como comprobar hipótesis sobre la validéz del uso del flujo vibrotáctil para la determinación del tiempo de contacto.

Aplicación de métodos numéricos en la simulación del flujo sanguíneo en el sistema arterial

Revisión de algunas aplicaciones de mecánica computacional de fluidos en la simulación del flujo sanguíneo.

Gestión del flujo de tráfico en la red española de carreteras desde el punto de vista energético y ambiental

Actualmente, son muchos los debates sobre la eficacia de las medidas de eficiencia energética que se han implantado en el sector transporte, y concretamente en el transporte por carretera. En el caso de España, en el último Plan de Acción 2011-2020 (IDAE, 2011b), se consideran prioritarias las medidas propuestas para los sectores difusos, como es el caso del transporte. En éste se pretende conseguir hasta un 33% de ahorro energético con medidas referidas al cambio modal, uso racional de medios, renovación de flotas, etc. Sin embargo, el transporte por carretera sigue creciendo y su principal consecuencia es el aumento de consumo energético y de emisiones de gases con efecto invernadero (GEIs). En este contexto, el proyecto OASIS (programa CENIT) ha definido una metodología de medición de huella energética de la autopista en su ciclo de vida completo. La metodología incluye un modelo de cálculo de la huella energética debido al flujo de tráfico por carretera en el que se obtiene consumo de combustible, consumo de energía y emisiones de CO2. En base a la metodología OASIS, y debido a que la red de alta capacidad de carreteras españolas está infrautilizada, la ponencia se centra en la proposición de estrategias de gestión del flujo de tráfico que pueden ser implementadas para conseguir una reducción de emisiones de CO2 y un aumento de la eficiencia energética y se aplica al caso concreto del paso de la cordillera Cantábrica, puerto de Pajares. Las estrategias incluyen la gestión de los flujos de vehículos ligeros y pesados, así como la gestión de la velocidad. Dichas estrategias se han diseñado mediante el análisis comparado de escenarios alternativos, lo que permite identificar las estrategias eficientes. En el caso de estudio se demuestra que la combinación de diferentes medidas de gestión: flujos de tráfico y velocidad son las más eficaces, con ahorros anuales de consumo energético entorno al 12%; además trazado con pendientes suaves pueden llegar a ahorrar entorno al 113% con respecto a trazados con grandes desniveles. Los resultados permiten comprobar la eficacia de cada estrategia desde el punto de vista energético y ambiental.

Evolución del flujo de residuos generados en obras de edificación residencial

En la actualidad, conocer los residuos de construcción y demolición (RCD) generados en una obra es indispensable para optimizar su gestión. Por tanto, cualquier herramienta que establezca una estimación de los residuos generados debe ser considerada como un instrumento más para plantear soluciones reales en el campo de la sostenibilidad. La presente comunicación da a conocer la evolución del flujo de residuos generados en la construcción de edificios residenciales, identificando la fase de obra que más RCD produce. Para ello, se han analizado varias obras reales construidas en la Comunidad de Madrid. Los resultados obtenidos permiten conocer, con la antelación suficiente, la cantidad de RCD que se va a generar, y por tanto planificar la cantidad y dimensión de los contenedores necesarios en un determinado momento de la obra, así como el espacio preciso en la obra para su correcta gestión.

Simulación computacional hidrdinámica de flujo incompresible a través de un codo

El objetivo del presente estudio es obtener los campos de presión, velocidad y patrones de líneas de corriente, caída de presión a través de un codo, mediante simulación computacional utilizando un programa comercial de CFD. Para ello se emplea como fluido de trabajo agua como líquido saturado a 25ºC y Número de Reynolds (Re) igual a 2x10 5, la cual fluye a través de un codo liso de 90º con relación radio de curvatura/radio de la tubería (R/r) igual a 1. Se utiliza una malla no estructurada de elementos tetraédricos y hexaédricos en las zonas cercanas a las paredes internas y el modelo de turbulencia k-E con leyes de pared escalables.

Microtermometría, direcciones de flujo hidrotermal y metalogenia del yacimiento de tungsteno de Pasto Bueno (Ancash, Perú)

El yacimiento de Pasto Bueno se localiza en el extremo nordeste del Batolito de la Cordillera Blanca, comprende diversas vetas, generalmente subverticales, asociadas al stock cuarzomonzonitico de Consuzo, datado como Terciario Superior, que intruye a las pizarras de la fm. Chicama y cuarcitas de la fm. Chimu. Las principales vetas discurren con direccion N-S cortando al stock, aunque tambien existen sistemas NE-SW asi como NW-SE encajados sobre las rocas metamorficas. La mineralogia de mena reconocida comprende wolframita (hubnerita), tetraedrita/tenantita, esfalerita y galena, en una ganga de cuarzo, fluorita, sericita, pirita y carbonatos, ademas de molibdenita, calcopirita, bornita, arsenopirita, enargita (luzonita), stolzita, scheelita, zinnwaldita, topacio, tungstita y arsenico nativo. Estudios previos han caracterizado Pasto Bueno como un yacimiento con una gran componente de greisen, con una evolucion de las vetas desde un episodio temprano esteril de 400 oC, depositando la mineralizacion economica en torno a los 220-250 oC y con un evento postumo de 175-220 oC rico en CO2. La precipitacion de la wolframita se produjo a partir de un fluido netamente hidrotermal, sin embargo, dicha precipitacion estuvo controlada por el aporte al sistema de aguas externas meteoricas y/o metamorficas. El trabajo llevado a cabo ha consistido en la realizacion de un estudio microtermometrico de las 3 principales estructuras del distrito: Consuelo, Alonso-Fenix y Chabuca, para caracterizar la evolucion del fluido mineralizador desde el stock (veta Consuelo) hacia las rocas metasedimentarias de las fm. Chicama y Chimu (manto Alonso-Fenix y veta Chabuca). Para ello se realizo un muestreo sobre el evento principal de mineralizacion. Dichas muestras se sometieron a un estudio petrografico de lamina gruesa para seleccionar las muestras optimas para el posterior estudio microtermometrico. Previamente a la obtencion de las medidas de temperatura de fusion del hielo (criotermometria) y de homogenizacion del fluido; se realizo un estudio de petrografia de inclusiones fluidas para caracterizarlas y seleccionar las representativas. La interpretacion de los resultados ha permitido confirmar la existencia de un episodio previo de alta temperatura, superior a 282 oC y un evento mineralizador con temperaturas en torno a los 200-240 oC. Sin embargo, las salinidades obtenidas son mucho menores que las previamente publicadas, en torno al 5 % peso eq. NaCl, frente a 11-17 % peso eq. NaCl. Tambien se ha observado un fluido postumo rico en CO2, pero de temperatura superior, en torno a los 270 oC. Los gradientes isotermicos muestran dos focos para dichos fluidos hidrotermales: el primero asociado al stock en la veta Consuelo, y el segundo en la veta Chabuca, asociado a la zona de cabalgamiento de las pizarras de la fm. Chicama sobre las cuarcitas de la fm. Chimu. Este segundo foco puede corresponder con los aportes externos de aguas metamorficas. Para finalizar, se dan una serie de pautas para guiar las futuras exploraciones en el yacimiento. ABSTRACT The Pasto Bueno deposit is located at the northeastern end of the Cordillera Blanca Batholith. It comprises several veins, generally subvertical, associated with the quartz-monzonite stock of Consuzo, dated as Tertiary, which intrudes the Chicama fm. slates and the Chimu fm. quartzites. The main veins trend N-S cutting the stock, although there are also NE-SW and NWSE systems, hosted by the metamorphic rocks. The ore mineralogy comprises wolframite (hubnerite), tetrahedrite/tennantite, sphalerite and galena in a gangue of quartz, fluorite, sericite, pyrite and carbonates, and minor molybdenite, chalcopyrite, bornite, arsenopyrite, enargite (luzonite), stolzite, scheelite, zinnwaldite, topaz, tungstite and native arsenic. Previous studies have characterized Pasto Bueno as a deposit with a large component of greisen, with an evolution of the veins from an early barren 400 oC event , followed by economic mineralization of about 220-250 °C and a late event of 175 -220 oC rich in CO2. Wolframite precipitation occurred from a purely hydrothermal fluid; however, this precipitation was controlled by an external flux of meteoric and/or metamorphic waters. Microthermometric studies of the 3 main structures of the district (Consuelo, Alonso-Fenix and Chabuca veins) have been carried out to depict the evolution of the mineralizing fluid coming from the stock (Consuelo vein) into the metasedimentary rocks of the Chimu and Chicama fm. (Alonso-Fenix and Chabuca veins). The sampling was performed over the main event of mineralization. These samples were subject to a quick plate petrography study in order to select the optimal samples for further microthermometry studies. Before the freezing/heating measures, a fluid inclusion petrography study was done to characterize and select the representative F.I. Interpretation of results has confirmed the existence of a previous episode of higher temperature, over 282 °C, and a mineralizing event with temperatures of about 200-240 °C. However, obtained salinities, about 5 wt% NaCl equivalents, are much lower than those previously reported, about 11-17 wt% NaCl equivalents. A last fluid, rich in CO2, but of higher temperature, about 270 oC, has been characterized. Isothermal gradients show two foci for the hydrothermal fluids: the first one associated to the Consuzo stock as shown in the Consuelo vein, and the second one related to the thrust fault which places the Chicama fm. slates over the Chimu fm. quartzites in the Chabuca vein area. This second focus may correspond to an external input of metamorphic waters. Finally, some guidelines have been given to guide future explorations.

La hidráulica kárstica como aplicación de la hidrodinámica general. Estudio del flujo en un terreno yesífero fisurado

Se pretende plantear y desarrollar algunas de las leyes clásicas de hidrodinámica introduciendo las características que permiten su aplicación al flujo subterráneo en general y a la hidráulica kárstica en particular. Se estudia la consideración de si el movimiento del agua subterránea en el karst se puede definir como flujo a través de conductos individuales o como un medio continuo con huecos saturados en una matriz sólida. Los trabajos de Hagen (1839) y Poiseuille (1846), junto con los de Darcy (1856) configuran la referencia básica para este estudio (in Crespo, 2006). A ellos puede añadirse también las aportaciones de Couette y Chezy (Rouse, 1951) para flujo libre sobre superficies rocosas. Se presentan varios casos en los que las leyes del movimiento laminar unidireccional ofrecen soluciones que pueden ser válidas para la definición de los parámetros fundamentales del flujo del agua en el subsuelo.

Modificaciones técnicas en fermentadores de flujo continuo: efectos sobre la población de protozoos y la fermentación

Los sistemas in vitro de simulación de la fermentación ruminal constituyen herramientas muy útiles para el estudio de la función ruminal. El mantenimiento de poblaciones microbianas representativas de las existentes en el rumen de los animales es uno de los requisitos que deberían cumplir los sistemas in vitro, pero la disminución del número de protozoos en los fermentadores es un hecho bien constatado (Moumen et al., 2009; Martínez et al., 2011). Dado el importante papel que los protozoos desempeñan en la función ruminal, el diseño de modificaciones técnicas que permitan mantener estas poblaciones en fermentadores supondría un importante avance para su utilización en el estudio del ecosistema ruminal. El objetivo del presente trabajo fue analizar dos modificaciones técnicas para mejorar la retención de los protozoos en fermentadores de flujo continuo.

Estudio numérico de flujo turbulento cargado con partículas sólidas a través de canales y tuberías de sección variable

La presente Tesis Doctoral tiene como objetivo el estudio de flujo turbulento cargado con partículas sólidas a través de canales y tuberías de sección constante usando un enfoque Euleriano-Lagrangiano. El campo de flujo de la fase de transporte (aire) se resuelve usando simulación de grandes escalas (LES), implementada en un programa de volúmenes finitos mientras que las ecuaciones gobernantes de la fase dispersa son resueltas por medio de un algoritmo de seguimiento Lagrangiano de partículas que ha sido desarrollado y acoplado al programa que resuelve el flujo. Se estudia de manera sistemática y progresiva la interacción fluido→partícula (one-way coupling), a través de diferentes configuraciones geométricas en coordenadas cartesianas (canales de sección constante y variable) y en coordenadas cilíndricas (tuberías de sección constante y sección variable) abarcando diferentes números de Reynolds y diferentes tamaños de partículas; todos los resultados obtenidos han sido comparados con datos publicados previamente. El estudio de flujo multifásico a través de, tuberías de sección variable, ha sido abordada en otras investigaciones mayoritariamente de forma experimental o mediante simulación usando modelos de turbulencia menos complejos y no mediante LES. El patrón de flujo que se verifica en una tubería con expansión es muy complejo y dicha configuración geométrica se halla en múltiples aplicaciones industriales que involucran el transporte de partículas sólidas, por ello es de gran interés su estudio. Como hecho innovador, en esta tesis no solo se resuelven las estadísticas de velocidad del fluido y las partículas en tuberías con diferentes tamaños de expansión y diferentes regímenes de flujo sino que se caracteriza, usando diversas formulaciones del número de Stokes y el parámetro de arrastre, el ingreso y acumulación de partículas dentro de la zona de recirculación, obteniéndose resultados coincidentes con datos experimentales. ABSTRACT The objective of this Thesis research is to study the turbulent flow laden with solid particles through channels and pipes with using Eulerian-Lagrangian approach. The flow field of the transport phase (air ) is solved using large eddy simulation ( LES ) implemented in a program of finite volume while the governing equations of the dispersed phase are resolved by means of a particle Lagrangian tracking algorithm which was developed and coupled to principal program flow solver . We studied systematically and progressively the fluid interaction → particle ( one- way coupling ) , through different geometric configurations in Cartesian coordinates ( channel with constant and variable section) and in cylindrical coordinates ( pipes with constant section and variable section ) covering different Reynolds numbers and different particle sizes, all results have been compared with previously published data . The study of multiphase flow through, pipes with variable section has been addressed in other investigations predominantly experimentally or by simulation using less complex models and no turbulence by LES. The flow pattern is verified in a pipe expansion is very complex and this geometry is found in many industrial applications involving the transport of solid particles, so it is of great interest to study. As an innovator fact , in this Thesis not only finds fluid velocity statistics and particles with different sizes of pipe expansion and different flow regimes but characterized, using various formulations of the Stokes number and the drag parameter are resolved, the entry and accumulation of particles within the recirculation zone , matching results obtained with experimental data.