Seawater carbonate chemistry and biological processes of foraminifera, Globigerinoides ruber and Globigerinella siphonifera during experiments, 2011

Specimens of two species of planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera, were grown under controlled laboratory conditions at a range of temperatures (18-31 °C), salinities (32-44 psu) and pH levels (7.9-8.4). The shells were examined for their calcium isotope compositions (d44/40Ca) and strontium to calcium ratios (Sr/Ca) using Thermal Ionization Mass Spectrometry and Inductively Coupled Plasma Mass Spectrometry. Although the total variation in d44/40Ca (~0.3 per mill) in the studied species is on the same order as the external reproducibility, the data set reveals some apparent trends that are controlled by more than one environmental parameter. There is a well-defined inverse linear relationship between d44/40Ca and Sr/Ca in all experiments, suggesting similar controls on these proxies in foraminiferal calcite independent of species. Analogous to recent results from inorganically precipitated calcite, we suggest that Ca isotope fractionation and Sr partitioning in planktic foraminifera are mainly controlled by precipitation kinetics. This postulation provides us with a unique tool to calculate precipitation rates and draws support from the observation that Sr/Ca ratios are positively correlated with average growth rates. At 25 °C water temperature, precipitation rates in G. siphonifera and G. ruber are calculated to be on the order of 2000 and 3000 µmol/m**2/h, respectively. The lower d44/40Ca observed at 29 °C in both species is consistent with increased precipitation rates at high water temperatures. Salinity response of d44/40Ca (and Sr/Ca) in G. siphonifera implies that this species has the highest precipitation rates at the salinity of its natural habitat, whereas increasing salinities appear to trigger higher precipitation rates in G. ruber. Isotope effects that cannot be explained by precipitation rate in planktic foraminifera can be explained by a biological control, related to a vacuolar pathway for supply of ions during biomineralization and a pH regulation mechanism in these vacuoles. In case of an additional pathway via cross-membrane transport, supplying light Ca for calcification, the d44/40Ca of the reservoir is constrained as -0.2 per mill relative to seawater. Using a Rayleigh distillation model, we calculate that calcification occurs in a semi-open system, where less than half of the Ca supplied by vacuolization is utilized for calcite precipitation. Our findings are relevant for interpreting paleo-proxy data on d44/40Ca and Sr/Ca in foraminifera as well as understanding their biomineralization processes.

Mass transfer effects on reactive separation processes - discovery of reactive arrheotropes

Membrane reactor, reactive membrane separation, arrheotrope, azeotrope, dusty gas model, esterification, residue curve map, distillation, kinetics, singular point, bifurcation

Carbon isotope fractionation during calcium carbonate precipitation induced by urease-catalysed hydrolysis of urea

Stable carbon isotopic fractionation during calcium carbonate precipitation induced by urease-catalysed hydrolysis of urea was experimentally investigated in artificial water at a constant temperature of 30 degrees C. Carbon isotope fractionation during urea hydrolysis follows a Rayleigh distillation trend characterized by a C-13-enrichment factor of -20 to -22 parts per thousand. CaCO3 precipitate is up to 17.9 parts per thousand C-13-depleted relative to the urea substrate (-48.9 +/- 0.07 parts per thousand). Initial CaCO3 precipitate forms close to isotopic equilibrium with dissolved inorganic carbon. Subsequent precipitation occurs at -2 to -3 parts per thousand offset from isotopic equilibrium, suggesting that the initial delta C-13 value of CaCO3 is reset through dissolution followed by reprecipitation with urease molecules playing a role in offsetting the delta C-13 value of CaCO3 from isotopic equilibrium. Potentially, this isotopic systematics may provide a tool for the diagnosis of ureolytically-formed carbonate cements used as sealing agent. Moreover, it may serve as a basis to develop a carbon isotope tool for the quantification of ureolytically-induced CO2 sequestration. Finally, it suggests carbon isotope disequilibrium as a hallmark of past enzymatic activity in ancient microbial carbonate formation. (C) 2012 Elsevier B.V. All rights reserved.

Improved water vapour spectroscopy in the 4174-4300 cm-1 region and its impact on SCIAMACHY HDO/H2O measurements

