Adaptive multiscale principal components analysis for online monitoring of wastewater treatment

Fault detection and isolation (FDI) are important steps in the monitoring and supervision of industrial processes. Biological wastewater treatment (WWT) plants are difficult to model, and hence to monitor, because of the complexity of the biological reactions and because plant influent and disturbances are highly variable and/or unmeasured. Multivariate statistical models have been developed for a wide variety of situations over the past few decades, proving successful in many applications. In this paper we develop a new monitoring algorithm based on Principal Components Analysis (PCA). It can be seen equivalently as making Multiscale PCA (MSPCA) adaptive, or as a multiscale decomposition of adaptive PCA. Adaptive Multiscale PCA (AdMSPCA) exploits the changing multivariate relationships between variables at different time-scales. Adaptation of scale PCA models over time permits them to follow the evolution of the process, inputs or disturbances. Performance of AdMSPCA and adaptive PCA on a real WWT data set is compared and contrasted. The most significant difference observed was the ability of AdMSPCA to adapt to a much wider range of changes. This was mainly due to the flexibility afforded by allowing each scale model to adapt whenever it did not signal an abnormal event at that scale. Relative detection speeds were examined only summarily, but seemed to depend on the characteristics of the faults/disturbances. The results of the algorithms were similar for sudden changes, but AdMSPCA appeared more sensitive to slower changes.

Ab initio modelling of basal plane oxidation of graphenes and implications for modelling char combustion

Ab initio calculations have been performed to determine the energetics of oxygen atoms adsorbed onto graphene planes and the possible reaction path extracting carbon atorns in the form of carbon monoxide. Front the energetics it is confirmed that this reaction path will not significantly contribute to the gasification of well ordered carbonaceous chars. Modelling results which explore this limit Lire presented. (C) 2002 Elsevier Science Ltd, All rights reserved.

Work domain analysis and sensors II: Pasteurizer II case study

In this paper we use sensor-annotated abstraction hierarchies (Reising & Sanderson, 1996, 2002a,b) to show that unless appropriately instrumented, configural displays designed according to the principles of ecological interface design (EID) might be vulnerable to misinterpretation when sensors become unreliable or are unavailable. Building on foundations established in Reising and Sanderson (2002a) we use a pasteurization process control example to show how sensor-annotated AHs help the analyst determine the impact of different instrumentation engineering policies on a configural display that is part of an ecological interface. Our analyses suggest that configural displays showing higher-order properties of a system are especially vulnerable under some conservative instrumentation configurations. However, sensor-annotated AHs can be used to indicate where corrective instrumentation might be placed. We argue that if EID is to be effectively employed in the design of displays for complex systems, then the information needs of the human operator need to be considered while instrumentation requirements are being formulated. Rasmussen's abstraction hierarchy-and particularly its extension to the analysis of information captured by sensors and derived from sensors-may therefore be a useful adjunct to up-stream instrumentation design. (C) 2002 Elsevier Science Ltd. All rights reserved.

Scale-up of mixer granulators for effective liquid distribution

There is considerable anecdotal evidence from industry that poor wetting and liquid distribution can lead to broad granule size distributions in mixer granulators. Current scale-up scenarios lead to poor liquid distribution and a wider product size distribution. There are two issues to consider when scaling up: the size and nature of the spray zone and the powder flow patterns as a function of granulator scale. Short, nucleation-only experiments in a 25L PMA Fielder mixer using lactose powder with water and HPC solutions demonstrated the existence of different nucleation regimes depending on the spray flux Psi(a)-from drop-controlled nucleation to caking. In the drop-controlled regime at low Psi(a) values. each drop forms a single nucleus and the nuclei distribution is controlled by the spray droplet size distribution. As Psi(a) increases, the distribution broadens rapidly as the droplets overlap and coalesce in the spray zone. The results are in excellent agreement with previous experiments and confirm that for drop-controlled nucleation. Psi(a) should be less than 0.1. Granulator flow studies showed that there are two powder flow regimes-bumping and roping. The powder flow goes through a transition from bumping to roping as impeller speed is increased. The roping regime gives good bed turn over and stable flow patterns. This regime is recommended for good liquid distribution and nucleation. Powder surface velocities as a function of impeller speed were measured using high-speed video equipment and MetaMorph image analysis software, Powder surface velocities were 0.2 to 1 ms(-1)-an order of magnitude lower than the impeller tip speed. Assuming geometrically similar granulators, impeller speed should be set to maintain constant Froude number during scale-up rather than constant tip speed to ensure operation in the roping regime. (C) 2002 Published by Elsevier Science B.V.

Growth and compaction behaviour of copper concentrate granules in a rotating drum

