Real time measurement and control of viscosity for extrusion processes using recycled materials

The aim of this paper is to develop a new generation of extruder control system for recycled materials which has ability to automatically maintain constant a polymer melt viscosity of mixed recycled polymers during extrusion, regardless of variations in the Melt Flow Index (MFI) of recycled mixed grade high density polyethylene (HDPE) feedstock. The variations in MFI are due to differences in the source of the recycled material used. The work describes how melt viscosity for specific extruder/die system is calculated in real time using the rheological properties of the materials, the pressure drop through the extruder die and the actual throughput measurements using a gravimetric loss-in-weight hopper feeder. A closed-loop controller is also developed to automatically regulate screw speed and barrel temperature profile to achieve constant viscosity and enable consistent processing of variable grade recycled HDPE materials. Such a system will improve processability of mixed MFI polymers may also reduce the risk of polymer melt degradation, reduce producing large volumes of scrap/waste and lead to improvement in product quality. The experimental results of real time viscosity measurement and control using a 38 mm single screw extruder with different recycled HDPEs with widely different MFIs are reported in this work.

The determination of optimal processing conditions for PLA and PLA/CaCO3 in twin screw extrusion using DoE and multivariate analysis

Biodegradable polymers, such as PLA (Polylactide), come from renewable resources like corn starch and if disposed of correctly, degrade and become harmless to the ecosystem making them attractive alternatives to petroleum based polymers. PLA in particular is used in a variety of applications including medical devices, food packaging and waste disposal packaging. However, the industry faces challenges in melt processing of PLA due to its poor thermal stability which is influenced by processing temperatures and shearing.
Identification and control of suitable processing conditions is extremely challenging, usually relying on trial and error, and often sensitive to batch to batch variations. Off-line assessment in a lab environment can result in high scrap rates, long lead times and lengthy and expensive process development. Scrap rates are typically in the region of 25-30% for medical grade PLA costing between €2000-€5000/kg.
Additives are used to enhance material properties such as mechanical properties and may also have a therapeutic role in the case of bioresorbable medical devices, for example the release of calcium from orthopaedic implants such as fixation screws promotes healing. Additives can also reduce the costs involved as less of the polymer resin is required.
This study investigates the scope for monitoring, modelling and optimising processing conditions for twin screw extrusion of PLA and PLA w/calcium carbonate to achieve desired material properties. A DAQ system has been constructed to gather data from a bespoke measurement die comprising melt temperature; pressure drop along the length of the die; and UV-Vis spectral data which is shown to correlate to filler dispersion. Trials were carried out under a range of processing conditions using a Design of Experiments approach and samples were tested for mechanical properties, degradation rate and the release rate of calcium. Relationships between recorded process data and material characterisation results are explored.

Caracterização de novas membranas de titanossilicatos microporosos por ensaios de permeabilidade

