Doctor of Philosophy / Department of Chemical Engineering / James H. Edgar / The sublimation recondensation crystal growth of aluminum nitride, titanium nitride, and yttrium nitride were explored experimentally and theoretically. Single crystals of these nitrides are potentially suitable as substrates for AlGaInN epitaxial layers, which are employed in ultraviolet optoelectronics including UV light-emitting diodes and laser diodes, and high power high frequency electronic device applications. A thermodynamic analysis was applied to the sublimation crystal growth of aluminum nitride to predict impurities transport (oxygen, carbon, and hydrogen) and to study the aspects of impurities incorporation for different growth conditions. A source purification procedure was established to minimize the impurity concentration and avoid degradation of the crystal’s properties. More than 98% of the oxygen, 99.9% of hydrogen and 90% of carbon originally in the source was removed. The AlN crystal growth process was explored in two ways: self-seeded growth with spontaneous nucleation directly on the crucible lid or foil, and seeded growth on SiC and AlN. The oxygen concentration was 2 ~ 4 x 1018cm-3, as measured by secondary ion mass spectroscopy in the crystals produced by self-seeded growth. Crystals grown from AlN seeds have visible grain size expansion. The initial AlN growth on SiC at a low temperature range (1400°C ~1600°C) was examined to understand the factors controlling nucleation. Crystals were obtained from c-plane on-axis and off-axis, Si-face and C-face, as well as m-plane SiC seeds. In all cases, crystal growth was fastest perpendicular to the c-axis. The growth rate dependence on temperature and pressure was determined for TiN and YN crystals, and their activation energies were 775.8±29.8kJ/mol and 467.1±21.7kJ/mol respectively. The orientation relationship of TiN (001) || W (001) with TiN [100] || W [110], a 45o angle between TiN [100] and W [100], was seen for TiN crystals deposited on both (001) textured tungsten and randomly orientated tungsten. Xray diffraction confirmed that the YN crystals were rock-salt structure, with a lattice constant of 4.88Å. Cubic yttria was detected in YN sample from the oxidation upon its exposed to air for limited time by XRD, while non-cubic yttria was detected in YN sample for exposures more than one hour by Raman spectra.

THEREDA (Thermodynamic Reference Database) is a collaborative project, which has been addressed this challenge. The partners are Helmholtz-Zentrum Dresden-Rossendorf, Karlsruhe Institute of Technology (KIT-INE), Gesellschaft für Anlagen- und Reaktorsicherheit Braunschweig mbH (GRS), TU Bergakademie Freiberg (TUBAF) and AF-Consult Switzerland AG (Baden, Switzerland). The aim of the project is the establishment of a consistent and quality assured database for all safety relevant elements, temperature and pressure ranges, with its focus on saline systems. This implied the use of the Pitzer approach to compute activity coefficients suitable for such conditions. Data access is possible via commonly available internet browsers under the address http://www.thereda.de. One part of the project - the data collection and evaluation for uranium – was a task of the Helmholtz-Zentrum Dresden-Rossendorf. The aquatic chemistry and thermodynamics of U(VI) and U(IV) is of great importance for geochemical modelling in repository-relevant systems. The OECD/NEA Thermochemical Database (NEA TDB) compilation is the major source for thermodynamic data of the aqueous and solid uranium species, even though this data selection does not utilize the Pitzer model for the ionic strength effect correction. As a result of the very stringent quality demands, the NEA TDB is rather restrictive and therefore incomplete for extensive modelling calculations of real systems. Therefore, the THEREDA compilation includes additional thermodynamic data of solid secondary phases formed in the waste material, the backfill and the host rock, though falling into quality assessment (QA) categories of lower accuracy. The data review process prefers log K values from solubility experiments (if available) to those calculated from thermochemical data.

