An unstructured CVFEM and moving interface algorithm for non-Newtonian Hele-Shaw flows in injection molding

Purpose - The purpose of this paper is to develop a novel unstructured simulation approach for injection molding processes described by the Hele-Shaw model. Design/methodology/approach - The scheme involves dual dynamic meshes with active and inactive cells determined from an initial background pointset. The quasi-static pressure solution in each timestep for this evolving unstructured mesh system is approximated using a control volume finite element method formulation coupled to a corresponding modified volume of fluid method. The flow is considered to be isothermal and non-Newtonian. Findings - Supporting numerical tests and performance studies for polystyrene described by Carreau, Cross, Ellis and Power-law fluid models are conducted. Results for the present method are shown to be comparable to those from other methods for both Newtonian fluid and polystyrene fluid injected in different mold geometries. Research limitations/implications - With respect to the methodology, the background pointset infers a mesh that is dynamically reconstructed here, and there are a number of efficiency issues and improvements that would be relevant to industrial applications. For instance, one can use the pointset to construct special bases and invoke a so-called ""meshless"" scheme using the basis. This would require some interesting strategies to deal with the dynamic point enrichment of the moving front that could benefit from the present front treatment strategy. There are also issues related to mass conservation and fill-time errors that might be addressed by introducing suitable projections. The general question of ""rate of convergence"" of these schemes requires analysis. Numerical results here suggest first-order accuracy and are consistent with the approximations made, but theoretical results are not available yet for these methods. Originality/value - This novel unstructured simulation approach involves dual meshes with active and inactive cells determined from an initial background pointset: local active dual patches are constructed ""on-the-fly"" for each ""active point"" to form a dynamic virtual mesh of active elements that evolves with the moving interface.

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.

Phantom accretion by black holes and the generalized second law of thermodynamics

The accretion of a phantom fluid with non-zero chemical potential by black holes is discussed with basis on the generalized second law of thermodynamics. For phantom fluids with positive temperature and negative chemical potential we demonstrate that the accretion process is possible, and that the condition guaranteeing the positiveness of the phantom fluid entropy coincides with the one required by the generalized second law. (C) 2010 Elsevier B.V. All rights reserved.

Rhyacian (2.23-2.20 Ga) juvenile accretion in the southern Sao Francisco craton, Brazil: Geochemical and isotopic evidence from the Serrinha magmatic suite, Mineiro belt

The Serrinha magmatic suite (Mineiro belt) crops out in the southern edge of the Sao Francisco craton, comprising the Brito quartz-diorite, Brumado de Cima and Brumado de Baixo granodiorites, granophyres and felsic sub-volcanic and volcanic rocks, part of which intruded into the Nazareno greenstone belt. The suite rocks have petrographic features that are consistent with magma supercooling due to the low water content combined with volatile loss, leading to crystallization of quartz and alkaline feldspar at the rims of plagioclase phenocrysts (granophyric intergrowth). The investigated rocks are sub-alkaline, calc-alkaline and show low content in rare earth elements. The U-Pb zircon crystallization ages for the Brumado de Cima granodiorite [2227 +/- 22 (23) Ma] and a coeval granophyre [2211 +/- 22 (23) Ma], coupled with available single-zircon Pb evaporation ages for the Brito and Brumado de Baixo plutons, are significantly older than the ""Minas orogeny"" (ca. 2100-2050 Ga) of Quadrilatero Ferrifero area, eastward from the Serrinha suite. Our data establish an early Rhyacian event tectonically linked with the evolution of the Mineiro belt. The bulk Nd isotopic signature [low negative to positive epsilon(Nd(t)) values] of the Serrinha samples are consistent with the important role of Paleoproterozoic mantle components in the magma genesis. The integrated geologic, geochemical and isotopic information suggests that Paleoproterozoic evolution of the Mineiro belt initiated in a passive continental margin basin with deposition of the Minas Supergroup at ca. 2500 Ma. This stage was succeeded by outboard rupture of the oceanic lithosphere with development and coalescence of progressively younger magmatic arcs during Rhyacian time. One of the earliest arcs formed the Serrinha suite. The tectonic collage of the Serrinha and Ritapolis (2190-2120 Ma) arcs produced the NE-SW Lenheiro shear zone, resulting in mylonitization and recrystallization of both the granitoid intrusions and host rocks. As a matter of fact juxtaposition of distinct magmatic units in age and origin took place along the Lenheiros structure in this sector of the Mineiro belt. (C) 2009 Elsevier Ltd. All rights reserved.

Sr and Nd isotopic geochemistry of the early ultramafic-mafic rocks of the Mako bimodal volcanic belt of the Kedougou-Kenieba inlier (Senegal)

