Numerical and experimental investigations of three-dimensional container filling with Newtonian viscous fluids

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

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Assessment of a high-order finite difference upwind scheme for the simulation of convection-diffusion problems

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

Small-angle X-ray scattering and nitrogen adsorption study of the nanoporosity elimination in TEOS sonohydrolysis-derived xerogels

Small-angle X-ray scattering (SAXS) and nitrogen adsorption techniques were used to study the temperature and time structural evolution of the nanoporosity in silica xerogels prepared from acid- and ultrasound-catalyzed hydrolysis of tetraetboxysilane (TEOS). Silica xerogels present a structure of nanopores of fully random shape, size, and distribution, which can be described by an exponential correlation function gamma(r) = exp (-r/a), where a is the correlation distance, as predicted by the Debye, Anderson, and Brumberger (DAB) model. The mean pore size was evaluated as about 1.25 nm from SAXS and about 1.9 nm from nitrogen adsorption. The nanopore elimination in TEOS sonohydrolysis-derived silica xerogels is readily accelerated at temperatures around 900 degrees C probably by the action of a viscous flow mechanism. The nanopore elimination process takes place in such a way that the pore volume fraction and the specific surface are reduced while the mean pore size remains constant. (c) 2005 WILEY-VCH Verlag GmbH S Co. KGaA, Weinheim.

Sintering mechanisms of ZrO2 center dot MgO with addition of TiO2 and CuO

Substitutions of Ti and Cu in ZrO2.MgO (Z), cause transformation from monoclinic (m) to cubic (c) and tetragonal (t). According to the vacancy model and solid Solution formation models, neither CuO nor TiO2 cause zirconia stabilization, which derives front other phenomena. Data analysis by TMA using the CRH (constant rate of heating) method shows a solid state reaction of ZrO2.MgO.TiO2 (Z.TiO2) demonstrating a dominant mechanism of volume diffusion (n = 1). However, the sintering of ZrO2.MgO.CuO (Z.CuO) shows a viscous flow mechanism (n = 0), a similar phenomena to that of by sintering of glass. Transformations, such as: CuO to Cu2O at 1000 degreesC, ZrO2 (m) to ZrO2 (t) at 1100 degreesC and Cu2O (s) to Cu2O (l) at 1230 degreesC cause successive rearrangements of microstructure inside of region I (sintering process) and lead to interpretation errors when the Bannister equation is used. (C) 2003 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Viscosity of fluoroindate glasses

The viscosity of two fluoroindate glasses was measured as a function of temperature in the range of 310 °C - 362 °C. In such interval, the viscosity values were found to be similar to those reported for fluorozirconate glasses. The log η - 1/T plots had an unexpected behavior: two viscosity regions that seem to obey Arrhenius equation were identified and the activation energy for viscous flow (EA) for the region near Tg is smaller than the value found above the transition range. This behavior is probably due to structural changes occurred around Tg. The low values of the activation energy for viscous flow obtained for the indium fluoride-based glasses studied, suggest a good resistance against the devitrification process, what can make them suitable for fiber preparation.

Effect of lithium doping on the evolution of rheological and structural properties during gelation of siloxane-poly(oxypropylene) nanocomposites

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

A C2-continuous high-resolution upwind convection scheme

A bounded upwinding scheme for numerical solution of hyperbolic conservation laws and Navier-Stokes equations is presented. The scheme is based on convection boundedness criterion and total variation diminishing stability criteria and developed by employing continuously differentiable functions. The accuracy of the scheme is verified by assessing the error and observed convergence rate on 1-D benchmark test cases. A comparative study between the new scheme and conventional total variation diminishing/convection boundedness criterion-based upwind schemes to solve standard nonlinear hyperbolic conservation laws is also accomplished. The scheme is then examined in the simulation of Newtonian and non-Newtonian fluid flows of increasing complexity; a satisfactory agreement has been observed in terms of the overall behavior. Finally, the scheme is used to study the hydrodynamics of a gas-solid flow in a bubbling fluidized bed. © 2013 John Wiley & Sons, Ltd.

