Hierarchical control and optimisation of the Tennessee Eastman process

Based on integrated system optimisation and parameter estimation a method is described for on-line steady state optimisation which compensates for model-plant mismatch and solves a non-linear optimisation problem by iterating on a linear - quadratic representation. The method requires real process derivatives which are estimated using a dynamic identification technique. The utility of the method is demonstrated using a simulation of the Tennessee Eastman benchmark chemical process.

An overview of interactive visual data mining techniques for knowledge discovery

n the past decade, the analysis of data has faced the challenge of dealing with very large and complex datasets and the real-time generation of data. Technologies to store and access these complex and large datasets are in place. However, robust and scalable analysis technologies are needed to extract meaningful information from these datasets. The research field of Information Visualization and Visual Data Analytics addresses this need. Information visualization and data mining are often used complementary to each other. Their common goal is the extraction of meaningful information from complex and possibly large data. However, though data mining focuses on the usage of silicon hardware, visualization techniques also aim to access the powerful image-processing capabilities of the human brain. This article highlights the research on data visualization and visual analytics techniques. Furthermore, we highlight existing visual analytics techniques, systems, and applications including a perspective on the field from the chemical process industry.

Chemical characterisation and determination of sensory attributes of hydrolysates produced by enzymatic hydrolysis of whey proteins following a novel integrative process

The overall aim of this work was to characterize the major angiotensin converting enzyme (ACE) inhibitory peptides produced by enzymatic hydrolysis of whey proteins, through the application of a novel integrative process. This process consisted of the combination of adsorption and microfiltration within a stirred cell unit for the selective immobilization of β-lactoglobulin and casein derived peptides (CDP) from whey. The adsorbed proteins were hydrolyzed in-situ which resulted in the separation of peptide products from the substrate and fractionation of peptides. Two different hydrolysates were produced: (i) from CDP (IC50 =287μg/mL) and (ii) from β-lactoglobulin (IC50=128μg/mL). IC50 is the concentration of inhibitor needed to inhibit ACE by half. The well known antihypertensive peptide IPP and several novel peptides that have structural similarities with reported ACE inhibitory peptides were identified and characterized in both hydrolysates. Furthermore, the hydrolysates were assessed for bitterness. No significant difference was found between the control (milk with no hydrolysate) and hydrolysate samples at different concentrations (at, below and above the IC50).

Relative proportions of polycyclic aromatic hydrocarbons differ between accumulation bioassays and chemical methods to predict bioavailability

Chemical methods to predict the bioavailable fraction of organic contaminants are usually validated in the literature by comparison with established bioassays. A soil spiked with polycyclic aromatic hydrocarbons (PAHs) was aged over six months and subjected to butanol, cyclodextrin and tenax extractions as well as an exhaustive extraction to determine total PAH concentrations at several time points. Earthworm (Eisenia fetida) and rye grass root (Lolium multiflorum) accumulation bioassays were conducted in parallel. Butanol extractions gave the best relationship with earthworm accumulation (r2 ≤ 0.54, p ≤ 0.01); cyclodextrin, butanol and acetone–hexane extractions all gave good predictions of accumulation in rye grass roots (r2 ≤ 0.86, p ≤ 0.01). However, the profile of the PAHs extracted by the different chemical methods was significantly different (p < 0.01) to that accumulated in the organisms. Biota accumulated a higher proportion of the heavier 4-ringed PAHs. It is concluded that bioaccumulation is a complex process that cannot be predicted by measuring the bioavailable fraction alone. The ability of chemical methods to predict PAH accumulation in Eisenia fetida and Lolium multiflorum was hindered by the varied metabolic fate of the different PAHs within the organisms.

Human exposure modelling for chemical risk assessment: a review of current approaches and research and policy implications

