Performance assessment of exponential Rosenbrock method for large systems of ODE

This paper studies time integration methods for large stiff systems of ordinary differential equations (ODEs) of the form u'(t) = g(u(t)). For such problems, implicit methods generally outperform explicit methods, since the time step is usually less restricted by stability constraints. Recently, however, explicit so-called exponential integrators have become popular for stiff problems due to their favourable stability properties. These methods use matrix-vector products involving exponential-like functions of the Jacobian matrix, which can be approximated using Krylov subspace methods that require only matrix-vector products with the Jacobian. In this paper, we implement exponential integrators of second, third and fourth order and demonstrate that they are competitive with well-established approaches based on the backward differentiation formulas and a preconditioned Newton-Krylov solution strategy.

Correction of step artefact associated with MRI scanning of long bones

3D models of long bones are being utilised for a number of fields including orthopaedic implant design. Accurate reconstruction of 3D models is of utmost importance to design accurate implants to allow achieving a good alignment between two bone fragments. Thus for this purpose, CT scanners are employed to acquire accurate bone data exposing an individual to a high amount of ionising radiation. Magnetic resonance imaging (MRI) has been shown to be a potential alternative to computed tomography (CT) for scanning of volunteers for 3D reconstruction of long bones, essentially avoiding the high radiation dose from CT. In MRI imaging of long bones, the artefacts due to random movements of the skeletal system create challenges for researchers as they generate inaccuracies in the 3D models generated by using data sets containing such artefacts. One of the defects that have been observed during an initial study is the lateral shift artefact occurring in the reconstructed 3D models. This artefact is believed to result from volunteers moving the leg during two successive scanning stages (the lower limb has to be scanned in at least five stages due to the limited scanning length of the scanner). As this artefact creates inaccuracies in the implants designed using these models, it needs to be corrected before the application of 3D models to implant design. Therefore, this study aimed to correct the lateral shift artefact using 3D modelling techniques. The femora of five ovine hind limbs were scanned with a 3T MRI scanner using a 3D vibe based protocol. The scanning was conducted in two halves, while maintaining a good overlap between them. A lateral shift was generated by moving the limb several millimetres between two scanning stages. The 3D models were reconstructed using a multi threshold segmentation method. The correction of the artefact was achieved by aligning the two halves using the robust iterative closest point (ICP) algorithm, with the help of the overlapping region between the two. The models with the corrected artefact were compared with the reference model generated by CT scanning of the same sample. The results indicate that the correction of the artefact was achieved with an average deviation of 0.32 ± 0.02 mm between the corrected model and the reference model. In comparison, the model obtained from a single MRI scan generated an average error of 0.25 ± 0.02 mm when compared with the reference model. An average deviation of 0.34 ± 0.04 mm was seen when the models generated after the table was moved were compared to the reference models; thus, the movement of the table is also a contributing factor to the motion artefacts.

A deletion and point mutation study of the human papillomavirus type 16 major capsid gene

Recombinant human papillomavirus (HPV) virus-like particles (VLPs) made from the major capsid protein L1 are promising vaccine candidates for use as vaccines against genital and other HPV infections, and particularly against HPV-16. However, HPV-16 genotype variants have different binding affinities for neutralising mouse Mabs raised against HPV-16 L1 VLPs. This paper analyses, using a panel of well-characterised Mabs, the effects on the antigenicity of various C- and N-terminal deletants of HPV-16 L1 made in insect cells via recombinant baculovirus, of an A → T mutation at residue 266 (A266T), and of a C → G mutation at conserved position 428 (C428G). The effects of these changes on assembly of the variant L1s were studied by electron microscopy. Binding of Mab H16:E70 to A266T was reduced by almost half in comparison to wild type L1. Retention of the C-terminal region 428-483 was critical for the binding of conformation-specific Mabs (H16:V5, H16:E70, H16:U4 and H16:9A) whereas deletion of the nuclear localisation signal (NLS) or the C428G mutation or an N-terminal deletion (residues 2-9) did not affect the antigenicity. The N-terminal deletion resulted in a mixed population of 30 and 55 nm VLPs, which differs from the same construct expressed in Escherichia coli, whereas pentamer aggregates resulted from deletion of the 428-465 region or the C428G mutation. The results have implications both for considering use of single-genotype HPV vaccines, and for design of novel second-generation vaccines. © 2006 Elsevier B.V. All rights reserved.

Duration-dependant response of mixed-method pre-cooling for intermittent-sprint exercise in the heat

