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.

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.

An analytical solution for diffusion and nonlinear uptake of oxygen in a spherical cell

We develop a new analytical solution for a reactive transport model that describes the steady-state distribution of oxygen subject to diffusive transport and nonlinear uptake in a sphere. This model was originally reported by Lin (Journal of Theoretical Biology, 1976 v60, pp449–457) to represent the distribution of oxygen inside a cell and has since been studied extensively by both the numerical analysis and formal analysis communities. Here we extend these previous studies by deriving an analytical solution to a generalized reaction-diffusion equation that encompasses Lin’s model as a particular case. We evaluate the solution for the parameter combinations presented by Lin and show that the new solutions are identical to a grid-independent numerical approximation.

An advanced meshless method for time fractional diffusion equation

Recently, because of the new developments in sustainable engineering and renewable energy, which are usually governed by a series of fractional partial differential equations (FPDEs), the numerical modelling and simulation for fractional calculus are attracting more and more attention from researchers. The current dominant numerical method for modeling FPDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings including difficulty in simulation of problems with the complex problem domain and in using irregularly distributed nodes. Because of its distinguished advantages, the meshless method has good potential in simulation of FPDEs. This paper aims to develop an implicit meshless collocation technique for FPDE. The discrete system of FPDEs is obtained by using the meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach are investigated theoretically and numerically. The numerical examples with regular and irregular 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 fractional partial differential equations.

Models of collective cell motion for cell populations with different aspect ratio : diffusion, proliferation and travelling waves

Continuum, partial differential equation models are often used to describe the collective motion of cell populations, with various types of motility represented by the choice of diffusion coefficient, and cell proliferation captured by the source terms. Previously, the choice of diffusion coefficient has been largely arbitrary, with the decision to choose a particular linear or nonlinear form generally based on calibration arguments rather than making any physical connection with the underlying individual-level properties of the cell motility mechanism. In this work we provide a new link between individual-level models, which account for important cell properties such as varying cell shape and volume exclusion, and population-level partial differential equation models. We work in an exclusion process framework, considering aligned, elongated cells that may occupy more than one lattice site, in order to represent populations of agents with different sizes. Three different idealizations of the individual-level mechanism are proposed, and these are connected to three different partial differential equations, each with a different diffusion coefficient; one linear, one nonlinear and degenerate and one nonlinear and nondegenerate. We test the ability of these three models to predict the population level response of a cell spreading problem for both proliferative and nonproliferative cases. We also explore the potential of our models to predict long time travelling wave invasion rates and extend our results to two dimensional spreading and invasion. Our results show that each model can accurately predict density data for nonproliferative systems, but that only one does so for proliferative systems. Hence great care must be taken to predict density data for with varying cell shape.

Use of a zeolite synthesised from alkali treated kaolin as a K fertiliser : glasshouse experiments on leaching and uptake of K by wheat plants in sandy soil

Zeolite N, a zeolite referred to in earlier publications as MesoLite, is made by caustic reaction of kaolin at temperatures between 80 °C and 95 °C. This material has a very high cation exchange capacity (CEC ≈ 500 meq/100 g). Soil column leaching experiments have shown that K-zeolite N additions greatly reduce leaching of NH4+ fertilisers but the agronomic effectiveness of the retained K+ and NH4+ is unknown. To measure the bioavailability of K in this zeolite, wheat was grown in a glasshouse with K-zeolite N as the K fertiliser in highly-leached and non-leached pots for four weeks and compared with a soluble K fertiliser (KCl). The plants grown in non-leached pots and fertilised with K-zeolite N were slightly larger than those grown with KCl. The elemental compositions in the plants were similar except for Si being significantly more concentrated in the plants supplied with K-zeolite N. Thus K-zeolite N may be an effective K-fertiliser. Plants grown in highly-leached pots were significantly smaller than those grown in non-leached pots. Plants grown in highly-leached pots were severely K deficient as half of the K from both KCl and K-zeolite N was leached from the pots within three days.

Phenological changes in six Australian subalpine plants in response to experimental warming and year-to-year variation

The likely phenological responses of plants to climate warming can be measured through experimental manipulation of field sites, but results are rarely validated against year-to-year changes in climate. Here, we describe the response of 1-5 years of experimental warming on phenology (budding, flowering and seed maturation) of six common subalpine plant species in the Australian Alps using the International Tundra Experiment (ITEX) protocol.2. Phenological changes in some species (particularly the forb Craspedia jamesii) were detected in experimental plots within a year of warming, whereas changes in most other species (the forb Erigeron bellidioides, the shrub Asterolasia trymalioides and the graminoids Carex breviculmis and Poa hiemata) did not develop until after 2-4 years; thus, there appears to be a cumulative effect of warming for some species across multiple years.3. There was evidence of changes in the length of the period between flowering and seed maturity in one species (P. hiemata) that led to a similar timing of seed maturation, suggesting compensation.4. Year-to-year variation in phenology was greater than variation between warmed and control plots and could be related to differences in thawing degree days (particularly, for E. bellidioides) due to earlier timing of budding and other events under warmer conditions. However, in Carex breviculmis, there was no association between phenology and temperature changes across years.5. These findings indicate that, although phenological changes occurred earlier in response to warming in all six species, some species showed buffered rather than immediate responses.6. Synthesis. Warming in ITEX open-top chambers in the Australian Alps produced earlier budding, flowering and seed set in several alpine species. Species also altered the timing of these events, particularly budding, in response to year-to-year temperature variation. Some species responded immediately, whereas in others the cumulative effects of warming across several years were required before a response was detected.

Diffusion determinants for passive building technologies

This summary is based on an international review of leading peer reviewed journals, in both technical and management fields. It draws on highly cited articles published between 2000 and 2009 to investigate the research question, "What are the diffusion determinants for passive building technologies in Australia?". Using a conceptual framework drawn from the innovation systems literature, this paper synthesises and interprets the literature to map the current state of passive building technologies in Australia and to analyse the drivers for, and obstacles to, their optimal diffusion. The paper concludes that the government has a key role to play through its influence over the specification of building codes.

Gaseous and particle emissions from an ethanol fumigated compression ignition engine

A 4-cylinder Ford 2701C test engine was used in this study to explore the impact of ethanol fumigation on gaseous and particle emission concentrations. The fumigation technique delivered vaporised ethanol into the intake manifold of the engine, using an injector, a pump and pressure regulator, a heat exchanger for vaporising ethanol and a separate fuel tank and lines. Gaseous (Nitric oxide (NO), Carbon monoxide (CO) and hydrocarbons (HC)) and particulate emissions (particle mass (PM2.5) and particle number) testing was conducted at intermediate speed (1700 rpm) using 4 load settings with ethanol substitution percentages ranging from 10-40 % (by energy). With ethanol fumigation, NO and PM2.5 emissions were reduced, whereas CO and HC emissions increased considerably and particle number emissions increased at most test settings. It was found that ethanol fumigation reduced the excess air factor for the engine and this led to increased emissions of CO and HC, but decreased emissions of NO. PM2.5 emissions were reduced with ethanol fumigation, as ethanol has a very low “sooting” tendency. This is due to the higher hydrogen-to-carbon ratio of this fuel, and also because ethanol does not contain aromatics, both of which are known soot precursors. The use of a diesel oxidation catalyst (as an after-treatment device) is recommended to achieve a reduction in the four pollutants that are currently regulated for compression ignition engines. The increase in particle number emissions with ethanol fumigation was due to the formation of volatile (organic) particles; consequently, using a diesel oxidation catalyst will also assist in reducing particle number emissions.