A wide variety of exposure models are currently employed for health risk assessments. Individual models have been developed to meet the chemical exposure assessment needs of Government, industry and academia. These existing exposure models can be broadly categorised according to the following types of exposure source: environmental, dietary, consumer product, occupational, and aggregate and cumulative. Aggregate exposure models consider multiple exposure pathways, while cumulative models consider multiple chemicals. In this paper each of these basic types of exposure model are briefly described, along with any inherent strengths or weaknesses, with the UK as a case study. Examples are given of specific exposure models that are currently used, or that have the potential for future use, and key differences in modelling approaches adopted are discussed. The use of exposure models is currently fragmentary in nature. Specific organisations with exposure assessment responsibilities tend to use a limited range of models. The modelling techniques adopted in current exposure models have evolved along distinct lines for the various types of source. In fact different organisations may be using different models for very similar exposure assessment situations. This lack of consistency between exposure modelling practices can make understanding the exposure assessment process more complex, can lead to inconsistency between organisations in how critical modelling issues are addressed (e.g. variability and uncertainty), and has the potential to communicate mixed messages to the general public. Further work should be conducted to integrate the various approaches and models, where possible and regulatory remits allow, to get a coherent and consistent exposure modelling process. We recommend the development of an overall framework for exposure and risk assessment with common approaches and methodology, a screening tool for exposure assessment, collection of better input data, probabilistic modelling, validation of model input and output and a closer working relationship between scientists and policy makers and staff from different Government departments. A much increased effort is required is required in the UK to address these issues. The result will be a more robust, transparent, valid and more comparable exposure and risk assessment process. (C) 2006 Elsevier Ltd. All rights reserved.

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with oxygen

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with O-2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second order rate constants at 10 Torr were fitted to the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.08 +/- 0.04) + (1.57 +/- 0.32 kJ mol(-1))/RT ln10 The decrease in rate constant values with increasing temperature, although systematic is very small. The rate constants showed slight increases in value with pressure at each temperature, but this was scarcely beyond experimental uncertainty. From estimates of Lennard-Jones collision rates, this reaction is occurring at ca. 1 in 20 collisions, almost independent of pressure and temperature. Ab initio calculations at the G3 level backed further by multi-configurational (MC) SCF calculations, augmented by second order perturbation theory (MRMP2), support a mechanism in which the initial adduct, H2SiOO, formed in the triplet state (T), undergoes intersystem crossing to the more stable singlet state (S) prior to further low energy isomerisation processes leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are H2O + SiO. The decomposition of the intermediate cyclo-siladioxirane, via O-O bond fission, plays an important role in the overall process. The bottleneck for the overall process appears to be the T -> S process in H2SiOO. This process has a small spin orbit coupling matrix element, consistent with an estimate of its rate constant of 1 x 10(9) s(-1) obtained with the aid of RRKM theory. This interpretation preserves the idea that, as in its reactions in general, SiH2 initially reacts at the encounter rate with O-2. The low values for the secondary reaction barriers on the potential energy surface account for the lack of an observed pressure dependence. Some comparisons are drawn with the reactions of CH2 + O-2 and SiCl2 + O-2.

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.

Investigation of the prototype silylene reaction, SiH2+H2O(and D2O): time-resolved gas-phase kinetic studies, isotope effects, RRKM calculations, and quantum chemical calculations of the reaction energy surface

Time-resolved kinetic studies of the reaction of silylene, SiH2, with H2O and with D2O have been carried out in the gas phase at 296 and at 339 K, using laser flash photolysis to generate and monitor SiH2. The reaction was studied over the pressure range 10-200 Torr with SF6 as bath gas. The second-order rate constants obtained were pressure dependent, indicating that the reaction is a third-body assisted association process. Rate constants at 339 K were about half those at 296 K. Isotope effects, k(H)/k(D), were small averaging 1.076 0.080, suggesting no involvement of H- (or D-) atom transfer in the rate determining step. RRKM modeling was undertaken based on a transition state appropriate to formation of the expected zwitterionic donoracceptor complex, H2Si...OH2. Because the reaction is close to the low pressure (third order) region, it is difficult to be definitive about the activated complex structure. Various structures were tried, both with and without the incorporation of rotational modes, leading to values for the high-pressure limiting (i.e., true secondorder) rate constant in the range 9.5 x 10(-11) to 5 x 10(-10) cm(3) molecule' s(-1). The RRKM modeling and mechanistic interpretation is supported by ab initio quantum calculations carried out at the G2 and G3 levels. The results are compared and contrasted with the previous studies.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 1. HCl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with HCL The reaction was studied in the gas phase at 10 Torr total pressure in SF6 bath gas, at five temperatures in the range of 296-611 K. The second-order rate constants fitted the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.51 +/- 0.06) + (1.92 +/- 0.47 kJ mol(-1))/RTIn10 Experiments at other pressures showed that these rate constants were unaffected by pressure in the range of 10-100 Torr, but showed small decreases in value of no more than 20% ( +/- 10%) at I Toff, at both the highest and lowest temperatures. The data are consistent with formation of an initial weakly bound donor-acceptor complex, which reacts by two parallel pathways. The first is by chlorine-to-silicon H-shift to make vibrationally excited chlorosilane, SiH3Cl*, which yields HSiCl by H-2 elimination from silicon. In the second pathway, the complex proceeds via H-2 elimination (4-center process) to make chlorosilylene, HSiCl, directly. This interpretation is supported by ab initio quantum calculations carried out at the G3 level which reveal the direct H-2 elimination route for the first time. RRKM modeling predicts the approximate magnitude of the pressure effect but is unable to determine the proportions of each pathway. The experimental data agree with the only previous measurements at room temperature. Comparisons with other reactions of SiH2 are also drawn.

