Quantification of magmatic and hydrothermal processes in a peralkaline syenite-alkali granite complex based on textures, phase equilibria, and stable and radiogenic isotopes

The Puklen complex of the Mid-Proterozoic Gardar Province, South Greenland, consists of various silica-saturated to quartz-bearing syenites, which are intruded by a peralkaline granite. The primary mafic minerals in the syenites are augite +/- olivine + Fe-Ti oxide + amphibole. Ternary feldspar thermometry and phase equilibria among mafic silicates yield T = 950-750degreesC, a(SiO2) = 0.7-1 and an f(O2) of 1-3 log units below the fayalite-magnetite-quartz (FMQ) buffer at 1 kbar. In the granites, the primary mafic minerals are ilmenite and Li-bearing arfvedsonite, which crystallized at temperatures below 750degreesC and at f(O2) values around the FMQ buffer. In both rock types, a secondary post-magmatic assemblage overprints the primary magmatic phases. In syenites, primary Ca-bearing minerals are replaced by Na-rich minerals such as aegirine-augite and albite, resulting in the release of Ca. Accordingly, secondary minerals include ferro-actinolite, (calcite-siderite)(ss), titanite and andradite in equilibrium with the Na-rich minerals. Phase equilibria indicate that formation of these minerals took place over a long temperature interval from near-magmatic temperatures down to similar to300degreesC. In the course of this cooling, oxygen fugacity rose in most samples. For example, late-stage aegirine in granites formed at the expense of arfvedsonite at temperatures below 300degreesC and at an oxygen fugacity above the haematite-magnetite (HM) buffer. The calculated delta(18)O(melt) value for the syenites (+5.9 to +6.3parts per thousand) implies a mantle origin, whereas the inferred delta(18)O(melt) value of <+5.1parts per thousand for the granitic melts is significantly lower. Thus, the granites require an additional low-delta(18)O contaminant, which was not involved in the genesis of the syenites. Rb/Sr data for minerals of both rock types indicate open-system behaviour for Rb and Sr during post-magmatic metasomatism. Neodymium isotope compositions (epsilonNd(1170 Ma) = -3.8 to -6.4) of primary minerals in syenites are highly variable, and suggest that assimilation of crustal rocks occurred to variable extents. Homogeneous epsilon(Nd) values of -5.9 and -6.0 for magmatic amphibole in the granites lie within the range of the syenites. Because of the very similar neodymium isotopic compositions of magmatic and late- to post-magmatic minerals from the same syenite samples a principally closed-system behaviour during cooling is implied. In contrast, for the granites an externally derived fluid phase is required to explain the extremely low epsilon(Nd) values of about -10 and low delta(18)O between +2.0 and +0.5parts per thousand for late-stage aegirine, indicating an open system in the late-stage history. In this study we show that the combination of phase equilibria constraints with stable and radiogenic isotope data on mineral separates can provide much better constraints on magma evolution during emplacement and crystallization than conventional whole-rock studies.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Epidote in hydrous mafic magmas : near solidus phase relations in the lower crust