O interesse crescente das membranas inorgânicas deve-se à potencial aplicação em novas áreas de investigação e da indústria, e em alternativa a operações mais convencionais. Em particular, as membranas de titanossilicatos oferecem vantagens importantes sobre as de zeólitos, pois podem ser sintetizadas sem agentes estruturantes orgânicos, para evitar a calcinação subsequente usualmente responsável por defeitos irreversíveis, exibem novas possibilidades de substituição isomórfica da matriz, permitindo um ajuste mais fino das propriedades catalíticas e de adsorção, e são capazes de separar misturas com base em diferenças de afinidade e tamanho molecular (efeito de peneiro). Os objectivos principais deste trabalho foram: i) a caracterização dinâmica de membranas do tipo zeolítico sintetizadas no Laboratório Associado CICECO, realizando-se experiências de permeação com gases puros e misturas; ii) o desenvolvimento e validação de novos modelos para a transferência de massa multicomponente através de membranas porosas pela abordagem de Maxwell-Stefan, tendo em conta os mecanismos específicos encontrados, particularmente a contribuição por difusão superficial; e iii) a modelação dos pontos experimentais medidos, bem como dados compilados da literatura. De forma a realizar os ensaios de permeação, desenhou-se, montou-se e testou-se uma instalação experimental. Para gases puros, os objectivos principais foram a medição de permeâncias a temperatura constante, por variação da pressão transmembranar r ( ΔP ), e de permeâncias a temperatura programada, conduzidas a ΔP constante. Seguidamente, calcularam-se as selectividades ideais. Em relação a misturas, a determinação de selectividades reais requer as fracções molares no permeado e no retido. Na globalidade, estudaram-se três suportes diferentes (aço inoxidável e α − alumina) e dezanove membranas de AM-3, ETS-10, ZSM-5 e zeólito 4A, utilizando-se H2, He, N2, CO2, e O2. A primeira avaliação exploratória da qualidade das membranas foi feita permeando azoto à temperatura ambiente. Assim, permeâncias superiores a 10−6 mol/m2s.Pa evidenciavam defeitos grosseiros, levando-nos a efectuar cristalizações adicionais sobre as primeiras camadas. Este procedimento foi implementado com oito membranas. Um trabalho experimental mais detalhado foi conduzido com cinco membranas. Membranas com curvas permeância-temperatura ( Π −T ) decrescentes indicam tipicamente transporte viscoso e de Knudsen, i.e. meso e macrodefeitos. Por exemplo, a membrana nº 3 de AM-3 exibiu este comportamento com H2, He, N2 e CO2 puros. A contribuição de Knudsen foi confirmada pela relação linear encontrada entre as permeâncias e o inverso da raiz quadrada da massa molar. O mecanismo viscoso foi também identificado, pois as permeâncias eram inversamente proporcionais à viscosidade do gás ou, atendendo a equações do tipo de Chapman-Enskog, directamente proporcionais a 2 0.5 k d M (onde k d é o diâmetro cinético e M a massa molar). Um comportamento de permeação distinto observou-se com a membrana nº 5 de AM-3. As permeâncias registadas a temperatura programada eram aproximadamente constantes para o N2, CO2 e O2, enquanto com o H2 cresciam significativamente. Conjuntamente elas evidenciam a ocorrência de macro, meso e microdefeitos intercristalinos. O transporte gasoso activado através dos microporos compensa o impacto diminuidor dos meso e macroporos. Ao contrário do N2, CO2 e O2, o pequeno diâmetro do hidrogénio torna-lhe possível permear através dos microporos intracristalinos, o que lhe adiciona um mecanismo de transferência responsável por esse crescimento. No que respeita à difusão superficial, o sistema CO2/ZSM-5 pode ser tomado como um exemplo paradigmático. Uma vez que este zeólito adsorve o CO2, as permeâncias diminuem com o crescimento de ΔP , em virtude de as concentrações no sólido aumentarem de forma não linear e tenderem para a saturação. Os resultados contrastantes obtidos com azoto realçam ainda mais o mecanismo superficial, pois o N2 não é adsorvido e as permeâncias medidas são constantes. Globalmente, as selectividades ideais calculadas ( α* ) variam de cerca de 1 a 4.2. Este parâmetro foi também utilizado para discriminar as melhores membranas, uma vez que baixos valores de α* denotam o escoamento viscoso não-selectivo típico de macrodefeitos. Por exemplo, o H2/CO2 na membrana nº 3 de AM-3 apresentou α* = 3.6 − 4.2 para 40–120ºC, enquanto que na membrana nº 5 de AM-3 originou α* = 2.6 − 3.1. Estes resultados corroboraram as observações anteriores, segundo as quais a membrana nº 5 era melhor do que a nº 3. Alguns ensaios foram realizados com membranas saturadas com água para aumentar a selectividade: as medições mostraram claramente uma melhoria inicial seguida de uma redução consistente de α* com o aumento da temperatura, devido à remoção das moléculas de água responsáveis pela obstrução de alguns poros. Em relação às selectividades reais de misturas contendo hidrogénio, devem ser realizadas mais experiências e a quantificação do hidrogénio deve ser melhorada. No que concerne à modelação, novos factores termodinâmicos de Maxwell- Stefan foram derivados para as isotérmicas mono e multicomponente de Nitta, Langmuir-Freundlich e Toth, tendo sido testadas com dados de equilíbrio e de permeação da literatura. (É importante realçar que só estão publicadas equações para Langmuir e Dual-Site Langmuir de componentes puros e misturas). O procedimento de validação adoptado foi exigente: i) as isotérmicas multicomponente foram previstas a partir das de gás puro; ii) os parâmetros de difusão dos componentes puros foram ajustados a dados de permeação de cada gás; iii) depois, as difusividades cruzadas de Maxwell- Stefan foram estimadas pela relação de Vignes; finalmente, v) as novas equações foram testadas usando-se estes parâmetros, tendo sido capazes de estimar com sucesso fluxos binários. Paralelamente ao enfoque principal do trabalho, derivou-se um novo modelo para permuta iónica em materiais microporosos baseado nas equações de Maxwell-Stefan. Este foi validado com dados experimentais de remoção de Hg2+ e Cd2+ de soluções aquosas usando ETS-4. A sua capacidade preditiva foi também avaliada, sendo possível concluir que se comporta muito bem. Com efeito, conseguiram-se boas previsões com parâmetros optimizados a partir de conjuntos de dados independentes. Este comportamento pode ser atribuído aos princípios físicos sólidos da teoria de Maxwell-Stefan.

