A comparison of finite difference and finite volume methods for solving the space fractional advection-dispersion equation with variable coefficients

Transport processes within heterogeneous media may exhibit non-classical diffusion or dispersion; that is, not adequately described by the classical theory of Brownian motion and Fick's law. We consider a space fractional advection-dispersion equation based on a fractional Fick's law. The equation involves the Riemann-Liouville fractional derivative which arises from assuming that particles may make large jumps. Finite difference methods for solving this equation have been proposed by Meerschaert and Tadjeran. In the variable coefficient case, the product rule is first applied, and then the Riemann-Liouville fractional derivatives are discretised using standard and shifted Grunwald formulas, depending on the fractional order. In this work, we consider a finite volume method that deals directly with the equation in conservative form. Fractionally-shifted Grunwald formulas are used to discretise the fractional derivatives at control volume faces. We compare the two methods for several case studies from the literature, highlighting the convenience of the finite volume approach.

Comparison between analytical and 3D finite element solutions for borehole stresses in anisotropic elastic rock

We present a rigorous validation of the analytical Amadei solution for the stress concentration around an arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients b11 and b55 are not equal. It is shown from theoretical considerations and published experimental data that the b11 and b55 are not equal for realistic rocks. Second, we develop a 3D finite element elastic model within a hybrid analytical–numerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic, transverse isotropic and orthorhombic symmetries. It is concluded that the analytical Amadei solution is valid with no restriction on the borehole orientation or the symmetry of the elastic anisotropy.

Mathematical modelling of controlled drug release from polymer micro-spheres : incorporating the effects of swelling, diffusion and dissolution via moving boundary problems

Controlled drug delivery is a key topic in modern pharmacotherapy, where controlled drug delivery devices are required to prolong the period of release, maintain a constant release rate, or release the drug with a predetermined release profile. In the pharmaceutical industry, the development process of a controlled drug delivery device may be facilitated enormously by the mathematical modelling of drug release mechanisms, directly decreasing the number of necessary experiments. Such mathematical modelling is difficult because several mechanisms are involved during the drug release process. The main drug release mechanisms of a controlled release device are based on the device’s physiochemical properties, and include diffusion, swelling and erosion. In this thesis, four controlled drug delivery models are investigated. These four models selectively involve the solvent penetration into the polymeric device, the swelling of the polymer, the polymer erosion and the drug diffusion out of the device but all share two common key features. The first is that the solvent penetration into the polymer causes the transition of the polymer from a glassy state into a rubbery state. The interface between the two states of the polymer is modelled as a moving boundary and the speed of this interface is governed by a kinetic law. The second feature is that drug diffusion only happens in the rubbery region of the polymer, with a nonlinear diffusion coefficient which is dependent on the concentration of solvent. These models are analysed by using both formal asymptotics and numerical computation, where front-fixing methods and the method of lines with finite difference approximations are used to solve these models numerically. This numerical scheme is conservative, accurate and easily implemented to the moving boundary problems and is thoroughly explained in Section 3.2. From the small time asymptotic analysis in Sections 5.3.1, 6.3.1 and 7.2.1, these models exhibit the non-Fickian behaviour referred to as Case II diffusion, and an initial constant rate of drug release which is appealing to the pharmaceutical industry because this indicates zeroorder release. The numerical results of the models qualitatively confirms the experimental behaviour identified in the literature. The knowledge obtained from investigating these models can help to develop more complex multi-layered drug delivery devices in order to achieve sophisticated drug release profiles. A multi-layer matrix tablet, which consists of a number of polymer layers designed to provide sustainable and constant drug release or bimodal drug release, is also discussed in this research. The moving boundary problem describing the solvent penetration into the polymer also arises in melting and freezing problems which have been modelled as the classical onephase Stefan problem. The classical one-phase Stefan problem has unrealistic singularities existed in the problem at the complete melting time. Hence we investigate the effect of including the kinetic undercooling to the melting problem and this problem is called the one-phase Stefan problem with kinetic undercooling. Interestingly we discover the unrealistic singularities existed in the classical one-phase Stefan problem at the complete melting time are regularised and also find out the small time behaviour of the one-phase Stefan problem with kinetic undercooling is different to the classical one-phase Stefan problem from the small time asymptotic analysis in Section 3.3. In the case of melting very small particles, it is known that surface tension effects are important. The effect of including the surface tension to the melting problem for nanoparticles (no kinetic undercooling) has been investigated in the past, however the one-phase Stefan problem with surface tension exhibits finite-time blow-up. Therefore we investigate the effect of including both the surface tension and kinetic undercooling to the melting problem for nanoparticles and find out the the solution continues to exist until complete melting. The investigation of including kinetic undercooling and surface tension to the melting problems reveals more insight into the regularisations of unphysical singularities in the classical one-phase Stefan problem. This investigation gives a better understanding of melting a particle, and contributes to the current body of knowledge related to melting and freezing due to heat conduction.

