A finite difference scheme for inviscid flows with non-eequilibrium chemistry and internal energy

A finite difference scheme is presented for the inviscid terms of the equations of compressible fluid dynamics with general non-equilibrium chemistry and internal energy.

An efficient Riemann solver for unsteady flows with non-ideal gases

A Riemann solver is presented for the Euler equations of gas dynamics with real gases. This represents a more efficient version of an algorithm originally presented by the author.

Effects of plant pathogens on population dynamics and community composition in grassland ecosystems: two case studies

Grassland ecosystems comprise a major portion of the earth’s terrestrial surface, ranging from high-input cultivated monocultures or simple species mixtures to relatively unmanaged but dynamic systems. Plant pathogens are a component of these systems with their impact dependent on many interacting factors, including grassland species population dynamics and community composition, the topics covered in this paper. Plant pathogens are affected by these interactions and also act reciprocally by modifying their nature. We review these features of disease in grasslands and then introduce the 150-year long-term Park Grass Experiment (PGE) at Rothamsted Research in the UK. We then consider in detail two plant-pathogen systems present in the PGE, Tragopogon pratensis-Puccinia hysterium and Holcus lanata-Puccinia coronata. These two systems have very different life history characteristics: the first, a biennial member of the Asteraceae infected by its host-specific, systemic rust; the second, a perennial grass infected by a host-non-specific rust. We illustrate how observational, experimental and modelling studies can contribute to a better understanding of population dynamics, competitive interactions and evolutionary outcomes. With Tragopogon pratensis-Puccinia hysterium, characterised as an “outbreak” species in the PGE, we show that pathogen-induced mortality is unlikely to be involved in host population regulation; and that the presence of even a short-lived seed-bank can affect the qualitative outcomes of the host-pathogen dynamics. With Holcus lanata-Puccinia coronata, we show how nutrient conditions can affect adaptation in terms of host defence mechanisms, and that co-existence of competing species affected by a common generalist pathogen is unlikely.

Development and validation of a parent report measure for detection of cognitive delay in infancy

Aim To develop a brief, parent-completed instrument (‘ERIC’) for detection of cognitive delay in 10-24 month-olds born preterm, or with low birth weight, or with perinatal complications, and to establish its diagnostic properties. Method Scores were collected from parents of 317 children meeting ≥1 inclusion criteria (birth weight <1500g; gestational age <34 completed weeks; 5-minute Apgar <7; presence of hypoxic-ischemic encephalopathy) and meeting no exclusion criteria. Children were assessed for cognitive delay using a criterion score on the Bayley Scales of Infant and Toddler Development Cognitive Scale III1 <80. Items were retained according to their individual associations with delay. Sensitivity, specificity, Positive and Negative Predictive Values were estimated and a truncated ERIC was developed for use <14 months. Results ERIC detected 17 out of 18 delayed children in the sample, with 94.4% sensitivity (95% CI [confidence interval] 83.9-100%), 76.9% specificity (72.1-81.7%), 19.8% positive predictive value (11.4-28.2%); 99.6% negative predictive value (98.7-100%); 4.09 likelihood ratio positive; and 0.07 likelihood ratio negative; the associated Area under the Curve was .909 (.829-.960). Interpretation ERIC has potential value as a quickly-administered diagnostic instrument for the absence of early cognitive delay in preterm or premature infants of 10-24 months, and as a screen for cognitive delay. Further research may be needed before ERIC can be recommended for wide-scale use.

Rerum cognoscere causas: part I - How do the ideas of system dynamics relate to traditional social theories and the voluntarism/determinism debate?

This is the first half of a two-part paper which deals with the social theoretic assumptions underlying system dynamics. The motivation is that clarification in this area can help mainstream social scientists to understand how our field relates to their literature, methods and concerns. Part I has two main sections. The aim of the first is to answer the question: How do the ideas of system dynamics relate to traditional social theories? The theoretic assumptions of the field are seldom explicit but rather are implicit in its practice. The range of system dynamics practice is therefore considered and related to a framework - widely used in both operational research (OR) and systems science - that organises the assumptions behind traditional social theoretic paradigms. Distinct and surprisingly varied groupings of practice are identified, making it difficult to place system dynamics in any one paradigm with any certainty. The difficulties of establishing a social theoretic home for system dynamics are exemplified in the second main section. This is done by considering the question: Is system dynamics deterministic? An analysis shows that attempts to relate system dynamics to strict notions of voluntarism or determinism quickly indicate that the field does not fit with either pole of this dichotomous, and strictly paradigmatic, view. Part I therefore concludes that definitively placing system dynamics with respect to traditional social theories is highly problematic. The scene is therefore set for Part II of the paper, which proposes an innovative and potentially fruitful resolution to this problem.