The relative abundance of the heavy water isotopologue HDO provides a deeper insight into the atmospheric hydrological cycle. The SCanning Imaging Absorption spectroMeter for Atmospheric CartograpHY (SCIAMACHY) allows for global retrievals of the ratio HDO/H2O in the 2.3 micron wavelength range. However, the spectroscopy of water lines in this region remains a large source of uncertainty for these retrievals. We therefore evaluate and improve the water spectroscopy in the range 4174–4300 cm−1 and test if this reduces systematic uncertainties in the SCIAMACHY retrievals of HDO/H2O. We use a laboratory spectrum of water vapour to fit line intensity, air broadening and wavelength shift parameters. The improved spectroscopy is tested on a series of ground-based high resolution FTS spectra as well as on SCIAMACHY retrievals of H2O and the ratio HDO/H2O. We find that the improved spectroscopy leads to lower residuals in the FTS spectra compared to HITRAN 2008 and Jenouvrier et al. (2007) spectroscopy, and the retrievals become more robust against changes in the retrieval window. For both the FTS and SCIAMACHY measurements, the retrieved total H2O columns decrease by 2–4% and we find a negative shift of the HDO/H2O ratio, which for SCIAMACHY is partly compensated by changes in the retrieval setup and calibration software. The updated SCIAMACHY HDO/H2O product shows somewhat steeper latitudinal and temporal gradients and a steeper Rayleigh distillation curve, strengthening previous conclusions that current isotope-enabled general circulation models underestimate the variability in the near-surface HDO/H2O ratio.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

(Table 1) Water column biogenic silica concentrations and Si isotopic composition during Marion Dufresne cruise MD166 off South Africa

Southern Ocean biogeochemical processes have an impact on global marine primary production and global elemental cycling, e.g. by likely controlling glacial-interglacial pCO2 variation. In this context, the natural silicon isotopic composition (d30Si) of sedimentary biogenic silica has been used to reconstruct past Si-consumption:supply ratios in the surface waters. We present a new dataset in the Southern Ocean from a IPY-GEOTRACES transect (Bonus-GoodHope) which includes for the first time summer d30Si signatures of suspended biogenic silica (i) for the whole water column at three stations and (ii) in the mixed layer at seven stations from the subtropical zone up to the Weddell Gyre. In general, the isotopic composition of biogenic opal exported to depth was comparable to the opal leaving the mixed layer and did not seem to be affected by any diagenetic processes during settling, even if an effect of biogenic silica dissolution cannot be ruled out in the northern part of the Weddell Gyre. We develop a mechanistic understanding of the processes involved in the modern Si-isotopic balance, by implementing a mixed layer model. We observe that the accumulated biogenic silica (sensu Rayleigh distillation) should satisfactorily describe the d30Si composition of biogenic silica exported out of the mixed layer, within the limit of the current analytical precision on the d30Si. The failures of previous models (Rayleigh and steady state) become apparent especially at the end of the productive period in the mixed layer, when biogenic silica production and export are low. This results from (1) a higher biogenic silica dissolution:production ratio imposing a lower net fractionation factor and (2) a higher Si-supply:Si-uptake ratio supplying light Si-isotopes into the mixed layer. The latter effect is especially expressed when the summer mixed layer becomes strongly Si-depleted, together with a large vertical silicic acid gradient, e.g. in the Polar Front Zone and at the Polar Front.

Geochemistry and soil characteristics of sediment core GeoB4234

Stable isotopes of sedimentary nitrogen and organic carbon are widely used as proxy variables for biogeochemical parameters and processes in the water column. In order to investigate alterations of the primary isotopic signal by sedimentary diagenetic processes, we determined concentrations and isotopic compositions of inorganic nitrogen (IN), organic nitrogen (ON), total nitrogen (TN), and total organic carbon (TOC) on one short core recovered from sediments of the eastern subtropical Atlantic, between the Canary Islands and the Moroccan coast. Changes with depth in concentration and isotopic composition of the different fractions were related to early diagenetic conditions indicated by pore water concentrations of oxygen, nitrate, and ammonium. Additionally, the nature of the organic matter was investigated by Rock-Eval pyrolysis and microscopic analysis. A decrease in ON during aerobic organic matter degradation is accompanied by an increase of the 15N/14N ratio. Changes in the isotopic composition of ON can be described by Rayleigh fractionation kinetics which are probably related to microbial metabolism. The influence of IN depleted in 15N on the bulk sedimentary (TN) isotope signal increases due to organic matter degradation, compensating partly the isotopic changes in ON. In anoxic sediments, fixation of ammonium between clay lattices results in a decrease of stable nitrogen isotope ratio of IN and TN. Changes in the carbon isotopic composition of TOC have to be explained by Rayleigh fractionation in combination with different remineralization kinetics of organic compounds with different isotopic composition. We have found no evidence for preferential preservation of terrestrial organic carbon. Instead, both TOC and refractory organic carbon are dominated by marine organic matter. Refractory organic carbon is depleted in 13C compared to TOC.

