In situ characterization of sediment oxygen demand across the spectrum of non-resuspension to full resuspension with varying bed shear stress

Sediment oxygen demand (SOD) can be a significant oxygen sink in various types of water bodies, particularly slow-moving waters with substantial organic sediment accumulation. In most settings where SOD is a concern, the prevailing hydraulic conditions are such that the impact of sediment resuspension on SOD is not considered. However, in the case of Bubbly Creek in Chicago, Illinois, the prevailing slack water conditions are interrupted by infrequent intervals of very high flow rates associated with pumped combined sewer overflow (CSO) during intense hydrologic events. These events can cause resuspension of the highly organic, nutrient-rich bottom sediments, resulting in precipitous drawdown of dissolved oxygen (DO) in the water column. While many past studies have addressed the dependence of SOD on near-bed velocity and bed shear stress prior to the point of sediment resuspension, there has been limited research that has attempted to characterize the complex and dynamic phenomenon of resuspended-sediment oxygen demand. To address this issue, a new in situ experimental apparatus referred to as the U of I Hydrodynamic SOD Sampler was designed to achieve a broad range of velocities and associated bed shear stresses. This allowed SOD to be analyzed across the spectrum of no sediment resuspension associated with low velocity/ bed shear stress through full sediment resuspension associated with high velocity / bed shear stress. The current study split SOD into two separate components: (1) SODNR is the sediment oxygen demand associated with non-resuspension conditions and is a surface sink calculated using traditional methods to yield a value with units (g/m2/day); and (2) SODR is the oxygen demand associated with resuspension conditions, which is a volumetric sink most accurately characterized using non-traditional methods and units that reflect suspension in the water column (mg/L/day). In the case of resuspension, the suspended sediment concentration was analyzed as a function of bed shear stress, and a formulation was developed to characterize SODR as a function of suspended sediment concentration in a form similar to first-order biochemical oxygen demand (BOD) kinetics with Monod DO term. The results obtained are intended to be implemented into a numerical model containing hydrodynamic, sediment transport, and water quality components to yield oxygen demand varying in both space and time for specific flow events. Such implementation will allow evaluation of proposed Bubbly Creek water quality improvement alternatives which take into account the impact of SOD under various flow conditions. Although the findings were based on experiments specific to the conditions in Bubbly Creek, the techniques and formulations developed in this study should be applicable to similar sites.

Measurement and Simulation of Suppression Effects in a Buoyant Turbulent Line Fire

An experimental and numerical study of turbulent fire suppression is presented. For this work, a novel and canonical facility has been developed, featuring a buoyant, turbulent, methane or propane-fueled diffusion flame suppressed via either nitrogen dilution of the oxidizer or application of a fine water mist. Flames are stabilized on a slot burner surrounded by a co-flowing oxidizer, which allows controlled delivery of either suppressant to achieve a range of conditions from complete combustion through partial and total flame quenching. A minimal supply of pure oxygen is optionally applied along the burner to provide a strengthened flame base that resists liftoff extinction and permits the study of substantially weakened turbulent flames. The carefully designed facility features well-characterized inlet and boundary conditions that are especially amenable to numerical simulation. Non-intrusive diagnostics provide detailed measurements of suppression behavior, yielding insight into the governing suppression processes, and aiding the development and validation of advanced suppression models. Diagnostics include oxidizer composition analysis to determine suppression potential, flame imaging to quantify visible flame structure, luminous and radiative emissions measurements to assess sooting propensity and heat losses, and species-based calorimetry to evaluate global heat release and combustion efficiency. The studied flames experience notable suppression effects, including transition in color from bright yellow to dim blue, expansion in flame height and structural intermittency, and reduction in radiative heat emissions. Still, measurements indicate that the combustion efficiency remains close to unity, and only near the extinction limit do the flames experience an abrupt transition from nearly complete combustion to total extinguishment. Measurements are compared with large eddy simulation results obtained using the Fire Dynamics Simulator, an open-source computational fluid dynamics software package. Comparisons of experimental and simulated results are used to evaluate the performance of available models in predicting fire suppression. Simulations in the present configuration highlight the issue of spurious reignition that is permitted by the classical eddy-dissipation concept for modeling turbulent combustion. To address this issue, simple treatments to prevent spurious reignition are developed and implemented. Simulations incorporating these treatments are shown to produce excellent agreement with the experimentally measured data, including the global combustion efficiency.

