Analysis of dynamic process models for structural insight and model reduction .1. Structural identification measures

An important consideration in the development of mathematical models for dynamic simulation, is the identification of the appropriate mathematical structure. By building models with an efficient structure which is devoid of redundancy, it is possible to create simple, accurate and functional models. This leads not only to efficient simulation, but to a deeper understanding of the important dynamic relationships within the process. In this paper, a method is proposed for systematic model development for startup and shutdown simulation which is based on the identification of the essential process structure. The key tool in this analysis is the method of nonlinear perturbations for structural identification and model reduction. Starting from a detailed mathematical process description both singular and regular structural perturbations are detected. These techniques are then used to give insight into the system structure and where appropriate to eliminate superfluous model equations or reduce them to other forms. This process retains the ability to interpret the reduced order model in terms of the physico-chemical phenomena. Using this model reduction technique it is possible to attribute observable dynamics to particular unit operations within the process. This relationship then highlights the unit operations which must be accurately modelled in order to develop a robust plant model. The technique generates detailed insight into the dynamic structure of the models providing a basis for system re-design and dynamic analysis. The technique is illustrated on the modelling for an evaporator startup. Copyright (C) 1996 Elsevier Science Ltd

A peptide-binding motif for I-A(g7), the class II major histocompatibility complex (MHC) molecule of NOD and Biozzi AB/H mice

The class II major histocompatibility complex molecule I-A(g7) is strongly linked to the development of spontaneous insulin-dependent diabetes mellitus (IDDM) in non obese diabetic mice and to the induction of experimental allergic encephalomyelitis in Biozzi AB/H mice. Structurally, it resembles the HLA-DQ molecules associated with human IDDM, in having a non-Asp residue at position 57 in its beta chain. To identify the requirements for peptide binding to I-A(g7) and thereby potentially pathogenic T cell epitopes, we analyzed a known I-A(g7)-restricted T cell epitope, hen egg white lysozyme (HEL) amino acids 9-27. NH2- and COOH-terminal truncations demonstrated that the minimal epitope for activation of the T cell hybridoma 2D12.1 was M12-R21 and the minimum sequence for direct binding to purified I-A(g7) M12-Y20/K13-R21. Alanine (A) scanning revealed two primary anchors for binding at relative positions (p) 6 (L) and 9 (Y) in the HEL epitope. The critical role of both anchors was demonstrated by incorporating L and Y in poly(A) backbones at the same relative positions as in the HEL epitope. Well-tolerated, weakly tolerated, and nontolerated residues were identified by analyzing the binding of peptides containing multiple substitutions at individual positions. Optimally, p6 was a large, hydrophobic residue (L, I, V, M), whereas p9 was aromatic and hydrophobic (Y or F) or positively charged (K, R). Specific residues were not tolerated at these and some other positions. A motif for binding to I-A(g7) deduced from analysis of the model HEL epitope was present in 27/30 (90%) of peptides reported to be I-A(g7)-restricted T cell epitopes or eluted from I-A(g7). Scanning a set of overlapping peptides encompassing human proinsulin revealed the motif in 6/6 good binders (sensitivity = 100%) and 4/13 weak or non-binders (specificity = 70%). This motif should facilitate identification of autoantigenic epitopes relevant to the pathogenesis and immunotherapy of IDDM.

Proton transport model in the ionosphere .1. Multistream approach of the transport equations

The suprathermal particles, electrons and protons, coming from the magnetosphere and precipitating into the high-latitude atmosphere are an energy source of the Earth's ionosphere. They interact with ambient thermal gas through inelastic and elastic collisions. The physical quantities perturbed by these precipitations, such as the heating rate, the electron production rate, or the emission intensities, can be provided in solving the kinetic stationary Boltzmann equation. This equation yields particle fluxes as a function of altitude, energy, and pitch angle. While this equation has been solved through different ways for the electron transport and fully tested, the proton transport is more complicated. Because of charge-changing reactions, the latter is a set of two-coupled transport equations that must be solved: one for protons and the other for H atoms. We present here a new approach that solves the multistream proton/hydrogen transport equations encompassing the collision angular redistributions and the magnetic mirroring effect. In order to validate our model we discuss the energy conservation and we compare with another model under the same inputs and with rocket observations. The influence of the angular redistributions is discussed in a forthcoming paper.

