Numerical simulation of a radial diffuser turbulent airflow

In the present work are presented results from numerical simulations performed with the ANSYS-CFX (R) code. We have studied a radial diffuser flow case, which is the main academic problem used to study the flow behavior on flat plate valves. The radial flow inside the diffuser has important behavior such as the turbulence decay downstream and recirculation regions inside the valve flow channel due to boundary layer detachment. These flow structures are present in compressor reed valve configurations, influencing to a greater extent the compressor efficiency. The main target of the present paper was finding the simulation set-up (computational domain, boundary conditions and turbulence model) that better fits with experimental data published by Tabatabai and Pollard. The local flow turbulence and velocity profiles were investigated using four different turbulence models, two different boundary conditions set-up, two different computational domains and three different flow conditions (Re-in - Reynolds number at the diffuser inlet). We used the Reynolds stress (BSL); the k-epsilon; the RNG k-epsilon; and the shear stress transport (SST) k-omega turbulence models. The performed analysis and comparison of the computational results with experimental data show that the choice of the turbulence model, as well as the choice of the other computational conditions, plays an important role in the results physical quality and accuracy. (c) 2007 Elsevier B.V. All rights reserved.

Approximate solution for non-Darcy transient film condensation in a porous medium

The problem of non-darcian transient film condensation adjacent to a vertical flat plate embedded in a porous medium has been considered. The governing equation for the boundary layer thickness was obtained by an integral method and solved approximately by the method of integral relations. It is shown that the results are in good agreement with those obtained exactly by the method of characteristics.

Approximate solution for non-darcy transient film condensation in a porous medium

The problem of non-darcian transient film condensation adjacent to a vertical flat plate embedded in a porous medium has been considered. The governing equation for the boundary layer thickness was obtained by an integral method and solved approximately by the method of integral relations. It is shown that the results are in good agreement with those obtained exactly by the method of characteristics.

Active flutter suppression in a 2-D airfoil using linear matrix inequalities techniques

Flutter is an in-flight vibration of flexible structures caused by energy in the airstream absorbed by the lifting surface. This aeroelastic phenomenon is a problem of considerable interest in the aeronautic industry, because flutter is a potentially destructive instability resulting from an interaction between aerodynamic, inertial, and elastic forces. To overcome this effect, it is possible to use passive or active methodologies, but passive control adds mass to the structure and it is, therefore, undesirable. Thus, in this paper, the goal is to use linear matrix inequalities (LMIs) techniques to design an active state-feedback control to suppress flutter. Due to unmeasurable aerodynamic-lag states, one needs to use a dynamic observer. So, LMIs also were applied to design a state-estimator. The simulated model, consists of a classical flat plate in a two-dimensional flow. Two regulators were designed, the first one is a non-robust design for parametric variation and the second one is a robust control design, both designed by using LMIs. The parametric uncertainties are modeled through polytopic uncertainties. The paper concludes with numerical simulations for each controller. The open-loop and closed-loop responses are also compared and the results show the flutter suppression. The perfomance for both controllers are compared and discussed. Copyright © 2006 by ABCM.

Simulação numérica da convecção natural no interior de um refrigerador doméstico

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Análise numérica de evaporadores do tipo Roll-Bond usados em refrigeradores domésticos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Determinação da eficiência de um coletor solar plano com aletas de vidro

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

Avaliação da eficiência térmica de um coletor solar tipo tubo evacuado modificado

Pós-graduação em Engenharia Mecânica - FEB

Análise da relação desempenho/custo de alguns coletores solares comerciais

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Degradation of [D-Leu]-Microcystin-LR by solar heterogeneous photocatalysis (TiO2)

