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.

Attitudes and practices of general practitioners in the diagnosis and management of attention-deficit/hyperactivity disorder

Objective: To assess understanding of, and actual and potential roles in management of attention-deficit/hyperactivity disorder (ADHD) among GPs. Methods: A cross-sectional questionnaire survey of Queensland GPs selected randomly from the Royal Australian College of General Practitioners directory of members was carried out. Main outcome measures were knowledge levels of ADHD, current management practices, referral patterns and self-perceived information and training needs. Results: Three hundred and ninety-nine GPs returned a completed questionnaire (response rate 76%). Roles identified by GPs were: the provisional diagnosis of ADHD and referral to specialist services for confirmation of the diagnosis and initiation of management; assistance with monitoring progress once a management plan was in place; education of the child and their family regarding the disorder; and liaison with the school where necessary. Perceived barriers to increased involvement of GPs were: time and resource constraints of general practice; concerns regarding abuse and addiction liability of prescription stimulants; complex diagnostic issues associated with childhood behavioural problems; and lack of training and education regarding ADHD. Conclusions: General practitioners identify a role for themselves in ADHD care that is largely supportive in nature and involves close liaison with specialist services.

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.

General practitioner attitudes to case conferences: how can we increase participation and effectiveness?

Objectives: To identify general practitioners' views on the barriers to using case conferencing (as outlined in the Medical Benefits Schedule (MBS) Enhanced Primary Care package) and to develop a set of principles to encourage greater GP participation in case conferences. Design: Qualitative study, involving semistructured questions administered to focus groups of GPs, conducted between April and July 2001 as part of a broader study of case coordination in palliative care. Participants: 29 GPs from urban, regional, and rural areas of Queensland. Principal findings: Many of the GPs' work practices militated against participation in traditionally structured case conferences. GPs thought the range of MBS item numbers should be expanded to cover alternative methods of liaison (eg, phone consultations with other service providers). The onerous bureaucratic processes required to claim reimbursement were an additional disincentive. Conclusions: GPs would probably be more likely to participate in case conferences if they were initiated by specialist services and arranged more flexibly to suit GP work schedules.

Infusing the use of seasonal climate forecasting into crop management practice in North East Australia using discussion support software

Seasonal climate forecasting offers potential for improving management of crop production risks in the cropping systems of NE Australia. But how is this capability best connected to management practice? Over the past decade, we have pursued participative systems approaches involving simulation-aided discussion with advisers and decision-makers. This has led to the development of discussion support software as a key vehicle for facilitating infusion of forecasting capability into practice. In this paper, we set out the basis of our approach, its implementation and preliminary evaluation. We outline the development of the discussion support software Whopper Cropper, which was designed for, and in close consultation with, public and private advisers. Whopper Cropper consists of a database of simulation output and a graphical user interface to generate analyses of risks associated with crop management options. The charts produced provide conversation pieces for advisers to use with their farmer clients in relation to the significant decisions they face. An example application, detail of the software development process and an initial survey of user needs are presented. We suggest that discussion support software is about moving beyond traditional notions of supply-driven decision support systems. Discussion support software is largely demand-driven and can compliment participatory action research programs by providing cost-effective general delivery of simulation-aided discussions about relevant management actions. The critical role of farm management advisers and dialogue among key players is highlighted. We argue that the discussion support concept, as exemplified by the software tool Whopper Cropper and the group processes surrounding it, provides an effective means to infuse innovations, like seasonal climate forecasting, into farming practice. Crown Copyright (C) 2002 Published by Elsevier Science Ltd. All rights reserved.

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.

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.