Time-resolved gas-phase kinetic and quantum chemical studies of reactions of silylene with chlorine-containing species. 2. CH3Cl

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of both silacyclopent-3-ene and phenylsilane, have been carried out to obtain second-order rate constants for its reaction with CH3Cl. The reaction was studied in the gas phase at six temperatures in the range 294-606 K. The second-order rate constants gave a curved Arrhenius plot with a minimum value at T approximate to 370 K. The reaction showed no pressure dependence in the presence of up to 100 Torr SF6. The rate constants, however, showed a weak dependence on laser pulse energy. This suggests an interpretation requiring more than one contributing reaction pathway to SiH2 removal. Apart from a direct reaction of SiH2 with CH3Cl, reaction of SiH2 with CH3 (formed by photodissociation of CH3Cl) seems probable, with contributions of up to 30% to the rates. Ab initio calculations (G3 level) show that the initial step of reaction of SiH2 with CH3Cl is formation of a zwitterionic complex (ylid), but a high-energy barrier rules out the subsequent insertion step. On the other hand, the Cl-abstraction reaction leading to CH3 + ClSiH2 has a low barrier, and therefore, this seems the most likely candidate for the main reaction pathway of SiH2 with CH3Cl. RRKM calculations on the abstraction pathway show that this process alone cannot account for the observed temperature dependence of the rate constants. The data are discussed in light of studies of other silylene reactions with haloalkanes.

Chemical, physical and electrical properties of aged dodecylbenzene: thermal ageing of mixed isomers in air

A commercial dodecylbenzene (DDB) cable oil was aged at temperatures between 90 and 135 degrees C in air and was analyzed using various analytical techniques including optical and infra-red spectroscopy and dielectric analysis. On ageing, the oil darkened, significant oxidation features were found by infra-red spectroscopy and the acid number, water content and dielectric loss all increased. Ageing in the presence of paper or aluminum did not affect the ageing process, whereas ageing was significantly modified by the presence of copper. An absorption at 680 nm ("red absorbers") was detected by ultra-violet/visible spectroscopy followed by the production of an opaque precipitate. A reaction between copper and the acid generated on ageing is thought to produce copper carboxylates, and X-ray fluorescence confirmed that copper was indeed present in both the aged oil and the precipitate. Significantly, once red absorbers were detected, the dielectric loss increased to catastrophically high values and, therefore, the appearance of these compounds may serve as a useful diagnostic indicator. The development of acidity on ageing appears to be key in initiating the destructive copper conversion reaction and hence the control of oil acidity may be key to prolonging the life of DDB cable oils.

Chemical, physical and electrical properties of aged dodecylbenzene 2: thermal ageing of single isomers in air

The linear isomer of dodecylbenzene (DDB), 1-phenyldodecane, was aged at temperatures of 105 and 135 degrees C in air and the resultant products were analyzed using a range of analytical techniques. On ageing, the 1-phenyldodecane darkened, the acid number, dielectric loss and water content increased and significant oxidation peaks were detected in the infrared spectrum. When aged in the presence of copper, a characteristic peak at 680 nm was also detected by UV/visible spectroscopy but, compared with previous studies of a cable-grade DDB, the strength of this peak was much increased and no appreciable precipitate formation occurred. At the same time, very high values of dielectric loss were recorded. On ageing in the absence of copper, an unusually strong infrared carbonyl band was seen, which correlates well with the detection of dodecanophenone by gas chromatography / mass spectrometry and nuclear magnetic resonance spectroscopy. It was therefore concluded that the ageing process proceeds via the initial production of aromatic ketones, which may then be further oxidized to carboxylic acids. In the presence of copper, these oxidation products are present in lower quantities, most of these oxidation products being combined with the copper present in the oil to give copper carboxylates. The behavior is described in terms of a complex autoxidation mechanism, in which copper acts as both an oxidizing and a reducing agent, depending on its oxidation state and, in particular, promotes elimination via the oxidation of intermediate alkyl radical species to carbocations.