Crystallisation of hydrous mafic magmas at high pressure is a subject of numerous petrologic and experimental studies since the last century and is mainly related to the process of continental crust formation and the possible link between mantle derived melts and low pressure granitoids. Albeit the sequence of crystallization is well constrained by experimental studies, the origin of exposed lower crustal rocks exposed on the earth surface is controversial. Ones line of argument is favouring high pressure crystallization of dry or wet mafic magmas, whereas others invoke partial melting of pre-existing crust. Therefore studies involving field, textural and chemical observations of exposed lower crust such as in Kohistan (Pakistan) or Talkeetna (Alaska) are crucial to understand the continental crust formation processes via arc magmatism. Epidote-bearing gabbros are very sparse and always associated with the deep part of continental crust (>30 km) as in the Kohistan Arc Complex (Pakistan) or in the Chelan Complex (western U.S.). Magmatic epidote is restricted to a small temperature interval above the water-saturated solidus of MORB and represent the last crystallizing liquids in lower crustal regions. However, epidote and melt stability at lower crustal pressures are not clearly established.¦The Chelan complex (western U.S.) at the base of the Cascadian Arc is composed mainly by peraluminous tonalité associated with gabbroic and ultramafic rocks and was traditionally interpreted as a migmatitic terrain. However field, chemical and mineralogical observations rather suggest a magmatic origin and point to a protracted crystallization at intermediate to high pressure ~ 1.0 GPa dominated by amphibole fractionation and followed by isobaric cooling down to 650°C. Crystal fractionation modelling using whole rock composition and field constraints is able to generate peraluminous tonalité. The modelled crystallisation sequence and the volume proportions are in agreement with experimental studies performed at these pressures. The Chelan complex was thus not formed during a partial melting event, but represent the sequence of crystallisation occurring at the base of the crust. Massive fractionation of hornblende is able to generate peraluminous tonalité without significant assimilation of crustal rock.¦Similarly to the Chelan complex, the base of the Kohistan arc is composed of cumulates derived by high pressure crystallization of hydrous magma. In garnet gabbros, epidote occurs as magmatic phase, crystallising from hydrous interstitial melt trapped between grain boundaries at lower crustal pressures (Ρ ~ 1.2 GPa) for temperature of (650-700 °C). Trace and REE signature in epidote indicate that epidote was formed through peritectic reaction involving garnet, clinopyroxene and plagioclase. At the beginning of the crystallisation epidote signature is dominated by REE content in the melt, whereas at the end the signature is dominated by reacting phases. Melt in equilibrium with epidote inferred from the partition coefficients available is similar to intrusive tonalité up the section indicating that hydrous melt was extracted from the garnet gabbros. In some gabbros epidote shows single homogeneous compositions, while in others coexisting epidote have different compositions indicating the presence of solvi along the Al-Fe3+ join. The overgowths are only observed in presence of paragonite in the assemblage, suggesting high water content. At high water content, the hydrous solidus is shift to lower temperature and probably intersects the solvi observed along the Al-Fe3+ join. Therefore, several compositions of epidote is stable at high water content.¦-¦La composition chimique de la croûte continentale est considérée comme similaire à celle du magmatisme calco-alcalin de marge continentale active (enrichissement en éléments mobiles dans les fluides, anomalies négatives en Nb, Ta et éléments à haut potentiel électronique, etc...). Cependant la nature andésitique de la croûte continentale (Si02 > 60 wt%), résultant des nombreuses intrusions de granitoïdes dans la croûte supérieure, est sujette à polémique et le lien entre les magmas dérivés du manteau et les roches évoluées de faible profondeur n'est pas clairement établi (fusion partielle de croûte basaltique, cristallisation fractionnée à haute pression, etc...).¦Les affleurements de croûte profonde sont rares mais précieux, car ils permettent d'observer les phénomènes se passant à grande profondeur. Le complexe de Chelan (Washington Cascades) en est un exemple. Formé à environ 30 km de profondeur, il est composé de roches gabbroïques et ultramafiques, ainsi que de tonalités, qui furent souvent interprétés comme le produit de la fusion partielle de la croûte. Cependant, les relations de terrain, la chimie des éléments majeurs et des éléments traces sont cohérentes avec l'évolution d'un complexe magmatique mafique dans la croûte profonde ou moyenne ( 1.0 GPa), dominée par le fractionnement de l'amphibole. Après son emplacement, le complexe a subi un refroidissement isobare jusqu'à des températures de l'ordre de 650 °C, déduit de la composition chimique des minéraux. Un bilan de masse contraint pax les observations de terrain permet de calculer la séquence et les volumes de fractionnement. Les faciès évolués légèrement hyperalumineux observés sur le terrain peuvent être générés par la cristallisation de 3 % de websterite à olivine, 12 % d'hornblendite à pyroxène 33 % d'hornblendite, 19 % de gabbros, 15 % de diorite et 2 % de tonalité. Nous montrons ainsi qu'une série de fractionnement contrôlée par l'amphibole permet de générer des tonalités sans assimilation de matériel crustal et l'exemple de Chelan illustre la viabilité de ce processus dans la formation de croûte continentale.¦Les réactions proches du solidus saturé en H20 dans les systèmes basaltiques à des pressions élevées restent énigmatiques. Diverses expériences tendent à montrer que l'épidote est stable dans ces conditions, mais rarement observée (décrite ?) comme phase primaire dans les systèmes naturels. Les épidotes trouvées dans les gabbros de Jijal (nord-Pakistan) montrent des textures de type .magmatique telles qu'observées dans les roches évoluées. Le contenu en terres rares de ces épidotes est très variable allant de signatures enrichies en terres rares légères impliquant la présence de liquide interstitiel à des signatures complètement déprimées en ces mêmes éléments, évoquant une cristallisation en coexistence avec du grenat. Ces diverses signatures reflètent un chemin de cristallisation en présence de liquide interstitiel et enregistrent des réactions péritectiques impliquant grenat, clinopyroxene et plagioclase à des pressions de ~ 1.2 GPa pour des températures de 650-700 °C. Cependant dans quelques échantillons deux ou trois compositions d'épidotes coexistent démontrant la présence de lacunes d'immiscibilité le long de la solution solide épidote-clinozoïsite. La forte teneur en H20 du liquide magmatique est certainement à l'origine de la coexistence de deux compositions distinctes.

