Application of a bounded upwinding scheme to complex fluid dynamics problems

This paper is concerned with an overview of upwinding schemes, and further nonlinear applications of a recently introduced high resolution upwind differencing scheme, namely the ADBQUICKEST [V.G. Ferreira, F.A. Kurokawa, R.A.B. Queiroz, M.K. Kaibara, C.M. Oishi, J.A.Cuminato, A.F. Castelo, M.F. Tomé, S. McKee, assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems, International Journal for Numerical Methods in Fluids 60 (2009) 1-26]. The ADBQUICKEST scheme is a new TVD version of the QUICKEST [B.P. Leonard, A stable and accurate convective modeling procedure based on quadratic upstream interpolation, Computer Methods in Applied Mechanics and Engineering 19 (1979) 59-98] for solving nonlinear balance laws. The scheme is based on the concept of NV and TVD formalisms and satisfies a convective boundedness criterion. The accuracy of the scheme is compared with other popularly used convective upwinding schemes (see, for example, Roe (1985) [19], Van Leer (1974) [18] and Arora & Roe (1997) [17]) for solving nonlinear conservation laws (for example, Buckley-Leverett, shallow water and Euler equations). The ADBQUICKEST scheme is then used to solve six types of fluid flow problems of increasing complexity: namely, 2D aerosol filtration by fibrous filters; axisymmetric flow in a tubular membrane; 2D two-phase flow in a fluidized bed; 2D compressible Orszag-Tang MHD vortex; axisymmetric jet onto a flat surface at low Reynolds number and full 3D incompressible flows involving moving free surfaces. The numerical simulations indicate that this convective upwinding scheme is a good generic alternative for solving complex fluid dynamics problems. © 2012.

Simulação do escoamento bifásico da mistura óleo-refrigerante através da folga radial de compressores rotativos de pistão rolante

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

Modelagem do escoamento transiente ao longo de trocadores de calor tubo capilar-linha de sucção

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

Modelo distribuído aplicado à análise de evaporadores do tipo tubo aletado

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

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

AN EXPERIMENTAL EVALUATION OF THE PERFORMANCE OF A VORTEX-DRIVEN AIR LIFT PUMP

A prototype vortex-driven air lift pump was developed and experimentally evaluated. It was designed to be easily manufactured and scalable for arbitrary riser diameters. The model tested fit in a 2 inch diameter riser with six air injection nozzles through which airwas injected helically around the perimeter of the riser at an angle of 70º from pure tangential injection. The pump was intended to transport both water and sediment over a large range of submergence ratios. A test apparatus was designed to be able to simulate deep water or oceanic environments. The resulting test setup had a finite reservoir; over the course of a test, the submergence ratio varied from 0.48 to 0.39. For air injection pressures ranging from 10 to 60 psig and for air flow rates of 6 to 15 scfm, the induced water discharge flow rates varied only slightly, due to the limited range of available submergence ratios. The anticipated simulation of deep water environment, with a corresponding equivalent increase in thesubmergence ratio, proved unattainable. The pump prototype successfully transported both water and sediment (sand). Thepercent volume yield of the sediment was in an acceptable range. The pump design has been subsequently used successfully in a 4 inch configuration in a follow-on project. A computer program was written in Matlab to simulate the pump characteristics. The program output water pressures at the location of air injection which were physicallycompatible with the experimental data.

Estudo fenomenológico e numérico do escoamento estratificado óleo-água ondulado e com mistura na interface

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura.

Polimerização em emulsão de estireno em microrreator.

A polimerização em emulsão de estireno em um microrreator Syrris de 250 µL com misturador estático junção \"T\" foi estudada em duas etapas. Primeiro somente a fluidodinâmica deste dispositivo não convencional foi avaliada, depois, foi desenvolvida a reação de polimerização de forma a observar como este fator influencia no sistema. Os experimentos foram realizados procurando se atingir maiores conversões, mas mantendo a estabilidade da emulsão. Foi um trabalho exploratório, portanto se assemelha mais a um processo de evolução (evolutionary process). Foram verificados a partir de qual relação das vazões dos dois fluidos ocorre a formação de gotas, e que com o aumento da vazão da fase contínua, aquosa (Qc), mantendo constante a vazão da fase dispersa (Qd), foi verificado uma diminuição do diâmetro das gotas e um regime de fluxo laminar. Posteriormente, realizou-se a polimerização em emulsão do estireno no microrreator, porém com restrições para altas vazões. Os parâmetros de processo testados foram a proporção Qc e Qd, a temperatura e a concentração do iniciador para então verificar o efeito que a variação destas ocasionam na conversão de monômero, no diâmetro e número de partículas e nas massas moleculares médias. A polimerização foi feita para soma das vazões Qc e Qd da ordem de 100 µ L/min, com 15% de monômero na formulação e com o maior tempo de residência possível de 2,5 minutos. Para maiores concentrações de monômero, acima de 15% foi verificado entupimento do canal do microrreator. A taxa de conversão de monômero aumentou com o aumento da temperatura e com o aumento da concentração do iniciador, mas o maior valor atingido foi de apenas 37% devido ao baixo tempo de residência. Nos casos de maiores taxas de conversão, as massas moleculares obtidas foram as menores conforme o esperado pela teoria. Finalmente, os índices de polidispersão (PDI), obtidos foram da ordem de 2,5 a 3,5.