Slot Film Cooling: A Comprehensive Experimental Characterization

When components of a propulsion system are exposed to elevated flow temperatures there is a risk for catastrophic failure if the components are not properly protected from the thermal loads. Among several strategies, slot film cooling is one of the most commonly used, yet poorly understood active cooling techniques. Tangential injection of a relatively cool fluid layer protects the surface(s) in question, but the turbulent mixing between the hot mainstream and cooler film along with the presence of the wall presents an inherently complex problem where kinematics, thermal transport and multimodal heat transfer are coupled. Furthermore, new propulsion designs rely heavily on CFD analysis to verify their viability. These CFD models require validation of their results, and the current literature does not provide a comprehensive data set for film cooling that meets all the demands for proper validation, namely a comprehensive (kinematic, thermal and boundary condition data) data set obtained over a wide range of conditions. This body of work aims at solving the fundamental issue of validation by providing high quality comprehensive film cooling data (kinematics, thermal mixing, heat transfer). 3 distinct velocity ratios (VR=uc/u∞) are examined corresponding to wall-wake (VR~0.5), min-shear (VR ~ 1.0), and wall-jet (VR~2.0) type flows at injection, while the temperature ratio TR= T∞/Tc is approximately 1.5 for all cases. Turbulence intensities at injection are 2-4% for the mainstream (urms/u∞, vrms/u∞,), and on the order of 8-10% for the coolant (urms/uc, vrms/uc,). A special emphasis is placed on inlet characterization, since inlet data in the literature is often incomplete or is of relatively low quality for CFD development. The data reveals that min-shear injection provides the best performance, followed by the wall-jet. The wall-wake case is comparably poor in performance. The comprehensive data suggests that this relative performance is due to the mixing strength of each case, as well as the location of regions of strong mixing with respect to the wall. Kinematic and thermal data show that strong mixing occurs in the wall-jet away from the wall (y/s>1), while strong mixing in the wall-wake occurs much closer to the wall (y/s<1). Min-shear cases exhibit noticeably weaker mixing confined to about y/s=1. Additionally to these general observations, the experimental data obtained in this work is analyzed to reveal scaling laws for the inlets, near-wall scaling, detecting and characterizing coherent structures in the flow as well as to provide data reduction strategies for comparison to CFD models (RANS and LES).

Simulação numérica do escoamento monofásico em uma bomba centrífuga radial

Centrifugal pumps are vastly used in many industrial applications. Knowledge of how these components behave in several circumstances is crucial for the development of more efficient and, therefore, less expensive pumping installations. The combination of multiple impellers, vaned diffusers and a volute might introduce several complex flow characteristics that largely deviate from regular inviscid pump flow theory. Computational Fluid Dynamics can be very helpful to extract information about which physical phenomena are involved in such flows. In this sense, this work performs a numerical study of the flow in a two-stage centrifugal pump (Imbil ITAP 65-330/2) with a vaned diffuser and a volute. The flow in the pump is modeled using the software Ansys CFX, by means of a multi-block, transient rotor-stator technique, with structured grids for all pump parts. The simulations were performed using water and a mixture of water and glycerin as work fluids. Several viscosities were considered, in a range between 87 and 720 cP. Comparisons between experimental data obtained by Amaral (2007) and numerical head curves showed a good agreement, with an average deviation of 6.8% for water. The behavior of velocity, pressure and turbulence kinetic energy fields was evaluated for several operational conditions. In general, the results obtained by this work achieved the proposed goals and are a significant contribution to the understanding of the flow studied.

