Computationally efficient methods for solving time-variable-order time-space fractional reaction-diffusion equation

Fractional differential equations are becoming more widely accepted as a powerful tool in modelling anomalous diffusion, which is exhibited by various materials and processes. Recently, researchers have suggested that rather than using constant order fractional operators, some processes are more accurately modelled using fractional orders that vary with time and/or space. In this paper we develop computationally efficient techniques for solving time-variable-order time-space fractional reaction-diffusion equations (tsfrde) using the finite difference scheme. We adopt the Coimbra variable order time fractional operator and variable order fractional Laplacian operator in space where both orders are functions of time. Because the fractional operator is nonlocal, it is challenging to efficiently deal with its long range dependence when using classical numerical techniques to solve such equations. The novelty of our method is that the numerical solution of the time-variable-order tsfrde is written in terms of a matrix function vector product at each time step. This product is approximated efficiently by the Lanczos method, which is a powerful iterative technique for approximating the action of a matrix function by projecting onto a Krylov subspace. Furthermore an adaptive preconditioner is constructed that dramatically reduces the size of the required Krylov subspaces and hence the overall computational cost. Numerical examples, including the variable-order fractional Fisher equation, are presented to demonstrate the accuracy and efficiency of the approach.

A research framework for assessing the effects of heat stress on construction workers

Hong Kong in summer (June - October) is hot and humid. Construction workers have to undertake physically demanding activities and often in confined spaces. They are vulnerable to heat stress in summer hence health and safety measures associated to heat stress measured by scientific and clinical parameters are urgently needed. This paper provides an initial report of a research project funded by the Research Grants Council (RGC) of the HKSAR. The aim of this study is to develop a set of indices measured by clinical and scientific methods to detect impending attacks of heat stress. These indices would be of tremendous value in better safeguarding workers’ health and safety by reducing the occurrences of heat stress on site. This paper firstly reports on the statistics of construction incidents arising from heat stress. Qualitative and quantitative research methods applied in conducting the research are discussed. It is believed that the construction industry and the government would benefit a lot as a result of this study.

Investigating the change in three dimensional deformity for idiopathic scoliosis using axially loaded MRI

Background: Adolescent idiopathic scoliosis is a complex three-dimensional deformity, involving a lateral deformity in the coronal plane and axial rotation of the vertebrae in the transverse plane. Gravitational loading plays an important biomechanical role in governing the coronal deformity, however, less is known about how they influence the axial deformity. This study investigates the change in three-dimensional deformity of a series of scoliosis patients due to compressive axial loading. Methods: Magnetic resonance imaging scans were obtained and coronal deformity (measured using the coronal Cobb angle) and axial rotations measured for a group of 18 scoliosis patients (Mean major Cobb angle was 43.4 o). Each patient was scanned in an unloaded and loaded condition while compressive loads equivalent to 50% body mass were applied using a custom developed compressive device. Findings: The mean increase in major Cobb angle due to compressive loading was 7.4 o (SD 3.5 o). The most axially rotated vertebra was observed at the apex of the structural curve and the largest average intravertebral rotations were observed toward the limits of the coronal deformity. A level-wise comparison showed no significant difference between the average loaded and unloaded vertebral axial rotations (intra-observer error = 2.56 o) or intravertebral rotations at each spinal level. Interpretation: This study suggests that the biomechanical effects of axial loading primarily influence the coronal deformity, with no significant change in vertebral axial rotation or intravertebral rotation observed between the unloaded and loaded condition. However, the magnitude of changes in vertebral rotation with compressive loading may have been too small to detect given the resolution of the current technique.

A new approach to automatically producing schedules for cane railways

The scheduling of locomotive movements on cane railways has proven to be a very complex task. Various optimisation methods have been used over the years to try and produce an optimised schedule that eliminates or minimises bin supply delays to harvesters and the factory, while minimising the number of locomotives, locomotive shifts and cane bins, and also the cane age. This paper reports on a new attempt to develop an automatic scheduler using a mathematical model solved using mixed integer programming and constraint programming approaches and blocking parallel job shop scheduling fundamentals. The model solution has been explored using conventional constraint programming search techniques and found to produce a reasonable schedule for small-scale problems with up to nine harvesters. While more effort is required to complete the development of the full model with metaheuristic search techniques, the work completed to date gives confidence that the metaheuristic techniques will provide near optimal solutions in reasonable time.