The non-conservation of potential vorticity by a dynamical core compared with the effects of parametrized physical processes

Numerical models of the atmosphere combine a dynamical core, which approximates solutions to the adiabatic, frictionless governing equations for fluid dynamics, with tendencies arising from the parametrization of other physical processes. Since potential vorticity (PV) is conserved following fluid flow in adiabatic, frictionless circumstances, it is possible to isolate the effects of non-conservative processes by accumulating PV changes in an air-mass relative framework. This “PV tracer technique” is used to accumulate separately the effects on PV of each of the different non-conservative processes represented in a numerical model of the atmosphere. Dynamical cores are not exactly conservative because they introduce, explicitly or implicitly, some level of dissipation and adjustment of prognostic model variables which acts to modify PV. Here, the PV tracers technique is extended to diagnose the cumulative effect of the non-conservation of PV by a dynamical core and its characteristics relative to the PV modification by parametrized physical processes. Quantification using the Met Office Unified Model reveals that the magnitude of the non-conservation of PV by the dynamical core is comparable to those from physical processes. Moreover, the residual of the PV budget, when tracing the effects of the dynamical core and physical processes, is at least an order of magnitude smaller than the PV tracers associated with the most active physical processes. The implication of this work is that the non-conservation of PV by a dynamical core can be assessed in case studies with a full suite of physics parametrizations and directly compared with the PV modification by parametrized physical processes. The nonconservation of PV by the dynamical core is shown to move the position of the extratropical tropopause while the parametrized physical processes have a lesser effect at the tropopause level.

Rossby wave propagation on potential vorticity fronts with finite width

The horizontal gradient of potential vorticity (PV) across the tropopause typically declines with lead time in global numerical weather forecasts and tends towards a steady value dependent on model resolution. This paper examines how spreading the tropopause PV contrast over a broader frontal zone affects the propagation of Rossby waves. The approach taken is to analyse Rossby waves on a PV front of finite width in a simple single-layer model. The dispersion relation for linear Rossby waves on a PV front of infinitesimal width is well known; here an approximate correction is derived for the case of a finite width front, valid in the limit that the front is narrow compared to the zonal wavelength. Broadening the front causes a decrease in both the jet speed and the ability of waves to propagate upstream. The contribution of these changes to Rossby wave phase speeds cancel at leading order. At second order the decrease in jet speed dominates, meaning phase speeds are slower on broader PV fronts. This asymptotic phase speed result is shown to hold for a wide class of single-layer dynamics with a varying range of PV inversion operators. The phase speed dependence on frontal width is verified by numerical simulations and also shown to be robust at finite wave amplitude, and estimates are made for the error in Rossby wave propagation speeds due to the PV gradient error present in numerical weather forecast models.

Monotone positive solutions for a fourth order equation with nonlinear boundary conditions

This work is concerned with the existence of monotone positive solutions for a class of beam equations with nonlinear boundary conditions. The results are obtained by using the monotone iteration method and they extend early works on beams with null boundary conditions. Numerical simulations are also presented. (C) 2009 Elsevier Ltd. All rights reserved.

Infrared analysis of ion beam irradiated polymers

Irradiation with heavy ions can produce several modifications in the chain structure of polymers. These modifications can be related to scissioning and cross-linking of chemical bonds. which depend on the ion fluence and the density of energy deposited in the material. Stacked thin film Makrofol-KG (R) samples were irradiated with 350 MeV Au(26+) ions and FTIR absorption spectroscopy was used to determine the bond changes in the samples. Data on the absorption bands as a function of the fluence indicated a higher probability for simple-bonds scissioning than for double-bonds scissioning and no dependence on the number of double bonds breaking with ion fluence. Since sample irradiation was done in a non-track-overlapping regime, a novel process for double bonds formation is suggested: the excitation of a site in the material by only one incident ion followed by a double bond formation during the de-excitation process. (C) 2009 Elsevier B.V. All rights reserved.

Near threshold vibrational excitation of molecules by positron impact: A projection operator approach

