Estudo do escoamento ao redor de cilindros elípticos

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

Implementação e teste do método da fronteira imersa para a simulação do escoamento em torno de cilindros estacionários e rotativos

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

Modelagem numérica do escoamento em válvulas automáticas de compressores pelo Método da Fronteira Imersa

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Análise numérica de condensadores do tipo arame-sobre-tubo usados em refrigeradores domésticos

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

Modelagem do escoamento ao longo de evaporadores de serpentina com tubos aletados

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

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)

Estabilidade e controle dinâmico de roll waves

Pós-graduação em Engenharia Elétrica - FEIS

Análise experimental do comportamento fluidodinâmico e da secagem de sementes de linhaça (Linum usitatissimum L.) em leito de jorro

A linhaça é a semente da planta do linho (Linum usitatissimum L.), uma espécie polimorfa originária do linho, sendo considerada uma das 6 plantas atualmente reconhecidas pelo Instituto Nacional do Câncer dos Estados Unidos (US National Cancer Institute - NCI) por suas propriedades específicas no combate ao câncer. Parte desse reconhecimento deve-se a notável característica de ser a fonte mais rica de precursores de lignina (esteróide vegetal de ação análoga ao estrógeno de mamíferos) na dieta humana. A variedade utilizada neste trabalho foi a “Linseed” de cor marrom e o objetivo deste trabalho foi analisar a fluidodinâmica dessa partícula em leito de jorro, estabelecida pelas medidas de tomadas de queda de pressão no leito a partir das deflexões em relação às velocidades de ar crescente e decrescente, obtendo assim informações para a determinação de parâmetros correlacionados ao processo, como: velocidade de mínimo jorro, queda de pressão máxima, queda de pressão no jorro estável e queda de pressão no mínimo jorro. Estes valores foram comparados aos correspondentes valores obtidos por equações empíricas citadas na literatura. Foi também avaliado o comportamento da secagem da matéria prima, mediante um planejamento estatístico 2², tendo com as variáveis de entrada temperatura do gás (T_g) e o tempo de operação (t), para a quantificação das variáveis de resposta razão de umidade (Xr, adim.), germinação (G, %) e o índice de velocidade de germinação (IVG, t^-1) e analisadas estatisticamente pelo planejamento fatorial completo com três repetições no ponto central. A cinética de secagem das sementes de linhaça, previamente umidificadas, foi realizada nas temperaturas de 45, 55 e 65 °C, e dentre os três modelos propostos, o modelo de Midilli et al., foi que melhor descreveu aos dados experimentais. Para os parâmetros fluidodinâmicos observou-se que a correlação de Gorshtein e Mukhlenov (1965) apresentou os menores desvios para queda e pressão de mínimo jorro e jorro estável, Abdelrazek (1969) apresentou o menor desvio para a velocidade no mínimo jorro e Pallai e Németh (1969) descreveu adequadamente a queda de pressão máxima. Foi observado que a carga de sementes e a temperatura exerceram influência significativa nos parâmetros fluidodinâmicos em leito de jorro. Com base na análise do planejamento estatístico proposto pode-se concluir que os parâmetros de entrada temperatura do gás e tempo de operação exerceram influência significativa sobre todas as variáveis de resposta, sendo observada influência quadrática das variáveis de entrada ao se observar a significância da curvatura sobre os parâmetros: Razão de umidade, Germinação e Índice de Velocidade de Germinação, propondo-se modelos representativos destes parâmetros com a presença da curvatura, apresentando um coeficiente de determinação (R²) superior a 99 %.

Análise do escoamento turbulento por um queimador industrial a gás utilizando dinâmica dos fluidos computacional

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

Theoretical and experimental study of agricultural spraying using CFD

This paper presents a computational fluid dynamics (CFD) application about the axial fan design used in an agricultural spraying system with a theoretical and experimental analysis of comparative results between the characteristic curves of a fan for several rotations and numerical results for the influence of blade attack angle variation and optimization of the spraying system, both for a same rotation. Flow was considered three-dimensional, turbulent, isothermal, viscous and non-compressible in a steady state, disregarding any influence of the gravity field. The average turbulent field was obtained from the application of time average where the turbulence model required for closing the set of equations was the k-E model. Resolution of all connected phenomena was achieved with the help of a fluid dynamics computer, CFX, which uses the finite volumes technique as a numerical method. In order to validate the theoretical analysis, an experiment was conducted in a circular section of a horizontal wind tunnel, using a Pitot tube for pressure readings. The main results demonstrate that the methodology used, based on CFD techniques, is able to reproduce the phenomenological behavior of an axial fan in a spraying system because results were very reliable and similar to experimentally measured ones.