Calcium isotopes in planktonic foraminifera - core tops and sediment trap samples

For paleoceanographic studies, it is important to understand the processes that influence the calcium (Ca) isotopic composition of foraminiferal calcite tests preserved in the sediment record. Seven species of planktonic foraminifera from coretop sediments collectively exhibited a Ca temperature dependent fractionation of 0.013 per mil per °C. This is in agreement with previously published estimates for most species of planktonic foraminifera as well as biogenic and inorganic calcite and aragonite. Four species of planktonic foraminifera collected from a sediment trap showed a considerable amount of scatter and no consistent temperature dependent fractionation. Analyzed size fractions of coretop samples show no significant relationship with d44/40Ca. However, preliminary results suggest that the symbiotic and spinose foraminifera G. sacculifer might exhibit a relationship between test size and d44/40Ca. A one-box model in which Ca isotopes are allowed to fractionate by Rayleigh distillation from a biomineralization reservoir (internal pool) was used to constrain the isotopic composition of the original biomineralization Ca reservoir, assuming around 85% of the Ca reservoir is precipitated and the fractionation factor during precipitation is 0.9985 + 0.00002(T ºC). To explain the foraminiferal Ca isotope data, this model indicates that the Ca isotopic composition of the biomineralization reservoir is offset from seawater (approximately -0.8per mil).

Sr/Ca and Cd/Ca ratios of benthic foraminifera in surface deep-sea sediments

Measurements of Sr/Ca of benthic foraminifera show a linear decrease with water depth which is superimposed upon significant variability identified by analyses of individual foraminifera. New data for Cd/Ca support previous work in defining a contrast between waters shallower and deeper than ~2500 m. Measured element partition coefficients in foraminiferal calcium carbonate relative to sea water (D) have been described by means of a one-box model in which elements are extracted by Rayleigh distillation from a biomineralization reservoir that serves for calcification with a constant fractionation factor (alpha), such that D = (1 - f**alpha)/(l - f), where f is the proportion of Ca remaining after precipitation. A modification to the model recognises differences in element speciation. The model is consistent with differences between D[Sr], D[Ba], and D[Cd] in benthic but not planktonic foraminifera. Depth variations in D for Sr and Ba are consistent with the model, as are differences in depth variation of D[Cd] in calcitic and aragonitic benthic foraminifera. The shallower depth variations may reflect increasing calcification rates with increasing water depth to an optimum of about 2500 m. Observations of unusually lower DCd for some deep waters, not accompanied by similar [Sr], or D[Ba] may be because of dissolution or a calcification response to a lower carbonate saturation state.

Ethyl carbamate kinetics in double distillation of sugar cane spirit

The aim of this study was to verify the effect of a double distillation on the reduction of the ethyl carbamate content in sugar cane spirit. Ethyl carbamate is a potentially carcinogenic compound normally present at critical levels in sugar cane spirit, constituting a public health problem and therefore hindering the export of this beverage. The ethanol, copper and ethyl carbamate contents were evaluated, using gas chromatography/mass spectroscopy, during a double distillation of the fermented sugar cane juice. The distillate fraction from the first distillation accumulated 30% of the ethyl carbamate formed. In the second distillation, the ethyl carbamate and the copper content increased during the process as the alcohol content decreased, and only 3% of the ethyl carbamate formed was collected in the spirit. Double distillation decreased the ethyl carbamate content in the sugar cane spirit by 97%. (C) Copyright 2012 The Institute of Brewing & Distilling

Metamorphism and kinetics in the Torres del Paine contact aureole

