A systematic approach to error isolation in computerized wastewater simulation models

Activated sludge models are used extensively in the study of wastewater treatment processes. While various commercial implementations of these models are available, there are many people who need to code models themselves using the simulation packages available to them, Quality assurance of such models is difficult. While benchmarking problems have been developed and are available, the comparison of simulation data with that of commercial models leads only to the detection, not the isolation of errors. To identify the errors in the code is time-consuming. In this paper, we address the problem by developing a systematic and largely automated approach to the isolation of coding errors. There are three steps: firstly, possible errors are classified according to their place in the model structure and a feature matrix is established for each class of errors. Secondly, an observer is designed to generate residuals, such that each class of errors imposes a subspace, spanned by its feature matrix, on the residuals. Finally. localising the residuals in a subspace isolates coding errors. The algorithm proved capable of rapidly and reliably isolating a variety of single and simultaneous errors in a case study using the ASM 1 activated sludge model. In this paper a newly coded model was verified against a known implementation. The method is also applicable to simultaneous verification of any two independent implementations, hence is useful in commercial model development.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Using clusters of computers for large QU-GENE simulation experiments

The QU-GENE Computing Cluster (QCC) is a hardware and software solution to the automation and speedup of large QU-GENE (QUantitative GENEtics) simulation experiments that are designed to examine the properties of genetic models, particularly those that involve factorial combinations of treatment levels. QCC automates the management of the distribution of components of the simulation experiments among the networked single-processor computers to achieve the speedup.

Fluorescence spectrum of a two-level atom driven by a multiple modulated field

We investigate the fluorescence spectrum of a two-level atom driven by a multiple amplitude-modulated field. The driving held is modeled as a polychromatic field composed of a strong central (resonant) component and a large number of symmetrically detuned sideband fields displaced from the central component by integer multiples of a constant detuning. Spectra obtained here differ qualitatively from those observed for a single pair of modulating fields [B. Blind, P.R. Fontana, and P. Thomann, J. Phys. B 13, 2717 (1980)]. In the case of a small number of the modulating fields, a multipeaked spectrum is obtained with the spectral features located at fixed frequencies that are independent of the number of modulating fields and their Rabi frequencies. As the number of the modulating fields increases, the spectrum ultimately evolves to the well-known Mellow triplet with the sidebands shifted from the central component by an effective Rabi frequency whose magnitude depends on the initial relative phases of the components of the driving held. For equal relative phases, the effective Rabi frequency of the driving field can be reduced to zero resulting in the disappearance of fluorescence spectrum, i.e., the atom can stop interacting with the field. When the central component and the modulating fields are 180 degrees out of phase, the spectrum retains its triplet structure with the sidebands located at frequencies equal to the sum of the Rabi frequencies of the component of the driving field. Moreover, we shaw that the frequency of spontaneous emission can be controlled and switched from one frequency to another when the Rabi frequency or initial phase of the modulating fields are varied.

Effect of quantum interference on a three-level atom driven by two laser fields

We study the effect of quantum interference on the population distribution and absorptive properties of a V-type three-level atom driven by two lasers of unequal intensities and different angular frequencies. Three coupling configurations of the lasers to the atom are analysed: (a) both lasers coupled to the same atomic transition, (b) each laser coupled to different atomic transition and (c) each laser coupled to both atomic transitions. Dressed stales for the three coupling configurations are identified, and the population distribution and absorptive properties of the weaker field are interpreted in terms of transition dipole moments and transition frequencies among these dressed states. In particular, we find that in the first two cases there is no population inversion between the bare atomic states, but the population can be trapped in a superposition of the dressed states induced by quantum interference and the stronger held. We show that the trapping of the population, which results from the cancellation of transition dipole moments, does not prevent the weaker field to be coupled to the cancelled (dark) transitions. As a result, the weaker field can be strongly amplified on transparent transitions. In the case of each laser coupled to both atomic transitions the population can be trapped in a linear superposition of the excited bare atomic states leaving the ground state unpopulated in the steady state. Moreover, we find that the absorption rate of the weaker field depends on the detuning of the strong field from the atomic resonances and the splitting between the atomic excited states. When the strong held is resonant to one of the atomic transitions a quasi-trapping effect appears in one of the dressed states. In the quasi-trapping situation all the transition dipole moments are different from zero, which allows the weaker field to be amplified on the inverted transitions. When the strong field is tuned halfway between the atomic excited states, the population is completely trapped in one of the dressed states and no amplification is found for the weaker field.

