Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.

The Granada ignimbrite: A compound pyroclastic unit and its relationship with Upper Miocene caldera volcanism in the northern Puna

The Granada ignimbrite, an Upper Miocene volcanic unit from the northern Puna, previously has been interpreted as an extensive ignimbrite (>2300 km(2)) associated with eruptions from the Vilama caldera (trap-door event). On the basis of new data, we revise its correlation and redefine the unit as a compound, high aspect ratio ignimbrite, erupted at approximately 9.8 Ma. Calculated volumes (similar to 100 km(3)) are only moderate in comparison with other large volume (>1000 km(3)) ignimbrites that erupted approximately 2-6 m.y. later in the region (e.g. Vilama, Panizos, Atana). Six new volcanic units are recognized from sequences previously correlated with Granada (only one sourced from the same center). Consequently, the area ascribed to the Granada ignimbrite is substantially reduced (630 km(2)), and links to the Vilama caldera are not supported. Transport directions suggest the volcanic source for the Granada ignimbrite corresponds to vents buried under younger (>= 7.9-5 Ma) volcanic rocks of the Abra Granada volcanic complex. Episodes of caldera collapse at some stage of eruption are likely, though their nature and timing cannot be defined from available data. The eruption of the Granada ignimbrite marks the onset of a phase of large volume (caldera-sourced) volcanism in the northern Puna. (C) 2007 Elsevier Ltd. All rights reserved.

Well-preserved Late Paleoproterozoic volcanic centers in the Sao Felix do Xingu region, Amazonian Craton, Brazil

The Amazonian craton in the Sao Felix do Xingu city, southeast region of the Para state, north of Brazil, hosts exceptionally well-preserved Paleoproterozoic bimodal magmatic units grouped in the Sobreiro and Santa Rosa formations. These formations are correlated to the Uatuma magmatic event, which is largely distributed in the Amazonian craton occupying more than 1,500,000 km(2). Geological mapping and petrographical observations reveal distinct spectra of volcanic facies in both formations. The basal calc-alkaline Sobreiro Formation is composed mainly of andesitic and dacitic lava flows and associated volcaniclastic facies of autoclastic origin, with subordinate pyroclastic flow deposits. This formation shows inferred eruption style that is similar to those in Flood Basalt Provinces, with rare scutulum-type lava shields. The upper A-type Santa Rosa Formation was generated by multicyclic explosive and effusive episodes predominantly associated with large fissures and is materialized by voluminous ignimbrites with subordinated ash-fall tuff, crystal tuff, lapilli-tuff, co-ignimbritic breccias, rhyolitic dikes and domes, and associated granitic porphyries and equigranular granitic intrusions. Ignimbrite and rhyolite dikes reveal conspicuous vertical flow pattern pointing to a fissure-controlled eruption, similar to Sierra Madre Occidental ignimbrite province. The proposed evolutionary model for the Sao Felix do Xingu units differs from those of other occurrences related to the Uatuma magmatic event in the Amazonian craton, characterized by predominance of A-type volcanism and contemporaneous granites. (C) 2010 Elsevier B.V. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 Elsevier B.V. All rights reserved.

Accurate velocity reconstruction for Discontinuous Galerkin approximations of two-phase porous media flows

We consider Discontinuous Galerkin approximations of two-phase, immiscible porous media flows in the global pressure/fractional flow formulation with capillary pressure. A sequential approach is used with a backward Euler step for the saturation equation, equal-order interpolation for the pressure and the saturation, and without any limiters. An accurate total velocity field is recovered from the global pressure equation to be used in the saturation equation. Numerical experiments show the advantages of the proposed reconstruction. To cite this article: A. Ern et al., C R. Acad. Sci. Paris, Ser. 1347 (2009). (C) 2009 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Forced convection to laminar flow of liquid egg yolk in circular and annular ducts

The steady-state heat transfer in laminar flow of liquid egg yolk - an important pseudoplastic fluid food - in circular and concentric annular ducts was experimentally investigated. The average convection heat transfer coefficients, determined by measuring temperatures before and after heating sections with constant temperatures at the tube wall, were used to obtain simple new empirical expressions to estimate the Nusselt numbers for fully established flows at the thermal entrance of the considered geometries. The comparisons with existing correlations for Newtonian and non-Newtonian fluids resulted in excellent agreement. The main contribution of this work is to supply practical and easily applicable correlations, which are, especially for the case of annulus, rather scarce and extensively required in the design of heat transfer operations dealing with similar shear-thinning products. In addition, the experimental results may support existing theoretical analyses.