Project, modelling and simulation of air stripping columns for the treatment of groundwater contaminated with volatile organic compounds

Volatile organic compounds are a common source of groundwater contamination that can be easily removed by air stripping in columns with random packing and using a counter-current flow between the phases. This work proposes a new methodology for the column design for any particular type of packing and contaminant avoiding the necessity of a pre-defined diameter used in the classical approach. It also renders unnecessary the employment of the graphical Eckert generalized correlation for pressure drop estimates. The hydraulic features are previously chosen as a project criterion and only afterwards the mass transfer phenomena are incorporated, in opposition to conventional approach. The design procedure was translated into a convenient algorithm using C++ as programming language. A column was built in order to test the models used either in the design or in the simulation of the column performance. The experiments were fulfilled using a solution of chloroform in distilled water. Another model was built to simulate the operational performance of the column, both in steady state and in transient conditions. It consists in a system of two partial non linear differential equations (distributed parameters). Nevertheless, when flows are steady, the system became linear, although there is not an evident solution in analytical terms. In steady state the resulting system of ODE can be solved, allowing for the calculation of the concentration profile in both phases inside the column. In transient state the system of PDE was numerically solved by finite differences, after a previous linearization.

Practical software improve design and performance of a randomly packed air stripping column for removing volatile organic compounds from water

STRIPPING is a software application developed for the automatic design of a randomly packing column where the transfer of volatile organic compounds (VOCs) from water to air can be performed and to simulate it’s behaviour in a steady-state. This software completely purges any need of experimental work for the selection of diameter of the column, and allows a choice, a priori, of the most convenient hydraulic regime for this type of operation. It also allows the operator to choose the model used for the calculation of some parameters, namely between the Eckert/Robbins model and the Billet model for estimating the pressure drop of the gaseous phase, and between the Billet and Onda/Djebbar’s models for the mass transfer. Illustrations of the graphical interface offered are presented.

A New Approach in the Design of Air Stripping Columns for the Treatment of Groundwater Contaminated With Volatile Organic Compounds

