A numerical method for the fractional Fitzhugh–Nagumo monodomain model

A fractional FitzHugh–Nagumo monodomain model with zero Dirichlet boundary conditions is presented, generalising the standard monodomain model that describes the propagation of the electrical potential in heterogeneous cardiac tissue. The model consists of a coupled fractional Riesz space nonlinear reaction-diffusion model and a system of ordinary differential equations, describing the ionic fluxes as a function of the membrane potential. We solve this model by decoupling the space-fractional partial differential equation and the system of ordinary differential equations at each time step. Thus, this means treating the fractional Riesz space nonlinear reaction-diffusion model as if the nonlinear source term is only locally Lipschitz. The fractional Riesz space nonlinear reaction-diffusion model is solved using an implicit numerical method with the shifted Grunwald–Letnikov approximation, and the stability and convergence are discussed in detail in the context of the local Lipschitz property. Some numerical examples are given to show the consistency of our computational approach.

A reduction of diffusion in PVA Fricke hydrogels

A modification to the PVA-FX hydrogel whereby the chelating agent, xylenol orange, was partially bonded to the gelling agent, poly-vinyl alcohol, resulted in an 8% reduction in the post irradiation Fe3+ diffusion, adding approximately 1 hour to the useful timespan between irradiation and readout. This xylenol orange functionalised poly-vinyl alcohol hydrogel had an OD dose sensitivity of 0.014 Gy−1 and a diffusion rate of 0.133 mm2 h−1. As this partial bond yields only incremental improvement, it is proposed that more efficient methods of bonding xylenol orange to poly-vinyl alcohol be investigated to further reduce the diffusion in Fricke gels.

Green urbanism and diffusion issues

This paper addresses the research question, ‘What are the diffusion determinants for green urbanism innovations in Australia?’ This is a significant topic given the global movement towards green urbanism. The study reported here is based on desktop research that provides new insights through (1) synthesis of the latest research findings on green urbanism innovations and (2) interpretation of diffusion issues through our innovation system model. Although innovation determinants have been studied extensively overseas and in Australia, there is presently a gap in the literature when it comes to these determinants for green urbanism in Australia. The current paper fills this gap. Using a conceptual framework drawn from the innovation systems literature, this paper synthesises and interprets the literature to map the current state of green urbanism innovations in Australia and to analyse the drivers for, and obstacles to, their optimal diffusion. The results point to the importance of collaboration between project-based actors in the implementation of green urbanism. Education, training and regulation across the product system is also required to improve the cultural and technical context for implementation. The results are limited by their exploratory nature and future research is planned to quantify barriers to green urbanism.

The chain reaction

Synopsis and review of the Australian feature film The Chain Reaction, directed by Ian Barry.

Further development of discrete computational techniques for calculation of restricted diffusion propagators in porous media

Magnetic resonance is a well-established tool for structural characterisation of porous media. Features of pore-space morphology can be inferred from NMR diffusion-diffraction plots or the time-dependence of the apparent diffusion coefficient. Diffusion NMR signal attenuation can be computed from the restricted diffusion propagator, which describes the distribution of diffusing particles for a given starting position and diffusion time. We present two techniques for efficient evaluation of restricted diffusion propagators for use in NMR porous-media characterisation. The first is the Lattice Path Count (LPC). Its physical essence is that the restricted diffusion propagator connecting points A and B in time t is proportional to the number of distinct length-t paths from A to B. By using a discrete lattice, the number of such paths can be counted exactly. The second technique is the Markov transition matrix (MTM). The matrix represents the probabilities of jumps between every pair of lattice nodes within a single timestep. The propagator for an arbitrary diffusion time can be calculated as the appropriate matrix power. For periodic geometries, the transition matrix needs to be defined only for a single unit cell. This makes MTM ideally suited for periodic systems. Both LPC and MTM are closely related to existing computational techniques: LPC, to combinatorial techniques; and MTM, to the Fokker-Planck master equation. The relationship between LPC, MTM and other computational techniques is briefly discussed in the paper. Both LPC and MTM perform favourably compared to Monte Carlo sampling, yielding highly accurate and almost noiseless restricted diffusion propagators. Initial tests indicate that their computational performance is comparable to that of finite element methods. Both LPC and MTM can be applied to complicated pore-space geometries with no analytic solution. We discuss the new methods in the context of diffusion propagator calculation in porous materials and model biological tissues.

Finite volume methods for simulating anomalous transport

In this thesis a new approach for solving a certain class of anomalous diffusion equations was developed. The theory and algorithms arising from this work will pave the way for more efficient and more accurate solutions of these equations, with applications to science, health and industry. The method of finite volumes was applied to discretise the spatial derivatives, and this was shown to outperform existing methods in several key respects. The stability and convergence of the new method were rigorously established.

Effect of nano-silica on the mechanical and transport properties of lightweight concrete

This paper investigated the influence of nano-silica (NS) on the mechanical and transport properties of lightweight concrete (LWC). The resistance of LWC to water and chloride ions penetration was enhanced despite strength marginally increased. Water penetration depth, moisture sorptivity, chloride migration and diffusion coefficient was reduced by 23% and 49%, 23% and 10%, 5% and 0%, 22% and 12% compared to the two reference LWC mixes (pure cement and 60% slag blended cement), respectively with 1% NS. Such improvements were attributed to more compact microstructures because the micropore system was refined and the interface between aggregates and paste was enhanced.

A meshless method based on Point Interpolation Method (PIM) for the space fractional diffusion equation

This paper aims to develop a meshless approach based on the Point Interpolation Method (PIM) for numerical simulation of a space fractional diffusion equation. Two fully-discrete schemes for the one-dimensional space fractional diffusion equation are obtained by using the PIM and the strong-forms of the space diffusion equation. Numerical examples with different nodal distributions are studied to validate and investigate the accuracy and efficiency of the newly developed meshless approach.

A semi-alternating direction method for a 2-D fractional FitzHugh–Nagumo monodomain model on an approximate irregular domain

A FitzHugh-Nagumo monodomain model has been used to describe the propagation of the electrical potential in heterogeneous cardiac tissue. In this paper, we consider a two-dimensional fractional FitzHugh-Nagumo monodomain model on an irregular domain. The model consists of a coupled Riesz space fractional nonlinear reaction-diffusion model and an ordinary differential equation, describing the ionic fluxes as a function of the membrane potential. Secondly, we use a decoupling technique and focus on solving the Riesz space fractional nonlinear reaction-diffusion model. A novel spatially second-order accurate semi-implicit alternating direction method (SIADM) for this model on an approximate irregular domain is proposed. Thirdly, stability and convergence of the SIADM are proved. Finally, some numerical examples are given to support our theoretical analysis and these numerical techniques are employed to simulate a two-dimensional fractional Fitzhugh-Nagumo model on both an approximate circular and an approximate irregular domain.

A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients

In this paper, we derive a new nonlinear two-sided space-fractional diffusion equation with variable coefficients from the fractional Fick’s law. A semi-implicit difference method (SIDM) for this equation is proposed. The stability and convergence of the SIDM are discussed. For the implementation, we develop a fast accurate iterative method for the SIDM by decomposing the dense coefficient matrix into a combination of Toeplitz-like matrices. This fast iterative method significantly reduces the storage requirement of O(n2)O(n2) and computational cost of O(n3)O(n3) down to n and O(nlogn)O(nlogn), where n is the number of grid points. The method retains the same accuracy as the underlying SIDM solved with Gaussian elimination. Finally, some numerical results are shown to verify the accuracy and efficiency of the new method.