Amplification on dark transitions in a driven three-level atom

We study a three-level atomic system of the vee type, but driven on only one transition by a monochromatic laser. It is shown that the gain of a probe beam, recently predicted for this system by Menon and Agarwal (Menon S and Agarwal G 2000 Phys. Rev. A 61 13 807), is due to an unexpected amplification on a completely inverted, nondecaying (dark) transition. This prediction violates the well known balance condition between the population inversion and the coupling strength of the probe field to the inverted transition, which requires that the coupling strength reduces with increasing population inversion. We show that the condition may be violated only if the probe field selectively couples to just one of the atomic transitions: when it couples to both transitions, the balance condition is satisfied and the system is transparent for the probe field coupled to the dark transitions. No amplification is possible in the latter case.

Optically driven micromachine elements

We report on a proof of principle demonstration of an optically driven micromachine element. Optical angular momentum is transferred from a circularly polarized laser beam to a birefringent particle confined in an optical tweezers trap. The optical torque causes the particle to spin at up to 350 Hz, and this torque is harnessed to drive an optically trapped microfabricated structure. We describe a photolithographic method for producing the microstructures and show how a light driven motor could be used in a micromachine system. (C) 2001 American Institute of Physics.

A simulation model for analysing patient activity in dermatology

We developed a general model to assess patient activity within the primary and secondary health-care sectors following a dermatology outpatient consultation. Based on observed variables from the UK teledermatology trial, the model showed that up to 11 doctor-patient interactions occurred before a patient was ultimately discharged from care. In a cohort of 1000 patients, the average number of health-care visits was 2.4 (range 1-11). Simulation analysis suggested that the most important parameter affecting the total number of doctor-patient Interactions is patient discharge from care following the initial consultation. This implies that resources should be concentrated in this area. The introduction of teledermatology (either realtime or store and forward) changes the values of the model parameters. The model provides a quantitative tool for planning the future provision of dermatology health-care.

The semantic interference effect in the picture-word paradigm: An event-related fMRI study employing overt responses

We used event-related functional magnetic resonance imaging (fMRI) to investigate neural responses associated with the semantic interference (SI) effect in the picture-word task. Independent stage models of word production assume that the locus of the SI effect is at the conceptual processing level (Levelt et al. [1999]: Behav Brain Sci 22:1-75), whereas interactive models postulate that it occurs at phonological retrieval (Starreveld and La Heij [1996]: J Exp Psychol Learn Mem Cogn 22:896-918). In both types of model resolution of the SI effect occurs as a result of competitive, spreading activation without the involvement of inhibitory links. These assumptions were tested by randomly presenting participants with trials from semantically-related and lexical control distractor conditions and acquiring image volumes coincident with the estimated peak hemodynamic response for each trial. Overt vocalization of picture names occurred in the absence of scanner noise, allowing reaction time (RT) data to be collected. Analysis of the RT data confirmed the SI effect. Regions showing differential hemodynamic responses during the SI effect included the left mid section of the middle temporal gyrus, left posterior superior temporal gyrus, left anterior cingulate cortex, and bilateral orbitomedial prefrontal cortex. Additional responses were observed in the frontal eye fields, left inferior parietal lobule, and right anterior temporal and occipital cortex. The results are interpreted as indirectly supporting interactive models that allow spreading activation between both conceptual processing and phonological retrieval levels of word production. In addition, the data confirm that selective attention/response suppression has a role in resolving the SI effect similar to the way in which Stroop interference is resolved. We conclude that neuroimaging studies can provide information about the neuroanatomical organization of the lexical system that may prove useful for constraining theoretical models of word production. (C) 2001 Wiley-Liss, Inc.

Initiation of cell locomotility is a morphogenetic checkpoint in thyroid epithelial cells regulated by ERK and PI3-kinase signals