Low rank Runge-Kutta methods, symplecticity and stochastic Hamiltonian problems with additive noise

In this paper we extend the ideas of Brugnano, Iavernaro and Trigiante in their development of HBVM($s,r$) methods to construct symplectic Runge-Kutta methods for all values of $s$ and $r$ with $s\geq r$. However, these methods do not see the dramatic performance improvement that HBVMs can attain. Nevertheless, in the case of additive stochastic Hamiltonian problems an extension of these ideas, which requires the simulation of an independent Wiener process at each stage of a Runge-Kutta method, leads to methods that have very favourable properties. These ideas are illustrated by some simple numerical tests for the modified midpoint rule.

Numerical investigation of mechanical properties of nanowires : a review

Nanowires (NWs) have attracted intensive researches owing to the broad applications that arise from their remarkable properties. Over the last decade, immense numerical studies have been conducted for the numerical investigation of mechanical properties of NWs. Among these numerical simulations, the molecular dynamics (MD) plays a key role. Herein we present a brief review on the current state of the MD investigation of nanowires. Emphasis will be placed on the FCC metal NWs, especially the Cu NWs. MD investigations of perfect NWs’ mechanical properties under different deformation conditions including tension, compression, torsion and bending are firstly revisited. Following in succession, the studies for defected NWs including the defects of twin boundaries (TBs) and pre-existing defects are discussed. The different deformation mechanism incurred by the presentation of defects is explored and discussed. This review reveals that the numerical simulation is an important tool to investigate the properties of NWs. However, the substantial gaps between the experimental measurements and MD results suggest the urgent need of multi-scale simulation technique.

3D-printing of highly uniform CaSiO3 ceramic scaffolds : preparation, characterization and in vivo osteogenesis

Calcium silicate (CaSiO3, CS) ceramics have received significant attention for application in bone regeneration due to their excellent in vitro apatite-mineralization ability; however, how to prepare porous CS scaffolds with a controllable pore structure for bone tissue engineering still remains a challenge. Conventional methods could not efficiently control the pore structure and mechanical strength of CS scaffolds, resulting in unstable in vivo osteogenesis. The aim of this study is to set out to solve these problems by applying a modified 3D-printing method to prepare highly uniform CS scaffolds with controllable pore structure and improved mechanical strength. The in vivo osteogenesis of the prepared 3D-printed CS scaffolds was further investigated by implanting them in the femur defects of rats. The results show that the CS scaffolds prepared by the modified 3D-printing method have uniform scaffold morphology. The pore size and pore structure of CS scaffolds can be efficiently adjusted. The compressive strength of 3D-printed CS scaffolds is around 120 times that of conventional polyurethane templated CS scaffolds. 3D-Printed CS scaffolds possess excellent apatite-mineralization ability in simulated body fluids. Micro-CT analysis has shown that 3D-printed CS scaffolds play an important role in assisting the regeneration of bone defects in vivo. The healing level of bone defects implanted by 3D-printed CS scaffolds is obviously higher than that of 3D-printed b-tricalcium phosphate (b-TCP) scaffolds at both 4 and 8 weeks. Hematoxylin and eosin (H&E) staining shows that 3D-printed CS scaffolds induce higher quality of the newly formed bone than 3D-printed b-TCP scaffolds. Immunohistochemical analyses have further shown that stronger expression of human type I collagen (COL1) and alkaline phosphate (ALP) in the bone matrix occurs in the 3D-printed CS scaffolds than in the 3D-printed b-TCP scaffolds. Considering these important advantages, such as controllable structure architecture, significant improvement in mechanical strength, excellent in vivo osteogenesis and since there is no need for second-time sintering, it is indicated that the prepared 3D-printed CS scaffolds are a promising material for application in bone regeneration.

Benign Epilepsy with Centro-temporal Spikes: Spike Triggered fMRI Shows Somato-sensory Cortex Activity