In order to understand the growth and compaction behaviour of chalcopyrite (copper concentrate), batch granulation tests were carried out using a rotating drum. The granule growth exhibited induction-type behaviour, as defined by Iveson and Litster [AIChE J. 44 (1998) 15 10]. There were two consecutive stages during granulation: the induction stage, during which the granules are gradually being compacted and little or no growth occurs, and the rapid growth stage, which starts when the granules have become surface wet and are rapidly growing. In agreement with earlier findings. an increased amount of binder liquid shortened the induction time. The compaction behaviour was also investigated. A displaced volume method was adopted to determine the porosity of the granules. It was shown that this technique had a limitation as it was unable to detect the reduction of the volumes of the granule pores after the granules had become surface wet. Due to this, some of the measurements were not suited for fitting a three-parameter empirical model. Attempts were made to determine whether the rapid growth stage started with the pore saturation exceeding a certain critical value, but due to the scatter in the porosity measurements and the fact that some of the measurements could not be used, it was not possible to determine a critical pore saturation, However, the porosity measurements clearly demonstrated that the porosity of the granules decreased during the induction stage of an experiment and that when rapid growth occurred, the granules had a pore saturation was around 0.85. This value was slightly lower than unity, which is most likely due to trapped air bubbles. (C) 2002 Published by Elsevier Science B.V.

Polymer product engineering utilising oscillatory baffled reactors

The Oscillatory baffled reactor (OBR) can be used to produce particles with controlled size and morphology, in batch or continuous flow. This is due to the effect of the superimposed oscillations that radially mixes fluid but still allows plug-flow (or close to plug flow) behaviour in a continuous system. This mixing, combined with a close to a constant level of turbulence intensity in the reactor, leads to tight droplet and subsequent product particle size distributions. By applying population balance equations together with experimental droplet size distributions, breakage rates of droplets can be determined and this is a useful tool for understanding the product engineering in OBRs. (C) 2002 Elsevier Science B.V All rights reserved.

On the effect of tracer density on axial dispersion in a batch oscillatory baffled column

In this paper, we report our modelling evaluation on the effect of tracer density on axial dispersion in a batch oscillatory baffled column (OBC). Tracer solution of potassium nitrite, its specific density ranged from 1.0 to 1.5, was used in the study, and was injected to the vertical column from either the top or bottom. Local concentration profiles are measured using conductivity probes at two locations along the height of the column. Using the experimental measured concentration profiles together with both 'Tank-in-Series' and 'Plug Flow with Axial Dispersion' models, axial dispersion coefficients were determined and used to describe the effect of specific tracer density on mixing in the OBC. The results showed that the axial dispersion coefficients evaluated by the two models are very similar in both magnitudes and trends, and the range of variations in such coefficients is generally larger for the bottom injection than for the top one. Empirical correlations linking the mechanical energy for mixing, the specific density of tracer and axial dispersion coefficient were established. Using these correlations, we identified the enhancements of up to 269% on axial dispersion for various specific tracer densities. (C) 2002 Elsevier Science B.V. All rights reserved.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Structural ordering of coal char during heat treatment and its impact on reactivity

The effect of heat treatment on the structure of an Australian semi-anthracite char was studied in detail in the 850-1150degreesC temperature range using XRD, HRTEM, and electrical resistivity techniques. It was found that the carbon crystallite size in the char does not change significantly during heat treatment in the temperature range studied, for both the raw coal and its ash-free derivative obtained by acid treatment. However, the fraction of the organized carbon in the raw coal chars, determined by XRD, increased with increase of heat treatment time and temperature, while that for the ash-free coal chars remained almost unchanged. This suggests the occurrence of catalytic ordering during heat treatment, supported by the observation that the electrical resistivity of the raw coal chars decreased with heat treatment, while that of the ash-free coal chars did not vary significantly. Further confirmatory evidence was provided by high resolution transmission electron micrographs depicting well-organized carbon layers surrounding iron particles. It is also found that the fraction of organized carbon does not reach unity, but attains an apparent equilibrium value that increases with increase in temperature, providing an apparent heat of ordering of 71.7 kJ mol(-1) in the temperature range studied. Good temperature-independent correlation was found between the electrical resistivity and the organized carbon fraction, indicating that electrical resistivity is indeed structure sensitive. Good correlation was also found between the electrical resistivity and the reactivity of coal char. All these results strongly suggest that the thermal deactivation is the result of a crystallite-perfecting process, which is effectively catalyzed by the inorganic matter in the coal char. Based on kinetic interpretation of the data it is concluded that the process is diffusion controlled, most likely involving transport of iron in the inter-crystallite nanospaces in the temperature range studied. The activation energy of this transport process is found to be very low, at about 11.8 kJ mol(-1), which is corroborated by model-free correlation of the temporal variation of organized carbon fraction as well as electrical resistivity data using the superposition method, and is suggestive of surface transport of iron. (C) 2002 Elsevier Science Ltd. All rights reserved.

