Analytical and numerical solutions of the space and time fractional bloch-torrey equation

Fractional order dynamics in physics, particularly when applied to diffusion, leads to an extension of the concept of Brown-ian motion through a generalization of the Gaussian probability function to what is termed anomalous diffusion. As MRI is applied with increasing temporal and spatial resolution, the spin dynamics are being examined more closely; such examinations extend our knowledge of biological materials through a detailed analysis of relaxation time distribution and water diffusion heterogeneity. Here the dynamic models become more complex as they attempt to correlate new data with a multiplicity of tissue compartments where processes are often anisotropic. Anomalous diffusion in the human brain using fractional order calculus has been investigated. Recently, a new diffusion model was proposed by solving the Bloch-Torrey equation using fractional order calculus with respect to time and space (see R.L. Magin et al., J. Magnetic Resonance, 190 (2008) 255-270). However effective numerical methods and supporting error analyses for the fractional Bloch-Torrey equation are still limited. In this paper, the space and time fractional Bloch-Torrey equation (ST-FBTE) is considered. The time and space derivatives in the ST-FBTE are replaced by the Caputo and the sequential Riesz fractional derivatives, respectively. Firstly, we derive an analytical solution for the ST-FBTE with initial and boundary conditions on a finite domain. Secondly, we propose an implicit numerical method (INM) for the ST-FBTE, and the stability and convergence of the INM are investigated. We prove that the implicit numerical method for the ST-FBTE is unconditionally stable and convergent. Finally, we present some numerical results that support our theoretical analysis.

Preparation, characterisation and in vivo osteogenesis of mesoporous bioactive glasses

Bone defects, especially large bone defects, remain a major challenge in orthopaedic surgery. Autologous bone transplantation is considered the most effective treatment, but insufficient donor tissue, coupled with concerns about donor site morbidity, has hindered this approach in large-scale applications. Alternative approaches include implanting biomaterials such as bioactive glass (BG), which has been widely used for bone defect healing, due to having generally good biocompatibility, and can be gradually biodegraded during the process of new bone formation. Mesoporous bioactive glass (MBG) is a newly developed bioactive glass which has been proven to have enhanced in-vitro bioactivity; however the in-vivo osteogenesis has not been studied. A critical problem in using the bone tissue engineering approach to restore large bone defects is that the nutrient supply and cell viability at the centre of the scaffold is severely hampered since the diffusion distance of nutrients and oxygen for cell survival is limited to 150-200µm. Cobalt ions has been shown to mimic hypoxia, which plays a pivotal role in coupling angiogenesis with osteogenesis in-vivo by activating hypoxia inducing factor-1α (HIF-1α) transcription factor, subsequently initiating the expression of genes associated with tissue regeneration. Therefore, one aim of this study is to investigate the in-vivo osteogenesis of MBG by comparison with BG and β-TCP, which are widely used clinically. The other aim is to explore hypoxia-mimicking biomaterials by incorporating Cobalt into MBG and β-TCP. MBG and β-TCP incorporated with 5% cobalt (5Co-MBG and 5CCP) have also been studied in-vivo to determine whether the hypoxic effect has a beneficial effect on the bone formation. The composition and microstructure of synthesised materials (BG, MBG, 5Co-MBG, 5CCP) were characterised, along with the mesopore properties of the MBG materials. Dissolution and cytotoxicity of the Co-containing materials were also investigated. Femoral samples with defects harvested at 4 and 8 weeks were scanned using micro-CT followed by processing for histology (H&E staining) to determine bone formation. Histology of MBG showed a slower rate of bone formation at 4 weeks than BG, however at 8 weeks it could be clearly seen that MBG had more bone formation. The in-vivo results show that the osteogenesis of MBG reciprocates the enhanced performance shown in-vitro compared to BG. Dissolution study showed that Co ions can be efficiently released from MBG and β-TCP in a controllable way. Low amounts of Co incorporated into the MBG and β-TCP showed no significant cytotoxicity and the Co-MBG powders maintained a mesopore structure although not as highly ordered as pure MBG. Preliminary study has shown that Co incorporated samples showed little to no bone formation, instead incurring high lymphocyte activity. Further studies need to be done on Co incorporated materials to determine the cause for high lymphocyte activity in-vivo, which appear to hinder bone formation. In conclusion, this study demonstrated the osteogenic activity of MBG and provided some valuable information of tissue reaction to Co-incorporated MBG and TCP materials.