Contact aureoles provide an excellent geologic environment to study the mechanisms of metamorphic reactions in a natural system. The Torres del Paine (TP) intrusion is one of the most spectacular natural laboratories because of its excellent outcrop conditions. It formed in a period from 12.59 to 12.43 Ma and consists of three large granite and four smaller mafic batches. The oldest granite is on top, the youngest at the bottom of the granitic complex, and the granites overly the mafic laccolith. The TP intruded at a depth of 2-3 km into regional metamorphic anchizone to greenschist facies pelites, sandstones, and conglomerates of the Cerro Toro and Punta Barrosa formations. It formed a thin contact aureole of 150-400 m width. This thesis focuses on the reaction kinetics of the mineral cordierite in the contact aureole using quantitative textural analysis methods. First cordierite was formed from chlorite break¬down (zone I, ca. 480 °C, 750 bar). The second cordierite forming reaction was the muscovite break-down, which is accompanied by a modal decrease in biotite and the appearance of k- feldspar (zone II, 540-550 °C, 750 bar). Crystal sizes of the roundish, poikiloblastic cordierites were determined from microscope thin section images by manually marking each crystal. Images were then automatically processed with Matlab. The correction for the intersection probability of each crystal radius yields the crystal size distribution in the rock. Samples from zone I below the laccolith have the largest crystals (0.09 mm). Cordierites from zone II are smaller, with a maximum crystal radius of 0.057 mm. Rocks from zone II have a larger number of small cordierite crystals than rocks from zone I. A combination of these quantitative analysis with numerical modeling of nucleation and growth, is used to infer nucleation and growth parameters which are responsible for the observed mineral textures. For this, the temperature-time paths of the samples need to be known. The thermal history is complex because the main body of the intrusion was formed by several intrusive batches. The emplacement mechanism and duration of each batch can influence the thermal structure in the aureole. A possible subdivision of batches in smaller increments, so called pulses, will focus heat at the side of the intrusion. Focusing all pulses on one side increases the contact aureole size on that side, but decreases it on the other side. It forms a strongly asymmetric contact aureole. Detailed modeling shows that the relative thicknesses of the TP contact aureole above and below the intrusion (150 and 400 m) are best explained by a rapid emplacement of at least the oldest granite batch. Nevertheless, temperatures are significantly too low in all models, compared to observed mineral assemblages in the hornfelses. Hence, an other important thermal mechanisms needs to take place in the host rock. Clastic minerals in the immature sediments outside the contact aureole are hydrated due to small amounts of expelled fluids during contact metamorphism. This leads to a temperature increase of up to 50 °C. The origin of fluids can be traced by stable isotopes. Whole rock stable isotope data (6D and δ180) and chlorine concentrations in biotite document that the TP intrusion induced only very small amounts of fluid flow. Oxygen whole rock data show δ180 values between 9.0 and 10.0 %o within the first 5 m of the contact. Values increase to 13.0 - 15.0 %o further away from the intrusion. Whole rock 6D values display a more complex zoning. First, host rock values (-90 to -70 %o) smoothly decrease towards the contact by ca. 20 %o, up to a distance of ca. 150 m. This is followed by an increase of ca. 20 %o within the innermost 150 m of the aureole (-97.0 to -78 %o at the contact). The initial decrease in 6D values is interpreted to be due to Rayleigh fractionation accompanying the dehydration reactions forming cordierite, while the final increase reflects infiltration of water-rich fluids from the intrusion. An over-estimate on the quantity and the corresponding thermal effect yields a temperature increase of less than 30 °C. This suggests that fluid flow might have contributed only for a small amount to the thermal evolution of the system. A combination of the numerical growth model with the thermal model, including the hydration reaction enthalpies but neglecting fluid flow and incremental growth, can be used to numerically reproduce the observed cordierite textures in the contact aureole. This yields kinetic parameters which indicate fast cordierite crystallization before the thermal peak in the inner aureole, and continued reaction after the thermal peak in the outermost aureole. Only small temperature dependencies of the kinetic parameters seem to be needed to explain the obtained crystal size data. - Les auréoles de contact offrent un cadre géologique privilégié pour l'étude des mécanismes de réactions métamorphiques associés à la mise en place de magmas dans la croûte terrestre. Par ses conditions d'affleurements excellentes, l'intrusion de Torres del Paine représente un site exceptionnel pour améliorer nos connaissances de ces processus. La formation de cette intrusion composée de trois injections granitiques principales et de quatre injections mafiques de volume inférieur couvre une période allant de 12.50 à 12.43 Ma. Le plus vieux granite forme la partie sommitale de l'intrusion alors que l'injection la plus jeune s'observe à la base du complexe granitique; les granites recouvrent la partie mafique du laccolite. L'intrusion du Torres del Paine s'est mise en place a 2-3 km de profondeur dans un encaissant métamorphique. Cet encaissant est caractérisé par un métamorphisme régional de faciès anchizonal à schiste vert et est composé de pélites, de grès, et des conglomérats des formations du Cerro Toro et Punta Barrosa. La mise en place des différentes injections granitiques a généré une auréole de contact de 150-400 m d'épaisseur autour de l'intrusion. Cette thèse se concentre sur la cinétique de réaction associée à la formation de la cordiérite dans les auréoles de contact en utilisant des méthodes quantitatives d'analyses de texture. On observe plusieurs générations de cordiérite dans l'auréole de contact. La première cordiérite est formée par la décomposition de la chlorite (zone I, environ 480 °C, 750 bar), alors qu'une seconde génération de cordiérite est associée à la décomposition de la muscovite, laquelle est accompagnée par une diminution modale de la teneur en biotite et l'apparition de feldspath potassique (zone II, 540-550 °C, 750 bar). Les tailles des cristaux de cordiérites arrondies et blastic ont été déterminées en utilisant des images digitalisées des lames minces et en marquant individuellement chaque cristal. Les images sont ensuite traitées automatiquement à l'aide du programme Matlab. La correction de la probabilité d'intersection en fonction du rayon des cristaux permet de déterminer la distribution de la taille des cristaux dans la roche. Les échantillons de la zone I, en dessous du lacolite, sont caractérisés par de relativement grands cristaux (0.09 mm). Les cristaux de cordiérite de la zone II sont plus petits, avec un rayon maximal de 0.057 mm. Les roches de la zone II présentent un plus grand nombre de petits cristaux de cordiérite que les roches de la zone I. Une combinaison de ces analyses quantitatives avec un modèle numérique de nucléation et croissance a été utilisée pour déduire les paramètres de nucléation et croissance contrôlant les différentes textures minérales observées. Pour développer le modèle de nucléation et de croissance, il est nécessaire de connaître le chemin température - temps des échantillons. L'histoire thermique est complexe parce que l'intrusion est produite par plusieurs injections successives. En effet, le mécanisme d'emplace¬ment et la durée de chaque injection peuvent influencer la structure thermique dans l'auréole. Une subdivision des injections en plus petits incréments, appelés puises, permet de concentrer la chaleur dans les bords de l'intrusion. Une mise en place préférentielle de ces puises sur un côté de l'intrusion modifie l'apport thermique et influence la taille de l'auréole de contact produite, auréole qui devient asymétrique. Dans le cas de la première injection de granite, une modélisation détaillée montre que l'épaisseur relative de l'auréole de contact de Torres del Paine au-dessus et en dessous de l'intrusion (150 et 400 m) est mieux expliquée par un emplacement rapide du granite. Néanmoins, les températures calculées dans l'auréole de con¬tact sont trop basses pour que les modèles thermiques soient cohérants par rapport à la taille de cette auréole. Ainsi, un autre mecanisme exothermique est nécessaire pour permettre à la roche encais¬sante de produire les assemblages observés. L'observation des roches encaissantes entourant les granites montre que les minéraux clastiques dans les sédiments immatures au-dehors de l'auréole sont hydratés suite à la petite quantité de fluide expulsée durant le métamorphisme de contact et/ou la mise en place des granites. Les réactions d'hydratation peuvent permettre une augmentation de la température jusqu'à 50 °C. Afin de déterminer l'origine des fluides, une étude isotopique de roches de l'auréole de contact a été entreprise. Les isotopes stables d'oxygène et d'hydrogène sur la roche totale ainsi que la concentration en chlore dans la biotite indiquent que la mise en place des granites du Torres del Paine n'induit qu'une circulation de fluide limitée. Les données d'oxygène sur roche totale montrent des valeurs δ180 entre 9.0 et 10.0%o au sein des cinq premiers mètres du contact. Les valeurs augmentent jusqu'à 13.0 - 15.0 plus on s'éloigne de l'intrusion. Les valeurs 5D sur roche totale montrent une zonation plus complexe. Les valeurs de la roche encaissante (-90 à -70%o) diminuent progressivement d'environ 20%o depuis l'extérieur de l'auréole jusqu'à une distance d'environ 150 m du granite. Cette diminution est suivie par une augmentation d'environ 20%o au sein des 150 mètres les plus proches du contact (-97.0 à -78%o au contact). La diminution initiale des valeurs de 6D est interprétée comme la conséquence du fractionnement de Rayleigh qui accompagne les réactions de déshydratation formant la cordiérite, alors que l'augmentation finale reflète l'infiltration de fluide riche en eau venant de l'intrusion. A partir de ces résultats, le volume du fluide issu du granite ainsi que son effet thermique a pu être estimé. Ces résultats montrent que l'augmentation de température associée à ces fluides est limitée à un maximum de 30 °C. La contribution de ces fluides dans le bilan thermique est donc faible. Ces différents résultats nous ont permis de créer un modèle thermique associé à la for¬mation de l'auréole de contact qui intègre la mise en place rapide du granite et les réactions d'hydratation lors du métamorphisme. L'intégration de ce modèle thermique dans le modèle numérique de croissance minérale nous permet de calculer les textures des cordiérites. Cepen¬dant, ce modèle est dépendant de la vitesse de croissance et de nucléation de ces cordiérites. Nous avons obtenu ces paramètres en comparant les textures prédites par le modèle et les textures observées dans les roches de l'auréole de contact du Torres del Paine. Les paramètres cinétiques extraits du modèle optimisé indiquent une cristallisation rapide de la cordiérite avant le pic thermique dans la partie interne de l'auréole, et une réaction continue après le pic thermique dans la partie la plus externe de l'auréole. Seules de petites dépendances de température des paramètres de cinétique semblent être nécessaires pour expliquer les don¬nées obtenues sur la distribution des tailles de cristaux. Ces résultats apportent un éclairage nouveau sur la cinétique qui contrôle les réactions métamorphiques.