Towards desired crystalline product properties: In-situ monitoring of batch crystallization

The objective of industrial crystallization is to obtain a crystalline product which has the desired crystal size distribution, mean crystal size, crystal shape, purity, polymorphic and pseudopolymorphic form. Effective control of the product quality requires an understanding of the thermodynamics of the crystallizing system and the effects of operation parameters on the crystalline product properties. Therefore, obtaining reliable in-line information about crystal properties and supersaturation, which is the driving force of crystallization, would be very advantageous. Advanced techniques, such asRaman spectroscopy, attenuated total reflection Fourier transform infrared (ATR FTIR) spectroscopy, and in-line imaging techniques, offer great potential for obtaining reliable information during crystallization, and thus giving a better understanding of the fundamental mechanisms (nucleation and crystal growth) involved. In the present work, the relative stability of anhydrate and dihydrate carbamazepine in mixed solvents containing water and ethanol were investigated. The kinetics of the solvent mediated phase transformation of the anhydrate to hydrate in the mixed solvents was studied using an in-line Raman immersion probe. The effects of the operation parameters in terms of solvent composition, temperature and the use of certain additives on the phase transformation kineticswere explored. Comparison of the off-line measured solute concentration and the solid-phase composition measured by in-line Raman spectroscopy allowedthe identification of the fundamental processes during the phase transformation. The effects of thermodynamic and kinetic factors on the anhydrate/hydrate phase of carbamazepine crystals during cooling crystallization were also investigated. The effect of certain additives on the batch cooling crystallization of potassium dihydrogen phosphate (KDP) wasinvestigated. The crystal growth rate of a certain crystal face was determined from images taken with an in-line video microscope. An in-line image processing method was developed to characterize the size and shape of thecrystals. An ATR FTIR and a laser reflection particle size analyzer were used to study the effects of cooling modes and seeding parameters onthe final crystal size distribution of an organic compound C15. Based on the obtained results, an operation condition was proposed which gives improved product property in terms of increased mean crystal size and narrowersize distribution.

Phase behavior and rheological analysis of reverse liquid crystals and W/I2, W/H2 gel emulsions using an amphiphilic block copolymer.

This article reports the phase behavior determi- nation of a system forming reverse liquid crystals and the formation of novel disperse systems in the two-phase region. The studied system is formed by water, cyclohexane, and Pluronic L-121, an amphiphilic block copolymer considered of special interest due to its aggregation and structural proper- ties. This system forms reverse cubic (I2) and reverse hexagonal (H2) phases at high polymer concentrations. These reverse phases are of particular interest since in the two-phase region, stable high internal phase reverse emulsions can be formed. The characterization of the I2 and H2 phases and of the derived gel emulsions was performed with small-angle X-ray scattering (SAXS) and rheometry, and the influence of temperature and water content was studied. TheH2 phase experimented a thermal transition to an I2 phase when temperature was increased, which presented an Fd3m structure. All samples showed a strong shear thinning behavior from low shear rates. The elasticmodulus (G0) in the I2 phase was around 1 order of magnitude higher than in theH2 phase. G0 was predominantly higher than the viscousmodulus (G00). In the gel emulsions,G0 was nearly frequency-independent, indicating their gel type nature. Contrarily to water-in-oil (W/O) normal emulsions, in W/I2 and W/H2 gel emulsions, G0, the complex viscosity (|η*|), and the yield stress (τ0) decreased with increasing water content, since the highly viscous microstructure of the con- tinuous phase was responsible for the high viscosity and elastic behavior of the emulsions, instead of the volumefraction of dispersed phase and droplet size. A rheological analysis, in which the cooperative flow theory, the soft glass rheology model, and the slip plane model were analyzed and compared, was performed to obtain one single model that could describe the non-Maxwellian behavior of both reverse phases and highly concentrated emulsions and to characterize their microstructure with the rheological properties.

