Effects of cation substitutions on the physical properties of M2M1T2O6 pyroxenes

In this PhD study, the effects of the cation substitutions on the physical properties of pyroxenes have been discussed. The results of this work extend the knowledge on pyroxenes with different chemical compositions. These properties might be used in the development of ceramic pigments, advanced materials and for the mineralogical phase identification. First of all, the crystallographic differences between Ge and Si pyroxenes have been examined. The structure of C2/c Ca rich Ge clinopyroxenes is very close to the low pressure C2/c structural configuration found in Ca-rich Si-pyroxenes. The shear of the unit cell is very similar, and the difference between a Ge end member and the corresponding Si-rich one is less than 1°. Instead, a remarkable difference exists between Ca-poor Si and Ge clinopyroxenes. First, Ca-poor Ge pyroxenes do not display a P21/c symmetry, but retain the C2/c symmetry; second, the observed C2/c structure shows, at room pressure, the configuration with highly kinked tetrahedral chains characteristic of the high pressure C2/c symmetry of Si Ca-poor pyroxenes. In orthopyroxenes, with Pbca symmetry, Ge-pyroxenes have volume larger than Si-pyroxenes. Samples along the system CaCoGe2O6 - CoCoGe2O6 have been synthesized at three different temperatures: 1050 °C, 1200 °C and 1250 °C. The aim of these solid state syntheses was to obtain a solid solution at ambient pressure, since the analogues Si-system needs high pressure. Unfortunately, very limited solution occurs because the structure forms of the two end member (high temperature for CaCoGe2O6 and high pressure CoCoGe2O6) are incompatible. The phase diagram of this system has been sketched and compared to that of Si. The cobalt end member (CoCoGe2O6) is stable at ambient pressure in two symmetries: at 1050 °C C2/c and 1200 °C Pbca. The impurity phase formed during these experiments is cobalt spinel. Raman spectroscopy has been used to investigate the vibrational properties of Ca-pyroxenes CaCoGe2O6, CaMgGe2O6, CaMgSi2O6 and CaCoSi2O6. A comparison between silicate and germanate pyroxenes shows significant changes in peak positions of the corresponding modes caused mainly by the difference of the Ge-Si atomic weight along with the distortion and compression of the coordination polyhedra. Red shift in Raman spectra of germanates has been calculated by a rough scale factor calculated by a simple harmonic oscillator model, considering the different bond lengths for 4-coordinated Si ~ 1.60- 1.65 Å vs Ge–O distance ~1.70 - 1.80 Å. The Raman spectra of CaMgGe2O6 and CaCoGe2O6 have been classified, in analogy with silicate (Wang et al., 2001) counterparts, in different ranges: - R1 (880-640 cm-1): strong T-O stretching modes of Ge and non-bridging O1 and O2 atoms within the GeO4 tetrahedron; - R2 (640-480 cm-1): stretching/bending modes of Ge-Obr-Ge bonds (chain stretching and chain bending); - R4 (480-360 cm-1): O-Ge-O vibrations; - R3 (360-240 cm-1): motions of the cations in M2 and M1 sites correlated with tetrahedral chain motion and tilting tetrahedra; - R5 (below 240 cm-1): lattice modes. The largest shift with respect to CaMgSi2O6 - CaCoSi2O6 is shown by the T-O stretching and chain modes. High-pressure Raman spectroscopy (up to about 8 GPa) on the same samples of Ca-pyroxenes using an ETH-type diamond anvil cell shows no phase transition within the P-ranges investigated, as all the peak positions vary linearly as a function of pressure. Our data confirm previous experimental findings on Si-diopside (Chopelas and Serghiou, 2000). In the investigated samples, all the Raman peaks shift upon compression, but the major changes in wavenumber with pressure are attributed to the chain bending (Ge-Obr-Ge bonds) and tetrahedra stretching modes (Ge-Onbr). Upon compression, the kinking angle, the bond lengths and T-T distances between tetrahedra decrease and consequently the wavenumber of the bending chain mode and tetrahedra stretching mode increases. Ge-pyroxenes show the higher P-induced peak-position shifts, being more compressible than corresponding silicates. The vibrational properties of CaM2+Ge2O6 (M2+ =Mg, Mn, Fe, Co, Ni, Zn) are reported for the first time. The wavenumber of Ge-Obr-Ge bending modes decreases linearly with increasing ionic radius of the M1 cation. No simple correlation has been found with M1 atomic mass or size or crystallographic parameters for the peak at ~850 cm-1 and in the low wavenumber regions. The magnetic properties of the system CaCoSi2O6 - CoCoSi2O6 have been investigated by magnetometry. The join is always characterized by 1 a.p.f.u. of cobalt in M1 site and this causes a pure collinear antiferromagnetic behaviour of the intra-chain superexchange interaction involving Co ions detected in all the measurements, while the magnetic order developed by the cobalt ions in M2 site (intra-chain) is affected by weak ferromagnetism, due to the non-collinearity of their antiferromagnetic interaction. In magnetically ordered systems, this non-collinearity effect promotes a spin canting of anti-parallel aligned magnetic moments and thus is a source of weak ferromagnetic behaviour in an antiferromagnetic. The weak ferromagnetism can be observed only for the samples with Co content higher than 0.5 a.p.f.u. in M2, when the concentration is sufficiently high to create a long range order along the M2 chain which is magnetically independent of M1 chain. The ferromagnetism was detected both in the M(T) at 10 Oe and M(H).