Poor aqueous solubilities of drug candidates limit the biopharmaceutical usefulness in either oral or parenteral dosage forms. Lipid assemblies, such as micelles, may provide a means of enhancing solubility. Despite their usefulness, little is known about the means by which micelles accomplish this result. The goal of the current dissertation is to provide the molecular level understanding of the mechanism by which simple micelle systems solubilize drugs. Specifically, the location, orientation and amount of the drug molecules in micelle systems are the focuses of the work. Three series of model drugs, steroids, benzodiazepines and parabens, in three surfactant systems with anionic, cationic and neutral hydrophilic headgroups were studied. Solubilization power of each micelle system for each model drug was determined by equilibrium solubility. The observed strong surface activities of model drug at hydrocarbon/water interface and the ability of the drugs to compete with surfactants for the model oil/water interface lend support to the hypothesis that drug molecules are mainly solubilized in the interfacial region of the micelles. A surface-localized thermodynamic model that considered the surfactant-drug competition at micelle surface was successfully applied to predict the micelle/water partitioning coefficients. The predictions were made without the use of adjustable parameters in the case of both dilute and concentrated solutions. The orientation of drug at micelle surface was determined by matching calculated occupied areas by solutes at oil/water interface using molecular modeling method to the experimental values. To look into the micro-structure of micelles, twodimensional and diffusion (or PGSE) NMR techniques were employed to detect the specific drug-surfactant interactions and the micelle sizes influenced by model drugs and electrolytes.

The measurement of vapour-liquid equilibrium data is extremely important as such data are crucial for the accurate design, simulation and optimization of the majority of separation processes, including distillation, extraction and absorption. This study involved the measurement of vapour-liquid equilibrium data, using a modified version of the static total pressure apparatus designed within the Thermodynamics Research Unit by J.D. Raal and commissioned by Motchelaho, (Motchelaho, 2006 and Raal et al., 2011). This apparatus provides a very simple and accurate means of obtaining P-x data using only isothermal total pressure and overall composition (z) measurements. Phase sampling is not required. Phase equilibrium measurement procedures using this type of apparatus are often tedious, protracted and repetitive. It is therefore useful and realizable in the rapidly advancing digital age, to incorporate computer-aided operation, to decrease the man hours required to perform such measurements. The central objective of this work was to develop and implement a control scheme, to fully automate the original static total pressure apparatus of Raal et al. (2011). The scheme incorporates several pressure feedback closed loops, to execute process step re-initialization, valve positioning and motion control in a stepwise fashion. High resolution stepper motors were used to engage the dispensers, as they provided a very accurate method of regulating the introduction of precise desired volumes of components into the cell. Once executed, the control scheme requires approximately two days to produce a single forty data points (P-x) isotherm, and minimizes human intervention to two to three hours. In addition to automation, the apparatus was modified to perform moderate pressure measurements up to 1.5 MPa. Vapour-liquid equilibrium test measurements were performed using both the manual and automated operating modes to validate the operability and reproducibility of the apparatus. The test systems measured include the water (1) + propan-1-ol (2) system at 313.15 K and the n-hexane (1) + butan- 2-ol system at 329.15 K. Phase equilibrium data of binary systems, containing the solvent morpholine-4-carbaldehyde (NFM) was then measured. The availability of vapour-liquid equilibrium data for binary systems containing NFM is limited in the literature. The new systems measured include: n-hexane (1) + NFM (2) at 343.15, 363.15 and 393.15 K, as well as n-heptane (1) + NFM (2) at 343.15, 363.15 and 393.15 K. The modified apparatus is quite efficient as combinations of the slightly volatile NFM with highly volatile alkane constituents were easily and accurately measured. The apparatus also allows for accurate vapour-liquid equilibrium measurements in the dilute composition regions. A standard uncertainty in the equilibrium pressure reading, within the 0 to 100 kPa range was calculated to be 0.106 kPa, and 1.06 kPa for the 100 to 1000 kPa pressure range. A standard uncertainty in the equilibrium temperature of 0.05 K was calculated. The isothermal data obtained were modelled using the combined (-) method described by Barker (1953). This involved the calculation of binary interaction parameters, by fitting the data to various thermodynamic models. The virial equation of state with the Hayden-O’Connell (1975) and modified Tsonopoulos (Long et al., 2004) second virial coefficient correlations were used in this work to account for vapour phase non-ideality. The Wilson (1964), NRTL (Renon and Prausnitz, 1968), Tsuboka-Katayama-Wilson (1975) and modified Universal Quasi-Chemical (Anderson and Prausnitz, 1978) activity coefficient models were used to account for the liquid phase non-ideality. A stability analysis was carried out on all the new systems measured to ensure that two-liquid phase formation did not occur in the measured temperature range. A model-free method based on the numerical integration of the coexistence equation was also used to determine the vapour phase compositions and activity coefficients from the measured P-z data. These results compare well with the results obtained by the model-dependent method. The infinite dilution activity coefficients for the systems under consideration were determined by the method of Maher and Smith (1979b), and by suitable extrapolation methods. Excess enthalpy and excess entropy data were calculated for the systems measured, using the Gibbs-Helmholtz equation in conjunction with the fundamental excess property relation. / Thesis (M.Sc.Eng.)-University of KwaZulu-Natal, Durban, 2012.