Summary:  Objective: We performed spike triggered functional MRI (fMRI) in a 12 year old girl with Benign Epilepsy with Centro-temporal Spikes (BECTS) and left-sided spikes. Our aim was to demonstrate the cerebral origin of her interictal spikes. Methods: EEG was recorded within the 3 Tesla MRI. Whole brain fMRI images were acquired, beginning 2–3 seconds after spikes. Baseline fMRI images were acquired when there were no spikes for 20 seconds. Image sets were compared with the Student's t-test. Results: Ten spike and 20 baseline brain volumes were analysed. Focal activiation was seen in the inferior left sensorimotor cortex near the face area. The anterior cingulate was more active during baseline than spikes. Conclusions: Left sided epileptiform activity in this patient with BECTS is associated with fMRI activation in the left face region of the somatosensory cortex, which would be consistent with the facial sensorimotor involvement in BECT seizures. The presence of BOLD signal change in other regions raises the possibility that the scalp recorded field of this patient with BECTs may reflect electrical change in more than one brain region.

Myoinositol Abnormalities in Temporal Lobe Epilepsy

Purpose: This study used magnetic resonance spectroscopy (MRS) to examine metabolite abnormalities in the temporal and frontal lobe of patients with temporal lobe epilepsy (TLE) of differing severity. Methods: We investigated myoinositol in TLE by using short-echo MRS in 34 TLE patients [26 late onset (LO-TLE), eight hippocampal sclerosis (HS-TLE)], and 16 controls. Single-voxel short-echo (35 ms) MR spectra of temporal and frontal lobes were acquired at 1.5 T and analyzed by using LCModel. Results: The temporal lobe ipsilateral to seizure origin in HS-TLE, but not LO-TLE, had reduced N-acetylaspartate (NA) and elevated myoinositol (MI; HS-TLE NA, 7.8 ± 1.9 mM, control NA, 9.2 ± 1.3 mM; p < 0.05; HS-TLE MI, 6.1 ± 1.6 mM, control mI 4.9 ± 0.8 mM, p< 0.05). Frontal lobe MI was low in both patient groups (LO-TLE, 4.3 ± 0.8 mM; p < 0.05; HS-TLE, 3.6 ±.05 mM; p < 0.001; controls, 4.8 ± 0.5 mM). Ipsilateral frontal lobes had lower MI (3.8 ± 0.7 mM; p < 0.01) than contralateral frontal lobes (4.3 ± 0.8 mM; p < 0.05). Conclusions: MI changes may distinguish between the seizure focus, where MI is increased, and areas of seizure spread where MI is decreased.

Implications of lipid biology for the pathogenesis of schizophrenia

Objective: Preclinical and clinical data suggest that lipid biology is integral to brain development and neurodegeneration. Both aspects are proposed as being important in the pathogenesis of schizophrenia. The purpose of this paper is to examine the implications of lipid biology, in particular the role of essential fatty acids (EFA), for schizophrenia. Methods: Medline databases were searched from 1966 to 2001 followed by the crosschecking of references. Results: Most studies investigating lipids in schizophrenia described reduced EFA, altered glycerophospholipids and an increased activity of a calcium-independent phospholipase A2 in blood cells and in post-mortem brain tissue. Additionally, in vivo brain phosphorus-31 Magnetic Resonance Spectroscopy (31P-MRS) demonstrated lower phosphomonoesters (implying reduced membrane precursors) in first- and multi-episode patients. In contrast, phosphodiesters were elevated mainly in first-episode patients (implying increased membrane breakdown products), whereas inconclusive results were found in chronic patients. EFA supplementation trials in chronic patient populations with residual symptoms have demonstrated conflicting results. More consistent results were observed in the early and symptomatic stages of illness, especially if EFA with a high proportion of eicosapentaenoic acid was used. Conclusion: Peripheral blood cell, brain necropsy and 31P-MRS analysis reveal a disturbed lipid biology, suggesting generalized membrane alterations in schizophrenia. 31P-MRS data suggest increased membrane turnover at illness onset and persisting membrane abnormalities in established schizophrenia. Cellular processes regulating membrane lipid metabolism are potential new targets for antipsychotic drugs and might explain the mechanism of action of treatments such as eicosapentaenoic acid.

Factors affecting 39K NMR detectability in rat tissue