Volatile organic compounds are a common source of groundwater contamination that can be easily removed by air stripping in columns with random packing and using a counter-current flow between the phases. This work proposes a new methodology for column design for any type of packing and contaminant which avoids the necessity of an arbitrary chosen diameter. It also avoids the employment of the usual graphical Eckert correlations for pressure drop. The hydraulic features are previously chosen as a project criterion. The design procedure was translated into a convenient algorithm in C++ language. A column was built in order to test the design, the theoretical steady-state and dynamic behaviour. The experiments were conducted using a solution of chloroform in distilled water. The results allowed for a correction in the theoretical global mass transfer coefficient previously estimated by the Onda correlations, which depend on several parameters that are not easy to control in experiments. For best describe the column behaviour in stationary and dynamic conditions, an original mathematical model was developed. It consists in a system of two partial non linear differential equations (distributed parameters). Nevertheless, when flows are steady, the system became linear, although there is not an evident solution in analytical terms. In steady state the resulting ODE can be solved by analytical methods, and in dynamic state the discretization of the PDE by finite differences allows for the overcoming of this difficulty. To estimate the contaminant concentrations in both phases in the column, a numerical algorithm was used. The high number of resulting algebraic equations and the impossibility of generating a recursive procedure did not allow the construction of a generalized programme. But an iterative procedure developed in an electronic worksheet allowed for the simulation. The solution is stable only for similar discretizations values. If different values for time/space discretization parameters are used, the solution easily becomes unstable. The system dynamic behaviour was simulated for the common liquid phase perturbations: step, impulse, rectangular pulse and sinusoidal. The final results do not configure strange or non-predictable behaviours.

Optimized venous return with a self-expanding cannula: from computational fluid dynamics to clinical application.

The Smart canula concept allows for collapsed cannula insertion, and self-expansion within a vein of the body. (A) Computational fluid dynamics, and (B) bovine experiments (76+/-3.8 kg) were performed for comparative analyses, prior to (C) the first clinical application. For an 18F access, a given flow of 4 l/min (A) resulted in a pressure drop of 49 mmHg for smart cannula versus 140 mmHg for control. The corresponding Reynolds numbers are 680 versus 1170, respectively. (B) For an access of 28F, the maximal flow for smart cannula was 5.8+/-0.5 l/min versus 4.0+/-0.1 l/min for standard (P<0.0001), for 24F 5.5+/-0.6 l/min versus 3.2+/-0.4 l/min (P<0.0001), and for 20F 4.1+/-0.3 l/min versus 1.6+/-0.3 l/min (P<0.0001). The flow obtained with the smart cannula was 270+/-45% (20F), 172+/-26% (24F), and 134+/-13% (28F) of standard (one-way ANOVA, P=0.014). (C) First clinical application (1.42 m2) with a smart cannula showed 3.55 l/min (100% predicted) without additional fluids. All three assessment steps confirm the superior performance of the smart cannula design.

In vivo testing of a novel adjustable glaucoma drainage device.

PURPOSE: We report on the in vivo testing of a novel noninvasively adjustable glaucoma drainage device (AGDD), which features an adjustable outflow resistance, and assess the safety and efficiency of this implant. METHODS: Under general anesthesia, the AGDD was implanted on seven white New Zealand rabbits for a duration of 4 months under a scleral flap in a way analogous to the Ex-PRESS device and set in an operationally closed position. The IOP was measured on a regular basis on the operated and control eyes using a rebound tonometer. Once a month the AGDD was adjusted noninvasively from its fully closed to its fully open position and the resulting pressure drop was measured. The contralateral eye was not operated and served as control. After euthanization, the eyes were collected for histology evaluation. RESULTS: The mean preoperative IOP was 11.1 ± 2.4 mm Hg. The IOP was significantly lower for the operated eye (6.8 ± 2 mm Hg) compared to the nonoperated eye (13.1 ± 1.6 mm Hg) during the first 8 days after surgery. When opening the AGDD from its fully closed to fully open position, the IOP dropped significantly from 11.2 ± 2.9 to 4.8 ± 0.9 mm Hg (P < 0.05). CONCLUSIONS: Implanting the AGDD is a safe and uncomplicated surgical procedure. The fluidic resistance was noninvasively adjustable during the postoperative period with the AGDD between its fully closed and fully open positions.

