Numerical simulation of abrupt contraction flows using the Double Convected Pom-Pom model

The Double Convected Pom-Pom model was recently introduced to circumvent some numerical and theological defects found in other formulations of the Pom-Pom concept. It is used here for the simulation of a benchmark problem: the flow in an abrupt planar contraction. The predictions are compared with birefringence measurements and show reasonable quantitative agreement with experimental data. A parametric study is also carried out with the aim of analysing the effect of the branching parameter on vortex dynamics and extrudate swell. The results show that the Double Convected Pom-Pom model (DCPP) model is able to discriminate between branched and linear macromolecular structures in accordance with experimental observations. In that respect, the role of the extensional properties in determining complex flow behaviour is stressed. Also, the ratio of the first normal stress difference to the shear stress appears to play a major role in die swell observation. For the time being, the role of the second normal stress difference appears to be less obvious to evaluate in this complex flow. (C) 2004 Elsevier B.V. All rights reserved.

The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Finite mixture regression model with random effects: application to neonatal hospital length of stay

A two-component mixture regression model that allows simultaneously for heterogeneity and dependency among observations is proposed. By specifying random effects explicitly in the linear predictor of the mixture probability and the mixture components, parameter estimation is achieved by maximising the corresponding best linear unbiased prediction type log-likelihood. Approximate residual maximum likelihood estimates are obtained via an EM algorithm in the manner of generalised linear mixed model (GLMM). The method can be extended to a g-component mixture regression model with the component density from the exponential family, leading to the development of the class of finite mixture GLMM. For illustration, the method is applied to analyse neonatal length of stay (LOS). It is shown that identification of pertinent factors that influence hospital LOS can provide important information for health care planning and resource allocation. (C) 2002 Elsevier Science B.V. All rights reserved.

Differences in ethical attitudes between registered nurses and medical students

In this study we compared the ethical attitudes of a group of experienced, predominantly female, registered nurses (n = 67) with those of a group of final year, mixed sex, medical students (n = 125). The purpose was to determine the basis of differences in attitudes that could lead to ethical disagreements between these two groups when they came to work together. A questionnaire developed to explore ethical attitudes was administered and the responses of the two groups were compared using t-tests. Because of the preponderance of females among the nurses an analysis of variance of the gender-adjusted scores for each group was also carried out. On comparing the responses, the nurses differed significantly from the medical students in a number of ethical domains. A potential source of conflict between these two groups is that the nurses were inclined to adopt the perspective of patients but the medical students identified with their profession. When corrected for the effects of gender, the differences persisted, indicating that it was discipline that determined the differences. We recommend that students of nursing and medicine receive ethics education together, and that more open dialogue between doctors and nurses with respect to their different ethical viewpoints is needed in the work setting. This article will be of interest to educators of students of medicine and nursing, as well as to doctors and nurses who are eager to improve their professional relations and thereby improve patient care.

Oxidation of mixed-valence CoIII/FeII complexes reversed at high pH: A kinetico-mechanistic study of water oxidation

