A preconditioned numerical solver for stiff nonlinear reaction-diffusion equations with fractional Laplacians that avoids dense matrices

The numerical solution of fractional partial differential equations poses significant computational challenges in regard to efficiency as a result of the spatial nonlocality of the fractional differential operators. The dense coefficient matrices that arise from spatial discretisation of these operators mean that even one-dimensional problems can be difficult to solve using standard methods on grids comprising thousands of nodes or more. In this work we address this issue of efficiency for one-dimensional, nonlinear space-fractional reaction–diffusion equations with fractional Laplacian operators. We apply variable-order, variable-stepsize backward differentiation formulas in a Jacobian-free Newton–Krylov framework to advance the solution in time. A key advantage of this approach is the elimination of any requirement to form the dense matrix representation of the fractional Laplacian operator. We show how a banded approximation to this matrix, which can be formed and factorised efficiently, can be used as part of an effective preconditioner that accelerates convergence of the Krylov subspace iterative solver. Our approach also captures the full contribution from the nonlinear reaction term in the preconditioner, which is crucial for problems that exhibit stiff reactions. Numerical examples are presented to illustrate the overall effectiveness of the solver.

Systematic review and meta-analysis of the effects of exercise for those with cancer-related lymphedema

Objective: To evaluate the effects of exercise on cancer-related lymphedema and related symptoms, and to determine the need for those with lymphedema to wear compression during exercise. Data Sources: CINAHL, Cochrane, Ebscohost, MEDLINE, Pubmed, ProQuest Health and Medical Complete, ProQuest Nursing and Allied Health Source, Science Direct and SPORTDiscus databases were searched for trials published prior to 1 January, 2015. Study Selection: Randomised and non-randomised, controlled trials, and single group pre-post studies published in English-language were included. Twenty-one (exercise) and four (compression and exercise) studies met inclusion criteria. Data Extraction: Data was extracted into tabular format using predefined data fields by one reviewer and assessed for accuracy by a second reviewer. Study quality was evaluated using the Effective Public Health Practice Project assessment tool. Data Synthesis: Data was pooled using a random effects model to assess the effects of acute and long-term exercise on lymphedema and lymphedema-associated symptoms, with subgroup analyses for exercise mode and intervention length. There was no effect of exercise (acute or intervention) on lymphedema or associated symptoms with standardised mean differences from all analyses ranging between −0.2 and 0.1 (p-values ≥0.22). Findings from subgroup analyses for exercise mode (aerobic, resistance, mixed, other) and intervention duration (>12 weeks or ≤12 weeks) were consistent with these findings; that is, no effect on lymphedema or associated symptoms. There were too few studies evaluating the effect of compression during regular exercise to conduct a meta-analysis. Conclusions: Individuals with secondary lymphedema can safely participate in progressive, regular exercise without experiencing a worsening of lymphedema or related-symptoms. However, the results also do not suggest any improvements will occur in lymphedema. At present, there is insufficient evidence to support or refute the current clinical recommendation to wear compression garments during regular exercise.

Common polygenic variation contributes to risk of migraine in the Norfolk Island population

Migraine has been defined as a common disabling primary headache disorder. Epidemiology studies have provided with the undeniable evidence of genetic components as active players in the development of the disease under a polygenic model in which multiple risk alleles exert modest individual effects. Our objective was to test the contribution of a polygenic effect to migraine risk in the Norfolk Island population using a panel of SNPs reported to be disease associated in published migraine GWAS. We also investigated whether individual SNPs were associated with gene expression levels measured in whole-blood. Polygenic scores were calculated in a total of 285 related individuals (74 cases, 211 controls) from the Norfolk Island using 51 SNPs previously reported to be associated with migraine in published GWAS. The association between polygenic score and migraine case-control status was tested using logistic regression. Results indicate that a migraine polygenic risk score was associated with migraine case-control status in this population (P=0.016). This supports the hypothesis that multiple SNPs with weak effects collectively contribute to migraine risk in this population. Amongst the SNPs included in the polygenic model, 4 were associated with the expression of the USMG5 gene, including rs171251 (P = 0.012). Results from this study provide evidence for a polygenic contribution to migraine risk in an isolated population and highlight specific SNPs that regulate the expression of USMG5, a gene critical for mitochondrial function.

Brief report: Intestinal dysbiosis in ankylosing spondylitis