Removal of toxic metals in a multi metal system using sorbents for potential application to urban stormwater treatment

In the context of increasing demand for potable water and the depletion of water resources, stormwater is a logical alternative. However, stormwater contains pollutants, among which metals are of particular interest due to their toxicity and persistence in the environment. Hence, it is imperative to remove toxic metals in stormwater to the levels prescribed by drinking water guidelines for potable use. Consequently, various techniques have been proposed, among which sorption using low cost sorbents is economically viable and environmentally benign in comparison to other techniques. However, sorbents show affinity towards certain toxic metals, which results in poor removal of other toxic metals. It was hypothesised in this study that a mixture of sorbents that have different metal affinity patterns can be used for the efficient removal of a range of toxic metals commonly found in stormwater. The performance of six sorbents in the sorption of Al, Cr, Cu, Pb, Ni, Zn and Cd, which are the toxic metals commonly found in urban stormwater, was investigated to select suitable sorbents for creating the mixtures. For this purpose, a multi criteria analytical protocol was developed using the decision making methods: PROMETHEE (Preference Ranking Organisation METHod for Enrichment Evaluations) and GAIA (Graphical Analysis for Interactive Assistance). Zeolite and seaweed were selected for the creation of trial mixtures based on their metal affinity pattern and the performance on predetermined selection criteria. The metal sorption mechanisms employed by seaweed and zeolite were defined using kinetics, isotherm and thermodynamics parameters, which were determined using the batch sorption experiments. Additionally, the kinetics rate-limiting steps were identified using an innovative approach using GAIA and Spearman correlation techniques developed as part of the study, to overcome the limitation in conventional graphical methods in predicting the degree of contribution of each kinetics step in limiting the overall metal removal rate. The sorption kinetics of zeolite was found to be primarily limited by intraparticle diffusion followed by the sorption reaction steps, which were governed mainly by the hydrated ionic diameter of metals. The isotherm study indicated that the metal sorption mechanism of zeolite was primarily of a physical nature. The thermodynamics study confirmed that the energetically favourable nature of sorption increased in the order of Zn < Cu < Cd < Ni < Pb < Cr < Al, which is in agreement with metal sorption affinity of zeolite. Hence, sorption thermodynamics has an influence on the metal sorption affinity of zeolite. On the other hand, the primary kinetics rate-limiting step of seaweed was the sorption reaction process followed by intraparticle diffusion. The boundary layer diffusion was also found to limit the metal sorption kinetics at low concentration. According to the sorption isotherm study, Cd, Pb, Cr and Al were sorbed by seaweed via ion exchange, whilst sorption of Ni occurred via physisorption. Furthermore, ionic bonding is responsible for the sorption of Zn. The thermodynamics study confirmed that sorption by seaweed was energetically favourable in the order of Zn < Cu < Cd < Cr . Al < Pb < Ni. However, this did not agree with the affinity series derived for seaweed suggesting a limited influence of sorption thermodynamics on metal affinity for seaweed. The investigation of zeolite-seaweed mixtures indicated that mixing sorbents have an effect on the kinetics rates and the sorption affinity. Additionally, the theoretical relationships were derived to predict the boundary layer diffusion rate, intraparticle diffusion rate, the sorption reaction rate and the enthalpy of mixtures based on that of individual sorbents. In general, low coefficient of determination (R2) for the relationships between theoretical and experimental data indicated that the relationships were not statistically significant. This was attributed to the heterogeneity of the properties of sorbents. Nevertheless, in relative terms, the intraparticle diffusion rate, sorption reaction rate and enthalpy of sorption had higher R2 values than the boundary layer diffusion rate suggesting that there was some relationship between the former set of parameters of mixtures and that of sorbents. The mixture, which contained 80% of zeolite and 20% of seaweed, showed similar affinity for the sorption of Cu, Ni, Cd, Cr and Al, which was attributed to approximately similar sorption enthalpy of the metal ions. Therefore, it was concluded that the seaweed-zeolite mixture can be used to obtain the same affinity for various metals present in a multi metal system provided the metal ions have similar enthalpy during sorption by the mixture.