The influence of crystallization conditions on the filtration characteristics of sulphathiazole suspensions

Cooling crystallization is one of the most important purification and separation techniques in the chemical and pharmaceutical industry. The product of the cooling crystallization process is always a suspension that contains both the mother liquor and the product crystals, and therefore the first process step following crystallization is usually solid-liquid separation. The properties of the produced crystals, such as their size and shape, can be affected by modifying the conditions during the crystallization process. The filtration characteristics of solid/liquid suspensions, on the other hand, are strongly influenced by the particle properties, as well as the properties of the liquid phase. It is thus obvious that the effect of the changes made to the crystallization parameters can also be seen in the course of the filtration process. Although the relationship between crystallization and filtration is widely recognized, the number of publications where these unit operations have been considered in the same context seems to be surprisingly small. This thesis explores the influence of different crystallization parameters in an unseeded batch cooling crystallization process on the external appearance of the product crystals and on the pressure filtration characteristics of the obtained product suspensions. Crystallization experiments are performed by crystallizing sulphathiazole (C9H9N3O2S2), which is a wellknown antibiotic agent, from different mixtures of water and n-propanol in an unseeded batch crystallizer. The different crystallization parameters that are studied are the composition of the solvent, the cooling rate during the crystallization experiments carried out by using a constant cooling rate throughout the whole batch, the cooling profile, as well as the mixing intensity during the batch. The obtained crystals are characterized by using an automated image analyzer and the crystals are separated from the solvent through constant pressure batch filtration experiments. Separation characteristics of the suspensions are described by means of average specific cake resistance and average filter cake porosity, and the compressibilities of the cakes are also determined. The results show that fairly large differences can be observed between the size and shape of the crystals, and it is also shown experimentally that the changes in the crystal size and shape have a direct impact on the pressure filtration characteristics of the crystal suspensions. The experimental results are utilized to create a procedure that can be used for estimating the filtration characteristics of solid-liquid suspensions according to the particle size and shape data obtained by image analysis. Multilinear partial least squares regression (N-PLS) models are created between the filtration parameters and the particle size and shape data, and the results presented in this thesis show that relatively obvious correlations can be detected with the obtained models.

Effect of hydrodynamics on modelling, monitoring and control of crystallization