Kinetics of adsorption on activated carbon: application of heterogeneous vacancy solution theory

The kinetics of single component adsorption on activated carbon is investigated here using a heterogeneous vacancy solution theory (VST) of adsorption. The adsorption isotherm is developed to account for the adsorbate non-ideality due to the size difference between the adsorbate molecule and the vacant site, while incorporating adsorbent heterogeneity through a pore-width-related potential energy. The transport process in the bidisperse carbon considers coupled mass transfer in both macropore and micropore phases simultaneously. Adsorbate diffusion in the micropore network is modeled through effective medium theory, thus considering pore network connectivity in the adsorbent, with the activation energy for adsorbate diffusion related to the adsorption energy, represented by the Steele 10-4-3 potential for carbons. Experimental data of five hydrocarbons, CO2 and SO2 on Ajax carbon at multiple temperatures, as well as three hydrocarbons on Norit carbon at three temperatures are first fitted by the heterogeneous VST model to obtain the isotherm parameters, followed by application of the kinetic model to uptake data on carbon particles of different sizes and geometry at various temperatures. For the hydrocarbons studied, the model can successfully correlate the experimental data for both adsorption equilibrium and kinetics. However, there is some deviation in the fit of the desorption kinetics for polar compounds such as CO2 and SO2, due to the inadequacy of the L-J potential model in this case. The significance of viscous transport in the micropores is also considered here and found to be negligible, consistent with recent molecular simulation studies. (C) 2002 Elsevier Science Ltd. All rights reserved.

A comprehensive insight into floc characteristics and their impact on compressibility and settleability of activated sludge

This paper presents a comprehensive study of sludge floc characteristics and their impact on compressibility and settleability of activated sludge in full scale wastewater treatment processes. The sludge flocs were characterised by morphological (floc size distribution, fractal dimension, filament index), physical (flocculating ability, viscosity, hydrophobicity and surface charge) and chemical (polymeric constituents and metal content) parameters. Compressibility and settleability were defined in terms of the sludge volume index (SVI) and zone settling velocity (ZSV). The floc morphological and physical properties had important influence on the sludge compressibility and settleability. Sludges containing large flocs and high quantities of filaments, corresponding to lower values of fractal dimension (D-f), demonstrated poor compressibility and settleability. Sludge flocs with high flocculating ability had lower SVI and higher ZSV, whereas high values of hydrophobicity, negative surface charge and viscosity of the sludge flocs correlated to high SVI and low ZSV. The quantity of the polymeric compounds protein. humic substances and carbohydrate in the sludge and the extracted extracellular polymeric substances (EPS) had significant positive correlations with SVI. The ZSV was quantitatively independent of the polymeric constituents. High concentrations of the extracted EPS were related to poor compressibility and settleability. The cationic ions Ca, Mg, Al and Fe in the sludge improved significantly the sludge compressibility and settleability. (C) 2003 Elsevier Science B.V. All rights reserved.

A new wavelet-based adaptive method for solving population balance equations

A new wavelet-based adaptive framework for solving population balance equations (PBEs) is proposed in this work. The technique is general, powerful and efficient without the need for prior assumptions about the characteristics of the processes. Because there are steeply varying number densities across a size range, a new strategy is developed to select the optimal order of resolution and the collocation points based on an interpolating wavelet transform (IWT). The proposed technique has been tested for size-independent agglomeration, agglomeration with a linear summation kernel and agglomeration with a nonlinear kernel. In all cases, the predicted and analytical particle size distributions (PSDs) are in excellent agreement. Further work on the solution of the general population balance equations with nucleation, growth and agglomeration and the solution of steady-state population balance equations will be presented in this framework. (C) 2002 Elsevier Science B.V. All rights reserved.

Prediction of multilayer adsorption and capillary condensation phenomena in cylindrical mesopores

MCM-41 periodic mesoporous silicates with a high degree of structural ordering are synthesized and used as model adsorbents to study the isotherm prediction of nitrogen adsorption. The nitrogen adsorption isotherm at 77 K for a macroporous silica is measured and used in high-resolution alpha(s)-plot comparative analysis to determine the external surface area, total surface area and primary mesopore volume of the MCM-41 materials. Adsorption equilibrium data of nitrogen on the different pore size MCM-41 samples (pore diameters from 2.40 to 4.92 nm) are also obtained. Based on the Broekhoff and de Boer' thermodynamic analysis, the nitrogen adsorption isotherms for the different pore size MCM-41 samples are interpreted using a novel strategy, in which the parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting only the multilayer region prior to capillary condensation for C-16 MCM-41. Subsequently the entire isotherm, including the phase transition, is predicted for all the different pore size MCM-41 samples without any fitting. The results show that the prediction of multilayer adsorption and total adsorbed amount are in good agreement with the experimental isotherms. The predictions of the relative pressure corresponding to capillary equilibrium (coexistence) transition agree remarkably with experimental data on the adsorption branch even for hysteretic isotherms, confirming that this is the branch appropriate for pore size distribution analysis. The impact of pore radius on the adsorption film thickness and capillary coexistence pressure is also investigated, and found to agree with the experimental data. (C) 2003 Elsevier Inc. All rights reserved.

