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.

Multicomponent kinetic spectrophotometric determination of pefloxacin and norfloxacin in pharmaceutical preparations and human plasma samples with the aid of chemometrics

A spectrophotometric method for the simultaneous determination of the important pharmaceuticals, pefloxacin and its structurally similar metabolite, norfloxacin, is described for the first time. The analysis is based on the monitoring of a kinetic spectrophotometric reaction of the two analytes with potassium permanganate as the oxidant. The measurement of the reaction process followed the absorbance decrease of potassium permanganate at 526 nm, and the accompanying increase of the product, potassium manganate, at 608 nm. It was essential to use multivariate calibrations to overcome severe spectral overlaps and similarities in reaction kinetics. Calibration curves for the individual analytes showed linear relationships over the concentration ranges of 1.0–11.5 mg L−1 at 526 and 608 nm for pefloxacin, and 0.15–1.8 mg L−1 at 526 and 608 nm for norfloxacin. Various multivariate calibration models were applied, at the two analytical wavelengths, for the simultaneous prediction of the two analytes including classical least squares (CLS), principal component regression (PCR), partial least squares (PLS), radial basis function-artificial neural network (RBF-ANN) and principal component-radial basis function-artificial neural network (PC-RBF-ANN). PLS and PC-RBF-ANN calibrations with the data collected at 526 nm, were the preferred methods—%RPET not, vert, similar 5, and LODs for pefloxacin and norfloxacin of 0.36 and 0.06 mg L−1, respectively. Then, the proposed method was applied successfully for the simultaneous determination of pefloxacin and norfloxacin present in pharmaceutical and human plasma samples. The results compared well with those from the alternative analysis by HPLC.

Look before you leap : a confidence-based method for selecting species criticality while avoiding negative populations in τ-leaping

The stochastic simulation algorithm was introduced by Gillespie and in a different form by Kurtz. There have been many attempts at accelerating the algorithm without deviating from the behavior of the simulated system. The crux of the explicit τ-leaping procedure is the use of Poisson random variables to approximate the number of occurrences of each type of reaction event during a carefully selected time period, τ. This method is acceptable providing the leap condition, that no propensity function changes “significantly” during any time-step, is met. Using this method there is a possibility that species numbers can, artificially, become negative. Several recent papers have demonstrated methods that avoid this situation. One such method classifies, as critical, those reactions in danger of sending species populations negative. At most, one of these critical reactions is allowed to occur in the next time-step. We argue that the criticality of a reactant species and its dependent reaction channels should be related to the probability of the species number becoming negative. This way only reactions that, if fired, produce a high probability of driving a reactant population negative are labeled critical. The number of firings of more reaction channels can be approximated using Poisson random variables thus speeding up the simulation while maintaining the accuracy. In implementing this revised method of criticality selection we make use of the probability distribution from which the random variable describing the change in species number is drawn. We give several numerical examples to demonstrate the effectiveness of our new method.

An experimental study of a reactive plume in grid turbulence

Experimental results for a reactive non-buoyant plume of nitric oxide (NO) in a turbulent grid flow doped with ozone (O3) are presented. The Damkohler number (Nd) for the experiment is of order unity indicating the turbulence and chemistry have similar timescales and both affect the chemical reaction rate. Continuous measurements of two components of velocity using hot-wire anemometry and the two reactants using chemiluminescent analysers have been made. A spatial resolution for the reactants of four Kolmogorov scales has been possible because of the novel design of the experiment. Measurements at this resolution for a reactive plume are not found in the literature. The experiment has been conducted relatively close to the grid in the region where self-similarity of the plume has not yet developed. Statistics of a conserved scalar, deduced from both reactive and non-reactive scalars by conserved scalar theory, are used to establish the mixing field of the plume, which is found to be consistent with theoretical considerations and with those found by other investigators in non-reative flows. Where appropriate the reactive species means and higher moments, probability density functions, joint statistics and spectra are compared with their respective frozen, equilibrium and reaction-dominated limits deduced from conserved scalar theory. The theoretical limits bracket reactive scalar statistics where this should be so according to conserved scalar theory. Both reactants approach their equilibrium limits with greater distance downstream. In the region of measurement, the plume reactant behaves as the reactant not in excess and the ambient reactant behaves as the reactant in excess. The reactant covariance lies outside its frozen and equilibrium limits for this value of Vd. The reaction rate closure of Toor (1969) is compared with the measured reaction rate. The gradient model is used to obtain turbulent diffusivities from turbulent fluxes. Diffusivity of a non-reactive scalar is found to be close to that measured in non-reactive flows by others.