Crystallization is a purification method used to obtain crystalline product of a certain crystal size. It is one of the oldest industrial unit processes and commonly used in modern industry due to its good purification capability from rather impure solutions with reasonably low energy consumption. However, the process is extremely challenging to model and control because it involves inhomogeneous mixing and many simultaneous phenomena such as nucleation, crystal growth and agglomeration. All these phenomena are dependent on supersaturation, i.e. the difference between actual liquid phase concentration and solubility. Homogeneous mass and heat transfer in the crystallizer would greatly simplify modelling and control of crystallization processes, such conditions are, however, not the reality, especially in industrial scale processes. Consequently, the hydrodynamics of crystallizers, i.e. the combination of mixing, feed and product removal flows, and recycling of the suspension, needs to be thoroughly investigated. Understanding of hydrodynamics is important in crystallization, especially inlargerscale equipment where uniform flow conditions are difficult to attain. It is also important to understand different size scales of mixing; micro-, meso- and macromixing. Fast processes, like nucleation and chemical reactions, are typically highly dependent on micro- and mesomixing but macromixing, which equalizes the concentrations of all the species within the entire crystallizer, cannot be disregarded. This study investigates the influence of hydrodynamics on crystallization processes. Modelling of crystallizers with the mixed suspension mixed product removal (MSMPR) theory (ideal mixing), computational fluid dynamics (CFD), and a compartmental multiblock model is compared. The importance of proper verification of CFD and multiblock models is demonstrated. In addition, the influence of different hydrodynamic conditions on reactive crystallization process control is studied. Finally, the effect of extreme local supersaturation is studied using power ultrasound to initiate nucleation. The present work shows that mixing and chemical feeding conditions clearly affect induction time and cluster formation, nucleation, growth kinetics, and agglomeration. Consequently, the properties of crystalline end products, e.g. crystal size and crystal habit, can be influenced by management of mixing and feeding conditions. Impurities may have varying impacts on crystallization processes. As an example, manganese ions were shown to replace magnesium ions in the crystal lattice of magnesium sulphate heptahydrate, increasing the crystal growth rate significantly, whereas sodium ions showed no interaction at all. Modelling of continuous crystallization based on MSMPR theory showed that the model is feasible in a small laboratoryscale crystallizer, whereas in larger pilot- and industrial-scale crystallizers hydrodynamic effects should be taken into account. For that reason, CFD and multiblock modelling are shown to be effective tools for modelling crystallization with inhomogeneous mixing. The present work shows also that selection of the measurement point, or points in the case of multiprobe systems, is crucial when process analytical technology (PAT) is used to control larger scale crystallization. The thesis concludes by describing how control of local supersaturation by highly localized ultrasound was successfully applied to induce nucleation and to control polymorphism in reactive crystallization of L-glutamic acid.

Crystallization characteristics of bioactive glasses

Bioactive glasses are excellent candidates for implant materials, because they can form a chemical bond to bone or guide bone growth, depending on the glass composition. Some compositions have even shown soft tissue attachment and antimicrobial effects. So far, most clinical applications are based on monoliths, plates and particulates of different grain sizes. There is a growing interest in special products such as porous implants sintered from microspheres and fibers drawn from preforms or glass melts. The viscosity range at which these are formed coincides with the crystallization temperature range for most bioactive glasses, thus complicating the manufacturing process. In this work, the crystallization tendency and its kinetics for a series of glasses with their compositions within the range of bioactivity were investigated. The factors affecting crystallization and how it is related to composition were studied by means of thermal analysis and hot stage microscopy. The crystal compositions formed during isothermal and non-isothermal heat treatments were analyzed with SEM-EDXA and X-ray diffraction analysis. The temperatures at which sintering and fiber drawing can take place without interfering with crystallization were determined and glass compositions which are suitable for these purposes were established. The bioactivity of glass fibers and partly crystallized glass plates was studied by soaking them in simulated body fluid (SBF). The thickness of silica, calcium and phosphate rich reaction layers on the glass surface after soaking was used as an indication of the bioactivity. The results indicated that the crystallization tendencies of the experimental glasses are strongly dependent on composition. The main factor affecting the crystallization was found to be the alkali oxide content: the higher the alkali oxide content the lower the crystallization temperature. The primary crystalline phase formed at low temperatures in these glasses was sodium calcium silicate. The crystals were found to form through internal nucleation, leading to bulk crystallization. These glasses had high bioactivity in vitro. Even when partially crystalline, they formed typical reaction layers, indicating bioactivity. In fact, sodium calcium silicate crystals were shown to transform in vitro into hydroxyapatite during soaking. However, crystallization should be avoided because it was shown to retard dissolution, bioactivity reactions and complicate fiber drawing process. Glass compositions having low alkali oxide content showed formation of wollastonite crystals on the surface, at about 300°C above the glass transition temperature. The wide range between glass transition and crystallization allowed viscous flow sintering of these compositions. These glasses also withstood the thermal treatments required for fiber drawing processing. Precipitation of calcium and phosphate on fibers of these glasses in SBF suggested that they were osteoconductive. Glasses showing bioactivity crystallize easily, making their hot working challenging. Undesired crystallization can be avoided by choosing suitable compositions and heat treatment parameters, allowing desired product forms to be attained. Small changes in the oxide composition of the glass can have large effects and therefore a thorough understanding of glass crystallization behavior is a necessity for a successful outcome, when designing and manufacturing implants containing bioactive glasses.