Epithelial locomotility is a fundamental determinant of tissue patterning that is subject to strict physiological regulation. The current, study sought to identify cellular signals that initiate cell migration in cultured thyroid epithelial cells. Porcine thyroid cells cultured as 3-dimensional follicles convert to 2-dimensional monolayers when deprived of agents that stimulate cAMP/PKA signaling. This morphogenetic event is driven by the activation of cell-on-substrate locomotility, providing a convenient assay for events that regulate the initiation of locomotion. In this system, the extracellular signal regulated kinase (ERK) pathway became activated as follicles converted to monolayer, as demonstrated by immunoblotting for activation-specific phosphorylation and nuclear accumulation of ERK. Inhibition of ERK activation using the drug PD98059 effectively prevented cells from beginning to migrate. PD98059 inhibited cell spreading, actin filament reorganization and the assembly of focal adhesions, cellular events that mediate the initiation of thyroid cell locomotility. Akt (PKB) signaling was also activated during follicle-to-monolayer conversion and the phosphoinositide 3-kinase (PI3-kinase) inhibitor, wortmannin, also blocked the initiation of cell movement. Wortmannin did not, however, block activation of ERK signaling. These findings, therefore, identify the ERK and PI3-kinase signaling pathways as important stimulators of thyroid cell locomotility. These findings are incorporated into a model where the initiation of thyroid cell motility constitutes a morphogenetic checkpoint regulated by coordinated changes in stimulatory (ERK, PI3-kinase) and tonic inhibitory (cAMP/PKA) signaling pathways. Cell Motil. Cytoskeleton 49:93-103, 2001. (C) 2001 Wiley-Liss, Inc.

Increasing profits and reducing risks in crop production using participatory systems simulation approaches

The development of cropping systems simulation capabilities world-wide combined with easy access to powerful computing has resulted in a plethora of agricultural models and consequently, model applications. Nonetheless, the scientific credibility of such applications and their relevance to farming practice is still being questioned. Our objective in this paper is to highlight some of the model applications from which benefits for farmers were or could be obtained via changed agricultural practice or policy. Changed on-farm practice due to the direct contribution of modelling, while keenly sought after, may in some cases be less achievable than a contribution via agricultural policies. This paper is intended to give some guidance for future model applications. It is not a comprehensive review of model applications, nor is it intended to discuss modelling in the context of social science or extension policy. Rather, we take snapshots around the globe to 'take stock' and to demonstrate that well-defined financial and environmental benefits can be obtained on-farm from the use of models. We highlight the importance of 'relevance' and hence the importance of true partnerships between all stakeholders (farmer, scientists, advisers) for the successful development and adoption of simulation approaches. Specifically, we address some key points that are essential for successful model applications such as: (1) issues to be addressed must be neither trivial nor obvious; (2) a modelling approach must reduce complexity rather than proliferate choices in order to aid the decision-making process (3) the cropping systems must be sufficiently flexible to allow management interventions based on insights gained from models. The pro and cons of normative approaches (e.g. decision support software that can reach a wide audience quickly but are often poorly contextualized for any individual client) versus model applications within the context of an individual client's situation will also be discussed. We suggest that a tandem approach is necessary whereby the latter is used in the early stages of model application for confidence building amongst client groups. This paper focuses on five specific regions that differ fundamentally in terms of environment and socio-economic structure and hence in their requirements for successful model applications. Specifically, we will give examples from Australia and South America (high climatic variability, large areas, low input, technologically advanced); Africa (high climatic variability, small areas, low input, subsistence agriculture); India (high climatic variability, small areas, medium level inputs, technologically progressing; and Europe (relatively low climatic variability, small areas, high input, technologically advanced). The contrast between Australia and Europe will further demonstrate how successful model applications are strongly influenced by the policy framework within which producers operate. We suggest that this might eventually lead to better adoption of fully integrated systems approaches and result in the development of resilient farming systems that are in tune with current climatic conditions and are adaptable to biophysical and socioeconomic variability and change. (C) 2001 Elsevier Science Ltd. All rights reserved.