Time-dependent fractional advection–diffusion equations by an implicit MLS meshless method

Recently, many new applications in engineering and science are governed by a series of fractional partial differential equations (FPDEs). Unlike the normal partial differential equations (PDEs), the differential order in a FPDE is with a fractional order, which will lead to new challenges for numerical simulation, because most existing numerical simulation techniques are developed for the PDE with an integer differential order. The current dominant numerical method for FPDEs is Finite Difference Method (FDM), which is usually difficult to handle a complex problem domain, and also hard to use irregular nodal distribution. This paper aims to develop an implicit meshless approach based on the moving least squares (MLS) approximation for numerical simulation of fractional advection-diffusion equations (FADE), which is a typical FPDE. The discrete system of equations is obtained by using the MLS meshless shape functions and the meshless strong-forms. The stability and convergence related to the time discretization of this approach are then discussed and theoretically proven. Several numerical examples with different problem domains and different nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of the FADE.

An implicit RBF meshless approach for time fractional diffusion equations

This paper aims to develop an implicit meshless approach based on the radial basis function (RBF) for numerical simulation of time fractional diffusion equations. The meshless RBF interpolation is firstly briefed. The discrete equations for two-dimensional time fractional diffusion equation (FDE) are obtained by using the meshless RBF shape functions and the strong-forms of the time FDE. The stability and convergence of this meshless approach are discussed and theoretically proven. Numerical examples with different problem domains and different nodal distributions are studied to validate and investigate accuracy and efficiency of the newly developed meshless approach. It has proven that the present meshless formulation is very effective for modeling and simulation of fractional differential equations.

Stochastic modelling of T cell homeostasis for two competing clonotypes via the master equation

Stochastic models for competing clonotypes of T cells by multivariate, continuous-time, discrete state, Markov processes have been proposed in the literature by Stirk, Molina-París and van den Berg (2008). A stochastic modelling framework is important because of rare events associated with small populations of some critical cell types. Usually, computational methods for these problems employ a trajectory-based approach, based on Monte Carlo simulation. This is partly because the complementary, probability density function (PDF) approaches can be expensive but here we describe some efficient PDF approaches by directly solving the governing equations, known as the Master Equation. These computations are made very efficient through an approximation of the state space by the Finite State Projection and through the use of Krylov subspace methods when evolving the matrix exponential. These computational methods allow us to explore the evolution of the PDFs associated with these stochastic models, and bimodal distributions arise in some parameter regimes. Time-dependent propensities naturally arise in immunological processes due to, for example, age-dependent effects. Incorporating time-dependent propensities into the framework of the Master Equation significantly complicates the corresponding computational methods but here we describe an efficient approach via Magnus formulas. Although this contribution focuses on the example of competing clonotypes, the general principles are relevant to multivariate Markov processes and provide fundamental techniques for computational immunology.

Supplement : Efficient weak second-order stochastic Runge-Kutta methods for non-commutative Stratonvich stochastic differential equations

This paper gives a modification of a class of stochastic Runge–Kutta methods proposed in a paper by Komori (2007). The slight modification can reduce the computational costs of the methods significantly.

Preparing for Research Assessment: Co-evolution and Gamesmanship

In this paper we adopt a complex systems perspective to examine the perturbations caused by the introduction of the Research Quality Framework (RQF) at a research-intensive Australian university. This case is instructive as it 1) presents a Federal policy initiative that attempted to fundamentally alter the recognition and reward mechanism within a regulated funding environment, 2) analyses the strategies of an institution and its research groups as they sought to not only comply with the implementation of the RQF but to maximise their outcome,and 3) it reveals the ways that some actors used this perturbation to advance their own interests. In short, this case represents an instrumental study into the dynamics of how information systems, organisations, and individuals co-evolve in practice as they seek to navigate a complex problem scenario.

Novel numerical methods for time-space fractional reaction diffusion equations in two dimensions

