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.

Numerical simulation of stochastic ordinary differential equations in biomathematical modelling

In this work we discuss the effects of white and coloured noise perturbations on the parameters of a mathematical model of bacteriophage infection introduced by Beretta and Kuang in [Math. Biosc. 149 (1998) 57]. We numerically simulate the strong solutions of the resulting systems of stochastic ordinary differential equations (SDEs), with respect to the global error, by means of numerical methods of both Euler-Taylor expansion and stochastic Runge-Kutta type. (C) 2003 IMACS. Published by Elsevier B.V. All rights reserved.

Numerical evaluation of three dimensional effects in planar flow birefringence

The influence of three dimensional effects on isochromatic birefringence is evaluated for planar flows by means of numerical simulation. Two fluid models are investigated in channel and abrupt contraction geometries. In practice, the flows are confined by viewing windows, which alter the stresses along the optical path. The observed optical properties differ therefore from their counterpart in an ideal two-dimensional flow. To investigate the influence of these effects, the stress optical rule and the differential propagation Mueller matrix are used. The material parameters are selected so that a retardation of multiple orders is achieved, as is typical for highly birefringent melts. Errors due to three dimensional effects are mainly found on the symmetry plane, and increase significantly with the flow rate. Increasing the geometric aspect ratio improve the accuracy provided that the error on the retardation is less than one order. (C) 2004 Elsevier B.V. All rights reserved.

Dynamic and steady-state rheology of Australian honeys at subzero temperatures

The rheology of 10 Australian honeys was investigated at temperatures -15C to 0C by a strain-controlled rheometer. The honeys exhibited Newtonian behavior irrespective of the temperature, and follow the Cox-Merz rule. G/G' and omega are quadratically related, and the crossover frequencies for liquid to solid transformation and relaxation times were obtained. The composition of the honeys correlates well (r(2) > 0.83) with the viscosity, and with 24 7 data sets (Australian and Greek honeys), the following equation was obtained: mu = 1.41 x 10(-17) exp [-1.20M + 0.01F - 0.0G + (18.6 X 10(3)/T)] The viscosity of the honeys showed a strong dependence on temperature, and four models were examined to describe this. The models gave good fits (r(2) > 0.95), but better fits were obtained for the WLF model using T-g of the honeys and mu(g) = 10(11) Pa.s. The WLF model with its universal values poorly predicted the viscosity, and the implications of the measured rheological behaviors of the honeys in their processing and handling are discussed.

A genetic estimation algorithm for parameters of stochastic ordinary differential equations

A generic method for the estimation of parameters for Stochastic Ordinary Differential Equations (SODEs) is introduced and developed. This algorithm, called the GePERs method, utilises a genetic optimisation algorithm to minimise a stochastic objective function based on the Kolmogorov-Smirnov statistic. Numerical simulations are utilised to form the KS statistic. Further, the examination of some of the factors that improve the precision of the estimates is conducted. This method is used to estimate parameters of diffusion equations and jump-diffusion equations. It is also applied to the problem of model selection for the Queensland electricity market. (C) 2003 Elsevier B.V. All rights reserved.

Investigation of the starch gelatinisation phenomena in water-glycerol systems: application of modulated temperature differential scanning calorimetry

The use of modulated temperature differential scanning calorimetry (MTDSC) has provided further insight into the gelatinisation process since it allows the detection of glass transition during gelatinisation process. It was found in this work that the glass transition overlapped with the gelatinisation peak temperature for all maize starch formulations studied. Systematic investigation on maize starch gelatinisation over a range of water-glycerol concentrations with MTDSC revealed that the addition of glycerol increased the gelatinisation onset temperature with an extent that depended on the water content in the system. Furthermore, the addition of glycerol promoted starch gelatinisation at low water content (0.4 g water/g dry starch) and the enthalpy of gelatinisation varied with glycerol concentration (0.73-19.61 J/g dry starch) depending on the water content and starch type. The validities of published gelatinisation models were explored. These models failed to explain the glass transition phenomena observed during the course of gelatinisation and failed to describe the gelatinisation behaviour observed over the water-glycerol concentrations range investigated. A hypothesis for the mechanisms involved during gelatinisation was proposed based on the side chain liquid crystalline polymer model for starch structure and the concept that the order-disorder transition in starch requires that the hydrogen bonds (the major structural element in the granule packing) to be broken before the collapse of order (helix-coil transition) can take place. (C) 2004 Elsevier Ltd. All rights reserved.

A cognitive model of binge drinking: The influence of alcohol expectancies and drinking refusal self-efficacy