The present study investigates the use of solar heterogeneous photocatalyis (TiO2) for the destruction of [D-Leu]-Microcystin-LR, powerful toxin of widespread occurrence within cyanobacteria blooms. We extracted [D-Leu]-Microcystin-LR from a culture of Microcystis spp. and used a flat plate glass reactor coated with TiO2 (Degussa, P25) for the degradation studies. The irradiance was measured during the experiments with the aid of a spectroradiometer. After the degradation experiments, toxin concentrations were determined by HPLC and mineralization by TOC analyses. Acute and chronic toxicities were, quantified using mice and phosphatase inhibition in vitro assays, respectively. According to the performed experiments, 150 min were necessary to reduce the toxin concentration to the WHO's guideline for drinking water (from 10 to 1 mu g L-1) and to mineralize 90% of the initial carbon content. Another important finding is that solar heterogeneous photocatalysis was a destructive process indeed, not only for the toxin, but also for the other extract components and degradation products generated. Moreover, toxicity tests using mice have shown that the acute effect caused by the initial sample was removed. However, tests using the phosphatase enzyme indicated that it may be formed products capable of inducing chronic effects on mammals. The performed experiments indicate the feasibility of using solar heterogeneous photocatalysis for treating contaminated water with [D-Leu]-Microcystin-LR, not only due to its destruction, but also to the significant removal of organic matter and acute toxicity that can be achieved. (C) 2012 Elsevier Ltd. All rights reserved.

Study of the Compressive Strength of Aluminum Curved Elements

Currently, the Specification for Aluminum Structures (Aluminum Association, 2010) shows thin-walled aluminum plate sections with radii greater than eight inches have a lower compressive strength capacity than a flat plate with the same width and thickness. This inconsistency with intuition, which suggests any degree of folding a plate should increase its elastic buckling strength, inspired this study. A wide range of curvatures are studied—from a nearly flat plate to semi-circular. To quantify the curvature, a single non-dimensional parameter is used to represent all combinations of width, thickness and radius. Using the finite strip method (CU-FSM), elastic local buckling stresses are investigated. Using the ratio of stress values of curved plates compared to flat plates of the same size, equivalent plate-buckling coefficients are calculated. Using this data, nonlinear regression analyses are performed to develop closed form equations for five different edge support conditions. These equations can be used to calculate the elastic critical buckling stress for any curved aluminum section when the geometric properties (width, thickness, and radius) and the material properties (elastic modulus and Poisson’s ratio) are known. This procedure is illustrated in examples, each showing the applicability of the derived equations to geometries other than those investigated in this study and also providing comparisons with theoretically exact numerical analysis results.

In vivo electroporation and ubiquitin promoter--a protocol for sustained gene expression in the lung

BACKGROUND: Gene therapy applications require safe and efficient methods for gene transfer. Present methods are restricted by low efficiency and short duration of transgene expression. In vivo electroporation, a physical method of gene transfer, has evolved as an efficient method in recent years. We present a protocol involving electroporation combined with a long-acting promoter system for gene transfer to the lung. METHODS: The study was designed to evaluate electroporation-mediated gene transfer to the lung and to analyze a promoter system that allows prolonged transgene expression. A volume of 250 microl of purified plasmid DNA suspended in water was instilled into the left lung of anesthetized rats, followed by left thoracotomy and electroporation of the exposed left lung. Plasmids pCiKlux and pUblux expressing luciferase under the control of the cytomegalovirus immediate-early promoter/enhancer (CMV-IEPE) or human polyubiquitin c (Ubc) promoter were used. Electroporation conditions were optimized with four pulses (200 V/cm, 20 ms at 1 Hz) using flat plate electrodes. The animals were sacrificed at different time points up to day 40, after gene transfer. Gene expression was detected and quantified by bioluminescent reporter imaging (BLI) and relative light units per milligram of protein (RLU/mg) was measured by luminometer for p.Pyralis luciferase and immunohistochemistry, using an anti-luciferase antibody. RESULTS: Gene expression with the CMV-IEPE promoter was highest 24 h after gene transfer (2932+/-249.4 relative light units (RLU)/mg of total lung protein) and returned to baseline by day 3 (382+/-318 RLU/mg of total lung protein); at day 5 no expression was detected, whereas gene expression under the Ubc promoter was detected up to day 40 (1989+/-710 RLU/mg of total lung protein) with a peak at day 20 (2821+/-2092 RLU/mg of total lung protein). Arterial blood gas (PaO2), histological assessment and cytokine measurements showed no significant toxicity neither at day 1 nor at day 40. CONCLUSIONS: These results provide evidence that in vivo electroporation is a safe and effective tool for non-viral gene delivery to the lungs. If this method is used in combination with a long-acting promoter system, sustained transgene expression can be achieved.

