On minimal models of process systems

Minimal representations are known to have no redundant elements, and are therefore of great importance. Based on the notions of performance and size indices and measures for process systems, the paper proposes conditions for a process model being minimal in a set of functionally equivalent models with respect to a size norm. Generalized versions of known procedures to obtain minimal process models for a given modelling goal, model reduction based on sensitivity analysis and incremental model building are proposed and discussed. The notions and procedures are illustrated and compared on a simple example, that of a simple nonlinear fermentation process with different modelling goals and on a case study of a heat exchanger modelling. (C) 2004 Elsevier Ltd. All rights reserved.

Modelling of entrainment in industrial flotation cells: the effect of solids suspension

Solids concentration and particle size distribution gradually change in the vertical dimension of industrial flotation cells, subject primarily to the flotation cell size and design and the cell operating conditions. As entrainment is a two-step process and involves only the suspended solids in the top pulp region near the pulp-froth interface, the solids suspension characteristics have a significant impact on the overall entrainment. In this paper, a classification function is proposed to describe the state of solids suspension in flotation cells, similar to the definition of degree of entrainment for classification in the froth phase found in the literature. A mathematical model for solids suspension is also developed, in which the classification function is expressed as an exponential function of the particle size. Experimental data collected from three different Outokumpu tank flotation cells in three different concentrators are well fitted by the proposed exponential model. Under the prevailing experimental conditions, it was found that the solids content in the top region was relatively independent of cell operating conditions such as froth height and air rate but dependent on the cell size. Moreover, the results obtained from the total solids tend to be similar to those from a particular gangue mineral and hence may be applied to all minerals in entrainment calculation. (C) 2004 Elsevier Ltd. All rights reserved.

Rice morphogenesis and plant architecture: Measurement, specification and the reconstruction of structural development by 3D architectural modelling

Background and Aims The morphogenesis and architecture of a rice plant, Oryza sativa, are critical factors in the yield equation, but they are not well studied because of the lack of appropriate tools for 3D measurement. The architecture of rice plants is characterized by a large number of tillers and leaves. The aims of this study were to specify rice plant architecture and to find appropriate functions to represent the 3D growth across all growth stages. Methods A japonica type rice, 'Namaga', was grown in pots under outdoor conditions. A 3D digitizer was used to measure the rice plant structure at intervals from the young seedling stage to maturity. The L-system formalism was applied to create '3D virtual rice' plants, incorporating models of phenological development and leaf emergence period as a function of temperature and photoperiod, which were used to determine the timing of tiller emergence. Key Results The relationships between the nodal positions and leaf lengths, leaf angles and tiller angles were analysed and used to determine growth functions for the models. The '3D virtual rice' reproduces the structural development of isolated plants and provides a good estimation of the fillering process, and of the accumulation of leaves. Conclusions The results indicated that the '3D virtual rice' has a possibility to demonstrate the differences in the structure and development between cultivars and under different environmental conditions. Future work, necessary to reflect both cultivar and environmental effects on the model performance, and to link with physiological models, is proposed in the discussion.

Functional-structural plant modelling using a combination of architectural analysis, L-systems and a canonical model of function

This paper presents a new method for producing a functional-structural plant model that simulates response to different growth conditions, yet does not require detailed knowledge of underlying physiology. The example used to present this method is the modelling of the mountain birch tree. This new functional-structural modelling approach is based on linking an L-system representation of the dynamic structure of the plant with a canonical mathematical model of plant function. Growth indicated by the canonical model is allocated to the structural model according to probabilistic growth rules, such as rules for the placement and length of new shoots, which were derived from an analysis of architectural data. The main advantage of the approach is that it is relatively simple compared to the prevalent process-based functional-structural plant models and does not require a detailed understanding of underlying physiological processes, yet it is able to capture important aspects of plant function and adaptability, unlike simple empirical models. This approach, combining canonical modelling, architectural analysis and L-systems, thus fills the important role of providing an intermediate level of abstraction between the two extremes of deeply mechanistic process-based modelling and purely empirical modelling. We also investigated the relative importance of various aspects of this integrated modelling approach by analysing the sensitivity of the standard birch model to a number of variations in its parameters, functions and algorithms. The results show that using light as the sole factor determining the structural location of new growth gives satisfactory results. Including the influence of additional regulating factors made little difference to global characteristics of the emergent architecture. Changing the form of the probability functions and using alternative methods for choosing the sites of new growth also had little effect. (c) 2004 Elsevier B.V. All rights reserved.

Subgeometric rates of convergence for a class of continuous-time Markov process

Let (Phi(t))(t is an element of R+) be a Harris ergodic continuous-time Markov process on a general state space, with invariant probability measure pi. We investigate the rates of convergence of the transition function P-t(x, (.)) to pi; specifically, we find conditions under which r(t) vertical bar vertical bar P-t (x, (.)) - pi vertical bar vertical bar -> 0 as t -> infinity, for suitable subgeometric rate functions r(t), where vertical bar vertical bar - vertical bar vertical bar denotes the usual total variation norm for a signed measure. We derive sufficient conditions for the convergence to hold, in terms of the existence of suitable points on which the first hitting time moments are bounded. In particular, for stochastically ordered Markov processes, explicit bounds on subgeometric rates of convergence are obtained. These results are illustrated in several examples.