A dynamic model with on-line identification for rotary sugar drying processes

A dynamic modelling methodology, which combines on-line variable estimation and parameter identification with physical laws to form an adaptive model for rotary sugar drying processes, is developed in this paper. In contrast to the conventional rate-based models using empirical transfer coefficients, the heat and mass transfer rates are estimated by using on-line measurements in the new model. Furthermore, a set of improved sectional solid transport equations with localized parameters is developed in this work to reidentified on-line using measurement data, the model is able to closely track the dynamic behaviour of rotary drying processes within a broad range of operational conditions. This adaptive model is validated against experimental data obtained from a pilot-scale rotary sugar dryer. The proposed modelling methodology can be easily incorporated into nonlinear model based control schemes to form a unified modelling and control framework.place the global correlation for the computation of solid retention time. Since a number of key model variables and parameters are identified on-line using measurement data, the model is able to closely track the dynamic behaviour of rotary drying processes within a broad range of operational conditions. This adaptive model is validated against experimental data obtained from a pilot-scale rotary sugar dryer. The proposed modelling methodology can be easily incorporated into nonlinear model based control schemes to form a unified modelling and control framework.

Set-valued Markov chains and negative semitrajectories of discretized dynamical systems

Computer simulation of dynamical systems involves a phase space which is the finite set of machine arithmetic. Rounding state values of the continuous system to this grid yields a spatially discrete dynamical system, often with different dynamical behaviour. Discretization of an invertible smooth system gives a system with set-valued negative semitrajectories. As the grid is refined, asymptotic behaviour of the semitrajectories follows probabilistic laws which correspond to a set-valued Markov chain, whose transition probabilities can be explicitly calculated. The results are illustrated for two-dimensional dynamical systems obtained by discretization of fractional linear transformations of the unit disc in the complex plane.

Viscoelasticity of radiation-formed PVA/PVP hydrogel

The dynamic theological behaviour of gamma-irradiated 12.8 wt% poly(vinyl alcohol) (PVA), 12.8 wt% poly(vinyl pyrrolidone) (PVP), and a blend of 8 wt% PVA and 4.8 wt% PVP aqueous solutions have been studied pre- and post-gelation. The non-irradiated solutions displayed theological behaviour typical of dilute to semi-dilute polymer solutions, with the complex viscosity being independent of the frequency and shear rate (i.e. Newtonian behaviour) over the range of frequencies tested and the loss modulus G(omega) and storage modulus G(omega) being nearly proportional to omega and omega(2) respectively. After a set of doses of gamma-radiation, the magnitudes of the dynamic moduli G'(omega) and G(omega) increased as the absorbed dose increased, with notable differences between the two homopolymers and the blend. The stages of gelation were effectively monitored by means of dynamic theological measurements, allowing the possible mechanisms of network formation to be elucidated. The doses required for gelation of the PVA, PVP, and blend samples, determined on the basis of the Winter and Chambon criteria for gelation, were found to be 12 kGy for the 12.8 wt% PVA, 4 kGy for the 12.8 wt% PVP, and 5 kGy for the 8 wt% PVA/4.8 wt% PVP solutions. The unexpected lower gelation dose demonstrated by the blend sample, compared with predictions based on the blend composition, and the associated gelation mechanism are also discussed.

Categorization and characterization of transcript-confirmed constitutively and alternatively spliced introns and exons from human

By spliced alignment of human DNA and transcript sequence data we constructed a data set of transcript-confirmed exons and introns from 2793 genes, 796 of which (28%) were seen to have multiple isoforms. We find that over one-third of human exons can translate in more than one frame, and that this is highly correlated with G+C content. Introns containing adenosine at donor site position +3 (A3), rather than guanosine (G3), are more common in low G+C regions, while the converse is true in high G+C regions. These two classes of introns are shown to have distinct lengths, consensus sequences and correlations among splice signals, leading to the hypothesis that A3 donor sites are associated with exon definition, and G3 donor sites with intron definition. Minor classes of introns, including GC-AG, U12-type GT-AG, weak, and putative AG-dependant introns are identified and characterized. Cassette exons are more prevalent in low G+C regions, while exon isoforms are more prevalent in high G+C regions. Cassette exon events outnumber other alternative events, while exon isoform events involve truncation twice as often as extension, and occur at acceptor sites twice as often as at donor sites. Alternative splicing is usually associated with weak splice signals, and in a majority of cases, preserves the coding frame. The reported characteristics of constitutive and alternative splice signals, and the hypotheses offered regarding alternative splicing and genome organization, have important implications for experimental research into RNA processing. The 'AltExtron' data sets are available at http://www.bit.uq.edu.au/altExtron/ and http://www.ebi.ac.uk/similar tothanaraj/altExtron/.