Hydrodynamic Study of Swirling Fluidized Bed and the Role of Distributor

The study aims to the hydrodynamic characteristics of swirling fluidized bed, using large particles (Geldart D-type) selected from locally available agricultural produce (coffee beans and black pepper). The important variables considered in the present study include percentage area of opening, angle of air injection and the percentage useful area of the distributor. A total of seven distributors have been designed and fabricated for a bed column of 300 mm, namely single row vane type distributors (15˚ and 20˚ vane angle), inclined hole type distributors (15˚ and 20˚ vane angle) and perforated plate distributors. The useful area of distributor of single row vane type, three now vane-type and inclined hole-type distributors are respectively 64%,91% and 94%. The hydrodynamic parameters considered in the present study include distributor pressure drop, air velocity, minimum fluidizing velocity, bed pressure drop, bed height and the bed behaviour. It has been observed that, in general, the distributor pressure drop decreases with an increase in the percentage area of opening, Further, and increase in the area of opening above 17% will not considerably reduce the distributor pressure drop. In the present study, for the distributor with an area of opening 17%, and corresponding to the maximum measured superficial velocity of 4.33 m/s, the distributor pressure drop obtained was 55.25mm of water. The study on the bed behavior revealed that, in a swirling fluidized bed, once swirl motion starts, the bed pressure drop increases with superficial velocity in the outer region and it decreases in the inner region. This means that, with higher superficial velocity, the air might get by-passed through the inner boundary of the bed (around the cone). So, depending on the process for which the bed is used, the maximum superficial velocity is to be limited to have an optimum bed performance.

Influence of Angle of Air Injection and Particles in Bed Hydrodynamics of Swirling Fluidized Bed

The thesis presented here unveils an experimental study of the hydrodynamic characteristics of swirling fluidized bed viz. pressure drop across the distributor and the bed, minimum fluidizing velocity, bed behaviour and angle of air injection. In swirling fluidized bed the air is admitted to the bed at an angle 'Ѳ' to the horizontal. The vertical component of the velocity v sin Ѳ causes fluidization and the horizontal component v cos Ѳ contributes to swirl motion of the bed material.The study was conducted using spherical particles having sizes 3.2 mm, 5.5 mm & 7.4 mm as the bed materials. Each of these particles was made from high density polyethylene, nylon and acetal having relative densities of 0.93, 1.05 and 1.47 respectively.The experiments were conducted using conidour type distributors having four rows of slits. Altogether four distributors having angles of air injection (Φ)- 0°, 5°, 10° & 15° were designed and fabricated for the study. The total number of slits in each distributor was 144. The area of opening was 6220 mm2 making the percentage area of opening to 9.17. But the percentage useful area of opening of the distributor was 96.The experiments on the variation of distributor pressure drop with superficial velocity revealed that the distributor pressure drop decreases with angle of air injection. Investigations related to bed hydrodynamics were conducted using 2.5 kg of bed material. The bed pressure drop measurements were made along the radial direction of the distributor at distances of 60 mm, 90 mm, 120 mm & 150 mm from the centre of the distributor. It was noticed that after attaining minimum fluidizing velocity, the bed pressure drop increases along the radial direction of the distributor. But at a radial distance of 90 mm from the distributor centre, after attaining minimum fluidizing velocity the bed pressure drop remains almost constant. It was also observed that the bed pressure drop varies inversely with particle size as well as particle density.An attempt was made to determine the effect of various parameters on minimum fluidizing velocity. It was noticed that the minimum fluidizing velocity varies directly with angle of air injection (Φ), particle size and particle density.The study on the bed behaviour showed that the superficial velocity required for initiating various bed phenomena (such as swirl motion and separation of particles from the cone at the centre) increase with increase in particle size as well as particle density. It was also observed that the particle size and particle density directly influence the superficial velocity required for various regimes of bed behaviour such as linear variation of bed pressure drop, constant bed pressure drop and sudden increase or decrease in bed pressure drop.Experiments were also performed to study the effect of angle of air injection (Φ). It was noticed that the bed pressure drop decreases with angle of air injection. It was also noticed that the angle of air injection directly influence the superficial velocity required for initiating various bed phenomena as well as the various regimes of bed behaviour.