Salts in solution should be removed by desalination techniques to prevent equipment fouling and corrosion. Common desalination technologies are energy intensive such as Multi Stage Flash (MSF) distillation which requires 14.5 J/m3 (Ribeiro. J, 1996) of energy. Desalination technologies produce purified water and a concentrated salt solution, where the salt concentration is dependent on the desalination technology used. This work investigates gas hydrate technology as a possible desalination technology. Hydrates are composed of guest molecules and host molecules. Guest molecules may be in the form of a liquid or gas. During hydrate formation, host molecules, water, form a cage enclosing the guest molecule. Common hydrate formers or guest molecules such as; methane, ethane, propane and carbon dioxide are currently being investigated in literature, for use in gas hydrate desalination technology. Common hydrate formers form hydrates at low temperatures; below 288 K and high pressures; above 2 MPa. To increase the temperature and reduce the pressure at which gas hydrates form, commercially available hydrofluorocarbon hydrate formers such as R14, R32, R116, R134a, R152a, R218, R404a, R407c, R410a and R507 are preliminarily investigated in this work. The criteria for choosing the most suitable fluorine-based formers require the former to be: environmentally acceptable where it is approved by the Montreal Protocol; non-toxic where it has a low acute toxicity; non-flammable; chemically stable; a structure II hydrate to simplify the washing process; available in commercial quantities; low cost in comparison to other hydrate formers; compatible with standard materials and contain a high critical point for a large heat of vaporisation (McCormack and Andersen, 1995). Taking all these criteria into account, R134a was chosen for further investigation as a possible hydrate former. In this work, hydrate-liquid-vapour phase equilibrium measurements are conducted using the isochoric method with a static high pressure stainless steel equilibrium cell. The Combined Standard Uncertainty for the 0-1 MPa pressure transducer, 0-10 MPa pressure transducer and the Pt100 temperature probes are ±0.64 MPa, ±5.00 MPa and ±0.09 K respectively. Vapour pressure measurements for Hydrofluoropropyleneoxide, CO2, R22 and R134a were measured to verify the pressure and temperature calibrations. Hydrate test systems for R22 (1) + water (2) and R134a (1) + water (2) were measured to verify calibrations, equipment and procedures. New systems measured included R134a (1) + water (2) + {5wt%, 10wt% or 15wt%} NaCl (3). For the system R134 (1) + water (2) at 281 K the dissociation pressure is 0.269 MPa. However, addition of NaCl to the system resulted in a shift of the HVL equilibrium phase boundary to lower temperatures or higher pressures. The average shift in temperature between the system R134a (1) + water (2) containing no salt and the systems containing {5, 10 and 15} wt% NaCl are -1.9K, -4.8K and -8.1K respectively. In this work, the measured systems were modelled using two methods of approach. The first method is where hydrofluorocarbon hydrate former solubility is included, (Parrish et al., 1972) and the second is where hydrofluorocarbon hydrate former solubility is ignored, (Eslamimanesh et al., 2011). From these models, it is found that hydrofluorocarbon solubility could not be neglected. In this work, the hydrate phase was modelled using modifications of the van der Waals and Platteeuw model, (Parrish et al., 1972). The liquid and vapour phases are modelled using the Peng- Robinson equation of state with classical mixing rules (Peng, 1976). The electrolyte component is modelled using the Aasberg-Peterson model (Aasberg-Petersen et al., 1991) modified by Tohidi (Tohidi et al., 1995). The percent absolute average deviation (%AAD) for the systems, which includes solubility, is 0.41 for R22 (1) + water (2) and 0.33 for R134a (1) + water (2). For the system R134a (1) + water (2) + {5 wt%, 10 wt% or 15 wt%} NaCl (3) the % AAD is 5.14. Using the hydrate former, R134a, is insufficient to ensure gas hydrate technology is competitive with other desalination technologies. Hydrate dissociation temperature should be increased and pressure decreased further to ambient conditions. As evident in literature, promoters, such as cyclopentane, are recommended to be added to the system to shift the HLV equilibrium phase boundary as close to ambient conditions as possible. / Thesis (M.Sc.Eng.)-University of KwaZulu-Natal, Durban, 2011.