Mathematical modelling of the drying of sol gel microspheres

This thesis presents a mathematical model of the evaporation of colloidal sol droplets suspended within an atmosphere consisting of water vapour and air. The main purpose of this work is to investigate the causes of the morphologies arising within the powder collected from a spray dryer into which the precursor sol for Synroc™ is sprayed. The morphology is of significant importance for the application to storage of High Level Liquid Nuclear Waste. We begin by developing a model describing the evaporation of pure liquid droplets in order to establish a framework. This model is developed through the use of continuum mechanics and thermodynamic theory, and we focus on the specific case of pure water droplets. We establish a model considering a pure water vapour atmosphere, and then expand this model to account for the presence of an atmospheric gas such as air. We model colloidal particle-particle interactions and interactions between colloid and electrolyte using DLVO Theory and reaction kinetics, then incorporate these interactions into an expression for net interaction energy of a single particle with all other particles within the droplet. We account for the flow of material due to diffusion, advection, and interaction between species, and expand the pure liquid droplet models to account for the presence of these species. In addition, the process of colloidal agglomeration is modelled. To obtain solutions for our models, we develop a numerical algorithm based on the Control Volume method. To promote numerical stability, we formulate a new method of convergence acceleration. The results of a MATLAB™ code developed from this algorithm are compared with experimental data collected for the purposes of validation, and further analysis is done on the sensitivity of the solution to various controlling parameters.

Modeling the co-precipitation of silica and calcium oxalate in sugar solutions

Solution chemistry plays a significant role in the rate and type of foulant formed on heated industrial surfaces. This paper describes the effect of sucrose, silica (SiO2), Ca2+ and Mg2+ ions, and trans-aconitic acid on the kinetics and solubility of SiO2 and calcium oxalate monohydrate (COM) in mixed salt solutions containing sucrose and refines models previously proposed. The developed SiO2 models show that sucrose and SiO2 concentrations are the main parameters that determine apparent order (n) and apparent rate of reaction (k) and SiO2 solubility over a 24 h period. The calcium oxalate solubility model shows that while increasing [Mg2+] increases COM solubility, the reverse is so with increasing sucrose concentrations. The role of solution species on COM crystal habit is discussed and the appearance of the uncommon (001) face is explained.

Understanding mild acid pretreatment of sugarcane bagasse through particle scale modeling

Sugarcane bagasse is an abundant and sustainable resource, generated as a by-product of sugarcane milling. The cellulosic material within bagasse can be broken down into glucose molecules and fermented to produce ethanol, making it a promising feedstock for biofuel production. Mild acid pretreatment hydrolyses the hemicellulosic component of biomass, thus allowing enzymes greater access to the cellulosic substrate during saccharification. A particle-scale mathematical model describing the mild acid pretreatment of sugarcane bagasse has been developed, using a volume averaged framework. Discrete population-balance equations are used to characterise the polymer degradation kinetics, and diffusive effects account for mass transport within the cell wall of the bagasse. As the fibrous material hydrolyses over time, variations in the porosity of the cell wall and the downstream effects on the reaction kinetics are accounted for using conservation of volume arguments. Non-dimensionalization of the model equations reduces the number of parameters in the system to a set of four dimensionless ratios that compare the timescales of different reaction and diffusion events. Theoretical yield curves are compared to macroscopic experimental observations from the literature and inferences are made as to constraints on these “unknown” parameters. These results enable connections to be made between experimental data and the underlying thermodynamics of acid pretreatment. Consequently, the results suggest that data-fitting techniques used to obtain kinetic parameters should be carefully applied, with prudent consideration given to the chemical and physiological processes being modeled.

Optimal catalyst thickness in titanium dioxide fixed film reactors: Mathematical modelling and experimental validation

