Aplicación de un tratamiento biológico a las aguas residuales provenientes de una destilería de alcohol de caña, utilizando un reactor UASB

En este trabajo se llevó a cabo el tratamiento de vinazas mediante dos tecnologías anaerobias. Se dividió en cuatro estudios técnicos. El primero fue el arranque y estabilización del reactor UASB (Upflow Anaerobic Sludge Blanket), en dónde se evaluó la estabilización mediante la eficiencia de remoción de DQO y la granulación del lodo. El segundo estudio evaluó el rendimiento del reactor UASB frente a diferentes Cva. El tercer estudio evaluó el efecto del TRH sobre la eficiencia del reactor UASB, y el cuarto de ellos fue evaluar el rendimiento del RABF (Reactor Anaerobio de Biomasa Fija). El reactor UASB de 2,6 L de capacidad, fue arrancado por lotes, con seis ensayos utilizando vinaza como sustrato. Se obtuvieron eficiencias de remoción en DQO en un rango de 79-91%, en los seis lotes. Se obtuvo formación de gránulos con diámetro (Ø) de 0,85-1,15 mm y un coeficiente de esfericidad (Є) de 0,7-0,77. Se logró la granulación de lodos tras 2 meses de operación. Alcanzada la estabilización del reactor UASB, se siguió una operación en flujo continuo. Las Cva probadas de 1, 2, 4 y 6 gDQO/L.d para el reactor UASB dan una respuesta bastante favorable con respecto al rendimiento del reactor, ya que presento eficiencias de remoción de DQOs del 51 hasta el 76%, eficiencias similares a los reportados por la literatura. En el estudio de TRH se operó con Cva de 6 gDQO/L.d y los TRH fueron de 24, 12 ,5 ,3 y 1 día. El % de eliminación de DQO fue de 51, 60, 57, 60 y 63 % remoción en DQOsoluble, respectivamente. Se alcanzó una producción de biogás máximo de 5.283 ml/d, pero al reducir el TRH se observó una reducción proporcional del volumen total de biogás. El %CH4 contenido en el biogás aumento al disminuir el TRH, reflejando valores de 80 al 92 % de CH4. El RABF con un volumen de 8,2 L, utilizo tubos de plástico corrugado como medio de soporte para las bacterias. Se aplicaron las siguientes Cva; 0,5, 1, 3 y 6 gDQO/L.d. El reactor RABF presento una excelente remoción de la materia orgánica (80% DQOs), una producción de biogás estable, y un contenido en CH4 del biogás muy interesante. Sin embargo, para una Cva superior a 3 gDQO/L.d empezó un comportamiento inesperado de reducción de capacidad. Las condiciones hidrodinámicas del reactor UASB son decisivas para la formación de los gránulos, condición previa para iniciar el flujo continuo. Al operar el reactor UASB en modo continuo, se pudo evaluar las mejores condiciones de operación para este tipo de residuo (vinaza). La Cva de 6 gDQO/L.d para el reactor UASB alimentado con vinaza bruta representa el límite de su capacidad. Sin embargo, al aumentar la Cva se genera una mayor producción de biogás y metano. La eficiencia de remoción de la DQO soluble es independiente del TRH, para una Cva de 6 g DQO/L•d y las condiciones de TRH probadas (24, 12, 5, 3 y 1 días). Los valores de remoción de DQO alcanzados son un poco superior a los valores de biodegradabilidad anaerobia de la vinaza observados de 50 %. De manera general, la reducción del TRH o bien la dilución de la vinaza no presenta un efecto significativo sobre la remoción de la materia orgánica soluble, pero si lo presenta en la remoción de sulfatos reduciendo indirectamente su toxicidad. El soporte termoplástico inoculado en el RABF y alimentado con vinaza bruta, actuó como un filtro, además de obtener buenos resultados en eliminación de DQO, pero dada las dimensiones y la altura del relleno se frena la evacuación del metano. This work was carried out by treatment vinasses with two anaerobic technologies. It was divided into four technical studies. The first was the start up and stabilization Upflow Anaerobic Sludge Blanket (UASB) reactor, where the stability was evaluated by the removal efficiency of COD and sludge granulation. The second study evaluated the performance of the UASB reactor against different OLR. The third study evaluated the effect of HRT on the efficiency of the UASB reactor, and the fourth of which was evaluate the performance Fixed Biomass Anaerobic (FBA) reactor. The UASB reactor of 2,6 L capacity, was started in batch, with six assays using vinasse as substrate. Were obtained removal efficiencies of COD in the range of 79- 91% in the six batches. Forming granules were obtained with a diameter (Ø) of 0,85- 1,15 mm and sphericity coefficient (Є) of 0,7 to 0,77. Sludge granulation was achieved after 2 months of operation. Once stabilization is achieved of the UASB reactor, it was followed by a continuous flow operation. The OLR tested 1, 2, 4 and 6 gCOD/L.d for UASB reactor gives a very favorable response regarding the performance of the reactor, as presented COD5 removal efficiencies of 51 to 76%, similar efficiencies those reported in the literature The HRT study was operated with an OLR of 6 gCOD/L.d and HRT were 24, 12, 5, 3 and 1 day. The removal efficiency was 51, 60, 57, 60 and 63% in soluble COD, respectively. It reached a maximum biogas production of 5.283 ml / d, but by reducing the HRT showed a proportional reduction in the total volume of biogas. The %CH4 content in the biogas increased with decreasing TRH, reflecting values of 80 to 92% of CH4. The FBA reactor with a volume of 8,2 L, used corrugated plastic tubes as carrier for bacteria transportation. The following OLR was applied, 0,5, 1, 3 and 6 gCOD/L.d. The FBA reactor showed an excellent removal of organic matter (80% CODS), a stable biogas production, and CH4 content very interesting. However, for more than 3 gCOD/L.d OLR began with unexpected behavior of capacity reduction. The UASB reactor hydrodynamic conditions are decisive for the formation of the granules, precondition to start the continuous flow. By operating the UASB reactor in continuous mode, it was possible to evaluate the best operating conditions for this type of waste (vinasse). The OLR of 6 gCOD/L.d for the UASB reactor fed with raw vinasse represents the limit of its capacity. However, with increasing OLR creates increased biogas production and methane. The removal efficiency of soluble COD is independent of HRT for OLR of 6 gCOD/L.d and HRT conditions tested (24, 12, 5, 3 and 1 day). COD Removal values achieved are slightly higher than the values of the vinasse anaerobic biodegradability of observed at 50%. Generally, reduction of HRT or vinasse dilution does not present a significant effect on the removal of the soluble organic matter; however if it occurs in the removal of sulfate reducing indirectly its toxicity. The thermoplastic support inoculated in FBA reactor and fed with raw vinasse, acted as a filter, in addition to obtaining good results in COD removal, but given the size and height of the filling slows evacuation of methane.