The first project in this work used classical molecular dynamics to study the ice recrystallization inhibition potential of a series of carbohydrates and alcochols, using the hydration index, partial molar volumes and isothermal compressibilities as parameters for measuring their cryogenic efficacy. Unfortunately, after 8 months of testing, this work demonstrates that the accuracy and precision of the density extracted from simulations is not sufficient in providing accurate partial molar volumes. As a result, this work clearly demonstrates that current classical molecular dynamics technology cannot probe the volumetric properties of interest with sufficient accuracy to aid in the research and development of novel cryoprotectants.The second project in this work used molecular dynamics simulations to evaluate the Gibbs free energy change of substituting CO2 in sI clathrate hydrates by N2,CH4, SO2 and H2S flue gas impurities under conditions proposed for CO2 capture (273 K, 10 bar). Our results demonstrate that CO2 substitutions by N2 in the small sI cages were thermodynamically favored. This substitution is problematic in terms of efficient CO2 capture, since the small cages make up 25% of the sI clathrate cages, therefore a significant amount of energy could be spent on removing N2 from the flue gas rather than CO2. The thermodynamics of CO2 substitution by CH4, SO2 and H2S in sI clathrate hydrates was also examined. The substitution of CO2 by these gases in both the small and large cages were determined to be favorable. This suggests that these gases may also disrupt the CO2 capture by sI clathrate hydrates if they are present in large concentrations in the combustion flue stream. Similar substitution thermodynamics at 200 K and 10 bar were also studied. With one exception, we found that the substitution free energies do not significantly change and do not alter the sign of thermodynamics. Thus, using a lower capture temperature does not significantly change the substitution free energies and their implications for CO2 capture by sI clathrate hydrates.

The principal theory concerning the origin of the elements heavier than the Fe-peak, such as Ba, strongly suggest that for old, metal-poor environments, the rapid (r-) process is the most likely path taken in their synthesis, while the slow (s-) process becomes more substantial in younger, more metal-rich stellar populations. In this work I test this theory by evaluating the isotope ratios of Ba. It is understood that Ba consists of seven stable isotopes, five of which are synthesised by the two neutron-capture processes. The two odd isotopes, 135,137Ba, as well as 138Ba are synthesised via both the r- and s-processes while two of the even isotopes, 134,136Ba are synthesised via the s-process only. The relative contribution of the r- and s-process to these isotopes can be understood via nucleosynthesis calculations and is described using the parameter fodd, where fodd = [N (135Ba) + N (137Ba)] /N (Ba). Low values of fodd (~0.11) indicate an s-process regime, while high values of fodd (~0.46) indicate an r-process regime. In the Ba II 4554 A line the even isotopes lie close to the line centre, while the odd isotopes, which are hyperfine split because of their non-zero nuclear spin, lie in the wings of the line. From an analysis of the line profile shape, one can determine whether Ba has been synthesised primarily through the r-process or s-process; a broad, asymmetric line would indicate a high r-process contribution, while a line with a deeper core and shallower wings would indicate a high s-process contribution. Using the radiative transfer code ATLAS, which assumes local thermodynamic equilibrium (LTE) and employs 1-dimensional (1D) KURUCZ06 model atmospheres, I synthesised line profiles for six metal-poor stars: HD140283, HD122563, HD88609, HD84937, BD-04 3208 and BD+26 3578 - for a range of isotope ratios. All six are of sufficiently low metallicity that Ba was expected to have an r-process origin. These were fit to high resolution (R\equiv \lamda/\Delta\lamda = 90 000 - 95 000), high signal-to-noise to the Ba II 4554 A line which has multiple components. In the first test, synthetic spectra were computed using the non local thermodynamic equilibrium (NLTE) radiative transfer code MULTI. The synthetic line profiles were fit to a number of lines in HD140283. Although this technique might have improved the fit in the line core, it was found that such a treatment did not improve upon fitting errors associated with the best fit 1D LTE synthetic profiles. The second test used a 3-dimensional (3D) radiative transfer code (LINFOR3D) that employed 3D, time-dependent atmospheres produced with CO5BOLD. The 3D synthetic pro les were fit to a selection of Fe lines and improvements over the poor fits produced by the 1D LTE synthesis were seen. It was found that the 3D synthesis could almost completely reproduce the line asymmetries seen in the observed stellar spectrum. This result suggests that further work to refine the 3D calculations and synthesis code would be valuable.