Simulação numérica da convecção mista em cavidade preenchida com meio poroso heterogêneo e homogêneo

In this work is presented mixed convection heat transfer inside a lid-driven cavity heated from below and filled with heterogeneous and homogeneous porous medium. In the heterogeneous approach, the solid domain is represented by heat conductive equally spaced blocks; the fluid phase surrounds the blocks being limited by the cavity walls. The homogeneous or pore-continuum approach is characterized by the cavity porosity and permeability. Generalized mass, momentum and energy conservation equations are obtained in dimensionless form to represent both the continuum and the pore-continuum models. The numerical solution is obtained via the finite volume method. QUICK interpolation scheme is set for numerical treatment of the advection terms and SIMPLE algorithm is applied for pressure-velocity coupling. Aiming the laminar regime, the flow parameters are kept in the range of 102≤Re≤103 and 103≤Ra≤106 for both the heterogeneous and homogeneous approaches. In the tested configurations for the continuous model, 9, 16, 36, and 64 blocks are considered for each combination of Re and Ra being the microscopic porosity set as constant φ=0,64 . For the pore-continuum model the Darcy number (Da) is set according to the number of blocks in the heterogeneous cavity and the φ. Numerical results of the comparative study between the microscopic and macroscopic approaches are presented. As a result, average Nusselt number equations for the continuum and the pore continuum models as a function of Ra and Re are obtained.

Quantifying the impacts of urban wind sheltering on the building energy consumption

Common building energy modeling approaches do not account for the influence of surrounding neighborhood on the energy consumption patterns. This thesis develops a framework to quantify the neighborhood impact on a building energy consumption based on the local wind flow. The airflow in the neighborhood is predicted using Computational Fluid Dynamics (CFD) in eight principal wind directions. The developed framework in this study benefits from wind multipliers to adjust the wind velocity encountering the target building. The input weather data transfers the adjusted wind velocities to the building energy model. In a case study, the CFD method is validated by comparing with on-site temperature measurements, and the building energy model is calibrated using utilities data. A comparison between using the adjusted and original weather data shows that the building energy consumption and air system heat gain decreased by 5% and 37%, respectively, while the cooling gain increased by 4% annually.

A mixer design for the pigtail braid

The stirring of a body of viscous fluid using multiple stirring rods is known to be particularly effective when the rods trace out a path corresponding to a nontrivial mathematical braid. The optimal braid is the so-called "pigtail braid", in which three stirring rods execute the usual "over-under" motion associated with braiding plaiting) hair. We show how to achieve this optimal braiding motion straightforwardly: one stirring rod is driven in a figure-of-eight motion, while the other two rods are baffles, which rotate episodically about their common centre. We also explore the extent to which the physical baffles may be replaced by flow structures (such as periodic islands).

Desenvolvimento de uma Metodologia para a Seletividade cinética aplicada ao processo de hidroconversão do petróleo, ocorrendo em reator gás-líquido

Este trabalho tem como objetivo desenvolver uma metodologia de seletividade cinética, para os pseudocomponentes do petróleo em escoamento gás-liquido em colunas de bolhas usando a Fluidodinâmica Computacional (CFD). Uma geometria cilíndrica de 2,5m de altura e 0,162m de diâmetro foi usada tanto na validação fluidodinâmica com base em dados experimentais da literatura, como na análise cinética do reator operando em dois modos distintos em relação a fase líquida: batelada e contínuo. Todos os casos de estudo operam em regime heterogêneo de escoamento, com velocidade superficial do gás igual a 8 cm/s e diâmetro médio de bolhas de 6 mm. O modelo fluidodinâmico validado apresentou boa concordância com os dados experimentais, sendo empregado como base para a implementação do modelo cinético de rede de Krishna e Saxena (1989). A análise da hidroconversão foi realizada a 371ºC, e os resultados mostraram o comportamento esperado para o processo reativo estudado, definindo-se os tempos (batelada) e posições axiais (contínuo) de coleta ideal para os pseudocomponentes leves. Em síntese, ressaltase o uso da ferramenta CFD no entendimento, desenvolvimento e otimização de processos.