Theoretical foundations of neural networks

Neural networks have often been motivated by superficial analogy with biological nervous systems. Recently, however, it has become widely recognised that the effective application of neural networks requires instead a deeper understanding of the theoretical foundations of these models. Insight into neural networks comes from a number of fields including statistical pattern recognition, computational learning theory, statistics, information geometry and statistical mechanics. As an illustration of the importance of understanding the theoretical basis for neural network models, we consider their application to the solution of multi-valued inverse problems. We show how a naive application of the standard least-squares approach can lead to very poor results, and how an appreciation of the underlying statistical goals of the modelling process allows the development of a more general and more powerful formalism which can tackle the problem of multi-modality.

Finite connectivity systems as error-correcting codes

We investigate the performance of parity check codes using the mapping onto spin glasses proposed by Sourlas. We study codes where each parity check comprises products of K bits selected from the original digital message with exactly C parity checks per message bit. We show, using the replica method, that these codes saturate Shannon's coding bound for K?8 when the code rate K/C is finite. We then examine the finite temperature case to asses the use of simulated annealing methods for decoding, study the performance of the finite K case and extend the analysis to accommodate different types of noisy channels. The analogy between statistical physics methods and decoding by belief propagation is also discussed.

Statistical mechanics of mutual information maximization

An unsupervised learning procedure based on maximizing the mutual information between the outputs of two networks receiving different but statistically dependent inputs is analyzed (Becker S. and Hinton G., Nature, 355 (1992) 161). By exploiting a formal analogy to supervised learning in parity machines, the theory of zero-temperature Gibbs learning for the unsupervised procedure is presented for the case that the networks are perceptrons and for the case of fully connected committees.

CDMA multiuser detection, neural networks, and statistical mechanics

A novel approach, based on statistical mechanics, to analyze typical performance of optimum code-division multiple-access (CDMA) multiuser detectors is reviewed. A `black-box' view ot the basic CDMA channel is introduced, based on which the CDMA multiuser detection problem is regarded as a `learning-from-examples' problem of the `binary linear perceptron' in the neural network literature. Adopting Bayes framework, analysis of the performance of the optimum CDMA multiuser detectors is reduced to evaluation of the average of the cumulant generating function of a relevant posterior distribution. The evaluation of the average cumulant generating function is done, based on formal analogy with a similar calculation appearing in the spin glass theory in statistical mechanics, by making use of the replica method, a method developed in the spin glass theory.