Oxigenoterapia domiciliar prolongada em crianças e adolescentes: uma análise do uso clínico e de custos de um programa assistencial

OBJETIVOS: Descrever as características clínicas e laboratoriais dos pacientes em oxigenoterapia domiciliar prolongada acompanhados pelo programa de atendimento domiciliar do Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo, durante um período de 8 anos, e comparar os grupos com e sem hipertensão pulmonar secundária. Estimar o custo do programa utilizando concentradores versus cilindros de oxigênio arcados pela instituição. MÉTODOS: Estudo descritivo retrospectivo e de coorte dos pacientes em oxigenoterapia domiciliar prolongada, em seguimento no período de 2002 a 2009, na Unidade de Pneumologia do Instituto da Criança do Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo. RESULTADOS: Foram estudados 165 pacientes, 53% do sexo masculino, com as medianas: idade de início da oxigenoterapia - 3,6 anos; tempo de oxigenoterapia - 7 anos; e tempo de sobrevida após início da oxigenoterapia - 3,4 anos. Os principais diagnósticos foram: fibrose cística (22%), displasia broncopulmonar (19%) e bronquiolite obliterante (15%). Dos 33 pacientes que realizaram espirometria, 70% apresentavam distúrbio ventilatório obstrutivo grave. O exame ecocardiográfico foi executado em 134 pacientes; 51% deles tinham hipertensão pulmonar secundária. Houve associação estatisticamente significante entre: presença de hipertensão pulmonar e necessidade de maiores fluxos de oxigênio (qui-quadrado, p = 0,011); e presença de hipertensão pulmonar e maior tempo de oxigenoterapia (Logrank, p = 0,0001). Não houve diferença estatisticamente significante entre tempo de sobrevida após início da oxigenoterapia e presença de hipertensão pulmonar. Os custos médios mensais do programa foram: US$ 7.392,93 para os concentradores e US$ 16.630,92 para cilindros. CONCLUSÕES: A oxigenoterapia domiciliar prolongada foi empregada em distintas doenças crônicas, predominantemente em lactentes e pré-escolares. Houve alta frequência de hipertensão pulmonar associada a maiores períodos de uso e fluxos de oxigênio, sem associação à sobrevida. A substituição dos cilindros por concentradores poderá reduzir custos significativamente.

An Eulerian Immersed Boundary Method for flow simulations over stationary and moving rigid bodies

The fluid flow over bodies with complex geometry has been the subject of research of many scientists and widely explored experimentally and numerically. The present study proposes an Eulerian Immersed Boundary Method for flows simulations over stationary or moving rigid bodies. The proposed method allows the use of Cartesians Meshes. Here, two-dimensional simulations of fluid flow over stationary and oscillating circular cylinders were used for verification and validation. Four different cases were explored: the flow over a stationary cylinder, the flow over a cylinder oscillating in the flow direction, the flow over a cylinder oscillating in the normal flow direction, and a cylinder with angular oscillation. The time integration was carried out by a classical 4th order Runge-Kutta scheme, with a time step of the same order of distance between two consecutive points in x direction. High-order compact finite difference schemes were used to calculate spatial derivatives. The drag and lift coefficients, the lock-in phenomenon and vorticity contour plots were used for the verification and validation of the proposed method. The extension of the current method allowing the study of a body with different geometry and three-dimensional simulations is straightforward. The results obtained show a good agreement with both numerical and experimental results, encouraging the use of the proposed method.

Numerical investigation of the three-dimensional secondary instabilities in the time-developing compressible mixing layer

Mixing layers are present in very different types of physical situations such as atmospheric flows, aerodynamics and combustion. It is, therefore, a well researched subject, but there are aspects that require further studies. Here the instability of two-and three-dimensional perturbations in the compressible mixing layer was investigated by numerical simulations. In the numerical code, the derivatives were discretized using high-order compact finite-difference schemes. A stretching in the normal direction was implemented with both the objective of reducing the sound waves generated by the shear region and improving the resolution near the center. The compact schemes were modified to work with non-uniform grids. Numerical tests started with an analysis of the growth rate in the linear regime to verify the code implementation. Tests were also performed in the non-linear regime and it was possible to reproduce the vortex roll-up and pairing, both in two-and three-dimensional situations. Amplification rate analysis was also performed for the secondary instability of this flow. It was found that, for essentially incompressible flow, maximum growth rates occurred for a spanwise wavelength of approximately 2/3 of the streamwise spacing of the vortices. The result demonstrated the applicability of the theory developed by Pierrehumbet and Widnall. Compressibility effects were then considered and the maximum growth rates obtained for relatively high Mach numbers (typically under 0.8) were also presented.

