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.

Globally networked public spheres? The Australian media reaction to WikiLeaks

The global release of 250,000 US Embassy diplomatic cables to selected media sites worldwide through the WikiLeaks website, was arguably the major global media event of 2010. As well as the implications of the content of the cables for international politics and diplomacy, the actions of WikiLeaks and its controversial editor-in-chief, the Australian Julian Assange, bring together a range of arguments about how the media, news and journalism are being transformed in the 21st century. This paper will focus on the reactions of Australian online news media sites to the release of the diplomatic cables by WikiLeaks, including both the online sites of established news outlets such as The Australian, Sydney Morning Herald and The Age, the ABC’s The Drum site, and online-only sites such as Crikey, New Matilda and On Line Opinion. The study focuses on opinion and commentary rather than straight news reportage, and analysis is framed around three issues: WikiLeaks and international diplomacy; implications of WikiLeaks for journalism; and WikiLeaks and democracy, including debates about the organisation and the ethics of its own practice. It also whether a “WikiLeaks Effect” has wider implications for how journalism is conducted in the future, particularly the method of ‘redaction’ of large amounts of computational data.

Investigation of phenol degradation : True reaction kinetics on fixed film titanium dioxide photocatalyst

This study of photocatalytic oxidation of phenol over titanium dioxide films presents a method for the evaluation of true reaction kinetics. A flat plate reactor was designed for the specific purpose of investigating the influence of various reaction parameters, specifically photocatalytic film thickness, solution flow rate (1–8 l min−1), phenol concentration (20, 40 and 80 ppm), and irradiation intensity (70.6, 57.9, 37.1and 20.4 W m−2), in order to further understand their impact on the reaction kinetics. Special attention was given to the mass transfer phenomena and the influence of film thickness. The kinetics of phenol degradation were investigated with different irradiation levels and initial pollutant concentration. Photocatalytic degradation experiments were performed to evaluate the influence of mass transfer on the reaction and, in addition, the benzoic acid method was applied for the evaluation of mass transfer coefficient. For this study the reactor was modelled as a batch-recycle reactor. A system of equations that accounts for irradiation, mass transfer and reaction rate was developed to describe the photocatalytic process, to fit the experimental data and to obtain kinetic parameters. The rate of phenol photocatalytic oxidation was described by a Langmuir–Hinshelwood type law that included competitive adsorption and degradation of phenol and its by-products. The by-products were modelled through their additive effect on the solution total organic carbon.

Thermal stability of crandallite CaAl3(PO4)2(OH)5•(H2O) – a ‘cave’ mineral from the Jenolan Caves

Thermogravimetry combined with evolved gas mass spectrometry has been used to characterise the mineral crandallite CaAl3(PO4)2(OH)5•(H2O) and to ascertain the thermal stability of this ‘cave’ mineral. X-ray diffraction proves the presence of the mineral and identifies the products after thermal decomposition. The mineral crandallite is formed through the reaction of calcite with bat guano. Thermal analysis shows that the mineral starts to decompose through dehydration at low temperatures at around 139°C while dehydroxylation occurs over the temperature range 200 to 700°C with loss of OH units. The critical temperature for OH loss is around 416°C and above this temperature the mineral structure is altered. Some minor loss of carbonate impurity occurs at 788°C. This study shows the mineral is unstable above 139°C. This temperature is well above the temperature in caves, which have a maximum temperature of 15°C. A chemical reaction for the synthesis of crandallite is offered and the mechanism for the thermal decomposition is given.

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.

Reaction to point and counterpoint : bold adventure or missed opportunity

The draft of the first stage of the national curriculum has now been published. Its final form to be presented in December 2010 should be the centrepiece of Labor’s Educational Revolution. All the other aspects – personal computers, new school buildings, rebates for uniforms and even the MySchool report card – are marginal to the prescription of what is to be taught and learnt in schools. The seven authors in this journal’s Point and Counterpoint (Curriculum Perspectives, 30(1) 2010, pp.53-74) raise a number of both large and small issues in education as a whole, and in science education more particularly. Two of them (Groves and McGarry) make brief reference to earlier attempts to achieve national curriculum in Australia. Those writing from New Zealand and USA will be unaware of just how ambitious this project is for Australia - a bold and overdue educational adventure or a foolish political decision destined to failure, as happened in the later 1970s and the 1990s.

Cyclic imides and an open-chain amide carboxylic acid from the facile reaction of cis-cyclohexane-1,2-dicarboxylic anhydride with the isomeric monofluoroanilines