This study examined the effects of pre-cooling duration on performance and neuromuscular function for self-paced intermittent-sprint shuttle running in the heat. Eight male, team-sport athletes completed two 35-min bouts of intermittent-sprint shuttle running separated by a 15-min recovery on three separate occasions (33°C, 34% relative humidity). Mixed-method pre-cooling was completed for 20 min (COOL20), 10-min (COOL10) or no cooling (CONT) and reapplied for 5-min mid-exercise. Performance was assessed via sprint times, percentage decline and shuttle-running distance covered. Maximal voluntary contractions (MVC), voluntary activation (VA) and evoked twitch properties were recorded pre- and post-intervention and mid- and post-exercise. Core temperature (T c), skin temperature, heart rate, capillary blood metabolites, sweat losses, perceptual exertion and thermal stress were monitored throughout. Venous blood draws pre- and post-exercise were analyzed for muscle damage and inflammation markers. Shuttle-running distances covered were increased 5.2 ± 3.3% following COOL20 (P < 0.05), with no differences observed between COOL10 and CONT (P > 0.05). COOL20 aided in the maintenance of mid- and post-exercise MVC (P < 0.05; d > 0.80), despite no conditional differences in VA (P > 0.05). Pre-exercise T c was reduced by 0.15 ± 0.13°C with COOL20 (P < 0.05; d > 1.10), and remained lower throughout both COOL20 and COOL10 compared to CONT (P < 0.05; d > 0.80). Pre-cooling reduced sweat losses by 0.4 ± 0.3 kg (P < 0.02; d > 1.15), with COOL20 0.2 ± 0.4 kg less than COOL10 (P = 0.19; d = 1.01). Increased pre-cooling duration lowered physiological demands during exercise heat stress and facilitated the maintenance of self-paced intermittent-sprint performance in the heat. Importantly, the dose-response interaction of pre-cooling and sustained neuromuscular responses may explain the improved exercise performance in hot conditions.

Mixed-method pre-cooling reduces physiological demand without improving performance of medium-fast bowling in the heat

This study examined physiological and performance effects of pre-cooling on medium-fast bowling in the heat. Ten, medium-fast bowlers completed two randomised trials involving either cooling (mixed-methods) or control (no cooling) interventions before a 6-over bowling spell in 31.9±2.1°C and 63.5±9.3% relative humidity. Measures included bowling performance (ball speed, accuracy and run-up speeds), physical characteristics (global positioning system monitoring and counter-movement jump height), physiological (heart rate, core temperature, skin temperature and sweat loss), biochemical (serum concentrations of damage, stress and inflammation) and perceptual variables (perceived exertion and thermal sensation). Mean ball speed (114.5±7.1 vs. 114.1±7.2 km · h−1; P = 0.63; d = 0.09), accuracy (43.1±10.6 vs. 44.2±12.5 AU; P = 0.76; d = 0.14) and total run-up speed (19.1±4.1 vs. 19.3±3.8 km · h−1; P = 0.66; d = 0.06) did not differ between pre-cooling and control respectively; however 20-m sprint speed between overs was 5.9±7.3% greater at Over 4 after pre-cooling (P = 0.03; d = 0.75). Pre-cooling reduced skin temperature after the intervention period (P = 0.006; d = 2.28), core temperature and pre-over heart rates throughout (P = 0.01−0.04; d = 0.96−1.74) and sweat loss by 0.4±0.3 kg (P = 0.01; d = 0.34). Mean rating of perceived exertion and thermal sensation were lower during pre-cooling trials (P = 0.004−0.03; d = 0.77−3.13). Despite no observed improvement in bowling performance, pre-cooling maintained between-over sprint speeds and blunted physiological and perceptual demands to ease the thermoregulatory demands of medium-fast bowling in hot conditions.

Effects of mixed-method cooling on recovery of medium-fast bowling performance in hot conditions on consecutive days