The outer-sphere oxidation of Fell in the mixed-valence complex trans-[(LCoNCFeII)-Co-14S-N-III(CN)(6)](-), being L-14S an N3S2 macrocylic donor set on the cobalt(III) center, has been studied. The comparison with the known processes of N-5 macrocycle complexes has been carried out in view of the important differences occurring on the redox potential of the cobalt center. The results indicate that the outer-sphere oxidation reactions with S2O82- and [Co(ox)(3)](3-) involve a great amount of solvent-assisted hydrogen bonding that, as a consequence from the change from two amines to sulfur donors, are more restricted. This is shown by the more positive values found for DeltaS(double dagger) and DeltaV(double dagger). The X-ray structure of the oxidized complex has been determined, and it is clearly indicative of the above-mentioned solvent-assisted hydrogen bonding between nitrogen and cyanide donors on the cobalt and iron centers, respectively. trans-[(LCoNCFeIII)-Co-14S-N-III(CN)(6)], as well as the analogous N-5 systems trans-[(LCoNCFeIII)-Co-14-N-III(CN)(6)], trans-[(LCoNCFeIII)-Co-15-N-III-(CN)(6)], and cis-[(LCoNCFeIII)-Co-n-N-III(CN)(6)], Oxidize water to hydrogen peroxide at pH > 10 with a rather simple stoichiometry, i.e., [(LCoNCFeIII)-Co-n-N-III(CN)(5)] + OH- - [(LCoNCFeII)-Co-n-N-III(CN)(5)](-) + 1/2H(2)O(2). In this way, the reversibility of the iron oxidation process is achieved. The determination of kinetic and thermal and pressure activation parameters for this water to hydrogen peroxide oxidation leads to the kinetic determination of a cyanide based OH- adduct of the complex. A second-order dependence on the base concentration is associated with deprotonation of this adduct to produce the final inner-sphere reduction process. The activation enthalpies are found to be extremely low (15 to 35 kJ mol(-1)) and responsible for the very fast reaction observed. The values of DeltaS(double dagger) and DeltaV(double dagger) (-76 to -113 J K-1 mol(-1) and -5.5 to -8.9 cm(3) mol(-1), respectively) indicate a highly organized but not very compressed transition state in agreement with the inner-sphere one-electron transfer from O2- to Fe-III.

Gelatinization of starch in mixed sugar systems

Sugars affect the gelatinization of starch, with the effect varying significantly between sugars. Since many food products contain a mixture of sugar sources, it is important to understand how their mixtures affect starch gelatinization. In a Rapid Visco Analyser study of maize starch gelatinization, changing proportions in binary mixtures of refined sugars saw a largely proportionate change in starch gelatinization properties. However, binary mixture of pure sugars and honey, or a model honey system (the main sugars in honey) and honey responded differently. Generally, replacing 25% or 50% of the refined sugar or model honey system with honey gave a large change in starch gelatinization properties, while further increases in honey level had little further effect. Differences between honey and buffered model honey system (either gluconic acid, or a mixture of citric acid and di-sodium phosphate) showed the sensitivity of starch gelatinization to the composition of the nonsaccharide component. (c) 2004 Swiss Society of Food Science and Technology. Published by Elsevier Ltd. All rights reserved.

Comparison of beam theory and finite-element analysis with in vivo bone strain data from the alligator cranium

The mechanical behavior of the vertebrate skull is often modeled using free-body analysis of simple geometric structures and, more recently, finite-element (FE) analysis. In this study, we compare experimentally collected in vivo bone strain orientations and magnitudes from the cranium of the American alligator with those extrapolated from a beam model and extracted from an FE model. The strain magnitudes predicted from beam and FE skull models bear little similarity to relative and absolute strain magnitudes recorded during in vivo biting experiments. However, quantitative differences between principal strain orientations extracted from the FE skull model and recorded during the in vivo experiments were smaller, and both generally matched expectations from the beam model. The differences in strain magnitude between the data sets may be attributable to the level of resolution of the models, the material properties used in the FE model, and the loading conditions (i.e., external forces and constraints). This study indicates that FE models and modeling of skulls as simple engineering structures may give a preliminary idea of how these structures are loaded, but whenever possible, modeling results should be verified with either in vitro or preferably in vivo testing, especially if precise knowledge of strain magnitudes is desired. (c) 2005 Wiley-Liss, Inc.

Anisotropy-based robust performance analysis of finite horizon linear discrete time varying systems

We consider a problem of robust performance analysis of linear discrete time varying systems on a bounded time interval. The system is represented in the state-space form. It is driven by a random input disturbance with imprecisely known probability distribution; this distributional uncertainty is described in terms of entropy. The worst-case performance of the system is quantified by its a-anisotropic norm. Computing the anisotropic norm is reduced to solving a set of difference Riccati and Lyapunov equations and a special form equation.