Tranversal motion and flow structure of fully nonlinear streaks in a laminar boundary layer

Typical streak computations present in the literature correspond to linear streaks or to small amplitude nonlinear streaks computed using DNS or nonlinear PSE. We use the Reduced Navier-Stokes (RNS) equations to compute the streamwise evolution of fully non-linear streaks with high amplitude in a laminar flat plate boundary layer. The RNS formulation provides Reynolds number independent solutions that are asymptotically exact in the limit $Re \gg 1$, it requires much less computational effort than DNS, and it does not have the consistency and convergence problems of the PSE. We present various streak computations to show that the flow configuration changes substantially when the amplitude of the streaks grows and the nonlinear effects come into play. The transversal motion (in the wall normal-streamwise plane) becomes more important and strongly distorts the streamwise velocity profiles, that end up being quite different from those of the linear case. We analyze in detail the resulting flow patterns for the nonlinearly saturated streaks and compare them with available experimental results.

Computation of nonlinear streaky structures in boundary layer flow

En esta tesis se integran numéricamente las ecuaciones reducidas de Navier Stokes (RNS), que describen el flujo en una capa límite tridimensional que presenta también una escala característica espacial corta en el sentido transversal. La formulación RNS se usa para el cálculo de “streaks” no lineales de amplitud finita, y los resultados conseguidos coinciden con los existentes en la literatura, obtenidos típicamente utilizando simulación numérica directa (DNS) o nonlinear parabolized stability equations (PSE). El cálculo de los “streaks” integrando las RNS es mucho menos costoso que usando DNS, y no presenta los problemas de estabilidad que aparecen en la formulación PSE cuando la amplitud del “streak” deja de ser pequeña. El código de integración RNS se utiliza también para el cálculo de los “streaks” que aparecen de manera natural en el borde de ataque de una placa plana en ausencia de perturbaciones en la corriente uniforme exterior. Los resultados existentes hasta ahora calculaban estos “streaks” únicamente en el límite lineal (amplitud pequeña), y en esta tesis se lleva a cabo el cálculo de los mismos en el régimen completamente no lineal (amplitud finita). En la segunda parte de la tesis se generaliza el código RNS para incluir la posibilidad de tener una placa no plana, con curvatura en el sentido transversal que varía lentamente en el sentido de la corriente. Esto se consigue aplicando un cambio de coordenadas, que transforma el dominio físico en uno rectangular. La formulación RNS se integra también expresada en las correspondientes coordenadas curvilíneas. Este código generalizado RNS se utiliza finalmente para estudiar el flujo de capa límite sobre una placa con surcos que varían lentamente en el sentido de la corriente, y es usado para simular el flujo sobre surcos que crecen en tal sentido. Abstract In this thesis, the reduced Navier Stokes (RNS) equations are numerically integrated. This formulation describes the flow in a three-dimensional boundary layer that also presents a short characteristic space scale in the spanwise direction. RNS equations are used to calculate nonlinear finite amplitude “streaks”, and the results agree with those reported in the literature, typically obtained using direct numerical simulation (DNS) or nonlinear parabolized stability equations (PSE). “Streaks” simulations through the RNS integration are much cheaper than using DNS, and avoid stability problems that appear in the PSE when the amplitude of the “streak” is not small. The RNS integration code is also used to calculate the “streaks” that naturally emerge at the leading edge of a flat plate boundary layer in the absence of any free stream perturbations. Up to now, the existing results for these “streaks” have been only calculated in the linear limit (small amplitude), and in this thesis their calculation is carried out in the fully nonlinear regime (finite amplitude). In the second part of the thesis, the RNS code is generalized to include the possibility of having a non-flat plate, curved in the spanwise direction and slowly varying in the streamwise direction. This is achieved by applying a change of coordinates, which transforms the physical domain into a rectangular one. The RNS formulation expressed in the corresponding curvilinear coordinates is also numerically integrated. This generalized RNS code is finally used to study the boundary layer flow over a plate with grooves which vary slowly in the streamwise direction; and this code is used to simulate the flow over grooves that grow in the streamwise direction.