We consider time-space fractional reaction diffusion equations in two dimensions. This equation is obtained from the standard reaction diffusion equation by replacing the first order time derivative with the Caputo fractional derivative, and the second order space derivatives with the fractional Laplacian. Using the matrix transfer technique proposed by Ilic, Liu, Turner and Anh [Fract. Calc. Appl. Anal., 9:333--349, 2006] and the numerical solution strategy used by Yang, Turner, Liu, and Ilic [SIAM J. Scientific Computing, 33:1159--1180, 2011], the solution of the time-space fractional reaction diffusion equations in two dimensions can be written in terms of a matrix function vector product $f(A)b$ at each time step, where $A$ is an approximate matrix representation of the standard Laplacian. We use the finite volume method over unstructured triangular meshes to generate the matrix $A$, which is therefore non-symmetric. However, the standard Lanczos method for approximating $f(A)b$ requires that $A$ is symmetric. We propose a simple and novel transformation in which the standard Lanczos method is still applicable to find $f(A)b$, despite the loss of symmetry. Numerical results are presented to verify the accuracy and efficiency of our newly proposed numerical solution strategy.

Cerclage, evolution and potential of a Cinderella Technology : an overview with reference to periprosthetic fractures

Periprosthetic fractures are increasingly frequent. The fracture may be located over the shaft of the prosthesis, at its tip or below (21). The treatment of explosion fractures is difficult because the shaft blocks the application of implants, like screws, which need to penetrate the medullary cavity. The cerclage, as a simple periosteal loop, made of wire or more recently cable, does not only avoid the medullary cavity. Its centripetal mode of action is well suited for reducing and maintaining radially displaced fractures. Furthermore, the cerclage lends itself well for minimally invasive internal fixation. New insight challenges the disrepute of which the cerclage technology suffered for decades. The outcome of cerclage fixation benefits from an improved understanding of its technology, mechano-biology and periosteal blood supply. Preconceived and generally accepted opinions like "strangulation of blood supply" need to be re-examined. Recent mechanical evaluations (22) demonstrate that the wire application may be improved but cable is superior in hand- ling, maintenance of tension and strength. Beside the classical concepts of absolute and relative stability a defined stability condition needs consideration. It is typical for cerclage. Called "loose-lock stability" it specifies the situation where a loosened implant allows first unimpeded displacement changing abruptly into a locked fixation preventing further dislocation.

Iron biogeochemistry and associated greenhouse gas evolution in a forested subtropical Australian coastal catchment : Poona Creek, Southeast Queensland

