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.

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.

Padrões espaciais na distribuição de abelhas Euglossina (Hymenoptera, Apidae) da região Neotropical

Abelhas das orquídeas (Apini, Euglossina) apresentam distribuição principalmente Neotropical, com cerca de 200 espécies e cinco gêneros descritos. Muitos levantamentos locais de fauna estão disponíveis na literatura, mas estudos comparativos sobre a composição e distribuição dos Euglossina são ainda escassos. O objetivo deste estudo é analisar os dados disponíveis de 29 assembleias a fim de entender os padrões gerais de distribuição espacial nas áreas amostradas ao longo do Neotrópico. Métodos de ordenação (DCA e NMDS) foram utilizados para descrever os agrupamentos de assembleias de acordo com as ocorrências de abelhas das orquídeas. As localidades de florestas da América Central e da Amazônia formaram grupos coesos em ambas as análises, enquanto as localidades de Mata Atlântica ficaram mais dispersas nos gráficos. Localidades na margem leste da Amazônia aparecem como áreas de transição características entre esta sub-região e a Mata Atlântica. As análises de variância entre o primeiro eixo da DCA e variáveis selecionadas apresentaram valores significantes quanto à influência dos gradientes de latitude, longitude e precipitação, bem como das sub-regiões biogeográficas nos agrupamentos das assembleias. O padrão geral encontrado é congruente com os padrões biogeográficos previamente propostos para a região Neotropical. Os resultados do DCA auxiliam ainda a identificar, de forma independente, os elementos das faunas de cada uma das formações vegetais estudadas.

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.

The Southern Brazilian shelf: general characteristics, quaternary evolution and sediment distribution

Extending from latitude 34ºS to 22ºS the Southern Brazilian shelf constitutes the only part of the Brazilian shelf with a subtropical to temperate environment. The studies on the different geological aspects of the area began in the 1960's and have recently been reassessed after studies related to the determination of the Economic Exclusive Zone. In terms of morphology, the Southern Brazilian shelf may be divided into three sectors, the São Paulo Bight, the Florianópolis-Mostardas Sector and the Rio Grande Cone, characterized by conspicuous differences in terms of geological determining factors, bathymetry, declivities and the presence of canyons and channels. Despite the existence of hundreds of radiocarbon datings the sea level changes curve of southern Brazil during the Last Glacial Cycle is still a matter of debate. A recent controversy on the Middle and late Holocene sea level changes curve raised the question of the amplitude of the oscillations which occurred in the period. Also, a few but relatively consistent radiocarbon datings suggest the occurrence of a high sea level during Isotope Stage 3. In terms of sedimentary cover the Southern Brazilian shelf exhibits a very strong hydrodynamic control, both latitudinal and bathymetrical. The sector southward from 25ºS is characterized by the influence of the plume of water carrying sediments originating from the Río de La Plata. Actually its presence is conspicuous up to 28ºS, with the area between this latitude and 25ºS constituting a transitional zone. In terms of bathymetry the outer shelf is marked by the "floor-polisher" effect of the Brazil Current, which is responsible for the maintenance of a relict facies in areas deeper than 100 meters.

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]

Immunomodulatory effects of recombinant BCG expressing pertussis toxin on TNF-alpha and IL-10 in a bladder cancer model

Background: Since successful treatment of superficial bladder cancer with BCG requires proper induction of Th1 immunity, we have developed a rBCG-S1PT strain that induced a stronger cellular immune response than BCG. This preclinical study was designed to compare the modulatory effects of BCG and rBCG-S1PT on bladder TNF-alpha and IL-10 expression and to evaluate antitumour activity. Methods: For Experiment I, the MB49 bladder cancer cell line was used in C57BL/6 mice. Chemical cauterization of the bladder was performed to promote intravesical tumor implantation. Mice were treated by intravesical instillation with BCG, rBCG-S1PT or PBS once a week for four weeks. After 35 days the bladders were removed and weighed. TNF-