Objective Ankylosing spondylitis (AS) is a common, highly heritable immune-mediated arthropathy that occurs in genetically susceptible individuals exposed to an unknown but likely ubiquitous environmental trigger. There is a close relationship between the gut and spondyloarthritis, as exemplified in patients with reactive arthritis, in whom a typically self-limiting arthropathy follows either a gastrointestinal or urogenital infection. Microbial involvement in AS has been suggested; however, no definitive link has been established. The aim of this study was to determine whether the gut in patients with AS carries a distinct microbial signature compared with that in the gut of healthy control subjects. Methods Microbial profiles for terminal ileum biopsy specimens obtained from patients with recent-onset tumor necrosis factor antagonist-naive AS and from healthy control subjects were generated using culture-independent 16S ribosomal RNA gene sequencing and analysis techniques. Results Our results showed that the terminal ileum microbial communities in patients with AS differ significantly (P < 0.001) from those in healthy control subjects, driven by a higher abundance of 5 families of bacteria (Lachnospiraceae [P = 0.001], Ruminococcaceae [P = 0.012], Rikenellaceae [P = 0.004], Porphyromonadaceae [P = 0.001], and Bacteroidaceae [P = 0.001]) and a decrease in the abundance of 2 families of bacteria (Veillonellaceae [P = 0.01] and Prevotellaceae [P = 0.004]). Conclusion We show evidence for a discrete microbial signature in the terminal ileum of patients with AS compared with healthy control subjects. The microbial composition was demonstrated to correlate with disease status, and greater differences were observed between disease groups than within disease groups. These results are consistent with the hypothesis that genes associated with AS act, at least in part, through effects on the gut microbiome.

Trans-scale granular modelling of cytoskeleton: A mini-review

Living cells are the functional unit of organs that controls reactions to their exterior. However, the mechanics of living cells can be difficult to characterize due to the crypticity of their microscale structures and associated dynamic cellular processes. Fortunately, multiscale modelling provides a powerful simulation tool that can be used to study the mechanical properties of these soft hierarchical, biological systems. This paper reviews recent developments in hierarchical multiscale modeling technique that aimed at understanding cytoskeleton mechanics. Discussions are expanded with respects to cytoskeletal components including: intermediate filaments, microtubules and microfilament networks. The mechanical performance of difference cytoskeleton components are discussed with respect to their structural and material properties. Explicit granular simulation methods are adopted with different coarse-grained strategies for these cytoskeleton components and the simulation details are introduced in this review.

On the effect of Reynolds number for flow past three side-by-side square cylinders for unequal gap spacings

Flow patterns and aerodynamic characteristics behind three side-by-side square cylinders has been found depending upon the unequal gap spacing (g1 = s1/d and g2 = s2/d) between the three cylinders and the Reynolds number (Re) using the Lattice Boltzmann method. The effect of Reynolds numbers on the flow behind three cylinders are numerically studied for 75 ≤ Re ≤ 175 and chosen unequal gap spacings such as (g1, g2) = (1.5, 1), (3, 4) and (7, 6). We also investigate the effect of g2 while keeping g1 fixed for Re = 150. It is found that a Reynolds number have a strong effect on the flow at small unequal gap spacing (g1, g2) = (1.5, 1.0). It is also found that the secondary cylinder interaction frequency significantly contributes for unequal gap spacing for all chosen Reynolds numbers. It is observed that at intermediate unequal gap spacing (g1, g2) = (3, 4) the primary vortex shedding frequency plays a major role and the effect of secondary cylinder interaction frequencies almost disappear. Some vortices merge near the exit and as a result small modulation found in drag and lift coefficients. This means that with the increase in the Reynolds numbers and unequal gap spacing shows weakens wakes interaction between the cylinders. At large unequal gap spacing (g1, g2) = (7, 6) the flow is fully periodic and no small modulation found in drag and lift coefficients signals. It is found that the jet flows for unequal gap spacing strongly influenced the wake interaction by varying the Reynolds number. These unequal gap spacing separate wake patterns for different Reynolds numbers: flip-flopping, in-phase and anti-phase modulation synchronized, in-phase and anti-phase synchronized. It is also observed that in case of equal gap spacing between the cylinders the effect of gap spacing is stronger than the Reynolds number. On the other hand, in case of unequal gap spacing between the cylinders the wake patterns strongly depends on both unequal gap spacing and Reynolds number. The vorticity contour visualization, time history analysis of drag and lift coefficients, power spectrum analysis of lift coefficient and force statistics are systematically discussed for all chosen unequal gap spacings and Reynolds numbers to fully understand this valuable and practical problem.

The impact of orbital errors on the estimation of satellite clock errors and PPP