Advances in global instability computations: from incompressible to hypersonic flow

Esta tesis constituye un gran avance en el conocimiento del estudio y análisis de inestabilidades hidrodinámicas desde un punto de vista físico y teórico, como consecuencia de haber desarrollado innovadoras técnicas para la resolución computacional eficiente y precisa de la parte principal del espectro correspondiente a los problemas de autovalores (EVP) multidimensionales que gobiernan la inestabilidad de flujos con dos o tres direcciones espaciales inhomogéneas, denominados problemas de estabilidad global lineal. En el contexto del trabajo de desarrollo de herramientas computacionales presentado en la tesis, la discretización mediante métodos de diferencias finitas estables de alto orden de los EVP bidimensionales y tridimensionales que se derivan de las ecuaciones de Navier-Stokes linealizadas sobre flujos con dos o tres direcciones espaciales inhomogéneas, ha permitido una aceleración de cuatro órdenes de magnitud en su resolución. Esta mejora de eficiencia numérica se ha conseguido gracias al hecho de que usando estos esquemas de diferencias finitas, técnicas eficientes de resolución de problemas lineales son utilizables, explotando el alto nivel de dispersión o alto número de elementos nulos en las matrices involucradas en los problemas tratados. Como más notable consecuencia cabe destacar que la resolución de EVPs multidimensionales de inestabilidad global, que hasta la fecha necesitaban de superordenadores, se ha podido realizar en ordenadores de sobremesa. Además de la solución de problemas de estabilidad global lineal, el mencionado desarrollo numérico facilitó la extensión de las ecuaciones de estabilidad parabolizadas (PSE) lineales y no lineales para analizar la inestabilidad de flujos que dependen fuertemente en dos direcciones espaciales y suavemente en la tercera con las ecuaciones de estabilidad parabolizadas tridimensionales (PSE-3D). Precisamente la capacidad de extensión del novedoso algoritmo PSE-3D para el estudio de interacciones no lineales de los modos de estabilidad, desarrollado íntegramente en esta tesis, permite la predicción de transición en flujos complejos de gran interés industrial y por lo tanto extiende el concepto clásico de PSE, el cuál ha sido empleado exitosamente durante las pasadas tres décadas en el mismo contexto para problemas de capa límite bidimensional. Típicos ejemplos de flujos incompresibles se han analizado en este trabajo sin la necesidad de recurrir a restrictivas presuposiciones usadas en el pasado. Se han estudiado problemas vorticales como es el caso de un vórtice aislado o sistemas de vórtices simulando la estela de alas, en los que la homogeneidad axial no se impone y así se puede considerar la difusión viscosa del flujo. Además, se ha estudiado el chorro giratorio turbulento, cuya inestabilidad se utiliza para mejorar las características de funcionamiento de combustores. En la tesis se abarcan adicionalmente problemas de flujos compresibles. Se presenta el estudio de inestabilidad de flujos de borde de ataque a diferentes velocidades de vuelo. También se analiza la estela formada por un elemento rugoso aislado en capa límite supersónica e hipersónica, mostrando excelentes comparaciones con resultados obtenidos mediante simulación numérica directa. Finalmente, nuevas inestabilidades se han identificado en el flujo hipersónico a Mach 7 alrededor de un cono elíptico que modela el vehículo de pruebas en vuelo HIFiRE-5. Los resultados comparan favorablemente con experimentos en vuelo, lo que subraya aún más el potencial de las metodologías de análisis de estabilidad desarrolladas en esta tesis. ABSTRACT The present thesis constitutes a step forward in advancing the frontiers of knowledge of fluid flow instability from a physical point of view, as a consequence of having been successful in developing groundbreaking methodologies for the efficient and accurate computation of the leading part of the spectrum pertinent to multi-dimensional eigenvalue problems (EVP) governing instability of flows with two or three inhomogeneous spatial directions. In the context of the numerical work presented in this thesis, the discretization of the spatial operator resulting from linearization of the Navier-Stokes equations around flows with two or three inhomogeneous spatial directions by variable-high-order stable finite-difference methods has permitted a speedup of four orders of magnitude in the solution of the corresponding two- and three-dimensional EVPs. This improvement of numerical performance has been achieved thanks to the high-sparsity level offered by the high-order finite-difference schemes employed for the discretization of the operators. This permitted use of efficient sparse linear algebra techniques without sacrificing accuracy and, consequently, solutions being obtained on typical workstations, as opposed to the previously employed supercomputers. Besides solution of the two- and three-dimensional EVPs of global linear instability, this development paved the way for the extension of the (linear and nonlinear) Parabolized Stability Equations (PSE) to analyze instability of flows which depend in a strongly-coupled inhomogeneous manner on two spatial directions and weakly on the third. Precisely the extensibility of the novel PSE-3D algorithm developed in the framework of the present thesis to study nonlinear flow instability permits transition prediction in flows of industrial interest, thus extending the classic PSE concept which has been successfully employed in the same context to boundary-layer type of flows over the last three decades. Typical examples of incompressible flows, the instability of which was analyzed in the present thesis without the need to resort to the restrictive assumptions used in the past, range from isolated vortices, and systems thereof, in which axial homogeneity is relaxed to consider viscous diffusion, as well as turbulent swirling jets, the instability of which is exploited in order to improve flame-holding properties of combustors. The instability of compressible subsonic and supersonic leading edge flows has been solved, and the wake of an isolated roughness element in a supersonic and hypersonic boundary-layer has also been analyzed with respect to its instability: excellent agreement with direct numerical simulation results has been obtained in all cases. Finally, instability analysis of Mach number 7 ow around an elliptic cone modeling the HIFiRE-5 flight test vehicle has unraveled flow instabilities near the minor-axis centerline, results comparing favorably with flight test predictions.