Numerical prediction of thermal and dynamic characteristics of a fireinduced ceiling-jet

The aim of this thesis is to test the ability of some correlative models such as Alpert correlations on 1972 and re-examined on 2011, the investigation of Heskestad and Delichatsios in 1978, the correlations produced by Cooper in 1982, to define both dynamic and thermal characteristics of a fire induced ceiling-jet flow. The flow occurs when the fire plume impinges the ceiling and develops in the radial direction of the fire axis. Both temperature and velocity predictions are decisive for sprinklers positioning, fire alarms positions, detectors (heat, smoke) positions and activation times and back-layering predictions. These correlative models will be compared with a 3D numerical simulation software CFAST. For the results comparison of temperature and velocity near the ceiling. These results are also compared with a Computational Fluid Dynamics (CFD) analysis, using ANSYS FLUENT.

Mesh Adaption for Tracking Vortex Structures in OVERTURNS Simulation of the S-76 Rotor in Hover

The constant need to improve helicopter performance requires the optimization of existing and future rotor designs. A crucial indicator of rotor capability is hover performance, which depends on the near-body flow as well as the structure and strength of the tip vortices formed at the trailing edge of the blades. Computational Fluid Dynamics (CFD) solvers must balance computational expenses with preservation of the flow, and to limit computational expenses the mesh is often coarsened in the outer regions of the computational domain. This can lead to degradation of the vortex structures which compose the rotor wake. The current work conducts three-dimensional simulations using OVERTURNS, a three-dimensional structured grid solver that models the flow field using the Reynolds-Averaged Navier-Stokes equations. The S-76 rotor in hover was chosen as the test case for evaluating the OVERTURNS solver, focusing on methods to better preserve the rotor wake. Using the hover condition, various computational domains, spatial schemes, and boundary conditions were tested. Furthermore, a mesh adaption routine was implemented, allowing for the increased refinement of the mesh in areas of turbulent flow without the need to add points to the mesh. The adapted mesh was employed to conduct a sweep of collective pitch angles, comparing the resolved wake and integrated forces to existing computational and experimental results. The integrated thrust values saw very close agreement across all tested pitch angles, while the power was slightly over predicted, resulting in under prediction of the Figure of Merit. Meanwhile, the tip vortices have been preserved for multiple blade passages, indicating an improvement in vortex preservation when compared with previous work. Finally, further results from a single collective pitch case were presented to provide a more complete picture of the solver results.

Error estimation and adaptive mesh refinement for aerodynamic flows

This lecture course covers the theory of so-called duality-based a posteriori error estimation of DG finite element methods. In particular, we formulate consistent and adjoint consistent DG methods for the numerical approximation of both the compressible Euler and Navier-Stokes equations; in the latter case, the viscous terms are discretized based on employing an interior penalty method. By exploiting a duality argument, adjoint-based a posteriori error indicators will be established. Moreover, application of these computable bounds within automatic adaptive finite element algorithms will be developed. Here, a variety of isotropic and anisotropic adaptive strategies, as well as $hp$-mesh refinement will be investigated.

Análisis teórico del contacto plasma-superficie y sus aplicaciones industriales