This investigation examined physiological and performance effects of cooling on recovery of medium-fast bowlers in the heat. Eight, medium-fast bowlers completed two randomised trials, involving two sessions completed on consecutive days (Session 1: 10-overs and Session 2: 4-overs) in 31 ± 3°C and 55 ± 17% relative humidity. Recovery interventions were administered for 20 min (mixed-method cooling vs. control) after Session 1. Measures included bowling performance (ball speed, accuracy, run-up speeds), physical demands (global positioning system, counter-movement jump), physiological (heart rate, core temperature, skin temperature, sweat loss), biochemical (creatine kinase, C-reactive protein) and perceptual variables (perceived exertion, thermal sensation, muscle soreness). Mean ball speed was higher after cooling in Session 2 (118.9 ± 8.1 vs. 115.5 ± 8.6 km · h−1; P = 0.001; d = 0.67), reducing declines in ball speed between sessions (0.24 vs. −3.18 km · h−1; P = 0.03; d = 1.80). Large effects indicated higher accuracy in Session 2 after cooling (46.0 ± 11.2 vs. 39.4 ± 8.6 arbitrary units [AU]; P = 0.13; d = 0.93) without affecting total run-up speed (19.0 ± 3.1 vs. 19.0 ± 2.5 km · h−1; P = 0.97; d = 0.01). Cooling reduced core temperature, skin temperature and thermal sensation throughout the intervention (P = 0.001–0.05; d = 1.31–5.78) and attenuated creatine kinase (P = 0.04; d = 0.56) and muscle soreness at 24-h (P = 0.03; d = 2.05). Accordingly, mixed-method cooling can reduce thermal strain after a 10-over spell and improve markers of muscular damage and discomfort alongside maintained medium-fast bowling performance on consecutive days in hot conditions.

Optimal planning of electric-vehicle charging stations in distribution systems

With the progressive exhaustion of fossil energy and the enhanced awareness of environmental protection, more attention is being paid to electric vehicles (EVs). Inappropriate siting and sizing of EV charging stations could have negative effects on the development of EVs, the layout of the city traffic network, and the convenience of EVs' drivers, and lead to an increase in network losses and a degradation in voltage profiles at some nodes. Given this background, the optimal sites of EV charging stations are first identified by a two-step screening method with environmental factors and service radius of EV charging stations considered. Then, a mathematical model for the optimal sizing of EV charging stations is developed with the minimization of total cost associated with EV charging stations to be planned as the objective function and solved by a modified primal-dual interior point algorithm (MPDIPA). Finally, simulation results of the IEEE 123-node test feeder have demonstrated that the developed model and method cannot only attain the reasonable planning scheme of EV charging stations, but also reduce the network loss and improve the voltage profile.

Evaluation of the Suchey-Brooks Method of age estimation in an Australian subpopulation using computed tomography of the pubic symphyseal surface

Despite the prominent use of the Suchey-Brooks (S-B) method of age estimation in forensic anthropological practice, it is subject to intrinsic limitations, with reports of differential inter-population error rates between geographical locations. This study assessed the accuracy of the S-B method to a contemporary adult population in Queensland, Australia and provides robust age parameters calibrated for our population. Three-dimensional surface reconstructions were generated from computed tomography scans of the pubic symphysis of male and female Caucasian individuals aged 15–70 years (n = 195) in Amira® and Rapidform®. Error was analyzed on the basis of bias, inaccuracy and percentage correct classification for left and right symphyseal surfaces. Application of transition analysis and Chi-square statistics demonstrated 63.9% and 69.7% correct age classification associated with the left symphyseal surface of Australian males and females, respectively, using the S-B method. Using Bayesian statistics, probability density distributions for each S-B phase were calculated, providing refined age parameters for our population. Mean inaccuracies of 6.77 (±2.76) and 8.28 (±4.41) years were reported for the left surfaces of males and females, respectively; with positive biases for younger individuals (<55 years) and negative biases in older individuals. Significant sexual dimorphism in the application of the S-B method was observed; and asymmetry in phase classification of the pubic symphysis was a frequent phenomenon. These results recommend that the S-B method should be applied with caution in medico-legal death investigations of Queensland skeletal remains and warrant further investigation of reliable age estimation techniques.

In vitro primary cell culture as a physiologically relevant method for preclinical testing of human oncolytic adenovirus

Ad[I/PPT-E1A] is an oncolytic adenovirus that specifically kills prostate cells via restricted replication by a prostate-specific regulatory element. Off-target replication of oncolytic adenoviruses would have serious clinical consequences. As a proposed ex vivo test, we describe the assessment of the specificity of Ad[I/PPT-E1A] viral cytotoxicity and replication in human nonprostate primary cells. Four primary nonprostate cell types were selected to mimic the effects of potential in vivo exposure to Ad[I/PPT-E1A] virus: bronchial epithelial cells, urothelial cells, vascular endothelial cells, and hepatocytes. Primary cells were analyzed for Ad[I/PPT-E1A] viral cytotoxicity in MTS assays, and viral replication was determined by hexon titer immunostaining assays to quantify viral hexon protein. The results revealed that at an extreme multiplicity of infection of 500, unlikely to be achieved in vivo, Ad[I/PPT-E1A] virus showed no significant cytotoxic effects in the nonprostate primary cell types apart from the hepatocytes. Transmission electron microscopy studies revealed high levels of Ad[I/PPT-E1A] sequestered in the cytoplasm of these cells. Adenoviral green fluorescent protein reporter studies showed no evidence for nuclear localization, suggesting that the cytotoxic effects of Ad[I/PPT-E1A] in human primary hepatocytes are related to viral sequestration. Also, hepatocytes had increased amounts of coxsackie adenovirus receptor surface protein. Active viral replication was only observed in the permissive primary prostate cells and LNCaP prostate cell line, and was not evident in any of the other nonprostate cells types tested, confirming the specificity of Ad[I/PPT-E1A]. Thus, using a relevant panel of primary human cells provides a convenient and alternative preclinical assay for examining the specificity of conditionally replicating oncolytic adenoviruses in vivo.