Cold-formed steel columns subject to local buckling at elevated temperatures

Fire safety design of building structures has received greater attention in recent times due to continuing losses of properties and lives in fires. However, the structural behaviour of thin-walled cold-formed steel columns under fire conditions is not well understood despite the increasing use of light gauge steels in building construction. Cold-formed steel columns are often subject to local buckling effects. Therefore a series of laboratory tests of lipped and unlipped channel columns made of varying steel thicknesses and grades was undertaken at uniform elevated temperatures up to 700°C under steady state conditions. Finite element models of the tested columns were also developed, and their elastic buckling and nonlinear analysis results were compared with test results at elevated temperatures. Effects of the degradation of mechanical properties of steel with temperature were included in the finite element analyses. The use of accurately measured yield stress, elasticity modulus and stress-strain curves at elevated temperatures provided a good comparison of the ultimate loads and load-deflection curves from tests and finite element analyses. The commonly used effective width design rules and the direct strength method at ambient temperature were then used to predict the ultimate loads at elevated temperatures by using the reduced mechanical properties. By comparing these predicted ultimate loads with those from tests and finite element analyses, the accuracy of using this design approach was evaluated.

Review of current fire design rules for cold-formed steel wall systems

Light gauge steel frame wall systems are commonly used in industrial and commercial buildings, and there is a need for simple fire design rules to predict their load capacities and fire resistance ratings. During fire events, the light gauge steel frame wall studs are subjected to non-uniform temperature distributions that cause thermal bowing, neutral axis shift and magnification effects and thus resulting in a combined axial compression and bending action on the studs. In this research, a series of full-scale fire tests was conducted first to evaluate the performance of light gauge steel frame wall systems with eight different wall configurations under standard fire conditions. Finite element models of light gauge steel frame walls were then developed, analysed under transient and steady-state conditions and validated using full-scale fire tests. Using the results from fire tests and finite element analyses, a detailed investigation was undertaken into the prediction of axial compression strength and failure times of light gauge steel frame wall studs in standard fires using the available fire design rules based on Australian, American and European standards. The results from both fire tests and finite element analyses were used to investigate the ability of these fire design rules to include the complex effects of non-uniform temperature distributions and their accuracy in predicting the axial compression strength of wall studs and the failure times. Suitable modifications were then proposed to the fire design rules. This article presents the details of this investigation on the fire design rules of light gauge steel frame walls and the results.

Blocking layer optimisation of poly(3-hexylthiopene)based solid state dye sensitized solar cells

The optimisation study of the fabrication of a compact TiO2 blocking layer (via Spray Pyrolysis Deposition) for poly (3-hexylthiopene) (P3HT) for Solid State Dye Sensitized Solar Cells (SDSCs) is reported. We used a novel spray TiO2 precursor solution composition obtained by adding acetylacetone to a conventional formulation (Diisopropoxytitanium bis (acetylacetonate) in ethanol). By Scanning Electron Microscopy a TiO2 layer with compact morphology and thickness of around 100 nmis shown. Through a Tafel plot analysis an enhancement of the device diode-like behaviour induced by the acetylacetone blocking layer respect to the conventional one is observed. Significantly, the device fabricatedwith the acetylacetone blocking layer shows an overall increment of the cell performance with respect to the cellwith the conventional one (DJsc/Jsc = +13.8%, DFF/FF = +39.7%, DPCE/PCE = +55.6%). A conversion efficiency optimumis found for 15 successive spray cycles where the diode-like behaviour of the acetylacetone blocking layer is more effective. Over three batches of cells (fabricated with P3HT and dye D35) an average conversion efficiency value of 3.9% (under a class A sun simulator with 1 sun A.M. 1.5 illumination conditions) was measured. From the best cell we fabricated a conversion efficiency value of 4.5% was extracted. This represents a significant increment with respect to previously reported values for P3HT/dye D35 based SDSCs.

Political networks on twitter : tweeting the Queensland state election

This paper examines patterns of political activity and campaigning on Twitter in the context of the 2012 election in the Australian state of Queensland. Social media have been a visible component of political campaigning in Australia at least since the 2007 federal election, with Twitter, in particular, rising to greater prominence in the 2010 federal election. At state level, however, they have remained comparatively less important thus far. In this paper, we track uses of Twitter in the Queensland campaign from its unofficial start in February through to the election day of 24 March 2012. We both examine the overall patterns of activity in the hash tag #qldvotes, and track specific interactions between politicians and other users by following some 80 Twitter accounts of sitting members of parliament and alternative candidates. Such analysis provides new insights into the different approaches to social media campaigning which were embraced by specific candidates and party organisations, as well as an indication of the relative importance of social media activities, at present, for state-level election campaigns.