Qualitative air flow modelling and analysis of data centre air conditioning as multiple jet array

This work identifies the importance of plenum pressure on the performance of the data centre. The present methodology followed in the industry considers the pressure drop across the tile as a dependant variable, but it is shown in this work that this is the only one independent variable that is responsible for the entire flow dynamics in the data centre, and any design or assessment procedure must consider the pressure difference across the tile as the primary independent variable. This concept is further explained by the studies on the effect of dampers on the flow characteristics. The dampers have found to introduce an additional pressure drop there by reducing the effective pressure drop across the tile. The effect of damper is to change the flow both in quantitative and qualitative aspects. But the effect of damper on the flow in the quantitative aspect is only considered while using the damper as an aid for capacity control. Results from the present study suggest that the use of dampers must be avoided in data centre and well designed tiles which give required flow rates must be used in the appropriate locations. In the present study the effect of hot air recirculation is studied with suitable assumptions. It identifies that, the pressure drop across the tile is a dominant parameter which governs the recirculation. The rack suction pressure of the hardware along with the pressure drop across the tile determines the point of recirculation in the cold aisle. The positioning of hardware in the racks play an important role in controlling the recirculation point. The present study is thus helpful in the design of data centre air flow, based on the theory of jets. The air flow can be modelled both quantitatively and qualitatively based on the results.

“Qualitative Air Flow Modelling And Analysis Of Data Centre Air Conditioning As Multiple Jet Array

In the present study the effect of hot air recirculation is studied with suitable assumptions. It identifies that, the pressure drop across the tile is a dominant parameter which governs the recirculation. The rack suction pressure of the hardware along with the pressure drop across the tile determines the point of recirculation in the cold aisle. The positioning of hardware in the racks play an important role in controlling the recirculation point. The present study is thus helpful in the design of data centre air flow, based on the theory of jets. The air flow can be modelled both quantitatively and qualitatively based on the results

The Inertio-Elastic Planar Entry Flow of Low-Viscosity Elastic Fluids in Micro-fabricated Geometries

The non-Newtonian flow of dilute aqueous polyethylene oxide (PEO) solutions through microfabricated planar abrupt contraction-expansions is investigated. The contraction geometries are fabricated from a high-resolution chrome mask and cross-linked PDMS gels using the tools of soft-lithography. The small length scales and high deformation rates in the contraction throat lead to significant extensional flow effects even with dilute polymer solutions having time constants on the order of milliseconds. The dimensionless extra pressure drop across the contraction increases by more than 200% and is accompanied by significant upstream vortex growth. Streak photography and videomicroscopy using epifluorescent particles shows that the flow ultimately becomes unstable and three-dimensional. The moderate Reynolds numbers (0.03 ≤ Re ≤ 44) associated with these high Deborah number (0 ≤ De ≤ 600) microfluidic flows results in the exploration of new regions of the Re-De parameter space in which the effects of both elasticity and inertia can be observed. Understanding such interactions will be increasingly important in microfluidic applications involving complex fluids and can best be interpreted in terms of the elasticity number, El = De/Re, which is independent of the flow kinematics and depends only on the fluid rheology and the characteristic size of the device.

The competing roles of extensional viscosity and normal stress differences in complex flows of elastic liquids