Modelling of two-stage anaerobic digestion using the IWA Anaerobic Digestion Model No. 1 (ADM1)

The aim of the study presented was to implement a process model to simulate the dynamic behaviour of a pilot-scale process for anaerobic two-stage digestion of sewage sludge. The model implemented was initiated to support experimental investigations of the anaerobic two-stage digestion process. The model concept implemented in the simulation software package MATLAB(TM)/Simulink(R) is a derivative of the IWA Anaerobic Digestion Model No.1 (ADM1) that has been developed by the IWA task group for mathematical modelling of anaerobic processes. In the present study the original model concept has been adapted and applied to replicate a two-stage digestion process. Testing procedures, including balance checks and 'benchmarking' tests were carried out to verify the accuracy of the implementation. These combined measures ensured a faultless model implementation without numerical inconsistencies. Parameters for both, the thermophilic and the mesophilic process stage, have been estimated successfully using data from lab-scale experiments described in literature. Due to the high number of parameters in the structured model, it was necessary to develop a customised procedure that limited the range of parameters to be estimated. The accuracy of the optimised parameter sets has been assessed against experimental data from pilot-scale experiments. Under these conditions, the model predicted reasonably well the dynamic behaviour of a two-stage digestion process in pilot scale. (C) 2004 Elsevier Ltd. All rights reserved.

Population-based continuous optimization, probabilistic modelling and mean shift

Evolutionary algorithms perform optimization using a population of sample solution points. An interesting development has been to view population-based optimization as the process of evolving an explicit, probabilistic model of the search space. This paper investigates a formal basis for continuous, population-based optimization in terms of a stochastic gradient descent on the Kullback-Leibler divergence between the model probability density and the objective function, represented as an unknown density of assumed form. This leads to an update rule that is related and compared with previous theoretical work, a continuous version of the population-based incremental learning algorithm, and the generalized mean shift clustering framework. Experimental results are presented that demonstrate the dynamics of the new algorithm on a set of simple test problems.

Mechanistic modelling of pilling. Part II: Individual-fibre computational model

A one-dimensional computational model of pilling of a fibre assembly has been created. The model follows a set of individual fibres, as free ends and loops appear as fuzz and arc progressively withdrawn from the body of the assembly, and entangle to form pills, which eventually break off or are pulled out. The time dependence of the computation is given by ticks, which correspond to cycles of a wear and laundering process. The movement of the fibres is treated as a reptation process. A set of standard values is used as inputs to the computation. Predictions arc given of the change with a number Of cycles of mass of fuzz, mass of pills, and mass removed from the assembly. Changes in the standard values allow sensitivity studies to be carried out.

Modelling comminution patterns within a pilot scale AG/SAG mill

The patterns of rock comminution within tumbling mills, as well as the nature of forces, are of significant practical importance. Discrete element modelling (DEM) has been used to analyse the pattern of specific energy applied to rock, in terms of spatial distribution within a pilot AG/SAG mill. We also analysed in some detail the nature of the forces, which may result in rock comminution. In order to examine the distribution of energy applied within the mill, the DEM models were compared with measured particle mass losses, in small scale AG and SAG mill experiments. The intensity of contact stresses was estimated using the Hertz theory of elastic contacts. The results indicate that in the case of the AG mill, the highest intensity stresses and strains are likely to occur deep within the charge, and close to the base. This effect is probably more pronounced for large AG mills. In the SAG mill case, the impacts of the steel balls on the surface of the charge are likely to be the most potent. In both cases, the spatial pattern of medium-to-high energy collisions is affected by the rotational speed of the mill. Based on an assumed damage threshold for rock, in terms of specific energy introduced per single collision, the spatial pattern of productive collisions within each charge was estimated and compared with rates of mass loss. We also investigated the nature of the comminution process within AG vs. SAG mill, in order to explain the observed differences in energy utilisation efficiency, between two types of milling. All experiments were performed using a laboratory scale mill of 1.19 m diameter and 0.31 m length, equipped with 14 square section lifters of height 40 mm. (C) 2006 Elsevier Ltd. All rights reserved.

Modelling of entrainment in industrial flotation cells: Water recovery and degree of entrainment

Entrainment in flotation can be considered as a two-step process, including the transfer of the suspended solids in the top of the pulp region just below the pulp-froth interface to the froth phase and the transfer of the entrained particles in the froth phase to the concentrate. Both steps have a strong classification characteristic. The degree of entrainment describes the classification effect of the drainage process in the froth phase. This paper briefly reviews two existing models of degree of entrainment. Experimental data were collected from an Outokumpu 3 m(3) tank cell in the Xstrata Mt. Isa Mines copper concentrator. The data are fitted to the models and the effect of cell operating conditions including air rate and froth height on the degree of entrainment is examined on a size-by-size basis. It is found that there is a strong correlation between the entrainment and the water recovery, which is close to lineal. for the fines. The degree of entrainment decreases with increase in particle size. Within the normal range of cell operating conditions, few particles coarser than 50 mu m are recovered by entrainment. In general, the degree of entrainment increases with increase in the ail rate and decreases with increase in the froth height. Air rate and froth height strongly interact with each other and affect the entrainment process mainly via changes in the froth retention time, the froth structure and froth properties. As a result, other mechanisms such as entrapment may become important in recovering the coarse entrained particles. (c) 2005 Elsevier Ltd. All rights reserved.