Aplicação da transformada integral e da transformação conforme na solução de uma classe de problemas difusivo-convectivos em domínios de geometrias não-convencionais

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

Vitrocerâmicas porosas via conformação e consolidação com amidos

The search for new methods of manufacture of glass-ceramics with controlled porosity and permeability is interesting of the industrial and commercial point of view, and a challenge of great complexity. Porous glass-ceramics produced by sintering and crystallization of glasses can find applications, for example, as filters, materials with bactericidal properties, bio-implants, as catalytic and enzymes supports, among others. An alternative and low cost method of ceramic manufacture reaching different levels of porosity, for diverse purposes, is the conformation assisted by addition of starch, known as consolidation with starch. The objective of this project is to study the process of conformation with starch for making porous glass-ceramics from a commercial glass in the system Na2O-CaO-SiO2, whose kinetics of sintering by viscous flow and surface crystallization are known. The method of conformation with starches is innovative for glass-ceramics and its development opens the way for obtaining a new class of materials. We found a possible route for the production of porous compacts of glass particles, from the powder preparation to the removal of starch. It was observed that a glass powder obtained by dry milling in a ball mill with alumina balls for 24 h, afterwards mixed with water in an eccentric ball mill for 2 h, without the addition of a deflocculant, and subsequently mixed with starch also in an eccentric ball mill for 5 min resulted in slurries stable against sedimentation

Modelagem numérica do escoamento em válvulas de compressores alternativos pelo método da fronteira imersa com refinamento adaptativo

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

Heat transfer in laminar flow of non-Newtonian fluids in ducts of elliptical section

Laminar forced convection inside tubes of various cross-section shapes is of interest in the design of a low Reynolds number heat exchanger apparatus. Heat transfer to thermally developing, hydrodynamically developed forced convection inside tubes of simple geometries such as a circular tube, parallel plate, or annular duct has been well studied in the literature and documented in various books, but for elliptical duct there are not much work done. The main assumption used in this work is a laminar flow of a power flow inside elliptical tube, under a boundary condition of first kind with constant physical properties and negligible axial heat diffusion (high Peclet number). To solve the thermally developing problem, we use the generalized integral transform technique (GITT), also known as Sturm-Liouville transform. Actually, such an integral transform is a generalization of the finite Fourier transform where the sine and cosine functions are replaced by more general sets of orthogonal functions. The axes are algebraically transformed from the Cartesian coordinate system to the elliptical coordinate system in order to avoid the irregular shape of the elliptical duct wall. The GITT is then applied to transform and solve the problem and to obtain the once unknown temperature field. Afterward, it is possible to compute and present the quantities of practical interest, such as the bulk fluid temperature, the local Nusselt number and the average Nusselt number for various cross-section aspect ratios. (C) 2006 Elsevier. SAS. All rights reserved.

Non-newtonian flow and pressure drop of pineapple juice in a plate heat exchanger

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

Control of instabilities in non-Newtonian free surface fluid flows

Free surface flows in inclined channels can develop periodic instabilities that are propagated downstream as shock waves with well-defined wavelengths and amplitudes. Such disturbances are called roll waves and are common in channels, torrential lava, landslides, and avalanches. The prediction and detection of such waves over certain types of structures and environments are useful for the prevention of natural risks. In this work, a mathematical model is established using a theoretical approach based on Cauchy's equations with the Herschel-Bulkley rheological model inserted into the viscous part of the stress tensor. This arrangement can adequately represent the behavior of muddy fluids, such as water-clay mixture. Then, taking into account the shallow water and the Rankine-Hugoniot's (shock wave) conditions, the equation of the roll wave and its properties, profile, and propagation velocity are determined. A linear stability analysis is performed with an emphasis on determining the condition that allows the generation of such instabilities, which depends on the minimum Froude number. A sensitivity analysis on the numerical parameters is performed, and numerical results including the influence of the Froude number, the index flow and dimensionless yield stress on the amplitude, the wavelength of roll waves and the propagation velocity of roll waves are shown. We show that our numerical results were in agreement with Coussot's experimental results (1994).