THE GENERALIZED MAXIMUM LIKELIHOOD METHOD APPLIED TO HIGH PRESSURE PHASE EQUILIBRIUM

The generalized maximum likelihood method was used to determine binary interaction parameters between carbon dioxide and components of orange essential oil. Vapor-liquid equilibrium was modeled with Peng-Robinson and Soave-Redlich-Kwong equations, using a methodology proposed in 1979 by Asselineau, Bogdanic and Vidal. Experimental vapor-liquid equilibrium data on binary mixtures formed with carbon dioxide and compounds usually found in orange essential oil were used to test the model. These systems were chosen to demonstrate that the maximum likelihood method produces binary interaction parameters for cubic equations of state capable of satisfactorily describing phase equilibrium, even for a binary such as ethanol/CO2. Results corroborate that the Peng-Robinson, as well as the Soave-Redlich-Kwong, equation can be used to describe phase equilibrium for the following systems: components of essential oil of orange/CO2.

HIGH PRESSURE PHASE EQUILIBRIUM: PREDICTION OF ESSENTIAL OIL SOLUBILITY

This work describes a method to predict the solubility of essential oils in supercritical carbon dioxide. The method is based on the formulation proposed in 1979 by Asselineau, Bogdanic and Vidal. The Peng-Robinson and Soave-Redlich-Kwong cubic equations of state were used with the van der Waals mixing rules with two interaction parameters. Method validation was accomplished calculating orange essential oil solubility in pressurized carbon dioxide. The solubility of orange essential oil in carbon dioxide calculated at 308.15 K for pressures of 50 to 70 bar varied from 1.7± 0.1 to 3.6± 0.1 mg/g. For same the range of conditions, experimental solubility varied from 1.7± 0.1 to 3.6± 0.1 mg/g. Predicted values were not very sensitive to initial oil composition.

Microstructural evolution of primary crystallization in diffusion-controlled grain growth

A model has been developed for evaluating grain size distributions in primary crystallizations where the grain growth is diffusion controlled. The body of the model is grounded in a recently presented mean-field integration of the nucleation and growth kinetic equations, modified conveniently in order to take into account a radius-dependent growth rate, as occurs in diffusion-controlled growth. The classical diffusion theory is considered, and a modification of this is proposed to take into account interference of the diffusion profiles between neighbor grains. The potentiality of the mean-field model to give detailed information on the grain size distribution and transformed volume fraction for transformations driven by nucleation and either interface- or diffusion-controlled growth processes is demonstrated. The model is evaluated for the primary crystallization of an amorphous alloy, giving an excellent agreement with experimental data. Grain size distributions are computed, and their properties are discussed.

Effect of sequence distribution on the morphology, crystallization, melting, and biodegradation of poly(epsilon-caprolactone-co-epsilon-caprolactam) copolymers

