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.

Catalytic activity evaluation of industrial Pd/C catalyst via gray-box dynamic modeling and simulation of hydropurification reactor

In this paper, dynamic modeling and simulation of the hydropurification reactor in a purified terephthalic acid production plant has been investigated by gray-box technique to evaluate the catalytic activity of palladium supported on carbon (0.5 wt.% Pd/C) catalyst. The reaction kinetics and catalyst deactivation trend have been modeled by employing artificial neural network (ANN). The network output has been incorporated with the reactor first principle model (FPM). The simulation results reveal that the gray-box model (FPM and ANN) is about 32 percent more accurate than FPM. The model demonstrates that the catalyst is deactivated after eleven months. Moreover, the catalyst lifetime decreases about two and half months in case of 7 percent increase of reactor feed flowrate. It is predicted that 10 percent enhancement of hydrogen flowrate promotes catalyst lifetime at the amount of one month. Additionally, the enhancement of 4-carboxybenzaldehyde concentration in the reactor feed improves CO and benzoic acid synthesis. CO is a poison to the catalyst, and benzoic acid might affect the product quality. The model can be applied into actual working plants to analyze the Pd/C catalyst efficient functioning and the catalytic reactor performance.

Limitations of a laboratory scale model in predicting optimal pilot scale conditions for dilute acid pretreatment of sugarcane bagasse

Pilot and industrial scale dilute acid pretreatment data can be difficult to obtain due to the significant infrastructure investment required. Consequently, models of dilute acid pretreatment by necessity use laboratory scale data to determine kinetic parameters and make predictions about optimal pretreatment conditions at larger scales. In order for these recommendations to be meaningful, the ability of laboratory scale models to predict pilot and industrial scale yields must be investigated. A mathematical model of the dilute acid pretreatment of sugarcane bagasse has previously been developed by the authors. This model was able to successfully reproduce the experimental yields of xylose and short chain xylooligomers obtained at the laboratory scale. In this paper, the ability of the model to reproduce pilot scale yield and composition data is examined. It was found that in general the model over predicted the pilot scale reactor yields by a significant margin. Models that appear very promising at the laboratory scale may have limitations when predicting yields on a pilot or industrial scale. It is difficult to comment whether there are any consistent trends in optimal operating conditions between reactor scale and laboratory scale hydrolysis due to the limited reactor datasets available. Further investigation is needed to determine whether the model has some efficacy when the kinetic parameters are re-evaluated by parameter fitting to reactor scale data, however, this requires the compilation of larger datasets. Alternatively, laboratory scale mathematical models may have enhanced utility for predicting larger scale reactor performance if bulk mass transport and fluid flow considerations are incorporated into the fibre scale equations. This work reinforces the need for appropriate attention to be paid to pilot scale experimental development when moving from laboratory to pilot and industrial scales for new technologies.

Solution chemistry impacts on the seawater neutralisation process: Benefits of nanofiltered seawater and reverse osmosis brine

It is well known that the neutralisation of Bayer liquor with seawater causes the precipitation of stable alkaline products and a reduction in pH and dissolved metal concentrations in the effluent. However, there is limited information available on solution chemistry effects on the stability and reaction kinetics of these precipitates. This investigation shows the influence of reactive species (magnesium and calcium) in seawater on precipitate stabilities and volumetric efficiencies during the neutralisation of bauxite refinery residues. Correlations between synthetic seawater solutions and real samples of seawater (filtered seawater, nanofiltered seawater and reverse osmosis brine) have been made. These investigations have been used to confirm that alternative seawater sources can be used to increase the productivity potential of the neutralisation process with minimal implications on the composition and stability of precipitates formed. The volume efficiency of the neutralisation process using synthetic analogues has been shown to be almost directly proportional with the concentration of magnesium. This was further confirmed in the nanofiltered seawater and reverse osmosis brine that showed increases in the efficiency of neutralisation by factors of 3 and 2 compared to seawater, which corresponds with relatively the same increase in the concentration of magnesium in these alternative seawater sources. An assessment of the chemical stability of the precipitates, volumetric efficiency, and discharge water quality have been determined using numerous techniques that include pH, conductivity, inductively coupled plasma optical emission spectroscopy, infrared spectroscopy, thermogravimetric analysis coupled to mass spectrometry and X-ray diffraction. Correlations between synthetic solution compositions and alternative seawater sources have been used to determine if alternative seawater sources are potential substitutes for seawater based on improvements in productivity, implementation costs, savings to operations and environmental benefits.

A kinetic study of the methanation of CO2 over Ni-Al2O3 catalyst

CO2 methanation has been studied over Ni-Al2O3 and the intrinsic kinetic data obtained are related by a powerlaw type rate expression, which established a good agreement between calculated and experimental values of conversion.