Precise satellite orbit and clocks are essential for providing high accuracy real-time PPP (Precise Point Positioning) service. However, by treating the predicted orbits as fixed, the orbital errors may be partially assimilated by the estimated satellite clock and hence impact the positioning solutions. This paper presents the impact analysis of errors in radial and tangential orbital components on the estimation of satellite clocks and PPP through theoretical study and experimental evaluation. The relationship between the compensation of the orbital errors by the satellite clocks and the satellite-station geometry is discussed in details. Based on the satellite clocks estimated with regional station networks of different sizes (∼100, ∼300, ∼500 and ∼700 km in radius), results indicated that the orbital errors compensated by the satellite clock estimates reduce as the size of the network increases. An interesting regional PPP mode based on the broadcast ephemeris and the corresponding estimated satellite clocks is proposed and evaluated through the numerical study. The impact of orbital errors in the broadcast ephemeris has shown to be negligible for PPP users in a regional network of a radius of ∼300 km, with positioning RMS of about 1.4, 1.4 and 3.7 cm for east, north and up component in the post-mission kinematic mode, comparable with 1.3, 1.3 and 3.6 cm using the precise orbits and the corresponding estimated clocks. Compared with the DGPS and RTK positioning, only the estimated satellite clocks are needed to be disseminated to PPP users for this approach. It can significantly alleviate the communication burdens and therefore can be beneficial to the real time applications.

The burden of child maltreatment in the East Asia and Pacific region

This study estimated the health and economic burden of child maltreatment in the East Asia and Pacific region, addressing a significant gap in the current evidence base. Systematic reviews and meta-analyses were conducted to estimate the prevalence of child physical abuse, sexual abuse, emotional abuse, neglect, and witnessing parental violence. Population Attributable Fractions were calculated and Disability-Adjusted Life Years (DALYs) lost from physical and mental health outcomes and health risk behaviors attributable to child maltreatment were estimated using the most recent comparable Global Burden of Disease data. DALY losses were converted into monetary value by assuming that one DALY is equal to the sub-region’s per capita GDP. The estimated economic value of DALYs lost to violence against children as a percentage of GDP ranged from 1.24% to 3.46% across sub-regions defined by the World Health Organization. The estimated economic value of DALYs (in constant 2000 US$) lost to child maltreatment in the EAP region totaled US $151 billion, accounting for 1.88% of the region’s GDP. Updated to 2012 dollars, the estimated economic burden totaled US $194 billion. In sensitivity analysis, the aggregate costs as a percentage of GDP range from 1.36% to 2.52%. The economic burden of child maltreatment in the East Asia and Pacific region is substantial, indicating the importance of preventing and responding to child maltreatment in this region. More comprehensive research into the impact of multiple types of childhood adversity on a wider range of putative health outcomes is needed to guide policy and programs for child protection in the region, and globally.

X-ray Structure of Gelonin at 1.8 Å Resolution

Gelonin is a single chain ribosome inactivating protein (RIP) with potential application in the treatment of cancer and AIDS. Diffraction quality crystals grown using PEG3350, belong to the space group P2(1), with it a = 49.4 Angstrom b = 44.9 Angstrom, c = 137.4 Angstrom and beta = 98.4 degrees, and contain two molecules in the asymmetric unit. Diffraction data collected to 1.8 Angstrom resolution has a R(m) value of 7.3%. Structure of gelonin has been solved by the molecular replacement method, using ricin A chain as the search model. Crystallographic refinement using X-PLOR resulted in a model for which the r.m.s deviations from ideal bond lengths and bond angles are 0.012 Angstrom and 2.7 degrees, respectively The final R-factor is 18.4% for 39,806 reflections for which I > 1.0 sigma(I).The C-alpha atoms of the two molecules in the asymmetric unit superpose to within 0.38 Angstrom for 247 atom pairs. The overall fold of gelonin is similar to that of other RIPs such as ricin A chain and alpha-momorcharin, the r.m.s.d. for C-alpha superpositions being 1.3 and 1.4 Angstrom, respectively The-catalytic residues (Glu166, Arg169 and Tyr113) in the active site form a hydrogen bond scheme similar to that observed in other RIPs. The conformation of Tyr74 in the active site, however, is significantly different from that in alpha-momorcharin. Three well defined water molecules are located in the active site cavity and one of them, X319, superposes to within 0.2 Angstrom of a corresponding water molecule in the structure of alpha-momorcharin. Any of the three could be the substrate water molecule in the hydrolysis reaction catalysed by gelonin.Difference electron density for a N-linked sugar moiety has been observed near only one of the two potential glycosylation sites in the sequence. The amino acid at position 239 has been established as Lys by calculation of omit electron density maps.The two cysteine residues in the sequence, Cys44 and Cys50, form a disulphide bond, and are therefore not available for disulphide conjugation with antibodies. Based on the structure, the region of the molecule that is involved in intradimer interactions is suggested to be suitable for introducing a Cys residue for purposes of conjugation with an antibody to produce useful immunotoxins.