In this study we have found that NMR detectability of 39K in rat thigh muscle may be substantially higher (up to 100% oftotal tissue potassium) than values previously reported of around 40%. The signal was found to consist of two superimposed components, one broad and one narrow, of approximately equal area. Investigations involving improvements in spectral parameters such as signal-to-noise ratio and baseline roll, together with computer simulations of spectra, show that the quality of the spectra has a major effect on the amount of signal detected, which is largely due to the loss of detectability of the broad signal component. In particular, lower-field spectrometers using conventional probes and detection methods generally have poorer signal-to-noise and worse baseline roll artifacts, which make detection of a broad component of the muscle signal difficult.

Influence of plate bending stiffness on mechanical properties of healing fracture, mouse femora

Mechanically well-defined stabilization systems have only recently become available, providing standardized conditions for studying the role of the mechanical environment on mouse bone fracture healing. The aim of this study was to characterize the time course of strength recovery and callus development of mouse femoral osteotomies stabilized with either low or high flexibility (in bending and torsion) internal fixation plates. Animals were euthanized and femora excised at 14, 21, and 28 days post-osteotomy for microCT analysis and torsional strength testing. While a larger mineralized callus was observed in osteotomies under more flexible conditions at all time points, the earlier bridging of the mineralized callus under less flexible conditions by 1 week resulted in an earlier recovery of torsional strength in mice stabilized with low flexibility fixation. Ultimate torque values for these bones were significantly higher at 14 and 21 days post-osteotomy compared to bones with the more flexible stabilization. Our study confirms the high reproducibility of the results that are achieved with this new implant system, therefore making it ideal for studying the influence of the mechanical environment on murine fracture healing under highly standardized conditions.

Evolutionary design of bone scaffolds with reference to material selection

The favourable scaffold for bone tissue engineering should have desired characteristic features, such as adequate mechanical strength and three-dimensional open porosity, which guarantee a suitable environment for tissue regeneration. In fact, the design of such complex structures like bone scaffolds is a challenge for investigators. One of the aims is to achieve the best possible mechanical strength-degradation rate ratio. In this paper we attempt to use numerical modelling to evaluate material properties for designing bone tissue engineering scaffold fabricated via the fused deposition modelling technique. For our studies the standard genetic algorithm was used, which is an efficient method of discrete optimization. For the fused deposition modelling scaffold, each individual strut is scrutinized for its role in the architecture and structural support it provides for the scaffold, and its contribution to the overall scaffold was studied. The goal of the study was to create a numerical tool that could help to acquire the desired behaviour of tissue engineered scaffolds and our results showed that this could be achieved efficiently by using different materials for individual struts. To represent a great number of ways in which scaffold mechanical function loss could proceed, the exemplary set of different desirable scaffold stiffness loss function was chosen. © 2012 John Wiley & Sons, Ltd.

Non-combinatorial library screening reveals subsite cooperativity and identifies new high-efficiency substrates for kallikrein-related peptidase 14

An array of substrates link the tryptic serine protease, kallikrein-related peptidase 14 (KLK14), to physiological functions including desquamation and activation of signaling molecules associated with inflammation and cancer. Recognition of protease cleavage sequences is driven by complementarity between exposed substrate motifs and the physicochemical signature of an enzyme's active site cleft. However, conventional substrate screening methods have generated conflicting subsite profiles for KLK14. This study utilizes a recently developed screening technique, the sparse matrix library, to identify five novel high-efficiency sequences for KLK14. The optimal sequence, YASR, was cleaved with higher efficiency (k(cat)/K(m)=3.81 ± 0.4 × 10(6) M(-1) s(-1)) than favored substrates from positional scanning and phage display by 2- and 10-fold, respectively. Binding site cooperativity was prominent among preferred sequences, which enabled optimal interaction at all subsites as indicated by predictive modeling of KLK14/substrate complexes. These simulations constitute the first molecular dynamics analysis of KLK14 and offer a structural rationale for the divergent subsite preferences evident between KLK14 and closely related KLKs, KLK4 and KLK5. Collectively, these findings highlight the importance of binding site cooperativity in protease substrate recognition, which has implications for discovery of optimal substrates and engineering highly effective protease inhibitors.