We report vibrational excitation (v(i) = 0 -> v(f) = 1) cross-sections for positron scattering by H(2) and model calculations for the (v(i) = 0 -> v(f) = 1) excitation of the C-C symmetric stretch mode of C(2)H(2). The Feshbach projection operator formalism was employed to vibrationally resolve the fixed-nuclei phase shifts obtained with the Schwinger multichannel method. The near threshold behavior of H(2) and C(2)H(2) significantly differ in the sense that no low lying singularity (either virtual or bound state) was found for the former, while a e(+)-acetylene virtual state was found at the equilibrium geometry (this virtual state becomes a bound state upon stretching the molecule). For C(2)H(2), we also performed model calculations comparing excitation cross-sections arising from virtual (-i kappa(0)) and bound (+i kappa(0)) states symmetrically located around the origin of the complex momentum plane (i.e. having the same kappa(0)). The virtual state is seen to significantly couple to vibrations, and similar cross-sections were obtained for shallow bound and virtual states. (c) 2007 Elsevier B.V. All rights reserved.

Cascades of Hopf bifurcations from boundary delay

We consider a 1-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We deal with non-negative solutions and analyze the stability behavior of its unique positive equilibrium solution, which is given by the constant function u equivalent to 1. We show that if the delay is small, this equilibrium solution is asymptotically stable, similar as in the case without delay. We also show that, as the delay goes to infinity, this equilibrium becomes unstable and undergoes a cascade of Hopf bifurcations. The structure of this cascade will depend on the parameters appearing in the equation. This equation shows some dynamical behavior that differs from the case where the nonlinearity with delay is in the interior of the domain. (C) 2009 Elsevier Inc. All rights reserved.

ON THE SUPERCRITICALITY OF THE FIRST HOPF BIFURCATION IN A DELAY BOUNDARY PROBLEM

The goal of this paper is to analyze the character of the first Hopf bifurcation (subcritical versus supercritical) that appears in a one-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We showed in the previous work [Arrieta et al., 2010] that if the delay is small, the unique non-negative equilibrium solution is asymptotically stable. We also showed that, as the delay increases and crosses certain critical value, this equilibrium becomes unstable and undergoes a Hopf bifurcation. This bifurcation is the first one of a cascade occurring as the delay goes to infinity. The structure of this cascade will depend on the parameters appearing in the equation. In this paper, we show that the first bifurcation that occurs is supercritical, that is, when the parameter is bigger than the delay bifurcation value, stable periodic orbits branch off from the constant equilibrium.

Tayloring Brazil: a system dynamics model for monetary policy feedback

The thesis introduces a system dynamics Taylor rule model of new Keynesian nature for monetary policy feedback in Brazil. The nonlinear Taylor rule for interest rate changes con-siders gaps and dynamics of GDP growth and inflation. The model closely tracks the 2004 to 2011 business cycle and outlines the endogenous feedback between the real interest rate, GDP growth and inflation. The model identifies a high degree of endogenous feedback for monetary policy and inflation, while GDP growth remains highly exposed to exogenous eco-nomic conditions. The results also show that the majority of the monetary policy moves during the sample period was related to GDP growth, despite higher coefficients of inflation parameters in the Taylor rule. This observation challenges the intuition that inflation target-ing leads to a dominance of monetary policy moves with respect to inflation. Furthermore, the results suggest that backward looking price-setting with respect to GDP growth has been the dominant driver of inflation. Moreover, simulation exercises highlight the effects of the new BCB strategy initiated in August 2011 and also consider recession and inflation avoid-ance versions of the Taylor rule. In methodological terms, the Taylor rule model highlights the advantages of system dynamics with respect to nonlinear policies and to the stock-and-flow approach. In total, the strong historical fit and some counterintuitive observations of the Taylor rule model call for an application of the model to other economies.

Dynamic Tracking with Zero Variation and Disturbance Rejection Applied to Discrete-Time Systems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Pair term in the electromagnetic current within the Front-Form dynamics: spin-0 case

The frame and scale dependence of the pair-term contribution to the electromagnetic form factor of a spin-zero composite system of two-fermions is studied within the Light Front. The form factor is evaluated from the plus-component of the current in the Breit frame, using for the first time a nonconstant, symmetric ansatz for the Bethe-Salpeter amplitude. The frame dependence is analyzed by allowing a nonvanishing plus component of the momentum transfer, while the dynamical scale is set by the masses of the constituents and by mass and size of the composite system. A transverse momentum distribution, associated with the Bethe-Salpeter amplitude, is introduced which allows to define strongly and weakly relativistic systems. In particular, for strongly relativistic systems, the pair term vanishes for the Drell-Yan condition, while is dominant for momentum transfer along the light-front direction. For a weakly relativistic system, fitted to the deuteron scale, the pair term is negligible up to momentum transfers of 1 (GeV/c)(2). A comparison with results obtained within the Front-Form Hamiltonian dynamics with a fixed number of constituents is also presented. (C) 2002 Elsevier B.V. B.V. All rights reserved.