The use of immobilised TiO2 for the purification of polluted water streams introduces the necessity to evaluate the effect of mechanisms such as the transport of pollutants from the bulk of the liquid to the catalyst surface and the transport phenomena inside the porous film. Experimental results of the effects of film thickness on the observed reaction rate for both liquid-side and support-side illumination are here compared with the predictions of a one-dimensional mathematical model of the porous photocatalytic slab. Good agreement was observed between the experimentally obtained photodegradation of phenol and its by-products, and the corresponding model predictions. The results have confirmed that an optimal catalyst thickness exists and, for the films employed here, is 5 μm. Furthermore, the modelling results have highlighted the fact that porosity, together with the intrinsic reaction kinetics are the parameters controlling the photocatalytic activity of the film. The former by influencing transport phenomena and light absorption characteristics, the latter by naturally dictating the rate of reaction.

Ultraviolet Photodissociation of the N-Methylpyridinium Ion: Action Spectroscopy and Product Characterization

The ultraviolet photodissociation of gas-phase N-methylpyridinium ions is studied at room temperature using laser photodissociation mass spectrometry and structurally diagnostic ion-molecule reaction kinetics. The C5H5N-CH3+ (m/z 94), C5H5N-CD3+ (m/z 97), and C5D5N-CH3+(m/z 99) isotopologues are investigated, and it is shown that the N-methylpyridinium ion photodissociates by the loss of methane in the 36 000 - 43 000 cm(-1) (280 - 230 nm) region. The dissociation likely occurs on the ground state surface following internal conversion from the SI state. For each isotopologue, by monitoring the photofragmentation yield as a function of photon wavenumber, a broad vibronically featured band is recorded with origin (0-0) transitions assigned at 38 130, 38 140 and 38 320 cm(-1) for C5H5N-CH3+ C5H5N-CD3+ and C5D5N-CH3+, respectively. With the aid of quantum chemical calculations (CASSCF(6,6)/aug-cc-pVDZ), most of the observed vibronic detail is assigned to two in-plane ring deformation modes. Finally, using ion-molecule reactions, the methane coproduct at m/z 78 is confirmed as a 2-pyridinylium ion.

Direct Observation of Photodissociation Products from Phenylperoxyl Radicals Isolated in the Gas Phase

Gas phase peroxyl radicals are central to our chemical understanding of combustion and atmospheric processes and are typically characterized by strong absorption in the UV (lambda(max) approximate to 240 nm). The analogous maximum absorption feature for arylperoxyl radicals is predicted to shift to the visible but has not previously been characterized nor have any photoproducts arising from this transition been identified. Here we describe the controlled synthesis and isolation in vacuo of an array of charge-substituted phenylperoxyl radicals at room temperature, including the 4-(N,N,N-trimethylammonium)methyl phenylperoxyl radical cation (4-Me3N[+]CH2-C6H4OO center dot), using linear ion-trap mass spectrometry. Photodissociation mass spectra obtained at wavelengths ranging from 310 to 500 nm reveal two major photoproduct channels corresponding to homolysis of aryl-OO and arylO-O bonds resulting in loss of O-2 and O, respectively. Combining the photodissociation yields across this spectral window produces a broad (FWHM approximate to 60 nm) but clearly resolved feature centered at lambda(max) = 403 nm (3.08 eV). The influence of the charge-tag identity and its proximity to the radical site are investigated and demonstrate no effect on the identity of the two dominant photoproduct channels. Electronic structure calculations have located the vertical (B) over tilde <- (X) over tilde transition of these substituted phenylperoxyl radicals within the experimental uncertainty and further predict the analogous transition for unsubstituted phenylperoxyl radical (C6H5OO center dot) to be 457 nm (2.71 eV), nearly 45 nm shorter than previous estimates and in good agreement with recent computational values.

Formation of the heterocumulene anion SCCCN- by a cyano migration from the radical anion of 1,2-dicyanoethylenedithiolate

The anionic heterocumulene SCCCN- was generated in the gas phase by collisional activation of the radical anion of 1,2-dicyanoethylenedithiolate. The mechanism of this reaction, as well as the structures of neutral and anionic products, was investigated by hybrid density functional theory (DFT) calculations. Dissociation to form SCCCN- and SCN is proposed to occur by a radical directed cyano migration reaction, with calculations suggesting this is the lowest energy fragmentation pathway available to the precursor anion. In contrast, the even-electron protonated 1,2-dicyanoethylenedithiolate anion fragmented by loss of HCN.

One-pot synthesis of imines from benzyl alcohol and amines on Au/ZrO2 catalyst