Actualmente, la física de plasmas constituye una parte importante de la investigación en física que está siendo desarrollada. Su campo de aplicación varía desde el estudio de plasmas interestelares y cósmicos, como las estrellas, las nebulosas, el medio intergaláctico, etc.; hasta aplicaciones más terrenales como la producción de microchips o los dispositivos de iluminación. Resulta particularmente interesante el estudio del contacto de una superficie metálica con un plasma. Siendo la razón que, la dinámica de la interfase formada entre un plasma imperturbado y una superficie metálica, resulta de gran importancia cuando se trata de estudiar problemas como: la implantación iónica en una oblea de silicio, el grabado por medio de plasmas, la carga de una aeronave cuando atraviesa la ionosfera y la diagnosis de plasmas mediante sondas de Langmuir. El uso de las sondas de Langmuir está extendido a través de multitud de aplicaciones tecnológicas e industriales como método de diagnosis de plasmas. Algunas de estas aplicaciones han sido mencionadas justo en el párrafo anterior. Es más, su uso también es muy popular en la investigación en física de plasmas, por ser una de las pocas técnicas de diagnosis que proporciona información local sobre el plasma. El equipamiento donde es habitualmente implementado varía desde plasmas de laboratorio de baja temperatura hasta plasmas de fusión en dispositivos como tokamaks o stellerators. La geometría más popular de este tipo de sondas es cilíndrica, y la principal magnitud que se usa para diagnosticar el plasma es la corriente recogida por la sonda cuando se encuentra polarizada a un cierto potencial. Existe un interes especial en diagnosticar por medio de la medida de la corriente iónica recogida por la sonda, puesto que produce una perturbación muy pequeña del plasma en comparación con el uso de la corriente electrónica. Dada esta popularidad, no es de extrañar que grandes esfuerzos se hayan realizado en la consecución de un modelo teórico que explique el comportamiento de una sonda de Langmuir inmersa en un plasma. Hay que remontarse a la primera mitad del siglo XX para encontrar las primeras teorías que permiten diagnosticar parámetros del plasma mediante la medida de la corriente iónica recogida por la sonda de Langmuir. Desde entonces, las mejoras en estos modelos y el desarrollo de otros nuevos ha sido una constante en la investigación en física de plasmas. No obstante, todavía no está claro como los iones se aproximan a la superficie de la sonda. Las dos principales, a la par que opuestas, aproximaciones al problema que están ampliamente aceptadas son: la radial y la orbital; siendo el problema que ambas predicen diferentes valores para la corriente iónica. Los experimentos han arrojado resultados de acuerdo con ambas teorías, la radial y la orbital; y lo que es más importante, una transición entre ambos ha sido recientemente observada. La mayoría de los logros conseguidos a la hora de comprender como los iones caen desde el plasma hacia la superficie de la sonda, han sido llevados a cabo en el campo de la dinámica de fluidos o la teoría cinética. Por otra parte, este problema puede ser abordado mediante el uso de simulaciones de partículas. La principal ventaja de las simulaciones de partículas sobre los modelos de fluidos o cinéticos es que proporcionan mucha más información sobre los detalles microscópicos del movimiento de las partículas, además es relativamente fácil introducir interacciones complejas entre las partículas. No obstante, estas ventajas no se obtienen gratuitamente, ya que las simulaciones de partículas requieren grandísimos recursos. Por esta razón, es prácticamente obligatorio el uso de técnicas de procesamiento paralelo en este tipo de simulaciones. El vacío en el conocimiento de las sondas de Langmuir, es el que motiva nuestro trabajo. Nuestra aproximación, y el principal objetivo de este trabajo, ha sido desarrollar una simulación de partículas que nos permita estudiar el problema de una sonda de Langmuir inmersa en un plasma y que está negativamente polarizada con respecto a éste. Dicha simulación nos permitiría estudiar el comportamiento de los iones en los alrededores de una sonda cilíndrica de Langmuir, así como arrojar luz sobre la transición entre las teorías radiales y orbitales que ha sido observada experimentalmente. Justo después de esta sección introductoria, el resto de la tesis está dividido en tres partes tal y como sigue: La primera parte está dedicada a establecer los fundamentos teóricos de las sondas de Langmuir. En primer lugar, se realiza una introducción general al problema y al uso de sondas de Langmuir como método de diagnosis de plasmas. A continuación, se incluye una extensiva revisión bibliográfica sobre las diferentes teorías que proporcionan la corriente iónica recogida por una sonda. La segunda parte está dedicada a explicar los detalles de las simulaciones de partículas que han sido desarrolladas a lo largo de nuestra investigación, así como los resultados obtenidos con las mismas. Esta parte incluye una introducción sobre la teoría que subyace el tipo de simulaciones de partículas y las técnicas de paralelización que han sido usadas en nuestros códigos. El resto de esta parte está dividido en dos capítulos, cada uno de los cuales se ocupa de una de las geometrías consideradas en nuestras simulaciones (plana y cilíndrica). En esta parte discutimos también los descubrimientos realizados relativos a la transición entre el comportamiento radial y orbital de los iones en los alrededores de una sonda cilíndrica de Langmuir. Finalmente, en la tercera parte de la tesis se presenta un resumen del trabajo realizado. En este resumen, se enumeran brevemente los resultados de nuestra investigación y se han incluido algunas conclusiones. Después de esto, se enumeran una serie de perspectivas futuras y extensiones para los códigos desarrollados.