Flow patterns and performance of distillation trays

It is known that distillation tray efficiency depends on the liquid flow pattern, particularly for large diameter trays. Scale·up failures due to liquid channelling have occurred, and it is known that fitting flow control devices to trays sometirr.es improves tray efficiency. Several theoretical models which explain these observations have been published. Further progress in understanding is at present blocked by lack of experimental measurements of the pattern of liquid concentration over the tray. Flow pattern effects are expected to be significant only on commercial size trays of a large diameter and the lack of data is a result of the costs, risks and difficulty of making these measurements on full scale production columns. This work presents a new experiment which simulates distillation by water cooling. and provides a means of testing commercial size trays in the laboratory. Hot water is fed on to the tray and cooled by air forced through the perforations. The analogy between heat and mass transfer shows that the water temperature at any point is analogous to liquid concentration and the enthalpy of the air is analogous to vapour concentration. The effect of the liquid flow pattern on mass transfer is revealed by the temperature field on the tray. The experiment was implemented and evaluated in a column of 1.2 m. dia. The water temperatures were measured by thennocouples interfaced to an electronic computerised data logging system. The "best surface" through the experimental temperature measurements was obtained by the mathematical technique of B. splines, and presented in tenos of lines of constant temperature. The results revealed that in general liquid channelling is more imponant in the bubbly "mixed" regime than in the spray regime. However, it was observed that severe channelling also occurred for intense spray at incipient flood conditions. This is an unexpected result. A computer program was written to calculate point efficiency as well as tray efficiency, and the results were compared with distillation efficiencies for similar loadings. The theoretical model of Porter and Lockett for predicting distillation was modified to predict water cooling and the theoretical predictions were shown to be similar to the experimental temperature profiles. A comparison of the repeatability of the experiments with an errors analysis revealed that accurate tray efficiency measurements require temperature measurements to better than ± 0.1 °c which is achievable with conventional techniques. This was not achieved in this work, and resulted in considerable scatter in the efficiency results. Nevertheless it is concluded that the new experiment is a valuable tool for investigating the effect of the liquid flow pattern on tray mass transfer.

Local shell-side coefficients in shell and tube exchanges the use of a mass transfer technique

A diffusion-controlled electrochemical mass transfer technique has been employed in making local measurements of shell-side coefficients in segmentally baffled shell and tube heat exchangers. Corresponding heat transfer data are predicted through the Chilton and Colburn heat and mass transfer analogy. Mass transfer coefficients were measured for baffle spacing lengths of individual tubes in an internal baffle compartment. Shell-side pressure measurements were also made. Baffle compartment average coefficients derived from individual tube coefficients are shown to be in good agreement with reported experimental bundle average heat transfer data for a heat exchanger model of similar geometry. Mass transfer coefficients of individual tubes compare favourably with those obtained previously by another mass transfer technique. Experimental data are reported for a variety of segmental baffle configurations over the shell-side Reynolds number range 100 to 42 000. Baffles with zero clearances were studied at three baffle cuts and two baffle spacings. Baffle geometry is shown to have a large effect on the distribution of tube coefficients within the baffle compartment. Fluid "jetting" is identified with some baffle configurations. No simple characteristic velocity is found to correlate zonal or baffle compartment average mass transfer data for the effect of both baffle cut and baffle spacing. Experiments with baffle clearances typical of commercial heat exchangers are also reported. The effect of leakage streams associated with these baffles is identified. Investigations were extended to double segmental baffles for which no data had previously been published. The similarity in the shell-side characteristics of this baffle arrangement and two parallel single segmental baffle arrangements is demonstrated. A general relationship between the shell-side mass transfer performance and pressure drop was indicated by the data for all the baffle configurations examined.

Computer simulation of liquid flow patterns on distillation trays