Iron (Fe) is the fourth most abundant element in the Earth’s crust. Excess Fe mobilization from terrestrial into aquatic systems is of concern for deterioration of water quality via biofouling and nuisance algal blooms in coastal and marine systems. Substantial Fe dissolution and transport involve alternate Fe(II) oxidation followed by Fe(III) reduction, with a diversity of Bacteria and Archaea acting as the key catalyst. Microbially-mediated Fe cycling is of global significance with regard to cycles of carbon (C), sulfur (S) and manganese (Mn). However, knowledge regarding microbial Fe cycling in circumneutral-pH habitats that prevail on Earth has been lacking until recently. In particular, little is known regarding microbial function in Fe cycling and associated Fe mobilization and greenhouse (CO2 and CH4, GHG) evolution in subtropical Australian coastal systems where microbial response to ambient variations such as seasonal flooding and land use changes is of concern. Using the plantation-forested Poona Creek catchment on the Fraser Coast of Southeast Queensland (SEQ), this research aimed to 1) study Fe cycling-associated bacterial populations in diverse terrestrial and aquatic habitats of a representative subtropical coastal circumneutral-pH (4–7) ecosystem; and 2) assess potential impacts of Pinus plantation forestry practices on microbially-mediated Fe mobilization, organic C mineralization and associated GHG evolution in coastal SEQ. A combination of wet-chemical extraction, undisturbed core microcosm, laboratory bacterial cultivation, microscopy and 16S rRNA-based molecular phylogenetic techniques were employed. The study area consisted primarily of loamy sands, with low organic C and dissolved nutrients. Total reactive Fe was abundant and evenly distributed within soil 0–30 cm profiles. Organic complexation primarily controlled Fe bioavailability and forms in well-drained plantation soils and water-logged, native riparian soils, whereas tidal flushing exerted a strong “seawater effect” in estuarine locations and formed a large proportion of inorganic Fe(III) complexes. There was a lack of Fe(II) sources across the catchment terrestrial system. Mature, first-rotation plantation clear-felling and second-rotation replanting significantly decreased organic matter and poorly crystalline Fe in well-drained soils, although variations in labile soil organic C fractions (dissolved organic C, DOC; and microbial biomass C, MBC) were minor. Both well-drained plantation soils and water-logged, native-vegetation soils were inhabited by a variety of cultivable, chemotrophic bacterial populations capable of C, Fe, S and Mn metabolism via lithotrophic or heterotrophic, (micro)aerobic or anaerobic pathways. Neutrophilic Fe(III)-reducing bacteria (FeRB) were most abundant, followed by aerobic, heterotrophic bacteria (heterotrophic plate count, HPC). Despite an abundance of FeRB, cultivable Fe(II)-oxidizing bacteria (FeOB) were absent in associated soils. A lack of links between cultivable Fe, S or Mn bacterial densities and relevant chemical measurements (except for HPC correlated with DOC) was likely due to complex biogeochemical interactions. Neither did variations in cultivable bacterial densities correlate with plantation forestry practices, despite total cultivable bacterial densities being significantly lower in estuarine soils when compared with well-drained plantation soils and water-logged, riparian native-vegetation soils. Given that bacterial Fe(III) reduction is the primary mechanism of Fe oxide dissolution in soils upon saturation, associated Fe mobilization involved several abiotic and biological processes. Abiotic oxidation of dissolved Fe(II) by Mn appeared to control Fe transport and inhibit Fe dissolution from mature, first-rotation plantation soils post-saturation. Such an effect was not observed in clear-felled and replanted soils associated with low SOM and potentially low Mn reactivity. Associated GHG evolution post-saturation mainly involved variable CO2 emissions, with low, but consistently increasing CH4 effluxes in mature, first-rotation plantation soil only. In comparison, water-logged soils in the riparian native-vegetation buffer zone functioned as an important GHG source, with high potentials for Fe mobilization and GHG, particularly CH4 emissions in riparian loam soils associated with high clay and crystalline Fe fractions. Active Fe–C cycling was unlikely to occur in lower-catchment estuarine soils associated with low cultivable bacterial densities and GHG effluxes. As a key component of bacterial Fe cycling, neutrophilic FeOB widely occurred in diverse aquatic, but not terrestrial, habitats of the catchment study area. Stalked and sheathed FeOB resembling Gallionella and Leptothrix were limited to microbial mat material deposited in surface fresh waters associated with a circumneutral-pH seep, and clay-rich soil within riparian buffer zones. Unicellular, Sideroxydans-related FeOB (96% sequence identity) were ubiquitous in surface and subsurface freshwater environments, with highest abundance in estuary-adjacent shallow coastal groundwater water associated with redox transition. The abundance of dissolved C and Fe in the groundwater-dependent system was associated with high numbers of cultivable anaerobic, heterotrophic FeRB, microaerophilic, putatively lithotrophic FeOB and aerobic, heterotrophic bacteria. This research represents the first study of microbial Fe cycling in diverse circumneutral-pH environments (terrestrial–aquatic, freshwater–estuarine, surface–subsurface) of a subtropical coastal ecosystem. It also represents the first study of its kind in the southern hemisphere. This work highlights the significance of bacterial Fe(III) reduction in terrestrial, and bacterial Fe(II) oxidation in aquatic catchment Fe cycling. Results indicate the risk of promotion of Fe mobilization due to plantation clear-felling and replanting, and GHG emissions associated with seasonal water-logging. Additional significant outcomes were also achieved. The first direct evidence for multiple biomineralization patterns of neutrophilic, microaerophilic, unicellular FeOB was presented. A putatively pure culture, which represents the first cultivable neutrophilic FeOB from the southern hemisphere, was obtained as representative FeOB ubiquitous in diverse catchment aquatic habitats.