Explaining internationalisation of small to medium sized enterprises within the Queensland food and beverage industry

The prevalence of International New Ventures (INVs) has increased during the past twenty years. Nevertheless, to date there has been no general consensus within the literature on an explanation as to the rapid internationalisation of some firms. Do they follow a similar process to other firms that internationalise based on a more ‘measured’ incremental sequential process of internationalisation. This paper proposes and tests an innovation diffusion model of the internationalisation of small firm INVs and others by drawing on key innovation diffusion models from the literature. The results of this analysis indicate that the synthesised model of export adoption is effective in explaining the internationalisation process of INVs and other firms within the Queensland Food and Beverage Industry. Significantly, the features of the original innovation diffusion models developed in the consumer behaviour literature, which had limited examination within the internationalisation literature, were confirmed. This includes the ability of firms, or specifically decision-makers, to skip stages based on previous experience.

Solution methods for advection-diffusion-reaction equations on growing domains and subdomains, with application to modelling skin substitutes

Problems involving the solution of advection-diffusion-reaction equations on domains and subdomains whose growth affects and is affected by these equations, commonly arise in developmental biology. Here, a mathematical framework for these situations, together with methods for obtaining spatio-temporal solutions and steady states of models built from this framework, is presented. The framework and methods are applied to a recently published model of epidermal skin substitutes. Despite the use of Eulerian schemes, excellent agreement is obtained between the numerical spatio-temporal, numerical steady state, and analytical solutions of the model.

Langevin dynamics modeling of the water diffusion tensor in partially aligned collagen networks

In this work, a Langevin dynamics model of the diffusion of water in articular cartilage was developed. Numerical simulations of the translational dynamics of water molecules and their interaction with collagen fibers were used to study the quantitative relationship between the organization of the collagen fiber network and the diffusion tensor of water in model cartilage. Langevin dynamics was used to simulate water diffusion in both ordered and partially disordered cartilage models. In addition, an analytical approach was developed to estimate the diffusion tensor for a network comprising a given distribution of fiber orientations. The key findings are that (1) an approximately linear relationship was observed between collagen volume fraction and the fractional anisotropy of the diffusion tensor in fiber networks of a given degree of alignment, (2) for any given fiber volume fraction, fractional anisotropy follows a fiber alignment dependency similar to the square of the second Legendre polynomial of cos(θ), with the minimum anisotropy occurring at approximately the magic angle (θMA), and (3) a decrease in the principal eigenvalue and an increase in the transverse eigenvalues is observed as the fiber orientation angle θ progresses from 0◦ to 90◦. The corresponding diffusion ellipsoids are prolate for θ < θMA, spherical for θ ≈ θMA, and oblate for θ > θMA. Expansion of the model to include discrimination between the combined effects of alignment disorder and collagen fiber volume fraction on the diffusion tensor is discussed.

Does occupational aggregation matter in determining the size of the gender wage gap?