Cognitive modelling and the behaviour genetics of reading

While it is well known that reading is highly heritable, less has been understood about the bases of these genetic influences. In this paper, we review the research that we have been conducting in recent years to examine genetic and environmental influences on the particular reading processes specified in the dual-route cognitive model of reading. We argue that a detailed understanding of the role of genetic factors in reading acquisition requires the delineation and measurement of precise phenotypes, derived from well-articulated models of the reading process. We report evidence for independent genetic influences on the lexical and nonlexical reading processes represented in the dual-route model, based on studies of children with particular subtypes of dyslexia, and on univariate and multivariate genetic modelling of reading performance in the normally reading population.

Cognitive change process during group cognitive behaviour therapy for depression

Background: This study extended that of Kwon and Oei [Kwon, S.M., Oei, T.P.S., 2003. Cognitive change processes in a group cognitive behavior therapy of depression. J. Behav. Ther. Exp. Psychiatry, 3, 73-85], which outlined a number of testable models based on Beck's cognitive theory of depression. Specifically, the current study tested the following four competing models: the causal, consequential, fully and partially interactive cognitive models in patients with major depressive disorder. Methods: A total of 168 clinically depressed outpatients were recruited into a 12-week group cognitive behaviour therapy program. Data was collected at three time points: baseline, mid- and at termination of therapy using the ATQ DAS and BD1. The data were analysed with Amos 4.01 (Arbuckle, J.L., 1999. Amos 4.1. Smallwaters, Chicago.) structural equation modelling. Results: Results indicated that dysfunctional attitudes, negative automatic thoughts and symptoms of depression reduced significantly during treatment. Both the causal and consequential models equally provided an adequate fit to the data. The fully interactive model provided the best fit. However, after removing non-significant pathways, it was found that reduced depressive symptom contributed to reduced depressogenic automatic thoughts and dysfunctional attitudes, not the reverse. Conclusion: These findings did not fully support Beck's cognitive theory of depression that cognitions are primary in the reduction of depressed mood. (c) 2006 Elsevier B.V. All rights reserved.

Extinction times for a birth-death process with two phases

Many populations have a negative impact on their habitat or upon other species in the environment if their numbers become too large. For this reason they are often subjected to some form of control. One common control regime is the reduction regime: when the population reaches a certain threshold it is controlled (for example culled) until it falls below a lower predefined level. The natural model for such a controlled population is a birth-death process with two phases, the phase determining which of two distinct sets of birth and death rates governs the process. We present formulae for the probability of extinction and the expected time to extinction, and discuss several applications. (c) 2006 Elsevier Inc. All rights reserved.

Multidimensional modelling to investigate interspecies hydrogen transfer in anaerobic biofilms

Anaerobic digestion is a multistep process, mediated by a functionally and phylogenetically diverse microbial population. One of the crucial steps is oxidation of organic acids, with electron transfer via hydrogen or formate from acetogenic bacteria to methanogens. This syntrophic microbiological process is strongly restricted by a thermodynamic limitation on the allowable hydrogen or formate concentration. In order to study this process in more detail, we developed an individual-based biofilm model which enables to describe the processes at a microbial resolution. The biochemical model is the ADM1, implemented in a multidimensional domain. With this model, we evaluated three important issues for the syntrophic relationship: (i) is there a fundamental difference in using hydrogen or formate as electron carrier? (ii) Does a thermodynamic-based inhibition function produced substantially different results from an empirical function? and; (iii) Does the physical colocation of acetogens and methanogens follow directly from a general model. Hydrogen or formate as electron carrier had no substantial impact on model results. Standard inhibition functions or thermodynamic inhibition function gave similar results at larger substrate field grid sizes (> 10 mu m), but at smaller grid sizes, the thermodynamic-based function reduced the number of cells with long interspecies distances (> 2.5 mu m). Therefore, a very fine grid resolution is needed to reflect differences between the thermodynamic function, and a more generic inhibition form. The co-location of syntrophic bacteria was well predicted without a need to assume a microbiological based mechanism (e.g., through chemotaxis) of biofilm formation.

Oscillatory regulation of hes1: Discrete stochastic delay modelling and simulation

Discrete stochastic simulations are a powerful tool for understanding the dynamics of chemical kinetics when there are small-to-moderate numbers of certain molecular species. In this paper we introduce delays into the stochastic simulation algorithm, thus mimicking delays associated with transcription and translation. We then show that this process may well explain more faithfully than continuous deterministic models the observed sustained oscillations in expression levels of hes1 mRNA and Hes1 protein.