Final-state effects in neutron Compton scattering measurements on zirconium deuteride and beryllium

We report inelastic neutron scattering measurements of the neutron Compton profile, J(y), for Be and for D in polycrystalline ZrD2 over a range of momentum transfers, q between 27 and 178 °A−1. The measurements were performed using the inverse geometry spectrometer eVS which is situated at the UK pulsed spallation neutron source ISIS. We have investigated deviations from impulse approximation (IA) scattering which are generically referred to as final state effects (FSEs) using a method described by Sears. This method allows both the magnitude and the q dependence of the FSE to be studied. Analysis of the measured data was compared with analysis of numerical simulations based on the harmonic approximation and good agreement was found for both ZrD2 and Be. Finally we have shown how (∇2V), where V is the interatomic potential, can be extracted from the antisymmetric component of J(y).

Half metallicity in a zigzag double-walled nanotube nanodot : an Ab initio prediction

Ab initio density functional calculations were performed to study finite-length zigzag (7, 0) @ (16, 0) double-walled carbon nanotubes (DWCNTs) with H-termination at the open ends. We find that such a DWCNT nanodot displays a very large magnetic moment at the zigzag edges and the ground state displays symmetric anti-ferromagnetic coupling. When an external electric field is applied along the direction of tube axis, a gap is opened for one spin channel, whereas another spin channel remains metallic, i.e. half metallicity occurs. Our results suggest an important new avenue for the development of CNT-based spintronic materials with enhanced properties.

Effects of free flow state highway traffic parameters on the performance of amplify and forward cooperative multi-hop vehicular networks

In this paper we analyse the effects of highway traffic flow parameters like vehicle arrival rate and density on the performance of Amplify and Forward (AF) cooperative vehicular networks along a multi-lane highway under free flow state. We derive analytical expressions for connectivity performance and verify them with Monte-Carlo simulations. When AF cooperative relaying is employed together with Maximum Ratio Combining (MRC) at the receivers the average route error rate shows 10-20 fold improvement compared to direct communication. A 4-8 fold increase in maximum number of traversable hops can also be observed at different vehicle densities when AF cooperative communication is used to strengthen communication routes. However the theorical upper bound of maximum number of hops promises higher performance gains.

Energy efficient virtual machine placement in data centers

Electricity cost has become a major expense for running data centers and server consolidation using virtualization technology has been used as an important technology to improve the energy efficiency of data centers. In this research, a genetic algorithm and a simulation-annealing algorithm are proposed for the static virtual machine placement problem that considers the energy consumption in both the servers and the communication network, and a trading algorithm is proposed for dynamic virtual machine placement. Experimental results have shown that the proposed methods are more energy efficient than existing solutions.

A bondgraph approach to formation control using relative state measurements

In this paper we apply port-Hamiltonian theory with the bondgraph modelling approach to the problem of formation control using partial measurements of relative positions. We present a control design that drives a group of vehicles to a desired formation without requiring inter-vehicle communications or global position and velocity measurements to be available. Our generic approach is applicable to any form of relative measurement between vehicles, but we specifically consider the important cases of relative bearings and relative distances. In the case of bearings, our theory closely relates to the field of image-based visual servo (IBVS) control. We present simulation results to support the developed theory.

Carbon dioxide reforming of methane to produce synthesis gas over metal-supported catalysts: State of the art

Carbon dioxide reforming of methane produces synthesis gas with a low hydrogen to carbon monoxide ratio, which is desirable for many industrial synthesis processes. This reaction also has very important environmental implications since both methane and carbon dioxide contribute to the greenhouse effect. Converting these gases into a valuable feedstock may significantly reduce the atmospheric emissions of CO2 and CH4. In this paper, we present a comprehensive review on the thermodynamics, catalyst selection and activity, reaction mechanism, and kinetics of this important reaction. Recently, research has centered on the development of catalysts and the feasible applications of this reaction in industry. Group VIII metals supported on oxides are found to be effective for this reason. However, carbon deposition causing catalyst deactivation is the major problem inhibiting the industrial application of the CO2/CH4 reaction. Ni-based catalysts impregnated on certain supports show carbon-free operation and thus attract much attention. To develop an effective catalyst for CO2 reforming of CH4 and accelerate the commercial application of the reaction, the following are identified to be the most important areas for future work: (1) selection of metal and support and studying the effect of their interaction on catalyst activity; (2) the effect of different promoter on catalyst activity; (3) the reaction mechanism and kinetics; and (4) pilot reactor performance and scale-up operation.