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.

Application of the maximum entropy method to the (2 + 1)D four-fermion model

We investigate spectral functions extracted using the maximum entropy method from correlators measured in lattice simulations of the (2+1)-dimensional four-fermion model. This model is particularly interesting because it has both a chirally broken phase with a rich spectrum of mesonic bound states and a symmetric phase where there are only resonances. In the broken phase we study the elementary fermion, pion, sigma, and massive pseudoscalar meson; our results confirm the Goldstone nature of the π and permit an estimate of the meson binding energy. We have, however, seen no signal of σ→ππ decay as the chiral limit is approached. In the symmetric phase we observe a resonance of nonzero width in qualitative agreement with analytic expectations; in addition the ultraviolet behavior of the spectral functions is consistent with the large nonperturbative anomalous dimension for fermion composite operators expected in this model.

Development of a generic crop model template in the cropping system model APSIM

The Agricultural Production Systems slMulator, APSIM, is a cropping system modelling environment that simulates the dynamics of soil-plant-management interactions within a single crop or a cropping system. Adaptation of previously developed crop models has resulted in multiple crop modules in APSIM, which have low scientific transparency and code efficiency. A generic crop model template (GCROP) has been developed to capture unifying physiological principles across crops (plant types) and to provide modular and efficient code for crop modelling. It comprises a standard crop interface to the APSIM engine, a generic crop model structure, a crop process library, and well-structured crop parameter files. The process library contains the major science underpinning the crop models and incorporates generic routines based on physiological principles for growth and development processes that are common across crops. It allows APSIM to simulate different crops using the same set of computer code. The generic model structure and parameter files provide an easy way to test, modify, exchange and compare modelling approaches at process level without necessitating changes in the code. The standard interface generalises the model inputs and outputs, and utilises a standard protocol to communicate with other APSIM modules through the APSIM engine. The crop template serves as a convenient means to test new insights and compare approaches to component modelling, while maintaining a focus on predictive capability. This paper describes and discusses the scientific basis, the design, implementation and future development of the crop template in APSIM. On this basis, we argue that the combination of good software engineering with sound crop science can enhance the rate of advance in crop modelling. Crown Copyright (C) 2002 Published by Elsevier Science B.V. All rights reserved.

Polyethylene flow prediction with a differential multi-mode Pom-Pom model

We report the first steps of a collaborative project between the University of Queensland, Polyflow, Michelin, SK Chemicals, and RMIT University; on simulation, validation and application of a recently introduced constitutive model designed to describe branched polymers. Whereas much progress has been made on predicting the complex flow behaviour of many - in particular linear - polymers, it sometimes appears difficult to predict simultaneously shear thinning and extensional strain hardening behaviour using traditional constitutive models. Recently a new viscoelastic model based on molecular topology, was proposed by McLeish and Larson (1998). We explore the predictive power of a differential multi-mode version of the pom-pom model for the flow behaviour of two commercial polymer melts: a (long-chain branched) low-density polyethylene (LDPE) and a (linear) high-density polyethylene (HDPE). The model responses are compared to elongational recovery experiments published by Langouche and Debbaut (1999), and start-up of simple shear flow, stress relaxation after simple and reverse step strain experiments carried out in our laboratory.

Application of the Williams-Landel-Ferry model to the viscosity-temperature relationship of Australian honeys

The rheological behaviour of nine unprocessed Australian honeys was investigated for the applicability of the Williams-Landel-Ferry (WLF) model. The viscosity of the honeys was obtained over a range of shear rates (0.01-40 s(-1)) from 2degrees to 40 degreesC, and all the honeys exhibited Newtonian behaviour with viscosity reducing as the temperature was increased. The honeys with high moisture were of lower viscosity, The glass transition temperatures of the honeys, as measured with a differential scanning calorimeter (DSC), ranged from -40degrees to -46 degreesC, and four models (WLF. Arrhenius, Vogel-Tammann-Fulcher (VTF), and power-law) were investigated to describe the temperature dependence of the viscosity. The WLF was the most suitable and the correlation coefficient averaged 0.999 +/- 0.0013 as against 0.996 +/- 0.0042 for the Arrhenius model while the mean relative deviation modulus was 0-12% for the WLF model and 10-40% for the Arrhenius one. With the universal values for the WLF constants, the temperature dependence of the viscosity was badly predicted. From non-linear regression analysis, the constants of the WLF models for the honeys were obtained (C-1 = 13.7-21.1: C-2 = 55.9-118.7) and are different from the universal values. These WLF constants will be valuable for adequate modeling of the rheology of the honeys, and they can be used to assess the temperature sensitivity of the honeys. (C) 2002 Elsevier Science Ltd. All rights reserved.

Distribution kinetics of solutes in the isolated in-situ perfused rat head using the multiple indicator dilution technique and a physiological two-barrier model