In various attempts to relate the behaviour of highly-elastic liquids in complex flows to their rheometrical behaviour, obvious candidates for study have been the variation of shear viscosity with shear rate, the two normal stress differences N(1) and N(2) especially N(1), and the extensional viscosity eta(E). In this paper, we shall be mainly interested in `constant-viscosity` Boger fluids, and, accordingly, we shall limit attention to N(1) and eta(E). We shall concentrate on two important flows - axisymmetric contraction flow and ""splashing"" (particularly that which arises when a liquid drop falls onto the free Surface of the same liquid). Modem numerical techniques are employed to provide the theoretical predictions. It is shown that the two obvious manifestations of viscoelastic rheometrical behaviour can sometimes be opposing influences in determining flow characteristics. Specifically, in an axisymmetric contraction flow, high eta(E) , can retard the flow, whereas high N(1) can have the opposite effect. In the splashing experiment, high eta(E) can certainly reduce the height of the so-called Worthington jet, thus confirming some early suggestions, but, again, other rheometrical influences can also have a role to play and the overall picture may not be as clear as it was once envisaged.

P-T-Fluid evolution and graphite deposition during retrograde metamorphism in Ribeira Fold Belt, SE Brazil: Oxygen fugacity, fluid inclusions and C-O-H isotopic evidence

Combined fluid inclusion (FI) microthermometry, Raman spectroscopy, X-ray diffraction, C-O-H isotopes and oxygen fugacities of granulites from central Ribeira Fold Belt, SE Brazil, provided the following results: i) Magnetite-Hematite fO(2) estimates range from 10(-11.5) bar (QFM + 1) to 10(-18.3) bar (QFM - 1) for the temperature range of 896 degrees C-656 degrees C, implying fO(2) decrease from metamorphic peak temperatures to retrograde conditions; ii) 5 main types of fluid inclusions were observed: a) CO(2) and CO(2)-N(2) (0-11 mol%) high to medium density (1.01-0.59 g/cm(3)) FI; b) CO(2) and CO(2)-N(2) (0-36 mol%) low density (0.19-0.29 g/cm(3)) FI; c) CO(2) (94-95 mol%)-N(2) (3 mol%)-CH(4) (2-3 mol%)-H(2)O (water phi(v) (25 degrees C) = 0.1) FI; d) low-salinity H(2)O-CO(2) FI; and e) late low-salinity H(2)O FI; iii) Raman analyses evidence two graphite types in khondalites: an early highly ordered graphite (T similar to 450 degrees C) overgrown by a disordered kind (T similar to 330 degrees C); iv) delta(18)O quartz results of 10.3-10.7 parts per thousand, imply high-temperature CO(2) delta(18)O values of 14.4-14.8 parts per thousand, suggesting the involvement of a metamorphic fluid, whereas lower temperature biotite delta(18)O and delta D results of 7.5-8.5 parts per thousand and -54 to -67 parts per thousand respectively imply H(2)O delta(18)O values of 10-11 parts per thousand and delta D(H2O) of -23 to -36 parts per thousand suggesting delta(18)O depletion and increasing fluid/rock ratio from metamorphic peak to retrograde conditions. Isotopic results are compatible with low-temperature H(2)O influx and fO(2) decrease that promoted graphite deposition in retrograde granulites, simultaneous with low density CO(2), CO(2)-N(2) and CO(2)-N(2)-CH(4)-H(2)O fluid inclusions at T = 450-330 degrees C. Graphite delta(13)C results of -10.9 to -11.4 parts per thousand imply CO(2) delta(13)C values of -0.8 to -1.3 parts per thousand suggesting decarbonation of Cambrian marine carbonates with small admixture of lighter biogenic or mantle derived fluids. Based on these results, it is suggested that metamorphic fluids from the central segment of Ribeira Fold Belt evolved to CO(2)-N(2) fluids during granulitic metamorphism at high fO(2), followed by rapid pressure drop at T similar to 400-450 degrees C during late exhumation that caused fO(2) reduction induced by temperature decrease and water influx, turning carbonic fluids into CO(2)-H(2)O (depleting biotite delta(18)O and delta D values), and progressively into H(2)O. When fO(2) decreased substantially by mixture of carbonic and aqueous fluids, graphite deposited forming khondalites. (C) 2010 Elsevier Ltd. All rights reserved.