Study of carotid arterial plaque stress for symptomatic and asymptomatic patients

Stroke is one of the leading causes of death in the world, resulting mostly from the sudden ruptures of atherosclerosis carotid plaques. Until now, the exact plaque rupture mechanism has not been fully understood, and also the plaque rupture risk stratification. The advanced multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components to be visualized in-vivo and reconstructed by computational modeling. In the study, plaque stress analysis using fully coupled fluid structure interaction was applied to 20 patients (12 symptomatic and 8 asymptomatic) reconstructed from in-vivo MRI, followed by a detailed biomechanics analysis, and morphological feature study. The locally extreme stress conditions can be found in the fibrous cap region, 85% at the plaque shoulder based on the present study cases. Local maximum stress values predicted in the plaque region were found to be significantly higher in symptomatic patients than that in asymptomatic patients (200±43. kPa vs. 127±37. kPa, p=0.001). Plaque stress level, defined by excluding 5% highest stress nodes in the fibrous cap region based on the accumulative histogram of stress experienced on the computational nodes in the fibrous cap, was also significantly higher in symptomatic patients than that in asymptomatic patients (154±32. kPa vs. 111±23. kPa, p<0.05). Although there was no significant difference in lipid core size between the two patient groups, symptomatic group normally had a larger lipid core and a significantly thinner fibrous cap based on the reconstructed plaques using 3D interpolation from stacks of 2D contours. Plaques with a higher stenosis were more likely to have extreme stress conditions upstream of plaque throat. The combined analyses of plaque MR image and plaque stress will advance our understanding of plaque rupture, and provide a useful tool on assessing plaque rupture risk.

Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected

Thermal conductivity of a new carbon nanotube analog: The diamond nanothread

Based on the non-equilibrium molecular dynamics simulations, we have studied the thermal conductivities of a novel ultra-thin one-dimensional carbon nanomaterial - diamond nanothread (DNT). Unlike single-wall carbon nanotube (CNT), the existence of the Stone-Wales transformations in DNT endows it with richer thermal transport characteristics. There is a transition from wave-dominated to particle-dominated transport region, which depends on the length of poly-benzene rings. However, independent of the transport region, strong length dependence in thermal conductivity is observed in DNTs with different lengths of poly-benzene ring. The distinctive SW characteristic in DNT provides more degrees of freedom to tune the thermal conductivity not found in the homogeneous structure of CNT. Therefore, DNT is an ideal platform to investigate various thermal transport mechanisms at the nanoscale. Its high tunability raises the potential to design DNTs for different applications, such as thermal connection and temperature management.

Depth profile composition studies of thin film CdS: Cu,S solar cells using XPS and AES

Surface composition and depth profile studies of hemiplated thin film CdS:CuzS solar cells have been carried out using x-ray photoelectron spectroscopy (XPS) and Auger electron spectroscopy (AES) techniques. These studies indicate that the junction is fairly diffused in the as-prepared cell. However, heat treatment of the cell at 210°C in air relatively sharpens the junction and improves the cell performance. Using the Cu(2p3p)/S(2p) ratio as well as the Cu(LVV)/(LMM) Auger intensity ratio, it can be inferred that the nominal valency of copper in the layers above the junction is Cut and it is essentially in the CUSS form. Copper signals are observed from layers deep down in the cell. These seem to appear mostly from the grain boundary region. From the observed concentration of Cd, Cu and S in these deeper layers and the Cu(LVV)/(LMM) ratio it appears that the signals from copper essentially originate partly from copper in CuS and partly from Cu2t trapped in the lattice. It is significant to note that the nominal valence state of copper changes rather abruptly from Cut to Cuz+ across the junction.

Odour emissions from anaerobic piggery ponds. 1. Results of a three season, 14-month survey

Odour emission rates were measured for seven different anaerobic ponds treating piggery wastes at six to nine discrete locations across the surface of each pond on each sampling occasion over a 14-month period. Emission rate values varied between ponds, between seasons for the same pond and even for the same pond on different days of a sampling week. Average seasonal emission rates ranged from 7.9 to 46.5 OU/m2 s, while average emission rates ranged from 16.0 to 29.0 OU/m2 s. Factors potentially responsible for the variability in emission rates were investigated, including air and pond liquor temperatures, time of day of sample collection, season and the impact of a prolonged drought.