While binge drinking-episodic or irregular consumption of excessive amounts of alcohol-is recognised as a serious problem affecting our youth, to date there has been a lack of psychological theory and thus theoretically driven research into this problem. The current paper develops a cognitive model using the key constructs of alcohol expectancies (AEs) and drinking refusal self-efficacy (DRSE) to explain the acquisition and maintenance of binge drinking. It is suggested that the four combinations of the AE and DRSE can explain the four drinking styles. These are normal/social drinkers, binge drinkers, regular heavy drinkers, and problem drinkers or alcoholics. Since AE and DRSE are cognitive constructs and therefore modifiable, the cognitive model can thus facilitate the design of intervention and-prevention strategies for binge drinking. (C) 2003 Elsevier Ltd. All rights reserved.

H-ras, K-ras, and inner plasma membrane raft proteins operate in nanoclusters with differential dependence on the actin cytoskeleton

Plasma membrane compartmentalization imposes lateral segregation on membrane proteins that is important for regulating signal transduction. We use computational modeling of immunogold spatial point patterns on intact plasma membrane sheets to test different models of inner plasma membrane organization. We find compartmentalization at the nanoscale level but show that a classical raft model of preexisting stable domains into which lipid raft proteins partition is incompatible with the spatial point patterns generated by the immunogold labeling of a palmitoylated raft marker protein. Rather, approximate to 30% of the raft protein exists in cholesterol-dependent nanoclusters, with approximate to 70% distributed as monomers. The cluster/monomer ratio (number of proteins in clusters/number of proteins outside clusters) is independent of expression level. H-rasG12V and K-rasG12V proteins also operate in nanoclusters with fixed cluster/monomer ratios that are independent of expression level. Detailed calibration of the immunogold imaging protocol suggests that radii of raft and RasG12V protein nanoclusters may be as small as 11 and 6 nm, respectively, and shows that the nanoclusters contain small numbers (6.0-7.7) of proteins. Raft nanoclusters do not form if the actin cytoskeleton is disassembled. The formation of K-rasG12V but not H-rasG12V nanoclusters also is actin-dependent. K-rasG12V but not H-rasG12V signaling is abrogated by actin cytoskeleton disassembly, which shows that nanoclustering is critical for Ras function. These findings argue against stable preexisting domains on the inner plasma membrane in favor of dynamic actively regulated nanoclusters similar to those proposed for the outer plasma membrane. RasG12V nanoclusters may facilitate the assembly of essential signal transduction complexes.

An experimental and finite element poroelastic creep response analysis of an intervertebral hydrogel disc model in axial compression

A hydrogel intervertebral disc (lVD) model consisting of an inner nucleus core and an outer anulus ring was manufactured from 30 and 35% by weight Poly(vinyl alcohol) hydrogel (PVA-H) concentrations and subjected to axial compression in between saturated porous endplates at 200 N for 11 h, 30 min. Repeat experiments (n = 4) on different samples (N = 2) show good reproducibility of fluid loss and axial deformation. An axisymmetric nonlinear poroelastic finite element model with variable permeability was developed using commercial finite element software to compare axial deformation and predicted fluid loss with experimental data. The FE predictions indicate differential fluid loss similar to that of biological IVDs, with the nucleus losing more water than the anulus, and there is overall good agreement between experimental and finite element predicted fluid loss. The stress distribution pattern indicates important similarities with the biological lVD that includes stress transference from the nucleus to the anulus upon sustained loading and renders it suitable as a model that can be used in future studies to better understand the role of fluid and stress in biological IVDs. (C) 2005 Springer Science + Business Media, Inc.

A simple implementation of a normal mixture approach to differential gene expression in multiclass microarrays

Motivation: An important problem in microarray experiments is the detection of genes that are differentially expressed in a given number of classes. We provide a straightforward and easily implemented method for estimating the posterior probability that an individual gene is null. The problem can be expressed in a two-component mixture framework, using an empirical Bayes approach. Current methods of implementing this approach either have some limitations due to the minimal assumptions made or with more specific assumptions are computationally intensive. Results: By converting to a z-score the value of the test statistic used to test the significance of each gene, we propose a simple two-component normal mixture that models adequately the distribution of this score. The usefulness of our approach is demonstrated on three real datasets.

Differential expression of the NMDA NR2B receptor subunit in motoneuron populations susceptible and resistant to amyotrophic lateral sclerosis

We have compared the expression pattern of NMDA receptor subunits (NR1 and NR2A-D)and NRI splice variants (NR1-1a/1b,-2a/2b,-3a/3b,4a/4b) in motor neuron populations from adult Wistar rats that are vulnerable (hypoglossal, XII) or resistant (oculomotor, III) to death in amyotrophic lateral sclerosis (ALS). The major finding was higher levels of expression of the NR2B subunit in the hypoglossal nucleus. Quantitative real-time PCR showed that NR1 was expressed at a greater level than any of the NR2 subunits (> 15 fold greater, P