The structures of the cyclic imides cis-2-(2-fluorophenyl)-3a,4,5,6,7,7a-hexahydroisoindole-1,3-dione, C14H14FNO2, (I), and cis-2-(4-fluorophenyl)-3a,4,5,6,7,7a-hexahydroisoindoline-1,3-dione, C14H14FNO2, (III), and the open-chain amide acid rac-cis-2-[(3-fluorophenyl)carbamoyl]cyclohexane-1-carboxylic acid, C14H16FNO3, (II), are reported. Cyclic imides (I) and (III) are conformationally similar, with comparable ring rotations about the imide N-Car bond [the dihedral angles between the benzene ring and the five-membered isoindole ring are 55.40 (8)° for (I) and 51.83 (7)° for (III)]. There are no formal intermolecular hydrogen bonds involved in the crystal packing of either (I) or (III). With the acid (II), in which the meta-related F-atom substituent is rotationally disordered (0.784:0.216), the amide group lies slightly out of the benzene plane [the interplanar dihedral angle is 39.7 (1)°]. Intermolecular amide-carboxyl N-HO hydrogen-bonding interactions between centrosymmetrically related molecules form stacks extending down b, and these are linked across c by carboxyl-amide O-HO hydrogen bonds, giving two-dimensional layered structures which lie in the (011) plane. The structures reported here represent examples of compounds analogous to the phthalimides or phthalanilic acids and have little precedence in the crystallographic literature.

The surface-structure sensitivity of dioxygen activation in the anatase-photocatalyzed oxidation reaction

The reaction pathways by which oxygen is incorporated into the substrate in the photocatalytic oxidation of terephthalic acid (TPTA) are vastly different on {001} and {101} facets of an anatase single crystal. This was established by controlling the percentage of {101} and {001} facets, isotopically tracing the origins of the hydroxy group, and studying dioxygen consumption and variance in the concentration of hydroxylation intermediate.

Hydrogen-bonding in cyclic imides and amide carboxylic acid derivatives from the facile reaction of cis-cyclohexane-1,2-carboxylic anhydride with o- and p-anisidine and m- and p-aminobenzoic acids

The structures of the open chain amide carboxylic acid rac-cis-[2-(2-methoxyphenyl)carbamoyl]cyclohexane-1-carboxylic acid, C15H19NO4, (I) and the cyclic imides rac-cis-2-(4-methoxyphenyl)-3a,4,5,6,7,7-hexahydroisoindole-1,3-dione,C15H17NO3, (II), chiral cis-2-(3-carboxyphenyl)-3a,4,5,6,7,7a-hexahydroisoindole-1,3-dione, C15H15NO4,(III) and rac-cis-2-(4-carboxyphenyl)- 3a,4,5,6,7,7a-hexahydroisoindole-1,3-dione monohydrate, C15H15NO4. H2O) (IV), are reported. In the amide acid (I), the phenylcarbamoyl group is essentially planar [maximum deviation from the least-squares plane = 0.060(1)Ang. for the amide O atom], the molecules form discrete centrosymmetric dimers through intermolecular cyclic carboxy-carboxy O-H...O hydrogen-bonding interactions [graph set notation R2/2(8)]. The cyclic imides (II)--(IV) are conformationally similar, with comparable phenyl ring rotations about the imide N-C(aromatic) bond [dihedral angles between the benzene and isoindole rings = 51.55(7)deg. in (II), 59.22(12)deg. in (III) and 51.99(14)deg. in (IV). Unlike (II) in which only weak intermolecular C-H...O(imide) hydrogen bonding is present, the crystal packing of imides (III) and (IV) shows strong intermolecular carboxylic acid O-H...O hydrogen-bonding associations. With (III), these involve imide O-atom acceptors, giving one-dimensional zigzag chains [graph set C(9)], while with the monohydrate (IV), the hydrogen bond involves the partially disordered water molecule which also bridges molecules through both imide and carboxyl O-atom acceptors in a cyclic R4/4(12) association, giving a two-dimensional sheet structure. The structures reported here expand the structural data base for compounds of this series formed from the facile reaction of cis-cyclohexane-1,2-dicarboxylic anhydride with substituted anilines, in which there is a much larger incidence of cyclic imides compared to amide carboxylic acids.

A least squares based finite volume method for the Cahn-Hilliard and Cahn-Hilliard-reaction equations

A vertex-centred finite volume method (FVM) for the Cahn-Hilliard (CH) and recently proposed Cahn-Hilliard-reaction (CHR) equations is presented. Information at control volume faces is computed using a high-order least-squares approach based on Taylor series approximations. This least-squares problem explicitly includes the variational boundary condition (VBC) that ensures that the discrete equations satisfy all of the boundary conditions. We use this approach to solve the CH and CHR equations in one and two dimensions and show that our scheme satisfies the VBC to at least second order. For the CH equation we show evidence of conservative, gradient stable solutions, however for the CHR equation, strict gradient-stability is more challenging to achieve.

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.

The persistence of phase-separation in LiFePO4 with two-dimensional Li+ transport : the Cahn-Hilliard-reaction equation and the role of defects

We examine the solution of the two-dimensional Cahn-Hilliard-reaction (CHR) equation in the xy plane as a model of Li+ intercalation into LiFePO4 material. We validate our numerical solution against the solution of the depth-averaged equation, which has been used to model intercalation in the limit of highly orthotropic diffusivity and gradient penalty tensors. We then examine the phase-change behaviour in the full CHR system as these parameters become more isotropic, and find that as the Li+ diffusivity is increased in the x direction, phase separation persists at high currents, even in small crystals with averaged coherency strain included. The resulting voltage curves decrease monotonically, which has previously been considered a hallmark of crystals that fill homogeneously.