Patterns in and kinetics of phase separation in binary mixtures

This work focused mainly on two aspects of kinetics of phase separation in binary mixtures. In the first part, we studied the interplay of hydrodynamics and the phase separation of binary mixtures. A considerably flat container (a laterally extended geometry), at an aspect ratio of 14:1 (diameter: height) was chosen, so that any hydrodynamic instabilities, if they arise, could be tracked. Two binary mixtures were studied. One was a mixture of methanol and hexane, doped with 5% ethanol, which phase separated under cooling. The second was a mixture of butoxyethanol and water, doped with 2% decane, which phase separated under heating. The dopants were added to bring down the phase transition temperature around room temperature.rnrnAlthough much work has been done already on classical hydrodynamic instabilities, not much has been done in the understanding of the coupling between phase separation and hydrodynamic instabilities. This work aimed at understanding the influence of phase separation in initiating any hydrodynamic instability, and also vice versa. Another aim was to understand the influence of the applied temperature protocol on the emergence of patterns characteristic to hydrodynamic instabilities. rnrnOn slowly cooling the system continuously, at specific cooling rates, patterns were observed in the first mixture, at the start of phase separation. They resembled the patterns observed in classical Rayleigh-Bénard instability, which arises when a liquid continuously is heated from below. To suppress this classical convection, the cooling setup was tuned such that the lower side of the sample always remained cooler by a few millikelvins, relative to the top. We found that the nature of patterns changed with different cooling rates, with stable patterns appearing for a specific cooling rate (1K/h). On the basis of the cooling protocol, we estimated a modified Rayleigh number for our system. We found that the estimated modified Rayleigh number is near the critical value for instability, for cooling rates between 0.5K/h and 1K/h. This is consistent with our experimental findings. rnrnThe origin of the patterns, in spite of the lower side being relatively colder with respect to the top, points to two possible reasons. 1) During phase separation droplets of either phases are formed, which releases a latent heat. Our microcalorimetry measurements show that the rise in temperature during the first phase separation is in the order of 10-20millikelvins, which in some cases is enough to reverse the applied temperature bias. Thus phase separation in itself initiates a hydrodynamic instability. 2) The second reason comes from the cooling protocol itself. The sample was cooled from above and below. At sufficiently high cooling rates, there are situations where the interior of the sample is relatively hotter than both top and bottom of the sample. This is sufficient to create an instability within the cell. Our experiments at higher cooling rates (5K/h and above) show complex patterns, which hints that there is enough convection even before phase separation occurs. Infact, theoretical work done by Dr.Hayase show that patterns could arise in a system without latent heat, with symmetrical cooling from top and bottom. The simulations also show that the patterns do not span the entire height of the sample cell. This is again consistent with the cell sizes measured in our experiment.rnrnThe second mixture also showed patterns at specific heating rates, when it was continuously heated inducing phase separation. In this case though, the sample was turbid for a long time until patterns appeared. A meniscus was most probably formed before the patterns emerged. We attribute the reason of patterns in this case to Marangoni convection, which is present in systems with an interface, where local differences in surface tension give rise to an instability. Our estimates for the Rayleigh number also show a significantly lower number than that's required for RB-type instability.rnrnIn the first part of the work, therefore, we identify two different kinds of hydrodynamic instabilities in two different mixtures. Both are observed during, or after the first phase separation. Our patterns compare with the classical convection patterns, but here the origins are from phase separation and the cooling protocol.rnrnIn the second part of the work, we focused on the kinetics of phase separation in a polymer solution (polystyrene and methylcyclohexane), which is cooled continuously far down into the two phase region. Oscillations in turbidity, denoting material exchange between the phases are seen. Three processes contribute to the phase separation: Nucleation of droplets, their growth and coalescence, and their subsequent sedimentation. Experiments in low molecular binary mixtures had led to models of oscillation [43] which considered sedimentation time scales much faster than the time scales of nucleation and growth. The size and shape of the sample therefore did not matter in such situations. The oscillations in turbidity were volume-dominated. The present work aimed at understanding the influence of sedimentation time scales for polymer mixtures. Three heights of the sample with same composition were studied side by side. We found that periods increased with the sample height, thus showing that sedimentation time determines the period of oscillations in the polymer solutions. We experimented with different cooling rates and different compositions of the mixture, and we found that periods are still determined by the sample height, and therefore by sedimentation time. rnrnWe also see that turbidity emerges in two ways; either from the interface, or throughout the sample. We suggest that oscillations starting from the interface are due to satellite droplets that are formed on droplet coalescence at the interface. These satellite droplets are then advected to the top of the sample, and they grow, coalesce and sediment. This type of an oscillation wouldn't require the system to pass the energy barrier required for homogenous nucleation throughout the sample. This mechanism would work best in sample where the droplets could be effectively advected throughout the sample. In our experiments, we see more interface dominated oscillations in the smaller cells and lower cooling rates, where droplet advection is favourable. In larger samples and higher cooling rates, we mostly see that the whole sample becomes turbid homogenously, which requires the system to pass the energy barrier for homogenous nucleation.rnrnOscillations, in principle, occur since the system needs to pass an energy barrier for nucleation. The height of the barrier decreases with increasing supersaturation, which in turn is from the temperature ramp applied. This gives rise to a period where the system is clear, in between the turbid periods. At certain specific cooling rates, the system can follow a path such that the start of a turbid period coincides with the vanishing of the last turbid period, thus eliminating the clear periods. This means suppressions of oscillations altogether. In fact we experimentally present a case where, at a certain cooling rate, oscillations indeed vanish. rnrnThus we find through this work that the kinetics of phase separation in polymer solution is different from that of a low molecular system; sedimentation time scales become relevant, and therefore so does the shape and size of the sample. The role of interface in initiating turbid periods also become much more prominent in this system compared to that in low molecular mixtures.rnrnIn summary, some fundamental properties in the kinetics of phase separation in binary mixtures were studied. While the first part of the work described the close interplay of the first phase separation with hydrodynamic instabilities, the second part investigated the nature and determining factors of oscillations, when the system was cooled deep into the two phase region. Both cases show how the geometry of the cell can affect the kinetics of phase separation. This study leads to further fundamental understandings of the factors contributing to the kinetics of phase separation, and to the understandings of what can be controlled and tuned in practical cases. rn