Hydrodynamics and suspended sediment transport in the Camboriú estuary - Brazil: pre jetty conditions

Estuarine hydrodynamics is a key factor in the definition of the filtering capacity of an estuary and results from the interaction of the processes that control the inlet morphodynamics and those that are acting in the mixing of the water in the estuary. The hydrodynamics and suspended sediment transport in the Camboriú estuary were assessed by two field campaigns conducted in 1998 that covered both neap and spring tide conditions. The period measured represents the estuarine hydrodynamics and sediment transport prior to the construction of the jetty in 2003 and provides important background information for the Camboriú estuary. Each field campaign covered two complete tidal cycles with hourly measurements of currents, salinity, suspended sediment concentration and water level. Results show that the Camboriú estuary is partially mixed with the vertical structure varying as a function of the tidal range and tidal phase. The dynamic estuarine structure can be balanced between the stabilizing effects generated by the vertical density gradient, which produces buoyancy and stratification flows, and the turbulent effects generated by the vertical velocity gradient that generates vertical mixing. The main sediment source for the water column are the bottom sediments, periodically resuspended by the tidal currents. The advective salt and suspended sediment transport was different between neap and spring tides, being more complex at spring tide. The river discharge term was important under both tidal conditions. The tidal correlation term was also important, being dominant in the suspended sediment transport during the spring tide. The gravitational circulation and Stokes drift played a secondary role in the estuarine transport processes.

Utilização de métodos acústicos em estudos de dinâmica costeira: exemplo na desembocadura lagunar de cananéia

Em função de suas condições de interface entre águas doces e salinas, desembocaduras estuarinas e lagunares constituem sistemas geomorfológicos altamente complexos e dinâmicos. Como conseqüência da variabilidade espacial e temporal dos fluxos de maré, o leito responde com uma grande variabilidade nas características morfológicas e sedimentares. Neste sentido, é possível relacionar diretamente a circulação de fundo e o transporte sedimentar com as feições submersas geradas. Perfis de ecossondagem, sonar de varredura lateral e sísmica de alta resolução, executados na desembocadura lagunar de Cananéia, revelaram a existência de uma dinâmica de fundo extremamente complexa, caracterizada por marcas onduladas e ondas de areia de alturas métricas. As maiores ondas de areia, localizadas em uma depressão na desembocadura lagunar, apresentam inversão de polaridade em sua assimetria, com a presença de ondas simétricas de grande tamanho no ponto de inversão. Este padrão morfológico não apresenta variação temporal em escala anual, sugerindo a persistência de um padrão de fluxos sobre o leito. Esta dinâmica revela, também, a constância de fluxos convergentes que aparentemente independem das condições de maré enchente ou vazante. Os resultados permitiram o estabelecimento de um primeiro modelo qualitativo de circulação de fundo na área, com aplicações potenciais na navegação e estudos de proteção da costa.

An Extension of the Invariance Principle for a Class of Differential Equations with Finite Delay

An extension of the uniform invariance principle for ordinary differential equations with finite delay is developed. The uniform invariance principle allows the derivative of the auxiliary scalar function V to be positive in some bounded sets of the state space while the classical invariance principle assumes that. V <= 0. As a consequence, the uniform invariance principle can deal with a larger class of problems. The main difficulty to prove an invariance principle for functional differential equations is the fact that flows are defined on an infinite dimensional space and, in such spaces, bounded solutions may not be precompact. This difficulty is overcome by imposing the vector field taking bounded sets into bounded sets.

Possible states in the flow around two circular cylinders in tandem with separations in the vicinity of the drag inversion spacing

The possible states in the flow around two identical circular cylinders in tandem arrangements are investigated for configurations in the vicinity of the drag inversion separation. By means of numerical simulations, the hysteresis in the transition between the shedding regimes is studied and the relationship between (three-dimensional) secondary instabilities and shedding regime determination is addressed. The differences observed in the behavior of two- and three-dimensional flows are analyzed, and the regions of bistable flow are delimited. Very good agreement is found between the proposed scenario and results available in the literature. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3420111]

The viscosity parameter's dependence on the profile of the disk scale height

It is shown that, for accretion disks, the height scale is a constant whenever hydrostatic equilibrium and the subsonic turbulence regime hold in the disk. In order to have a variable height scale, processes are needed that contribute an extra term to the continuity equation. This contribution makes the viscosity parameter much greater in the outer region and much smaller in the inner region. Under these circumstances, turbulence is the presumable source of viscosity in the disk.