<p>Thermodynamic stability measurements on proteins and protein-ligand complexes can offer insights not only into the fundamental properties of protein folding reactions and protein functions, but also into the development of protein-directed therapeutic agents to combat disease. Conventional calorimetric or spectroscopic approaches for measuring protein stability typically require large amounts of purified protein. This requirement has precluded their use in proteomic applications. Stability of Proteins from Rates of Oxidation (SPROX) is a recently developed mass spectrometry-based approach for proteome-wide thermodynamic stability analysis. Since the proteomic coverage of SPROX is fundamentally limited by the detection of methionine-containing peptides, the use of tryptophan-containing peptides was investigated in this dissertation. A new SPROX-like protocol was developed that measured protein folding free energies using the denaturant dependence of the rate at which globally protected tryptophan and methionine residues are modified with dimethyl (2-hydroxyl-5-nitrobenzyl) sulfonium bromide and hydrogen peroxide, respectively. This so-called Hybrid protocol was applied to proteins in yeast and MCF-7 cell lysates and achieved a ~50% increase in proteomic coverage compared to probing only methionine-containing peptides. Subsequently, the Hybrid protocol was successfully utilized to identify and quantify both known and novel protein-ligand interactions in cell lysates. The ligands under study included the well-known Hsp90 inhibitor geldanamycin and the less well-understood omeprazole sulfide that inhibits liver-stage malaria. In addition to protein-small molecule interactions, protein-protein interactions involving Puf6 were investigated using the SPROX technique in comparative thermodynamic analyses performed on wild-type and Puf6-deletion yeast strains. A total of 39 proteins were detected as Puf6 targets and 36 of these targets were previously unknown to interact with Puf6. Finally, to facilitate the SPROX/Hybrid data analysis process and minimize human errors, a Bayesian algorithm was developed for transition midpoint assignment. In summary, the work in this dissertation expanded the scope of SPROX and evaluated the use of SPROX/Hybrid protocols for characterizing protein-ligand interactions in complex biological mixtures.</p> / Dissertation

Ce travail a pour but l’acquisition de connaissances fondamentales dans les équilibres de phases et les transformations de phases des alliages ternaires contenant du fer, du nickel et du cuivre.Cette étude est composée de trois parties : la première est consacrée aux équilibres de phase entre600°C et 1000°C dans le coin riche en fer du système ternaire Cu-Fe-Ni, la seconde concerne l’analyse des transformations de phases dans des alliages modèles suivant deux vitesses de refroidissement et la troisième étudie les alliages industriels.Cette approche expérimentale est appuyée par des calculs thermodynamiques. Une comparaison est effectuée avec les données issues de la littérature. Les sections isothermes à 600°C, 800°C et1000°C ont été reconfirmées. Les domaines triphasés αFe/γFe/γCu ont été déterminés à 600°C et800°C et semblent décalés vers de plus fortes teneurs en Ni. La lacune de miscibilité est moins étendue que celle calculée et ce particulièrement du côté du binaire Cu-Ni.Les transformations de phase γ→α apparaissant lors du refroidissement d’alliages Fe-xCu-10Ni et Fe-10Cu-xNi (0<x<15 % en masse) ont été étudiées. Elles sont observées dans tous les alliages refroidis lentement par dilatométrie sauf dans l’alliage binaire Fe-Cu. Les températures de transformation sont systématiquement inférieures aux températures d’équilibre et cet écart augmente lorsque les teneurs en Ni et Cu augmentent. La formation de la ferrite bainitique est favorisée par l’addition de Cu dans les alliages trempés Fe-xCu-10Ni et par l’addition de Ni dans les alliages Fe-10Cu-xNi tandis que la ferrite massive se forme préférentiellement lors du refroidissement lent. Une diminution des températures de transformation se produit lorsque la taille de grains est affinée et est reliée au processus d’accommodation plastique liée à la transformation.La composition et la vitesse de refroidissement jouent un rôle sur la dureté de ces alliages via la présence de Ni en solution solide ainsi que le nombre et la taille des précipités de Cu. / This work aims in getting fundamental knowledge of phase equilibria and microstructures of ternary alloys containing copper, nickel and iron. The thesis is composed of three parts: a first part is devoted to phase equilibria in the Fe-rich corner of the ternary Cu-Fe-Ni system between 600°C and 1000°C, a second part is devoted to the microstructures in these materials for model alloys, for two cooling rates and a third part is devoted to industrial alloys.Phase equilibria of the system are investigated in the range 600-1000°C using diffusion multiples in conjunction with selected equilibrated alloys. This experimental approach is supplemented by thermodynamic calculations. A comparison is drawn with data reported in the literature. The isothermal sections at 600°C, 800°C and 1000°C have been reconfirmed. The three-phase regions αFe/γFe/γCu are determined at 600°C and 800°C and appear shifted to a higher Ni content. The miscibility gap is narrower than the calculated, particularly on the Cu-Ni binary sides.Austenite/ferrite phase transformations occurring on cooling in Fe-xCu-10Ni and Fe-10Cu-xNi alloys,0<x<15 (mass%), have been studied. The influence of copper and nickel additions and of the cooling rate on the microstructure is discussed. Metastable transformations in slowly cooled alloys have been detected by dilatometry in all alloys except in the binary Fe-10Cu alloy: all the cooling transformation temperatures are systematically lower than the equilibrium temperature and the Ni and Cu addition decrease this transformation temperature. The formation of bainitic ferrite is favoured by copper additions in quenched Fe-xCu-10Ni alloys and by Ni addition in Fe-10Cu-xNi alloys while massive ferrite form preferentially during slow cooling. A decrease of the transformation temperatures is recorded when the alloys have a finer grain size likely related to plastic accommodation processes.The composition and cooling rate influence the hardness of the alloys mostly dependant on the amount of Ni in solid solution and on the rate and size of Cu precipitates.