Experimental studies on transitional separated boundary layers

Separated transitional boundary layers appear on key aeronautical processes such as the flow around wings or turbomachinery blades. The aim of this thesis is the study of these flows in representative scenarios of technological applications, gaining knowledge about phenomenology and physical processes that occur there and, developing a simple model for scaling them. To achieve this goal, experimental measurements have been carried out in a low speed facility, ensuring the flow homogeneity and a low disturbances level such that unwanted transitional mechanisms are avoided. The studied boundary layers have been developed on a flat plate, by imposing a pressure gradient by means of contoured walls. They generate an initial acceleration region followed by a deceleration zone. The initial region is designed to obtain at the beginning of the deceleration the Blasius profile, characterized by its momentum thickness, and an edge boundary layer velocity, defining the problem characteristic velocity. The deceleration region is designed to obtain a linear evolution of the edge velocity, thereby defining the characteristic length of the problem. Several experimental techniques, both intrusive (hot wire anemometry, total pressure probes) as nonintrusive (PIV and LDV anemometry, high-speed filming), have been used in order to take advantage of each of them and allow cross-validation of the results. Once the boundary layer at the deceleration beginning has been characterized, ensuring the desired integral parameters and level of disturbance, the evolution of the laminar boundary layer up to the point of separation is studied. It has been compared with integral methods, and numerical simulations. In view of the results a new model for this evolution is proposed. Downstream from the separation, the flow near to the wall is configured as a shear layer that encloses low momentum recirculating fluid. The region where the shear layer remains laminar tends to be positioned to compensate the adverse pressure gradient associated with the imposed deceleration. Under these conditions, the momentum thickness remains almost constant. This laminar shear layer region extends up to where transitional phenomena appear, extension that scales with the momentum thickness at separation. These transitional phenomena are of inviscid type, similar to those found in free shear layers. The transitional region analysis begins with a study of the disturbances evolution in the linear growth region and the comparison of experimental results with a numerical model based on Linear Stability Theory for parallel flows and with data from other authors. The results’ coalescence for both the disturbances growth and the excited frequencies is stated. For the transition final stages the vorticity concentration into vortex blobs is found, analogously to what happens in free shear layers. Unlike these, the presence of the wall and the pressure gradient make the large scale structures to move towards the wall and quickly disappear under certain circumstances. In these cases, the recirculating flow is confined into a closed region saying the bubble is closed or the boundary layer reattaches. From the reattachment point, the fluid shows a configuration in the vicinity of the wall traditionally considered as turbulent. It has been observed that existing integral methods for turbulent boundary layers do not fit well to the experimental results, due to these methods being valid only for fully developed turbulent flow. Nevertheless, it has been found that downstream from the reattachment point the velocity profiles are self-similar, and a model has been proposed for the evolution of the integral parameters of the boundary layer in this region. Finally, the phenomenon known as bubble burst is analyzed. It has been checked the validity of existing models in literature and a new one is proposed. This phenomenon is blamed to the inability of the large scale structures formed after the transition to overcome with the adverse pressure gradient, move towards the wall and close the bubble. El estudio de capas límites transicionales con separación es de gran relevancia en distintas aplicaciones tecnológicas. Particularmente, en tecnología aeronáutica, aparecen en procesos claves, tales como el flujo alrededor de alas o álabes de turbomaquinaria. El objetivo de esta tesis es el estudio de estos flujos en situaciones representativas de las aplicaciones tecnológicas, ganando por un lado conocimiento sobre la fenomenología y los procesos físicos que aparecen y, por otra parte, desarrollando un modelo sencillo para el escalado de los mismos. Para conseguir este objetivo se han realizado ensayos en una instalación experimental de baja velocidad específicamente diseñada para asegurar un flujo homogéneo y con bajo nivel de perturbaciones, de modo que se evita el disparo de mecanismos transicionales no deseados. La capa límite bajo estudio se ha desarrollado sobre una placa plana, imponiendo un gradiente de presión a la misma por medio de paredes de geometría especificada. éstas generan una región inicial de aceleración seguida de una zona de deceleración. La región inicial se diseña para tener en al inicio de la deceleración un perfil de capa límite de Blasius, caracterizado por su espesor de cantidad de movimiento, y una cierta velocidad externa a la capa límite que se considera la velocidad característica del problema. La región de deceleración está concebida para que la variación de la velocidad externa a la capa límite sea lineal, definiendo de esta forma una longitud característica del problema. Los ensayos se han realizado explotando varias técnicas experimentales, tanto intrusivas (anemometría de hilo caliente, sondas de presión total) como no intrusivas (anemometrías láser y PIV, filmación de alta velocidad), de cara a aprovechar las ventajas de cada una de ellas y permitir validación cruzada de resultados entre las mismas. Caracterizada la capa límite al comienzo de la deceleración, y garantizados los parámetros integrales y niveles de perturbación deseados se procede al estudio de la zona de deceleración. Se presenta en la tesis un análisis de la evolución de la capa límite laminar desde el inicio de la misma hasta el punto de separación, comparando con métodos integrales, simulaciones numéricas, y proponiendo un nuevo modelo para esta evolución. Aguas abajo de la separación, el flujo en las proximidades de la pared se configura como una capa de cortadura que encierra una región de fluido recirculatorio de baja cantidad de movimiento. Se ha caracterizado la región en que dicha capa de cortadura permanece laminar, encontrando que se posiciona de modo que compensa el gradiente adverso de presión asociado a la deceleración de la corriente. En estas condiciones, el espesor de cantidad de movimiento permanece prácticamente constante y esta capa de cortadura laminar se extiende hasta que los fenómenos transicionales aparecen. Estos fenómenos son de tipo no viscoso, similares a los que aparecen en una capa de cortadura libre. El análisis de la región transicional comienza con un estudio de la evolución de las vii viii RESUMEN perturbaciones en la zona de crecimiento lineal de las mismas y la comparación de los resultados experimentales con un modelo numérico y con datos de otros autores. La coalescencia de los resultados tanto para el crecimiento de las perturbaciones como para las frecuencias excitadas queda demostrada. Para los estadios finales de la transición se observa la concentración de la vorticidad en torbellinos, de modo análogo a lo que ocurre en capas de cortadura libres. A diferencia de estas, la presencia de la pared y del gradiente de presión hace que, bajo ciertas condiciones, la gran escala se desplace hacia la pared y desaparezca rápidamente. En este caso el flujo recirculatorio queda confinado en una región cerrada y se habla de cierre de la burbuja o readherencia de la capa límite. A partir del punto de readherencia se tiene una configuración fluida en las proximidades de la pared que tradicionalmente se ha considerado turbulenta. Se ha observado que los métodos integrales existentes para capas límites turbulentas no ajustan bien a las medidas experimentales realizadas, hecho imputable a que no se obtiene en dicha región un flujo turbulento plenamente desarrollado. Se ha encontrado, sin embargo, que pasado el punto de readherencia los perfiles de velocidad próximos a la pared son autosemejantes entre sí y se ha propuesto un modelo para la evolución de los parámetros integrales de la capa límite en esta región. Finalmente, el fenómeno conocido como “estallido” de la burbuja se ha analizado. Se ha comprobado la validez de los modelos existentes en la literatura y se propone uno nuevo. Este fenómeno se achaca a la incapacidad de la gran estructura formada tras la transición para vencer el gradiente adverso de presión, desplazarse hacia la pared y cerrar la burbuja.