Within-subject intra- and inter-method consistency of two experimental risk attitude elicitation methods

We compare the consistency of choices in two methods to used elicit risk preferences on an aggregate as well as on an individual level. We asked subjects to choose twice from a list of nine decision between two lotteries, as introduced by Holt and Laury (2002, 2005) alternating with nine decisions using the budget approach introduced by Andreoni and Harbaugh (2009). We find that while on an aggregate(subject pool) level the results are (roughly) consistent, on an individual(within-subject) level,behavior is far from consistent. Within each method as well as across methods we observe low correlations. This again questions the reliability of experimental risk elicitation measures and the ability to use results from such methods to control for the risk aversion of subjects when explaining e�ects in other experimental games.

An analytical method for assessing stage-specific drug activity in Plasmodium vivax malaria : implications for ex vivo drug susceptibility testing

The emergence of highly chloroquine (CQ) resistant P. vivax in Southeast Asia has created an urgent need for an improved understanding of the mechanisms of drug resistance in these parasites, the development of robust tools for defining the spread of resistance, and the discovery of new antimalarial agents. The ex vivo Schizont Maturation Test (SMT), originally developed for the study of P. falciparum, has been modified for P. vivax. We retrospectively analysed the results from 760 parasite isolates assessed by the modified SMT to investigate the relationship between parasite growth dynamics and parasite susceptibility to antimalarial drugs. Previous observations of the stage-specific activity of CQ against P. vivax were confirmed, and shown to have profound consequences for interpretation of the assay. Using a nonlinear model we show increased duration of the assay and a higher proportion of ring stages in the initial blood sample were associated with decreased effective concentration (EC50) values of CQ, and identify a threshold where these associations no longer hold. Thus, starting composition of parasites in the SMT and duration of the assay can have a profound effect on the calculated EC50 for CQ. Our findings indicate that EC50 values from assays with a duration less than 34 hours do not truly reflect the sensitivity of the parasite to CQ, nor an assay where the proportion of ring stage parasites at the start of the assay does not exceed 66%. Application of this threshold modelling approach suggests that similar issues may occur for susceptibility testing of amodiaquine and mefloquine. The statistical methodology which has been developed also provides a novel means of detecting stage-specific drug activity for new antimalarials.

Suspected child abuse and neglect reporting behaviour