Submitted in fulfillment of the requirements of the degree of Master of Engineering, Durban University of Technology, Durban, South Africa, 2016. / Separation processes in the chemical process industries is dependent on the science of chemical thermodynamics. In the field of chemical separation process engineering, phase equilibrium is a primary area of interest. This is due to separation processes such as distillation and extraction which involves the contacting of different phases for effective separation. The focal point of this research project is the measurement and modeling of binary vapour-liquid equilibrium (VLE) phase data of systems containing ether-alcohol organic compounds. The VLE data were measured with the use of the modified apparatus of Raal and Mühlbauer, (1998). The systems of interest for this research arose from an industrial demand for VLE data for systems containing ether-alcohol organic compounds. This gave rise to the experimental VLE data isotherms being measured for the following binary systems: a) Methyl tert-butyl ether (1) + 1-pentanol (2) at 317.15 and 327.15 K b) Methyl tert-butyl ether (1) + 2, 2, 4-trimethylpentane (2) at 307.15, 317.15 and 327.15K c) 2, 2, 4-Trimethylpentane (1) + 1-pentanol (2) at 350.15, 360.15 and 370.15K d) Diisopropyl ether (1) + 2,2,4-trimethylpentane (2) at 320.15, 330.15 and 340.15K e) Diisopropyl ether (1) + 1-propanol (2) at 320.15, 330.15 and 340.15K f) Diisopropyl ether (1) + 2-butanol (2) at 320.15, 330.15 and 340.15K The data for all the measured binary systems investigated at these temperatures are currently not available in the open source literature found on the internet and in library text resources. The systems were not measured at the same temperatures because certain system isotherm temperatures correlate to a pressures above 1 bar. This pressure of 1 bar is the maximum operating pressure specification of the VLE apparatus used in this project. The experimental VLE data were correlated for model parameters for both the  and methods. For the method, the fugacity coefficients (vapour-phase non-idealities) were tabulated using the virial equation of state and the Hayden-O’Connell correlation (1975); chemical theory and the Nothnagel et al. (1973) correlation method. The activity coefficients (liquid phase non-idealities) were calculated using three local-composition based activity coefficients models: the Wilson (1964) model, the NRTL model (Renon and Prausnitz, 1968); and the UNIQUAC model (Abrams and Prausnitz, 1975). Regarding the direct method, the Soave-Redlich-Kwong (Redlich and Kwong, 1949) and Peng-Robinson (1976) equations of state ii were used with the temperature dependent alpha-function (α) of Mathias and Copeman (1983) with the Wong-Sandler (1992) mixing rule. Thermodynamic consistency testing, which presents an indication of the quality and reliability of the data, was also performed for all the experimental VLE data. All the systems measured showed good thermodynamic consistency for the point test of Van Ness et al. (1973) - the consistency test of choice for this research. This however, was based on the model chosen for the data regression of a particular system. Therefore, the combined method of VLE reduction produced the most favourable results for the NRTL and Wilson models. / M