DYNAMICS AND CONTROL OF REACTIVE DISTILLATION PROCESS FOR MONOMER SYNTHESIS OF POLYCARBONATE PLANTS

Polycarbonate (PC) is an important engineering thermoplastic that is currently produced in large industrial scale using bisphenol A and monomers such as phosgene. Since phosgene is highly toxic, a non-phosgene approach using diphenyl carbonate (DPC) as an alternative monomer, as developed by Asahi Corporation of Japan, is a significantly more environmentally friendly alternative. Other advantages include the use of CO2 instead of CO as raw material and the elimination of major waste water production. However, for the production of DPC to be economically viable, reactive-distillation units are needed to obtain the necessary yields by shifting the reaction-equilibrium to the desired products and separating the products at the point where the equilibrium reaction occurs. In the field of chemical reaction engineering, there are many reactions that are suffering from the low equilibrium constant. The main goal of this research is to determine the optimal process needed to shift the reactions by using appropriate control strategies of the reactive distillation system. An extensive dynamic mathematical model has been developed to help us investigate different control and processing strategies of the reactive distillation units to increase the production of DPC. The high-fidelity dynamic models include extensive thermodynamic and reaction-kinetics models while incorporating the necessary mass and energy balance of the various stages of the reactive distillation units. The study presented in this document shows the possibility of producing DPC via one reactive distillation instead of the conventional two-column, with a production rate of 16.75 tons/h corresponding to start reactants materials of 74.69 tons/h of Phenol and 35.75 tons/h of Dimethyl Carbonate. This represents a threefold increase over the projected production rate given in the literature based on a two-column configuration. In addition, the purity of the DPC produced could reach levels as high as 99.5% with the effective use of controls. These studies are based on simulation done using high-fidelity dynamic models.