Magnetic and Acoustic Investigations of Turbulent Spherical Couette Flow

Title of dissertation: MAGNETIC AND ACOUSTIC INVESTIGATIONS OF TURBULENT SPHERICAL COUETTE FLOW Matthew M. Adams, Doctor of Philosophy, 2016 Dissertation directed by: Professor Daniel Lathrop Department of Physics This dissertation describes experiments in spherical Couette devices, using both gas and liquid sodium. The experimental geometry is motivated by the Earth's outer core, the seat of the geodynamo, and consists of an outer spherical shell and an inner sphere, both of which can be rotated independently to drive a shear flow in the fluid lying between them. In the case of experiments with liquid sodium, we apply DC axial magnetic fields, with a dominant dipole or quadrupole component, to the system. We measure the magnetic field induced by the flow of liquid sodium using an external array of Hall effect magnetic field probes, as well as two probes inserted into the fluid volume. This gives information about possible velocity patterns present, and we extend previous work categorizing flow states, noting further information that can be extracted from the induced field measurements. The limitations due to a lack of direct velocity measurements prompted us to work on developing the technique of using acoustic modes to measure zonal flows. Using gas as the working fluid in our 60~cm diameter spherical Couette experiment, we identified acoustic modes of the container, and obtained excellent agreement with theoretical predictions. For the case of uniform rotation of the system, we compared the acoustic mode frequency splittings with theoretical predictions for solid body flow, and obtained excellent agreement. This gave us confidence in extending this work to the case of differential rotation, with a turbulent flow state. Using the measured splittings for this case, our colleagues performed an inversion to infer the pattern of zonal velocities within the flow, the first such inversion in a rotating laboratory experiment. This technique holds promise for use in liquid sodium experiments, for which zonal flow measurements have historically been challenging.

Adaptive coatings based on polyaniline for direct 2D observation of diffusion processes in microfluidic systems

[EN] Therefore the understanding and proper evaluation of the flow and mixing behaviour at microscale becomes a very important issue. In this study, the diffusion behaviour of two reacting solutions of HCI and NaOH were directly observed in a glass/polydimethylsiloxane microfluidic device using adaptive coatings based on the conductive polymer polyaniline that are covalently attached to the microchannel walls. The two liquid streams were combined at the junction of a Y-shaped microchannel, and allowed to diffuse into each other and react. The results showed excellent correlation between optical observation of the diffusion process and the numerical results. A numerical model which is based on finite volume method (FVM) discretisation of steady Navier-Stokes (fluid flow) equations and mass transport equations without reactions was used to calculate the flow variables at discrete points in the finite volume mesh element. The high correlation between theory and practical data indicates the potential of such coatings to monitor diffusion processes and mixing behaviour inside microfluidic channels in a dye free environment.