This thesis describes work carried out to improve the fundamental modelling of liquid flows on distillation trays. A mathematical model is presented based on the principles of computerised fluid dynamics. It models the liquid flow in the horizontal directions allowing for the effects of the vapour through the use of an increased liquid turbulence, modelled by an eddy viscosity, and a resistance to liquid flow caused by the vapour being accelerated horizontally by the liquid. The resultant equations are similar to the Navier-Stokes equations with the addition of a resistance term.A mass-transfer model is used to calculate liquid concentration profiles and tray efficiencies. A heat and mass transfer analogy is used to compare theoretical concentration profiles to experimental water-cooling data obtained from a 2.44 metre diameter air-water distillation simulation rig. The ratios of air to water flow rates are varied in order to simulate three pressures: vacuum, atmospheric pressure and moderate pressure.For simulated atmospheric and moderate pressure distillation, the fluid mechanical model constantly over-predicts tray efficiencies with an accuracy of between +1.7% and +11.3%. This compares to -1.8% to -10.9% for the stagnant regions model (Porter et al. 1972) and +12.8% to +34.7% for the plug flow plus back-mixing model (Gerster et al. 1958). The model fails to predict the flow patterns and tray efficiencies for vacuum simulation due to the change in the mechanism of liquid transport, from a liquid continuous layer to a spray as the liquid flow-rate is reduced. This spray is not taken into account in the development of the fluid mechanical model. A sensitivity analysis carried out has shown that the fluid mechanical model is relatively insensitive to the prediction of the average height of clear liquid, and a reduction in the resistance term results in a slight loss of tray efficiency. But these effects are not great. The model is quite sensitive to the prediction of the eddy viscosity term. Variations can produce up to a 15% decrease in tray efficiency. The fluid mechanical model has been incorporated into a column model so that statistical optimisation techniques can be employed to fit a theoretical column concentration profile to experimental data. Through the use of this work mass-transfer data can be obtained.

Use of computational fluid dynamics in the design of bubble column reactors

Investigations into the modelling techniques that depict the transport of discrete phases (gas bubbles or solid particles) and model biochemical reactions in a bubble column reactor are discussed here. The mixture model was used to calculate gas-liquid, solid-liquid and gasliquid-solid interactions. Multiphase flow is a difficult phenomenon to capture, particularly in bubble columns where the major driving force is caused by the injection of gas bubbles. The gas bubbles cause a large density difference to occur that results in transient multi-dimensional fluid motion. Standard design procedures do not account for the transient motion, due to the simplifying assumptions of steady plug flow. Computational fluid dynamics (CFD) can assist in expanding the understanding of complex flows in bubble columns by characterising the flow phenomena for many geometrical configurations. Therefore, CFD has a role in the education of chemical and biochemical engineers, providing the examples of flow phenomena that many engineers may not experience, even through experimentation. The performance of the mixture model was investigated for three domains (plane, rectangular and cylindrical) and three flow models (laminar, k-e turbulence and the Reynolds stresses). mThis investigation raised many questions about how gas-liquid interactions are captured numerically. To answer some of these questions the analogy between thermal convection in a cavity and gas-liquid flow in bubble columns was invoked. This involved modelling the buoyant motion of air in a narrow cavity for a number of turbulence schemes. The difference in density was caused by a temperature gradient that acted across the width of the cavity. Multiple vortices were obtained when the Reynolds stresses were utilised with the addition of a basic flow profile after each time step. To implement the three-phase models an alternative mixture model was developed and compared against a commercially available mixture model for three turbulence schemes. The scheme where just the Reynolds stresses model was employed, predicted the transient motion of the fluids quite well for both mixture models. Solid-liquid and then alternative formulations of gas-liquid-solid model were compared against one another. The alternative form of the mixture model was found to perform particularly well for both gas and solid phase transport when calculating two and three-phase flow. The improvement in the solutions obtained was a result of the inclusion of the Reynolds stresses model and differences in the mixture models employed. The differences between the alternative mixture models were found in the volume fraction equation (flux and deviatoric stress tensor terms) and the viscosity formulation for the mixture phase.