Flow around a hemisphere-cylinder at high angle of attack and low Reynolds number. Part II: POD and DMD applied to reduced domains

Three-dimensional direct numerical simulations (DNS) have been performed on a finite-size hemispherecylinder model at angle of attack AoA = 20◦ and Reynolds numbers Re = 350 and 1000. Under these conditions, massive separation exists on the nose and lee-side of the cylinder, and at both Reynolds numbers the flow is found to be unsteady. Proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are employed in order to study the primary instability that triggers unsteadiness at Re = 350. The dominant coherent flow structures identified at the lower Reynolds number are also found to exist at Re = 1000; the question is then posed whether the flow oscillations and structures found at the two Reynolds numbers are related. POD and DMD computations are performed using different subdomains of the DNS computational domain. Besides reducing the computational cost of the analyses, this also permits to isolate spatially localized oscillatory structures from other, more energetic structures present in the flow. It is found that POD and DMD are in general sensitive to domain truncation and noneducated choices of the subdomain may lead to inconsistent results. Analyses at Re = 350 show that the primary instability is related to the counter rotating vortex pair conforming the three-dimensional afterbody wake, and characterized by the frequency St ≈ 0.11, in line with results in the literature. At Re = 1000, vortex-shedding is present in the wake with an associated broadband spectrum centered around the same frequency. The horn/leeward vortices at the cylinder lee-side, upstream of the cylinder base, also present finite amplitude oscillations at the higher Reynolds number. The spatial structure of these oscillations, described by the POD modes, is easily differentiated from that of the wake oscillations. Additionally, the frequency spectra associated with the lee-side vortices presents well defined peaks, corresponding to St ≈ 0.11 and its few harmonics, as opposed to the broadband spectrum found at the wake.