Numerical analysis of a modified evacuated tubes solar collector

The need for renewable energy sources, facing the consequences of Climate Change, results in growing investment for solar collectors’ use. Research in this field has accompanied this expansion and evacuated tube solar collector stands as an important study focus. Thus, several works have been published for representing the stratification of the fluid inside the tubes and the reservoir, as well as analytical modeling for the heat flow problem. Based on recent publications, this paper proposes the study of solar water heating with evacuated tubes, their operation characteristics and operating parameters. To develop this work, a computational tool will be used - in this case, the application of computational fluid dynamics (CFD) software. In possession of the implemented model, a numerical simulation will be performed to evaluate the behavior of the fluid within this solar collector and possible improvements to be applied in the model.

Numerical analysis of control surface effects on AUV

Computational fluid dynamics, CFD, is becoming an essential tool in the prediction of the hydrodynamic efforts and flow characteristics of underwater vehicles for manoeuvring studies. However, when applied to the manoeuvrability of autonomous underwater vehicles, AUVs, most studies have focused on the de- termination of static coefficients without considering the effects of the vehicle control surface deflection. This paper analyses the hydrodynamic efforts generated on an AUV considering the combined effects of the control surface deflection and the angle of attack using CFD software based on the Reynolds-averaged Navier–Stokes formulations. The CFD simulations are also independently conducted for the AUV bare hull and control surface to better identify their individual and interference efforts and to validate the simulations by comparing the experimental results obtained in a towing tank. Several simulations of the bare hull case were conducted to select the k –ω SST turbulent model with the viscosity approach that best predicts its hydrodynamic efforts. Mesh sensitivity analyses were conducted for all simulations. For the flow around the control surfaces, the CFD results were analysed according to two different methodologies, standard and nonlinear. The nonlinear regression methodology provides better results than the standard methodology does for predicting the stall at the control surface. The flow simulations have shown that the occurrence of the control surface stall depends on a linear relationship between the angle of attack and the control surface deflection. This type of information can be used in designing the vehicle’s autopilot system.

Simulación URANS 2D del flujo alrededor de un cilindro para 104 ≤ Re ≤ 106 con tratamiento de paredes

Máster Universitario en Sistemas Inteligentes y Aplicaciones Numéricas en Ingeniería (SIANI)

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Microfluidic device and interfacial transport: application to biomolecules and nanostructures

The aim of my dissertation is to provide new knowledge and applications of microfluidics in a variety of problems, from materials science, devices, and biomedicine, where the control on the fluid dynamics and the local concentration of the solutions containing the relevant molecules (either materials, precursors, or biomolecules) is crucial. The control of interfacial phenomena occurring in solutions at dierent length scales is compelling in nanotechnology for devising new sensors, molecular electronics devices, memories. Microfluidic devices were fabricated and integrated with organic electronics devices. The transduction involves the species in the solution which infills the transistor channel and confined by the microfluidic device. This device measures what happens on the surface, at few nanometers from the semiconductor channel. Soft-lithography was adopted to fabricate platinum electrodes, starting from platinum carbonyl precursor. I proposed a simple method to assemble these nanostructures in periodic arrays of microstripes, and form conductive electrodes with characteristic dimension of 600 nm. The conductivity of these sub-microwires is compared with the values reported in literature and bulk platinum. The process is suitable for fabricating thin conductive patterns for electronic devices or electrochemical cells, where the periodicity of the conductive pattern is comparable with the diusion length of the molecules in solution. The ordering induced among artificial nanostructures is of particular interest in science. I show that large building blocks, like carbon nanotubes or core-shell nanoparticles, can be ordered and self-organised on a surface in patterns due to capillary forces. The eective probability of inducing order with microfluidic flow is modeled with finite element calculation on the real geometry of the microcapillaries, in soft-lithographic process. The oligomerization of A40 peptide in microconfined environment represents a new investigation of the extensively studied peptide aggregation. The added value of the approach I devised is the precise control on the local concentration of peptides together with the possibility to mimick cellular crowding. Four populations of oligomers where distinguished, with diameters ranging from 15 to 200 nm. These aggregates could not be addresses separately in fluorescence. The statistical analysis on the atomic force microscopy images together with a model of growth reveal new insights on the kinetics of amyloidogenesis as well as allows me to identify the minimum stable nucleus size. This is an important result owing to its implications in the understanding and early diagnosis and therapy of the Alzheimer’s disease