Further studies on flow distribution on distillation sieve trays

The thesis describes experimental work on the possibility of using deflection baffles in conventional distillation trays as flow straightening devices, with the view of enhancing tray efficiency. The mode of operation is based on deflecting part of the liquid momentum from the centre of the tray to the segment regions in order to drive stagnating liquid at the edges forward. The first part of the work was a detailed investigation into the two-phase flow patterns produced on a conventional sieve tray having 1 mm hole size perforations. The data provide a check on some earlier work and extend the range of the existing databank, particularly to conditions more typical of industrial operation. A critical survey of data collected on trays with different hole sizes (Hine, 1990; Chambers, 1993; Fenwick, 1996; this work) showed that the hole diameter has a significant influence on the flow regime, the size of the stagnant regions and the hydraulic and mass transfer performance. Five modified tray topologies were created with different configurations of baffles and tested extensively in the 2.44 m diameter air-water pilot distillation simulator for their efficacy in achieving uniform flow across the tray and for their impact on tray loading capacity and mass transfer efficiency. Special attention was given to the calibration of the over 100 temperature probes used in measuring the water temperature across the tray on which the heat and mass transfer analogy is based. In addition to normal tray capacity experiments, higher weir load experiments were conducted using a 'half-tray' mode in order to extend the range of data to conditions more typical of industrial operation. The modified trays show superior flow characteristics compared to the conventional tray in terms of the ability to replenish the zones of exceptionally low temperatures and high residence times at the edges of the tray, to lower the bulk liquid gradient and to achieve a more uniform flow across the tray. These superior flow abilities, however, tend to diminish with increasing weir load because of the increasing tendency for the liquid to jump over the barriers instead of flowing over them. The modified tray topologies showed no tendency to cause undue limitation to tray loading capacity. Although the improvement in the efficiency of a single tray over that of the conventional tray was moderate and in some cases marginal, the multiplier effect in a multiple tray column situation would be significant (Porter et al., 1972). These results are in good agreement with an associated CFD studies (Fischer, 1999) carried out by partners in the Advanced Studies in Distillation consortium. It is concluded that deflection baffles can be used in a conventional distillation sieve tray to achieve better liquid flow distribution and obtain enhanced mass transfer efficiency, without undermining the tray loading capacity. Unlike any other controlled-flow tray whose mechanical complexity impose stringent manufacturing and installation tolerances, the baffled-tray models are simple to design, manufacture and install and thus provide an economic method of retrofitting badly performing sieve trays both in terms of downtime and fabrication. NOTE APPENDICES 2-5 ARE ON A SEPARATE FLOPPY DISK ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY WITH PRIOR ARRANGEMENT

Novel hydrogel polymers

Hydrogels may be conveniently described as hydrophilic polymers that are swollen by, but do not dissolve in water. In this work a series of copolymer hydrogels and semi-interpenetrating polymer networks based on the monomers 2-hydroxyethyl methacrylate, N-vinyl pyrrolidone and N'N' dimethyl acrylamide, together with some less hydrophilic hydroxyalkyl acrylates and methacrylates have been synthesised. Variations in structure and composition have been correlated both with the total equilibrium water content of the resultant hydrogel and with the more detailed water binding behaviour, as revealed by differential scanning calorimetry studies. The water binding characteristics of the hydrogels were found to be primarily a function of the water structuring groups present in gel. The water binding abilities of these groups were, however, modified by steric effects. The mechanical properties of the hydrogels were also investigated. These were found to be dependent on both the polymer composition and the amount and nature of the water present in the gels. In biological systems, composite formation provides a means of producing strong, high water content materials. As an analogy with these systems hydrogel composites were prepared. In an initial study of these materials the water binding and mechanical properties of semi-interpenetrating polymer networks of N'N'dimethyl acrylamide with cellulosic type materials, with polyurethanes and with ester containing polymers were examined. A preliminary investigation of surface properties of both the copolymers and semi-interpenetrating polymer networks has been completed, using both contact angle measurements and anchorage dependent fibroblast cells. Measurable differences in surface properties attributable to structural variations in the polymers were detected by droplet techniques in the dehydrated state. However, in the hydrated state these differences were masked by the water in the gels. The use of cells enabled the underlying differences to be probed and the nature of the water structuring group was again found to be the dominant factor.