Flow around a hemisphere-cylinder at high angle of attack and low Reynolds number. Part I: Experimental and numerical investigation

Three-dimensional Direct Numerical Simulations combined with Particle Image Velocimetry experiments have been performed on a hemisphere-cylinder at Reynolds number 1000 and angle of attack 20◦. At these flow conditions, a pair of vortices, so-called “horn” vortices, are found to be associated with flow separation. In order to understand the highly complex phenomena associated with this fully threedimensional massively separated flow, different structural analysis techniques have been employed: Proper Orthogonal and Dynamic Mode Decompositions, POD and DMD, respectively, as well as criticalpoint theory. A single dominant frequency associated with the von Karman vortex shedding has been identified in both the experimental and the numerical results. POD and DMD modes associated with this frequency were recovered in the analysis. Flow separation was also found to be intrinsically linked to the observed modes. On the other hand, critical-point theory has been applied in order to highlight possible links of the topology patterns over the surface of the body with the computed modes. Critical points and separation lines on the body surface show in detail the presence of different flow patterns in the base flow: a three-dimensional separation bubble and two pairs of unsteady vortices systems, the horn vortices, mentioned before, and the so-called “leeward” vortices. The horn vortices emerge perpendicularly from the body surface at the separation region. On the other hand, the leeward vortices are originated downstream of the separation bubble, as a result of the boundary layer separation. The frequencies associated with these vortical structures have been quantified.