The traditional decomposition of the gender wage gap distinguishes between a component attributable to gender differences in productivity-related characteristics and a residual component that is often taken as a measure of discrimination. This study of data from the 1989 Canadian Labour Market Activity Survey shows that when occupation is treated as a productivity-related characteristic, the proportion of the gender wage gap labeled explained increases with the number of occupational classifications distinguished. However, on the basis of evidence that occupational differences reflect the presence of barriers faced by women attempting to enter male-dominated occupations, the authors conclude that occupation should not be treated as a productivity-related characteristic; and in a decomposition of the gender wage gap that treats occupation as endogenously determined, they find that the level of occupational aggregation has little effect on the size of the "explained" component of the gap.

Effect of body position on measurements of diffusion capacity after exercise

Background--Pulmonary diffusing capacity for carbon monoxide (Dlco), alveolar capillary membrane diffusing capacity (Dm), and pulmonary capillary blood volume (Vc) are all significantly reduced after exercise. Objective--To investigate whether measurement position affects this impaired gas transfer. Methods--Before and one, two, and four hours after incremental cycle ergometer exercise to fatigue, single breath Dlco, Dm, and Vc measurements were obtained in 10 healthy men in a randomly assigned supine and upright seated position. Results--After exercise, Dlco, Dm, and Vc were significantly depressed compared with baseline in both positions. The supine position produced significantly higher values over time for Dlco (5.22 (0.13) v 4.66 (0.15) ml/min/mm Hg/l, p = 0.022) and Dm (6.78 (0.19) v 6.03 (0.19) ml/min/mm Hg/l, p = 0.016), but there was no significant position effect for Vc. There was a similar pattern of change over time for Dlco, Dm, and Vc in the two positions. Conclusions--The change in Dlco after exercise appears to be primarily due to a decrease in Vc. Although the mechanism for the reduction in Vc cannot be determined from these data, passive relocation of blood to the periphery as the result of gravity can be discounted, suggesting that active vasoconstriction of the pulmonary vasculature and/or peripheral vasodilatation is occurring after exercise.

Measurement constructs to explore innovation diffusion in construction

Although the drivers of innovation have been studied extensively in construction, greater attention is required on how innovation diffusion can be effectively assessed within this complex and interdependent project-based industry. The authors draw on a highly cited innovation diffusion model by Rogers (2006) and develop a tailored conceptual framework to guide future empirical work aimed at assessing innovation diffusion in construction. The conceptual framework developed and discussed in this paper supports a five-stage process model of innovation diffusion namely: 1) knowledge and idea generation, 2) persuasion and evaluation; 3) decision to adopt, 4) integration and implementation, and 5) confirmation. As its theoretical contribution, this paper proposes three critical measurements constructs which can be used to assess the effectiveness of the diffusion process. These measurement constructs comprise: 1) nature and introduction of an innovative idea, 2) organizational capacity to acquire, assimilate, transform and exploit an innovation, and 3) rates of innovation facilitation and adoption. The constructs are interpreted in the project-based context of the construction industry, extending the contribution of general management theorists. Research planned by the authors will test the validity and reliability of the constructs developed in this paper.

A note on a reaction-diffusion model describing the bone morphogen protein gradient in Drosophila embryonic patterning

In this article, we consider the Eldar model [3] from embryology in which a bone morphogenic protein, a short gastrulation protein, and their compound react and diffuse. We carry out a perturbation analysis in the limit of small diffusivity of the bone morphogenic protein. This analysis establishes conditions under which some elementary results of [3] are valid.

Near-Wall Turbulent Pressure Diffusion Modeling and Influence in Three-Dimensional Secondary Flows

The purpose of this paper is to develop a second-moment closure with a near-wall turbulent pressure diffusion model for three-dimensional complex flows, and to evaluate the influence of the turbulent diffusion term on the prediction of detached and secondary flows. A complete turbulent diffusion model including a near-wall turbulent pressure diffusion closure for the slow part was developed based on the tensorial form of Lumley and included in a re-calibrated wall-normal-free Reynolds-stress model developed by Gerolymos and Vallet. The proposed model was validated against several one-, two, and three-dimensional complex flows.