The separation of hydrogen and carbon monoxide using polymer membranes

This work studies the development of polymer membranes for the separation of hydrogen and carbon monoxide from a syngas produced by the partial oxidation of natural gas. The CO product is then used for the large scale manufacture of acetic acid by reaction with methanol. A method of economic evaluation has been developed for the process as a whole and a comparison is made between separation of the H2/CO mixture by a membrane system and the conventional method of cryogenic distillation. Costs are based on bids obtained from suppliers for several different specifications for the purity of the CO fed to the acetic acid reactor. When the purity of the CO is set at that obtained by cryogenic distillation it is shown that the membrane separator offers only a marginal cost advantage. Cost parameters for the membrane separation systems have been defined in terms of effective selectivity and cost permeability. These new parameters, obtained from an analysis of the bids, are then used in a procedure which defines the optimum degree of separation and recovery of carbon monoxide for a minimum cost of manufacture of acetic acid. It is shown that a significant cost reduction is achieved with a membrane separator at the optimum process conditions. A method of "targeting" the properties of new membranes has been developed. This involves defining the properties for new (hypothetical -yet to be developed) membranes such that their use for the hydrogen/carbon monoxide separation will produce a reduced cost of acetic acid manufacture. The use of the targeting method is illustrated in the development of new membranes for the separation of hydrogen and carbon monoxide. The selection of polymeric materials for new membranes is based on molecular design methods which predict the polymer properties from the molecular groups making up the polymer molecule. Two approaches have been used. One method develops the analogy between gas solubility in liquids and that in polymers. The UNIFAC group contribution method is then used to predict gas solubility in liquids. In the second method the polymer Permachor number, developed by Salame, has been correlated with hydrogen and carbon monoxide permeabilities. These correlations are used to predict the permeabilities of gases through polymers. Materials have been tested for hydrogen and carbon monoxide permeabilities and improvements in expected economic performance have been achieved.

Effect of liquid flow patterns on distillation trays

The thesis describes experimental work on sieve trays in an air-water simulator, 2.44 m in diameter. The liquid flow pattern, for flowrates similar to those used in commercial scale distillation, was observed experimentally by water cooling experiments, in which the temperature of the water is measured at over 100 positions over the tray area. The water is cooled by the rising air which is forced through the tray. A heat and mass transfer analogy is drawn whereby the water temperature is mapped to liquid concentration in mass transfer, and the water temperature profiles reveal how liquid channelling may reduce the tray efficiency. The first experiment was to observe the flow of water only over an unperforated tray. With the exception of very low weir loads, the flow separated at the ends of the inlet downcomer. This caused liquid to flow straight across the tray between the downcomers and large circulating regions to be formed in the side regions of the tray. The effect of the air crossflow on the flow pattern was then observed on a sieve tray of 10% free area with 1 mm diameter holes (such as is used in cryogenic distillation). The flow patterns developed on the tray were similar to those produced with water only on the unperforated tray, but at low weir loads the air crossflow prevented separation of the water flow and the associated circulating regions. At higher weir loads, liquid channelling down the centre of the tray and circulation in the side regions occurred. The percentage of the tray occupied by circulating liquid depended upon the velocity of the liquid entering the tray, which was set by the weir load and size of the gap under the inlet downcomer. The water cooling experiments showed that the temperature of the water in a circulating region is much lower than in other parts of the tray, indicating that the driving force for heat transfer is reduced. In a column section where trays (and circulating areas) are mounted on top of each other, the circulating regions will cause air (or vapour) passing through them to have a reduced change in temperature or concentration leading a loss in tray efficiency.