The Limits of Rational Choice: New Institutionalism in the Test Bed of Central Banking Politics in Australia

This paper tests the explanatory capacities of different versions of new institutionalism by examining the Australian case of a general transition in central banking practice and monetary politics: namely, the increased emphasis on low inflation and central bank independence. Standard versions of rational choice institutionalism largely dominate the literature on the politics of central banking, but this approach (here termed RC1) fails to account for Australian empirics. RC1 has a tendency to establish actor preferences exogenously to the analysis; actors' motives are also assumed a priori; actor's preferences are depicted in relatively static, ahistorical terms. And there is the tendency, even a methodological requirement, to assume relatively simple motives and preference sets among actors, in part because of the game theoretic nature of RC1 reasoning. It is possible to build a more accurate rational choice model by re-specifying and essentially updating the context, incentives and choice sets that have driven rational choice in this case. Enter RC2. However, this move subtly introduces methodological shifts and new theoretical challenges. By contrast, historical institutionalism uses an inductive methodology. Compared with deduction, it is arguably better able to deal with complexity and nuance. It also utilises a dynamic, historical approach, and specifies (dynamically) endogenous preference formation by interpretive actors. Historical institutionalism is also able to more easily incorporate a wider set of key explanatory variables and incorporate wider social aggregates. Hence, it is argued that historical institutionalism is the preferred explanatory theory and methodology in this case.

A note on quasi-stationary distributions of birth-death processes and the SIS logistic epidemic

For Markov processes on the positive integers with the origin as an absorbing state, Ferrari, Kesten, Martinez and Picco studied the existence of quasi-stationary and limiting conditional distributions by characterizing quasi-stationary distributions as fixed points of a transformation Phi on the space of probability distributions on {1, 2,.. }. In the case of a birth-death process, the components of Phi(nu) can be written down explicitly for any given distribution nu. Using this explicit representation, we will show that Phi preserves likelihood ratio ordering between distributions. A conjecture of Kryscio and Lefevre concerning the quasi-stationary distribution of the SIS logistic epidemic follows as a corollary.

Intracellular sorting and transport of proteins

The secretory and endocytic pathways of eukaryotic organelles consist of multiple compartments, each with a unique set of proteins and lipids. Specific transport mechanisms are required to direct molecules to defined locations and to ensure that the identity, and hence function, of individual compartments are maintained. The localisation of proteins to specific membranes is complex and involves multiple interactions. The recent dramatic advances in understanding the molecular mechanisms of membrane transport has been due to the application of a multi-disciplinary approach, intergrating membrane biology, genetics, imaging, protein and lipid biochemistry and structural biology. The aim of this review is to summarise the general principles of protein sorting in the secretory and endocytic pathways and to highlight the dynamic nature of these processes. The molecular mechanisms involved in this transport along the secretory and endocytic pathways are discussed along with the signals responsible for targeting proteins to different intracellular locations. (C) 2003 Elsevier Science Ltd. All rights reserved.

Economies of violence: An autobiographical account

This paper uses autobiographical stories to explore various aspects of institutional violence. The stories are taken from a range of institutional settings including schools, academia and commercial enterprises. Each story is used in a dynamic interplay with extant theories and explanations of violence in a mutually explorative and informing manner. Each tale intersects with a different set of issues relating to organisations and violence. The stories and the analysis cohere around the central notion that there exist economies of violence in which violence is exchanged, transacted and within which it circulates. Such economies are constituted and legitimated by discourses, the economy shifts and alters as the surrounding discourses change.

Orientation of the genetic variance-covariance matrix and the fitness surface for multiple male sexually selected traits