Coalescence of secondary dispersions in fibrous beds

SUMMARY A study has been made of the coalescence of secondary dispersions in a fibrous bed. The literature pertaining to the formation, hydrodynamic behaviour and methods of separation of droplets less than one hundred micrometres in diameter has been reviewed with particular reference to fibrous bed coalescers. The main operating parameters were identified as inlet drop size distribution, phase ratio, superficial velocity and the thickness and voidage of the bed . A recirculatory rig with interchangeable fibrous bed pads was designed and operated with toluene-water dispersions generated by a combination of centrifugal pumps . Inlet drop sizes were analysed using a Coulter Counter and outlet drops were sized photographically. A novel technique, involving conductivity measur ements at different planes in the bed, was developed to measure hold up distribution. Single phase flow and two phase flow pressure drops were correlated by a Blake-Kozeny type equation. Exit drop size was independent of inlet drop size distribution and phase ratio but a function of superficialvelocity and packing thickness. Average bed hold up was independent of inlet drop size distribution and phase ratio, but decreased with increase in superficial velocity. Hold up was not evenly distributed in the bed, the highest value occurred at the inlet followed by a sharp -2 drop at approximately 1.2 x 10 m. Hold up remained constant throughout the rest of the bed until the exit plane, where it increased. From the results, a mechanism is postulated involving: (a) Capture of the inlet drops followed by interdrop coalescence until an equilibrium value is reached. (b) Equilibrium size droplets flowing as rivulets through the intermediate portion of the bed, and (c) Each rivulet forms droplets at the exit face, which detach by a 'drip point' mechanism.

Levantamento do Tamanho Médio de Bolhas em Escoamento Bifásico Vertical por Técnica Fotográfica

Para confirmar e validar códigos computacionais de escoamento de fluidos, importantes na segurança de reatores nucleares, é preciso que sejam realizadas medidas experimentais. Como a geração de vapor no fluido refrigerante diminui a eficiência de troca térmica a fração de vazio passa a ser um dos parâmetros mais importantes a ser obtido, e que, num escoamento vertical, pode ser conhecido através da quantidade de bolha e seu diâmetro médio.

Automatic dynamic mask extraction for PIV images containing an unsteady interface, bubbles, and a moving structure

When performing Particle Image Velocimetry (PIV) measurements in complex fluid flows with moving interfaces and a two-phase flow, it is necessary to develop a mask to remove non-physical measurements. This is the case when studying, for example, the complex bubble sweep-down phenomenon observed in oceanographic research vessels. Indeed, in such a configuration, the presence of an unsteady free surface, of a solid–liquid interface and of bubbles in the PIV frame, leads to generate numerous laser reflections and therefore spurious velocity vectors. In this note, an image masking process is developed to successively identify the boundaries of the ship and the free surface interface. As the presence of the solid hull surface induces laser reflections, the hull edge contours are simply detected in the first PIV frame and dynamically estimated for consecutive ones. As for the unsteady surface determination, a specific process is implemented like the following: i) the edge detection of the gradient magnitude in the PIV frame, ii) the extraction of the particles by filtering high-intensity large areas related to the bubbles and/or hull reflections, iii) the extraction of the rough region containing these particles and their reflections, iv) the removal of these reflections. The unsteady surface is finally obtained with a fifth-order polynomial interpolation. The resulted free surface is successfully validated from the Fourier analysis and by visualizing selected PIV images containing numerous spurious high intensity areas. This paper demonstrates how this data analysis process leads to PIV images database without reflections and an automatic detection of both the free surface and the rigid body. An application of this new mask is finally detailed, allowing a preliminary analysis of the hydrodynamic flow.