Imines were synthesized from benzyl alcohol and amines by using catalysts of gold nanoparticles supported on ZrO2 (Au/ZrO2). The effects of reaction time, temperature, gold loadings and base were investigated. High yields were achieved under moderate conditions (60 °C) in the presence of KOCH3. For instance, the yield of N-benzylidenebenzylamine produced from benzyl alcohol and benzylamine on 3 wt% Au/ZrO2 is 87 %. The synthesis of imine involves two reaction steps: selective oxidation of benzyl alcohol to benzaldehyde and the coupling reaction of amines with benzaldehyde. In the first step, the base promotes the selective oxidation. The reactions of benzyl alcohol with three different amines, aniline, n-butylamine and benzylamine, were conducted to produce corresponding imines. The results show that the amine with stronger nucleophilicity has better ability to react with benzaldehyde in the second step, resulting in higher yield of the corresponding imine. We proposed a tentative mechanism for the synthesis process.

Hydrolysis of genotoxic methyl-substituted oxiranes : Experimental kinetic and semiempirical studies

The kinetics of acid-catalyzed hydrolysis of seven methylated aliphatic epoxides - R1R2C(O)CR3R4 (A: R1=R2=R3=R4=H; B: R1=R2=R3=H, R4=Me; C: R1=R2=H, R3=R4=Me; D: R1=R3=H, R2=R4=Me(trans); E: R1=R3=H, R2=R4=Me(cis); F: R1=R3=R4=Me, R2=H; G: R1=R2=R3=R4=Me) - has been studied at 36 ± 1.5°C. Compounds with two methyl groups at the same carbon atom of the oxirane ring exhibit highest rate constants (k(eff) in reciprocal molar concentration per second: 11.0 ± 1.3 for C, 10.7 ± 2.1 for F, and 8.7 ± 0.7 for G as opposed to 0.124 ± 0.003 for B, 0.305 ± 0.003 for D, and 0.635 ± 0.036 for E). Ethylene oxide (A) displays the lowest rate of hydrolysis (0.027 M-1 s-1). The results are consistent with literature data available for compounds A, B, and C. To model the reactivities we have employed quantum chemical calculations (MNDO, AM1, PM3, and MINDO/3) of the main reaction species. There is a correlation of the logarithm k(eff) with the total energy of epoxide ring opening. The best correlation coefficients (r) were obtained using the AM1 and MNDO methods (0.966 and 0.957, respectively). However, unlike MNDO, AM1 predicts approximately zero energy barriers for the oxirane ring opening of compounds B, C, E and G, which is not consistent with published kinetic data. Thus, the MNDO method provides a preferential means of modeling the acidic hydrolysis of the series of methylated oxiranes. The general ranking of mutagenicity in vitro, A > B > C, is in line with the concept that this sequence also gradually leaves the expoxide reactivity optimal for genotoxicity toward reactivities leading to higher biological detoxifications.

Development of a log-quadratic model to describe microbial inactivation, illustrated by thermal inactivation of Clostridium botulinum

In the commercial food industry, demonstration of microbiological safety and thermal process equivalence often involves a mathematical framework that assumes log-linear inactivation kinetics and invokes concepts of decimal reduction time (DT), z values, and accumulated lethality. However, many microbes, particularly spores, exhibit inactivation kinetics that are not log linear. This has led to alternative modeling approaches, such as the biphasic and Weibull models, that relax strong log-linear assumptions. Using a statistical framework, we developed a novel log-quadratic model, which approximates the biphasic and Weibull models and provides additional physiological interpretability. As a statistical linear model, the log-quadratic model is relatively simple to fit and straightforwardly provides confidence intervals for its fitted values. It allows a DT-like value to be derived, even from data that exhibit obvious "tailing." We also showed how existing models of non-log-linear microbial inactivation, such as the Weibull model, can fit into a statistical linear model framework that dramatically simplifies their solution. We applied the log-quadratic model to thermal inactivation data for the spore-forming bacterium Clostridium botulinum and evaluated its merits compared with those of popular previously described approaches. The log-quadratic model was used as the basis of a secondary model that can capture the dependence of microbial inactivation kinetics on temperature. This model, in turn, was linked to models of spore inactivation of Sapru et al. and Rodriguez et al. that posit different physiological states for spores within a population. We believe that the log-quadratic model provides a useful framework in which to test vitalistic and mechanistic hypotheses of inactivation by thermal and other processes. Copyright © 2009, American Society for Microbiology. All Rights Reserved.