Some new strategies for the chemical separation of actinides and lanthanides

A future goal in nuclear fuel reprocessing is the conversion or transmutation of the long-lived radioisotopes of minor actinides, such as americium, into short-lived isotopes by irradiation with neutrons. In order to achieve this transmutation, it is necessary to separate the minor actinides(III), [An(Ill)], from the lanthanides(III), [Ln(Ill)], by solvent extraction (partitioning), because the lanthanides absorb neutrons too effectively and hence limit neutron capture by the transmutable actinides. Partitioning using ligands containing only carbon, hydrogen, nitrogen and oxygen atoms is desirable because they are completely incinerable and thus the final volume of waste is minimised [1]. Nitric acid media will be used in the extraction experiments because it is envisaged that the An(III)/Ln(III) separation process could take place after the PUREX process. There is no doubt that the correct design of a molecule that is capable of acting as a ligand or extraction reagent is required for the effective separation of metal ions such as actinides(III) from lanthanides. Recent attention has been directed towards heterocyclic ligands with for the preferential separation of the minor actinides. Although such molecules have a rich chemistry, this is only now becoming sufficiently well understood in relation to the partitioning process [2]. The molecules shown in Figures I and 2 will be the principal focus of this study. Although the examples chosen here are used rather specific, the guidelines can be extended to other areas such as the separation of precious metals [3].

Ultrasonic chemical oxidative degradations of 1,3-dialkylimidazolium ionic liquids and their mechanistic elucidations

A highly efficient process for oxidative degradation of 1,3-dialkylimidazolium ionic liquids in hydrogen peroxide/acetic acid aqueous medium assisted by ultrasonic chemical irradiation is, for the first time, described. It is shown that more than 93% of the 1,3-dialkylimidazolium cation with the corresponding Cl-, Br-, BF4- and PF6- counter-anions at a concentration of 2.5 mM can be degraded at 50 degrees C within 12 h while at 72 h the conversions approach 99%. A tentative mechanism for the degradation of these ILs is for the first time proposed through a detailed kinetic analysis of several characteristic transients and/or immediate products, which are identified during the ILs degradation using GC-MS. The results clearly indicate that three hydrogen atoms in the imidazolium ring are the first sites preferably oxidized, followed by cleavage of the alkyl groups attached to the N atoms from the ring. The nature of the alkyl chain length on the imidazolium ring and the type of counter anion do not seem to affect the degradation process. Further, selective fragmentations of C-N bonds of the imidazolium or derived ring lead to ring opening, forming degraded intermediates. It is also shown that acetoxyacetic acid and biurea are the final kinetically stable degraded products from the ILs under the degradation conditions.

Seasonal changes in olive fruit characteristics and oil accumulation during ripening process

Olive fruits of three of the most important Spanish and Italian cultivars, 'Picual', `Hojiblanca' and 'Frantoio', were harvested at bi-weekly periods during three crop seasons to study their development and ripening process. Fresh and dry weights and ripening index were determined for fruit, while dry matter, oil and moisture contents were determined in both fruit and pulp (flesh). Fruit growth rate and oil accumulation were calculated. Each olive cultivar showed a different ripening pattern, 'Hojiblanca' being the last one to maturate. Fruit weight increased, decreasing its growth rate from the middle of November. Dry matter and moisture contents decreased during ripening in pulp and fruit, 'Hojiblanca' showing the highest values for both. Oil content, when expressed on a fresh weight basis, increased in all cultivars, although for the last time period showed variations due to climatic conditions. During ripening, oil content on a dry weight basis increased in fruit, but oil biosynthesis in flesh ceased from November. Olive fruits presented lower oil and higher dry matter contents in the year of lowest rainfall. Therefore fruit harvesting should be carried out from the middle of November in order to obtain the highest oil yield and avoid natural fruit drop. (C) 2004 Society of Chemical Industry.