The purpose of this study was to determine the pharmacokinetics of [C-14]diclofenac, [C-14]salicylate and [H-3]clonidine using a single pass rat head perfusion preparation. The head was perfused with 3-[N-morpholino] propane-sulfonic acid-buffered Ringer's solution. Tc-99m-red blood cells and a drug were injected in a bolus into the internal carotid artery and collected from the posterior facial vein over 28 min. A two-barrier stochastic organ model was used to estimate the statistical moments of the solutes. Plasma, interstitial and cellular distribution volumes for the solutes ranged from 1.0 mL (diclofenac) to 1.6 mL (salicylate), 2.0 mL (diclofenac) to 4.2 mL (water) and 3.9 mL (salicylate) to 20.9 mL (diclofenac), respectively. A comparison of these volumes to water indicated some exclusion of the drugs from the interstitial space and salicylate from the cellular space. Permeability-surface area (PS) products calculated from plasma to interstitial fluid permeation clearances (CLPI) (range 0.02-0.40 mL s(-1)) and fractions of solute unbound in the perfusate were in the order: diclofenac>salicylate >clonidine>sucrose (from 41.8 to 0.10 mL s(-1)). The slow efflux of diclofenac, compared with clonidine and salicylate, may be related to its low average unbound fraction in the cells. This work accounts for the tail of disposition curves in describing pharmacokinetics in the head.

D2 dopamine receptor gene polymorphism discriminates two kinds of novelty seeking

Cloninger's psychobiological model of personality as applied to substance misuse has received mixed support. Contrary to the model, recent data suggest that a combination of high novelty seeking (NS) and high harm avoidance (HA) represents a significant risk for the development of severe substance misuse. A genetic polymorphism previously implicated in severe substance dependence, the A1 allele of the D2 dopamine receptor (DRD2) gene, was examined in relation to NS and HA amongst 203 adolescent boys. Specifically, we hypothesized that subjects with the A1 + allele (A1/A1 and A1/A2 genotypes) would report stronger NS and would exhibit a more positive relationship between NS and HA than those with the A1-allele (A2/A2 genotypes). These predictions were supported. The correlation between NS and HA in 81 A1 + allelic boys (r = 0.27, P = 0.02), and that in the 122 A1- allelic boys (r = -0.15, P = 0.09), indicated that this relationship differed according to allelic status (F = 8.52, P < 0:004). Among those with the A1-allele, the present results are consistent with the traditional view that novelty seeking provides positive reinforcement, or the fulfillment of appetitive drives. In contrast, novelty seeking in those with the A1 + allele appears to include a negative reinforcement or self-medicating function. (C) 2002 Elsevier Science Ltd. All rights reserved.

The tourism-leisure behavioural continuum

Tourism and leisure have historically been studied in isolation, but recent researchers attempt to understand the relationship between them. However, there have been few comparisons of how people behave in their environments. This paper provides a conceptual model of the relationship based oil a critical review of the current literature oil pleasure-oriented environments. The model indicates that despite potential differences in observable behavior, there is a commonality between the underlying influences that define how people behave during their leisure and tourism experiences. Therefore, the paper concludes that these should not be regarded as separate phenomenon. (C) 2002 Elsevier Science Ltd. All Lights reserved.

Cyclosporin A pretreatment in a rat model of warm ischaemia/reperfusion injury

Background/Aims: These studies investigated the role of apoptosis following ischaemia/reperfusion (I/R) injury to the liver and the effect of pretreatment with Cyclosporin A. Methods: Male Sprague-Dawley rats received 30 min of warm ischaemia followed by a period of reperfusion of 6 h. Rats were given olive oil or Cyclosporin A (30 mg/kg p.o.) the day before surgery. Neutrophil numbers were assessed in haematoxylin-eosin-stained sections of liver. In situ staining of sections using TdT-mediated dUTP-fluoreseein nick-end labelling was carried out to determine the extent of apoptosis, followed by electron microscopy. Semi-quantitative polymerase chain reaction (PCR) analysis of the transcript for Fas antigen was performed. Results and Conclusions: High levels of apoptosis were observed in I/R injury, which were greatly ameliorated in Cyclosporin A-pretreated groups. PCR analysis indicated a reduction in the level of expression of Fas transcript in Cyclosporin A-treated rats. Histological analysis showed a significant increase in the number of neutrophils infiltrating I/R-injured tissue (62 +/- 10.69, it = 16), which was markedly reduced by Cyclosporin A pretreatment (16 +/- 7, n = 6, P < 0.05). These results indicate a role of parenchymal apoptosis in the pathogenesis of I/R injury, which occurs in association with neutrophil infiltration, both of which can be significantly reduced by Cyclosporin A pretreatment. (C) 2002 European Association for the Study of the Liver. Published by Elsevier Science B.V. All rights reserved.

The E(NK) model: extending the NK model to incorporate gene-by-environment interactions and epistasis for diploid genomes

The haploid NK model developed by Kauffman can be extended to diploid genomes and to incorporate gene-by-environment interaction effects in combination with epistasis. To provide the flexibility to include a wide range of forms of gene-by-environment interactions, a target population of environment types (TPE) is defined. The TPE consists of a set of E different environment types, each with their own frequency of occurrence. Each environment type conditions a different NK gene network structure or series of gene effects for a given network structure, providing the framework for defining gene-by-environment interactions. Thus, different NK models can be partially or completely nested within the E environment types of a TPE, giving rise to the E(NK) model for a biological system. With this model it is possible to examine how populations of genotypes evolve in context with properties of the environment that influence the contributions of genes to the fitness values of genotypes. We are using the E(NK) model to investigate how both epistasis and gene-by-environment interactions influence the genetic improvement of quantitative traits by plant breeding strategies applied to agricultural systems. © 2002 Wiley Periodicals, Inc.