Executuve Summary Background and Aims Child abuse and neglect is a tragedy within our community, with over 10,000 substantiated reports of abuse and neglect in Queensland in the past year. The considerable consequences of child abuse and neglect are far-reaching, substantial and can be fatal. The reporting of suspicions of child abuse or neglect is often the first step in preventing further abuse or neglect. In the State of Queensland, medical practitioners are mandated by law to report their suspicions of child abuse and neglect. However, despite this mandate many still do not report their suspicions. A 1998 study indicated that 43% of medical practitioners had, at some time, made a conscious decision to not report suspected abuse or neglect (Van Haeringen, Dadds & Armstrong, 1998). The aim of this study was to gain a better understanding of beliefs about reporting suspected child abuse and neglect and the barriers to reporting suspected abuse and neglect by medical practitioners and parents and students. The findings have the potential to inform the training and education of members of the community who have a shared responsibility to protect the wellbeing of its most vulnerable members. Method In one of the largest studies of reporting behaviour in relation to suspected child abuse and neglect in Australia, we examined and compared medical practitioners’ responses with members of the community, namely parents and students. We surveyed 91 medical practitioners and 214 members of the community (102 parents and 112 students) regarding their beliefs and reporting behaviour related to suspected child abuse and neglect. We also examined reasons for not reporting suspected abuse or neglect, as well as awareness of responsibilities and the appropriate reporting procedures. To obtain such information, participants anonymously completed a comprehensive questionnaire using items from previous studies of reporting attitudes and behaviour. Executive Summary Abused Child Trust Report August 2003 5 Findings Key findings include: • The majority of medical practitioners (97%) were aware of their duty to report suspected abuse and neglect and believed they had a professional and ethical duty to do so. • A majority of parents (82%) and students (68%) also believed that they had a professional and ethical duty to report suspected abuse and neglect. • In accord with their statutory duty to report suspected abuse and neglect, 69% of medical practitioners had made a report at some point. • Sixteen percent of parents and 9% of students surveyed indicated that they had reported their suspicions of neglect and abuse. • The most endorsed belief associated with not reporting suspected child abuse and neglect was that, ‘unpleasant events would follow reporting’. • Over a quarter of medical practitioners (26%) admitted to making a decision not to report their suspicions of child abuse or neglect on at least one occasion. • Compared with previous research, there has been a decline in the number of medical practitioners who decided not to report suspected abuse or neglect from 43% (Van Haeringen et al., 1998) to 26% in the current study. • Fourteen percent of parents and 15% of students surveyed had also chosen not to report a case of suspected abuse or neglect. • Attitudes that most strongly influenced the decision to report or not report suspected abuse or neglect differed between groups (medical practitioners, parents, or students). A belief that, ‘the abuse was a single incident’ was the best predictor of non-reporting by medical practitioners, while having ‘no time to follow-up the report’ or failing to be ‘convinced of evidence of abuse’ best predicted failure to report abuse by students. A range of beliefs predicted non-reporting by parents, including the beliefs that reporting suspected abuse was ‘not their responsibility’ and ‘knowing the child had retracted their statement’. Conclusions Of major concern is that approximately 25% of medical practitioners with a mandated responsibility to report, as well as some members of the general public, revealed that they have suspected child neglect or abuse but have made the decision not to report their suspicions. Parents and students perceived the general community as having responsibility for reporting suspicions of abuse or neglect. Despite this perception, they felt that lodging a report may be overly demanding in terms of time and they had the confidence in their ability to identify child abuse and neglect. An explanation for medical practitioners deciding not to report may be based upon their optimistic belief that suspected abuse or neglect was a single incident. Our findings may best be understood from the ‘inflation of optimism’ hypothesis put forward by the Nobel Laureate, Daniel Kahneman. He suggests that in spite of rational evidence, human beings tend to make judgements based on an optimistic view rather than engaging in a rational decision-making process. In this case, despite past behaviour of abuse or neglect being the best predictor of future behaviour, medical practitioners have taken an optimistic view, choosing to believe that their suspicion of child abuse or neglect represents a single incident. The clear implication of findings in the current research is the need for the members of the general community and medical practitioners to be better appraised of the consequences of their decision-making in relation to suspicionsof child abuse and neglect. Finally findings from parents and students relating to their reporting behaviour suggest that members of the larger community represent an untapped resourcewho might, with appropriate awareness, play a more significant role in theidentification and reporting of suspected child abuse and neglect.

Contribution to single-point closure Reynolds-stress modelling of inhomogeneous flow

This paper is concerned with recent advances in the development of near wall-normal-free Reynolds-stress models, whose single point closure formulation, based on the inhomogeneity direction concept, is completely independent of the distance from the wall, and of the normal to the wall direction. In the present approach the direction of the inhomogeneity unit vector is decoupled from the coefficient functions of the inhomogeneous terms. A study of the relative influence of the particular closures used for the rapid redistribution terms and for the turbulent diffusion is undertaken, through comparison with measurements, and with a baseline Reynolds-stress model (RSM) using geometric wall normals. It is shown that wall-normal-free rsms can be reformulated as a projection on a tensorial basis that includes the inhomogeneity direction unit vector, suggesting that the theory of the redistribution tensor closure should be revised by taking into account inhomogeneity effects in the tensorial integrity basis used for its representation.

Numerical simulation of red blood cells' deformation using SPH method

A numerical simulation method for the Red Blood Cells’ (RBC) deformation is presented in this study. The two-dimensional RBC membrane is modeled by the spring network, where the elastic stretch/compression energy and the bending energy are considered with the constraint of constant RBC surface area. Smoothed Particle Hydrodynamics (SPH) method is used to solve the Navier-Stokes equation coupled with the Plasma-RBC membrane and Cytoplasm- RBC membrane interaction. To verify the method, the motion of a single RBC is simulated in Poiseuille flow and compared with the results reported earlier. Typical motion and deformation mechanism of the RBC is observed.