The Mako bimodal volcanic belt of the Kedougou-Kenieba inlier is composed of volcanic basalts and peridotites interbedded by quartzites and limestones intruded by different generations of granitoids. The early volcanic episode of the belt is constituted of submarine basalts with peridotite similar to those of the oceanic abyssal plains. It is intruded by the Badon Kakadian TTG-granitic batholite dated around 2200 Ma. The second volcanic phase is constituted of basaltic, andesitic, and felsitic flows exhibit structures of aerial volcanic rocks. It is intruded by granites dated between 2160 and 2070 Ma. The general pattern of trace element variation of submarine volcanic rocks is consistent with those of basalts from oceanic plateaus which are the modern equivalent of the Archean greenstones belts. The Nd and Sr isotopic systematics typical of juvenile material indicates that the source of these igneous rocks is derived from a depleted mantle source. These results are consistent with the idea of a major accretion within the West African Craton occurring at about 2.1 Ga and corresponding to an important process of mantle-oceanic lithosphere differentiation.

Evolution of a rift basin dominated by subaerial deposits: The Guaritas Rift, Early Cambrian, Southern Brazil

Most existing models for the evolution of rift basins predict the development of deep-water depositional systems during the stage of greatest tectonic subsidence, when accommodation generation potentially outpaces sedimentation. Despite this, some rift basins do not present deep-water systems, instead being dominated by subaerial deposits. This paper focuses on one of these particular rift basins, the Cambrian Guaritas Rift, Southern Brazil, characterized by more than 1500 m of alluvial and aeolian strata deposited in a 50-km-wide basin. The deposits of the Guaritas Rift can be ascribed to four depositional systems: basin-border alluvial fans, bedload-dominated ephemeral rivers, mixed-load ephemeral rivers and aeolian dune fields. These four systems are in part coeval and in part succeed each other, forming three stages of basin evolution: (i) Rift Initiation to Early Rift Climax stage, (ii) Mid to Late Rift Climax stage, and (iii) Early Post-Rift stage. The first stage comprises most of the Guaritas Group and is characterized by homogeneous bed-load-dominated river deposits, which do not clearly record the evolution of subsidence rates. The onset of sedimentation of finer-grained deposits occurred as a consequence of a reactivation event that changed the outline of the basin and the distribution of the nearby highlands. This strongly suggests that the capture of the main river system to another depression decreased the sediment supply to the basin. The study of the Guaritas Rift indicates that rift basins in which the sediment supply exceeds the accommodation generation occur as a consequence of moderate subsidence combined with the capture of a major river system to the basin during the initial stages of basin evolution. In these basins, changes in the average discharge of the river system or tectonic modification of the drainage network may be the major control on the stratigraphic architecture. (c) 2009 Published by Elsevier B.V.

Tectono-sedimentary evolution of the Neoproterozoic BIF-bearing Jacadigo Group, SW-Brazil

The Jacadigo Group contains one of the largest sedimentary iron and associated manganese deposits of the Neoproterozoic. Despite its great relevance, no detailed sedimentological study concerning the unit has been carried out to date. Here we present detailed sedimentological data and interpretation on depositional systems, system tracts, external controls on basin evolution, basin configuration and regional tectonic setting of the Jacadigo Basin. Six depositional systems were recognized: (I) an alluvial fan system; (II) a siliciclastic lacustrine system; (III) a fan-delta system; (IV) a bedload-dominated river system; (V) an iron formation-dominated lacustrine or marine gulf system; and (VI) a rimmed carbonate platform system. The interpreted depositional systems are related to three tectonic system tracts. The first four depositional systems are mainly made of continental siliciclastics and refer to the rift initiation to early rift climax stage; the lake/gulf system corresponds to the mid to late rift climax stage and the carbonate platform represents the immediate to late post rift stage (Bocaina Formation deposits of the Ediacaran fossil-bearing Corumba Group). The spatial distribution of the depositional systems and associated paleocurrent patterns indicate a WNW-ESE orientation of the master fault zone related to the formation of the Jacadigo Basin. Thus, the iron formations of the Jacadigo Group were deposited in a starved waterbody related to maximum fault displacement and accommodation rates in a restricted continental rift basin. The Fe-Si-Mn source was probably related to hydrothermal plume activity that reached the basin through the fault system during maximum fault displacement phases. Our results also suggest a restricted tectono-sedimentary setting for the type section of the Puga Formation. The Jacadigo Group and the Puga Formation, usually interpreted as glacial deposits, are readdressed here as basin margin gravitational deposits with no necessary relation to glacial processes. (C) 2011 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.