Stabilizing selection has been predicted to change genetic variances and covariances so that the orientation of the genetic variance-covariance matrix (G) becomes aligned with the orientation of the fitness surface, but it is less clear how directional selection may change G. Here we develop statistical approaches to the comparison of G with vectors of linear and nonlinear selection. We apply these approaches to a set of male sexually selected cuticular hydrocarbons (CHCs) of Drosophila serrata. Even though male CHCs displayed substantial additive genetic variance, more than 99% of the genetic variance was orientated 74.9degrees away from the vector of linear sexual selection, suggesting that open-ended female preferences may greatly reduce genetic variation in male display traits. Although the orientation of G and the fitness surface were found to differ significantly, the similarity present in eigenstructure was a consequence of traits under weak linear selection and strong nonlinear ( convex) selection. Associating the eigenstructure of G with vectors of linear and nonlinear selection may provide a way of determining what long-term changes in G may be generated by the processes of natural and sexual selection.

Quantum doubles from a class of noncocommutative weak Hopf algebras

The concept of biperfect (noncocommutative) weak Hopf algebras is introduced and their properties are discussed. A new type of quasi-bicrossed products is constructed by means of weak Hopf skew-pairs of the weak Hopf algebras which are generalizations of the Hopf pairs introduced by Takeuchi. As a special case, the quantum double of a finite dimensional biperfect (noncocommutative) weak Hopf algebra is built. Examples of quantum doubles from a Clifford monoid as well as a noncommutative and noncocommutative weak Hopf algebra are given, generalizing quantum doubles from a group and a noncommutative and noncocommutative Hopf algebra, respectively. Moreover, some characterizations of quantum doubles of finite dimensional biperfect weak Hopf algebras are obtained. (C) 2004 American Institute of Physics.

Accuracy of SRM and power tap power monitoring systems for bicycling

Purpose: Although manufacturers of bicycle power monitoring devices SRM and Power Tap (PT) claim accuracy to within 2.5%, there are limited scientific data available in support. The purpose of this investigation was to assess the accuracy of SRM and PT under different conditions. Methods: First, 19 SRM were calibrated, raced for 11 months, and retested using a dynamic CALRIG (50-1000 W at 100 rpm). Second, using the same procedure, five PT were repeat tested on alternate days. Third, the most accurate SRM and PT were tested for the influence of cadence (60, 80, 100, 120 rpm), temperature (8 and 21degreesC) and time (1 h at similar to300 W) on accuracy. Finally, the same SRM and PT were downloaded and compared after random cadence and gear surges using the CALRIG and on a training ride. Results: The mean error scores for SRM and PT factory calibration over a range of 50-1000 W were 2.3 +/- 4.9% and -2.5 +/- 0.5%, respectively. A second set of trials provided stable results for 15 calibrated SRM after 11 months (-0.8 +/- 1.7%), and follow-up testing of all PT units confirmed these findings (-2.7 +/- 0.1%). Accuracy for SRM and PT was not largely influenced by time and cadence; however. power output readings were noticeably influenced by temperature (5.2% for SRM and 8.4% for PT). During field trials, SRM average and max power were 4.8% and 7.3% lower, respectively, compared with PT. Conclusions: When operated according to manufacturers instructions, both SRM and PT offer the coach, athlete, and sport scientist the ability to accurately monitor power output in the lab and the field. Calibration procedures matching performance tests (duration, power, cadence, and temperature) are, however, advised as the error associated with each unit may vary.

Initiation of evasive manoeuvres during self-motion: a test of three hypotheses

To understand performance of evasive and interceptive actions it is important to know how people decide when to initiate a movement - initiating at the 'right' moment is often essential for successful performance. It has been proposed that initiation is triggered when a perceptually derived quantity reaches an invariant criterion value. Candidate quantities include time-to-collision (TTC), distance, and rate of image expansion ( ROE), all of which have received empirical support. We studied initiation of an evasive manoeuvre in a computer-simulated steering task in which the observer was required to steer through a stationary visual environment and avoid colliding with an obstacle in their path. The results could not be explained by hypotheses which propose that evasive manoeuvre initiation is based on a fixed criterion value of TTC or distance. The overall pattern was, however, consistent with the use of a criterion ROE value. This was further tested by analyses designed to directly evaluate whether the ROE value used to initiate the response was the same across experimental conditions. Only two of the six participants showed evidence for using the ROE strategy.