Cork industry wastewater characterization assessment of the biodegradability, reuse and of the relationship between BOD, COD and tannins with TOC

Cork processing involves a boiling step to make the cork softer, which consumes a high volume of water and generates a wastewater with a high organic content, rich in tannins. An assessment of the final wastewater characteristics and of the boiling water composition along the boiling process was performed. The parameters studied were pH, color, total organic carbon (TOC), chemical and biochemical oxygen demands (COD, BOD5, BOD20), total suspended solids (TSS), total phenols and tannins (TP, TT). It was observed that the water solutes extraction power is significantly reduced for higher quantities of cork processed. Valid relationships between parameters were established not only envisaging wastewater characterization but also to provide an important tool for wastewater monitoring and for process control/optimization. Boiling water biodegradability presented decreasing values with the increase of cork processed and for the final wastewater its value is always lower than 0.5, indicating that these wastewaters are very difficult to treat by biological processes. The biodegradability was associated with the increase of tannin content that can rise up to 0.7 g/L. These compounds can be used by other industries when concentrated and the clarified wastewater can be reused, which is a potential asset in this wastewater treatment.

Optimização da estação de tratamento de águas residuais de Crestuma

O presente trabalho, elaborado entre Fevereiro e Julho de 2011, pretendeu analisar o impacto de algumas estratégias de optimização, adoptadas no âmbito da gestão de processo. As soluções estudadas tiveram como objectivo fazer face a alguns problemas de funcionamento da estação de tratamento de águas residuais (ETAR) de Crestuma, com recurso aos meios existentes. Para esse efeito foi realizada uma análise ao histórico dos registos de exploração referentes ao período de Janeiro de 2008 a Julho de 2011, pois foi considerado que este seria o período necessário para a análise e apuramento de resultados conclusivos. A ETAR fica situada na freguesia de Crestuma, actualmente com 2619 habitantes, localizada na confluência do regato de Vessadas com o rio Douro. Destina-se a tratar as águas residuais provenientes do sistema de drenagem de três localidades, S. Fioso, S. Picoto e S. Vessada, contribuindo para a despoluição do rio Douro. Genericamente, a ETAR tem apresentado um desempenho bastante satisfatório, respeitando os valores limites de emissão previstos na legislação actual (Decreto-Lei 152/97). A permanência de grandes quantidades de lama no processo, associada a uma elevada idade das lamas, está na origem da sua má sedimentabilidade, dificultando, por sua vez, a sua remoção do sistema. Com objectivo de se conseguir remover as lamas biológicas de uma forma mais eficiente, foi adoptada uma estratégia que compreendeu duas etapas. Com a finalidade de melhorar a sedimentabilidade, numa primeira etapa foi introduzido no processo um produto biológico, composto essencialmente por fungos, comercialmente conhecido por Optibiom 7450L. Com vista a obter-se uma maior eficiência na desidratação das lamas, a segunda etapa consistiu na substituição do polímero usado na lama espessada, por outro, de composição diferente, comercialmente conhecido por Easy 6040. Com as alterações introduzidas, foi possível obter um melhor controlo do processo e convergir para parâmetros de funcionamento mais próximos dos valores referidos na bibliografia. As alterações mais relevantes foram verificadas nos valores de sólidos suspensos totais (SST) no reactor biológico e na recirculação, tendo-se obtido valores médios de 3846 e 9716 mg/L, respectivamente (correspondendo a valores anteriores ao tratamento de 5301 e 7430 mg/L). Assim, a razão entre as concentrações de SST na corrente de recirculação e no reactor biológico passou de 1,4 para 2,6. A idade das lamas diminuiu de 42 para 36 dias com a adição do polímero Easy 6040 e o índice volumétrico de lamas apresentou uma redução, de 134 mL/g para 100 mL/g, o que se traduziu numa melhor sedimentabilidade das lamas. Desta forma, obteve-se uma diminuição de 33,3% no volume de lama produzida por m3 de afluente tratado, bem como uma redução de 27% no número de horas de funcionamento da centrífuga e um decréscimo ligeiro no consumo global de energia da ETAR (0,89%). Conclui-se que as alterações introduzidas no processo, adição de um produto biológico composto por fungos (Optibiom 7450L) e substituição do polímero usado na lama espessada por outro de composição diferente (Easy 6040), conduziram a uma melhoria no funcionamento do tratamento biológico da ETAR de Crestuma.