Mass transfer and chemical reaction on a distillation plate

A multistage distillation column in which mass transfer and a reversible chemical reaction occurred simultaneously, has been investigated to formulate a technique by which this process can be analysed or predicted. A transesterification reaction between ethyl alcohol and butyl acetate, catalysed by concentrated sulphuric acid, was selected for the investigation and all the components were analysed on a gas liquid chromatograph. The transesterification reaction kinetics have been studied in a batch reactor for catalyst concentrations of 0.1 - 1.0 weight percent and temperatures between 21.4 and 85.0 °C. The reaction was found to be second order and dependent on the catalyst concentration at a given temperature. The vapour liquid equilibrium data for six binary, four ternary and one quaternary systems are measured at atmospheric pressure using a modified Cathala dynamic equilibrium still. The systems with the exception of ethyl alcohol - butyl alcohol mixtures, were found to be non-ideal. Multicomponent vapour liquid equilibrium compositions were predicted by a computer programme which utilised the Van Laar constants obtained from the binary data sets. Good agreement was obtained between the predicted and experimental quaternary equilibrium vapour compositions. Continuous transesterification experiments were carried out in a six stage sieve plate distillation column. The column was 3" in internal diameter and of unit construction in glass. The plates were 8" apart and had a free area of 7.7%. Both the liquid and vapour streams were analysed. The component conversion was dependent on the boilup rate and the reflux ratio. Because of the presence of the reaction, the concentration of one of the lighter components increased below the feed plate. In the same region a highly developed foam was formed due to the presence of the catalyst. The experimental results were analysed by the solution of a series of simultaneous enthalpy and mass equations. Good agreement was obtained between the experimental and calculated results.

Effects Of High Pressure Homogenization On The Activity, Stability, Kinetics And Three-dimensional Conformation Of A Glucose Oxidase Produced By Aspergillus Niger.

High pressure homogenization (HPH) is a non-thermal method, which has been employed to change the activity and stability of biotechnologically relevant enzymes. This work investigated how HPH affects the structural and functional characteristics of a glucose oxidase (GO) from Aspergillus niger. The enzyme was homogenized at 75 and 150 MPa and the effects were evaluated with respect to the enzyme activity, stability, kinetic parameters and molecular structure. The enzyme showed a pH-dependent response to the HPH treatment, with reduction or maintenance of activity at pH 4.5-6.0 and a remarkable activity increase (30-300%) at pH 6.5 in all tested temperatures (15, 50 and 75°C). The enzyme thermal tolerance was reduced due to HPH treatment and the storage for 24 h at high temperatures (50 and 75°C) also caused a reduction of activity. Interestingly, at lower temperatures (15°C) the activity levels were slightly higher than that observed for native enzyme or at least maintained. These effects of HPH treatment on function and stability of GO were further investigated by spectroscopic methods. Both fluorescence and circular dichroism revealed conformational changes in the molecular structure of the enzyme that might be associated with the distinct functional and stability behavior of GO.