Two types of poly(epsilon-caprolactone (CLo)-co-poly(epsilon-caprolactam (CLa)) copolymers were prepared by catalyzed hydrolytic ring-opening polymerization. Both cyclic comonomers were added simultaneously in the reaction medium for the First type or materials where copolymers have a random distribution of counits, as evidenced by H-1 and C-13 NMR. For the second type of copolymers, the cyclic comonomers were added sequentially, yielding diblock poly(ester-amides). The materials were characterized by differential scanning calorimetry (DSC), wide- and small-angle X-ray scattering (WAXS and SAXS), and transmission and scanning electron microscopies (TEM and SEM). Their biodegradation in compost was also studied. All copolymers were found to be miscible by the absence of structure in the melt. TEM revealed that all samples exhibited a crystalline lamellar morphology. DSC and WAXS showed that in a wide composition range (CLo contents from 6 to 55%) only the CLa units were capable of crystallization in the random copolymers. The block copolymer samples only experience a small reduction of crystallization and melting temperature with composition, and this was attributed to a dilution effect caused by the miscible noncrystalline CLo units. The comparison between block and random copolymers provided a unique opportunity to distinguish the dilution effect of the CLo units on the crystallization and melting of the polyamide phase from the chemical composition effect in the random copolymers case, where the CLa sequences are interrupted statistically by the CLo units, making the crystallization of the polyamide strongly composition dependent. Finally, the enzymatic degradation of the copolymers in composted soil indicate a synergistic behavior where much faster degradation was obtained for random copolymers witha CLo content larger than 30% than for neat PCL.

Quantitative model for the kinetics of lyotropic phase transitions involving changes in monolayer curvature

We present a kinetic model for transformations between different self-assembled lipid structures. The model shows how data on the rates of phase transitions between mesophases of different geometries can be used to provide information on the mechanisms of the transformations and the transition states involved. This can be used, for example, to gain an insight into intermediate structures in cell membrane fission or fusion. In cases where the monolayer curvature changes on going from the initial to the final mesophase, we consider the phase transition to be driven primarily by the change in the relaxed curvature with pressure or temperature, which alters the relative curvature elastic energies of the two mesophase structures. Using this model, we have analyzed previously published kinetic data on the inter-conversion of inverse bicontinuous cubic phases in the 1-monoolein-30 wt% water system. The data are for a transition between QII(G) and QII(D) phases, and our analysis indicates that the transition state more closely resembles the QII(D) than the QII(G) phase. Using estimated values for the monolayer mean curvatures of the QII(G) and QII(D) phases of -0.123 nm(-1) and -0.133 nm(-1), respectively, gives values for the monolayer mean curvature of the transition state of between -0.131 nm(-1) and -0.132 nm(-1). Furthermore, we estimate that several thousand molecules undergo the phase transition cooperatively within one "cooperative unit", equivalent to 1-2 unit cells of QII(G) or 4-10 unit cells of QII(D).

An oblivious shortest-path routing algorithm for fully connected cubic networks

Fully connected cubic networks (FCCNs) are a class of newly proposed hierarchical interconnection networks for multicomputer systems, which enjoy the strengths of constant node degree and good expandability. The shortest path routing in FCCNs is an open problem. In this paper, we present an oblivious routing algorithm for n-level FCCN with N = 8(n) nodes, and prove that this algorithm creates a shortest path from the source to the destination. At the costs of both an O(N)-parallel-step off-line preprocessing phase and a list of size N stored at each node, the proposed algorithm is carried out at each related node in O(n) time. In some cases the proposed algorithm is superior to the one proposed by Chang and Wang in terms of the length of the routing path. This justifies the utility of our routing strategy. (C) 2006 Elsevier Inc. All rights reserved.

Kinetics and mechanism of the interconversion of inverse bicontinuous cubic mesophases

This paper describes time-resolved x-ray diffraction data monitoring the transformation of one inverse bicontinuous cubic mesophase into another, in a hydrated lipid system. The first section of the paper describes a mechanism for the transformation that conserves the topology of the bilayer, based on the work of Charvolin and Sadoc, Fogden and Hyde, and Benedicto and O'Brien in this area. We show a pictorial representation of this mechanism, in terms of both the water channels and the lipid bilayer. The second section describes the experimental results obtained. The system under investigation was 2:1 lauric acid: dilauroylphosphatidylcholine at a hydration of 50% water by weight. A pressure-jump was used to induce a phase transition from the gyroid (Q(II)(G)) to the diamond (Q(II)(D)) bicontinuous cubic mesophase, which was monitored by time-resolved x-ray diffraction. The lattice parameter of both mesophases was found to decrease slightly throughout the transformation, but at the stage where the Q(II)(D) phase first appeared, the ratio of lattice parameters of the two phases was found to be approximately constant for all pressure-jump experiments. The value is consistent with a topology-preserving mechanism. However, the polydomain nature of our sample prevents us from confirming that the specific pathway is that described in the first section of the paper. Our data also reveal signals from two different intermediate structures, one of which we have identified as the inverse hexagonal (H-II) mesophase. We suggest that it plays a role in the transfer of water during the transformation. The rate of the phase transition was found to increase with both temperature and pressure-jump amplitude, and its time scale varied from the order of seconds to minutes, depending on the conditions employed.