Near-Wall Turbulent Pressure Diffusion Modeling and Influence in Three-Dimensional Secondary Flows

The purpose of this paper is to develop a second-moment closure with a near-wall turbulent pressure diffusion model for three-dimensional complex flows, and to evaluate the influence of the turbulent diffusion term on the prediction of detached and secondary flows. A complete turbulent diffusion model including a near-wall turbulent pressure diffusion closure for the slow part was developed based on the tensorial form of Lumley and included in a re-calibrated wall-normal-free Reynolds-stress model developed by Gerolymos and Vallet. The proposed model was validated against several one-, two, and three-dimensional complex flows.

Air crimes and atmospheric justice

This chapter explores the political economy of air pollution. It draws on discourses of power, harm and violence to analyse air pollution within emerging frameworks of 'eco-crime' and atmospheric justice' (see Vanderheiden 2008; Walters 2010). In doing so, it identifies how green criminology continues to push new boundaries by engaging with issues of both global and local concern.

Efficient solution of two-sided nonlinear space-fractional diffusion equations using fast Poisson preconditioners

We develop a fast Poisson preconditioner for the efficient numerical solution of a class of two-sided nonlinear space fractional diffusion equations in one and two dimensions using the method of lines. Using the shifted Gr¨unwald finite difference formulas to approximate the two-sided(i.e. the left and right Riemann-Liouville) fractional derivatives, the resulting semi-discrete nonlinear systems have dense Jacobian matrices owing to the non-local property of fractional derivatives. We employ a modern initial value problem solver utilising backward differentiation formulas and Jacobian-free Newton-Krylov methods to solve these systems. For efficient performance of the Jacobianfree Newton-Krylov method it is essential to apply an effective preconditioner to accelerate the convergence of the linear iterative solver. The key contribution of our work is to generalise the fast Poisson preconditioner, widely used for integer-order diffusion equations, so that it applies to the two-sided space fractional diffusion equation. A number of numerical experiments are presented to demonstrate the effectiveness of the preconditioner and the overall solution strategy.

A banded preconditioner for the two-sided, nonlinear space-fractional diffusion equation

The method of lines is a standard method for advancing the solution of partial differential equations (PDEs) in time. In one sense, the method applies equally well to space-fractional PDEs as it does to integer-order PDEs. However, there is a significant challenge when solving space-fractional PDEs in this way, owing to the non-local nature of the fractional derivatives. Each equation in the resulting semi-discrete system involves contributions from every spatial node in the domain. This has important consequences for the efficiency of the numerical solver, especially when the system is large. First, the Jacobian matrix of the system is dense, and hence methods that avoid the need to form and factorise this matrix are preferred. Second, since the cost of evaluating the discrete equations is high, it is essential to minimise the number of evaluations required to advance the solution in time. In this paper, we show how an effective preconditioner is essential for improving the efficiency of the method of lines for solving a quite general two-sided, nonlinear space-fractional diffusion equation. A key contribution is to show, how to construct suitable banded approximations to the system Jacobian for preconditioning purposes that permit high orders and large stepsizes to be used in the temporal integration, without requiring dense matrices to be formed. The results of numerical experiments are presented that demonstrate the effectiveness of this approach.

A preconditioned Lanczos method for space-fractional reaction-diffusion equations on two-dimensional unstructured meshes

We consider a two-dimensional space-fractional reaction diffusion equation with a fractional Laplacian operator and homogeneous Neumann boundary conditions. The finite volume method is used with the matrix transfer technique of Ilić et al. (2006) to discretise in space, yielding a system of equations that requires the action of a matrix function to solve at each timestep. Rather than form this matrix function explicitly, we use Krylov subspace techniques to approximate the action of this matrix function. Specifically, we apply the Lanczos method, after a suitable transformation of the problem to recover symmetry. To improve the convergence of this method, we utilise a preconditioner that deflates the smallest eigenvalues from the spectrum. We demonstrate the efficiency of our approach for a fractional Fisher’s equation on the unit disk.

Stratospheric dust collections : valuable resources for space and atmospheric scientists

Collections of solid particles from the Earths' stratosphere have been a significant part of atmospheric research programs since 1965 [1], but it has only been in the past decade that space-related disciplines have provided the impetus for a continued interest in these collections. Early research on specific particle types collected from the stratosphere established that interplanetary dust particles (IDP's) can be collected efficiently and in reasonable abundance using flat-plate collectors [2-4]. The tenacity of Brownlee and co-workers in this subfield of cosmochemistry has led to the establishment of a successful IDP collection and analysis program (using flat-plate collectors on high-flying aircraft) based on samples available for distribution from Johnson Space Center [5]. Other stratospheric collections are made, but the program at JSC offers a unique opportunity to study well-documented, individual particles (or groups of particles) from a wide variety of sources [6]. The nature of the collection and curation process, as well as the timeliness of some sampling periods [7], ensures that all data obtained from stratospheric particles is a valuable resource for scientists from a wide range of disciplines. A few examples of the uses of these stratospheric dust collections are outlined below.

The analytical solution and numerical solution of the fractional diffusion-wave equation with damping

Fractional partial differential equations have been applied to many problems in physics, finance, and engineering. Numerical methods and error estimates of these equations are currently a very active area of research. In this paper we consider a fractional diffusionwave equation with damping. We derive the analytical solution for the equation using the method of separation of variables. An implicit difference approximation is constructed. Stability and convergence are proved by the energy method. Finally, two numerical examples are presented to show the effectiveness of this approximation.

Deformation with coupled chemical diffusion

The deformation of rocks is commonly intimately associated with metamorphic reactions. This paper is a step towards understanding the behaviour of fully coupled, deforming, chemically reacting systems by considering a simple example of the problem comprising a single layer system with elastic-power law viscous constitutive behaviour where the deformation is controlled by the diffusion of a single chemical component that is produced during a metamorphic reaction. Analysis of the problem using the principles of non-equilibrium thermodynamics allows the energy dissipated by the chemical reaction-diffusion processes to be coupled with the energy dissipated during deformation of the layers. This leads to strain-rate softening behaviour and the resultant development of localised deformation which in turn nucleates buckles in the layer. All such diffusion processes, in leading to Herring-Nabarro, Coble or “pressure solution” behaviour, are capable of producing mechanical weakening through the development of a “chemical viscosity”, with the potential for instability in the deformation. For geologically realistic strain rates these chemical feed-back instabilities occur at the centimetre to micron scales, and so produce structures at these scales, as opposed to thermal feed-back instabilities that become important at the 100–1000 m scales.

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.

Characterisation of atmospheric deposited particles during a dust storm in urban areas of Eastern Australia

The characteristics of dust particles deposited during the 2009 dust storm in the Gold Coast and Brisbane regions of Australia are discussed in this paper. The study outcomes provide important knowledge in relation to the potential impacts of dust storm related pollution on ecosystem health in the context that the frequency of dust storms is predicted to increase due to anthropogenic desert surface modifications and climate change impacts. The investigated dust storm contributed a large fraction of fine particles to the environment with an increased amount of total suspended solids, compared to dry deposition under ambient conditions. Although the dust storm passed over forested areas, the organic carbon content in the dust was relatively low. The primary metals present in the dust storm deposition were aluminium, iron and manganese, which are common soil minerals in Australia. The dust storm deposition did not contain significant loads of nickel, cadmium, copper and lead, which are commonly present in the urban environment. Furthermore, the comparison between the ambient and dust storm chromium and zinc loads suggested that these metals were contributed to the dust storm by local anthropogenic sources. The potential ecosystem health impacts of the 2009 dust storm include, increased fine solids deposition on ground surfaces resulting in an enhanced capacity to adsorb toxic pollutants as well as increased aluminium, iron and manganese loads. In contrast, the ecosystem health impacts related to organic carbon and other metals from dust storm atmospheric deposition are not considered to be significant.

Comment on "Local accumulation times for source, diffusion, and degradation models in two and three dimensions" [J. Chem. Phys. 138, 104121 (2013)]

In a recent paper, Gordon, Muratov, and Shvartsman studied a partial differential equation (PDE) model describing radially symmetric diffusion and degradation in two and three dimensions. They paid particular attention to the local accumulation time (LAT), also known in the literature as the mean action time, which is a spatially dependent timescale that can be used to provide an estimate of the time required for the transient solution to effectively reach steady state. They presented exact results for three-dimensional applications and gave approximate results for the two-dimensional analogue. Here we make two generalizations of Gordon, Muratov, and Shvartsman’s work: (i) we present an exact expression for the LAT in any dimension and (ii) we present an exact expression for the variance of the distribution. The variance provides useful information regarding the spread about the mean that is not captured by the LAT. We conclude by describing further extensions of the model that were not considered by Gordon,Muratov, and Shvartsman. We have found that exact expressions for the LAT can also be derived for these important extensions...

First-principle studies of the formation and diffusion of hydrogen vacancies in magnesium hydride

Ab initio density functional theory (DFT) calculations are performed to study the formation and diffusion of hydrogen vacancies on MgH2(110) surface and in bulk. We find that the formation energies for a single H-vacancy increase slightly from the surface to deep layers. The energies for creating adjacent surface divancacies at two inplane sites and at an inplane and a bridge site are even smaller than that for the formation of a single H-vacancy, a fact that is attributed to the strong vacancy−vacancy interactions. The diffusion of an H-vacancy from an in-plane site to a bridge site on the surface has the smallest activation barrier calculated at 0.15 eV and should be fast at room temperature. The activation barriers computed for H-vacancy diffusion from the surface into sublayers are all less than 0.70 eV, which is much smaller than the activation energy for desorption of hydrogen on the MgH2(110) surface (1.78−2.80 eV/H2). This suggests that surface desorption is more likely than vacancy diffusion to be rate determining, such that finding effective catalyst on the MgH2 surface to facilitate desorption will be very important for improving overall dehydrogenation performance.

Atomic hydrogen diffusion in novel magnesium nanostructures : the impact of incorporated subsurface carbon atoms

Ab initio Density Functional Theory (DFT) calculations are performed to study the diffusion of atomic hydrogen on a Mg(0001) surface and their migration into the subsurface layers. A carbon atom located initially on a Mg(0001) surface can migrate into the sub-surface layer and occupy a fcc site, with charge transfer to the C atom from neighboring Mg atoms. The cluster of postively charged Mg atoms surrounding a sub-surface C is then shown to facilitate the dissociative chemisorption of molecular hydrogen on the Mg(0001) surface, and the surface migration and subsequent diffusion into the subsurface of atomic hydrogen. This helps rationalize the experimentally-observed improvement in absorption kinetics of H2 when graphite or single walled carbon nanotubes (SWCNT) are introduced into the Mg powder during ball milling.

Computer modeling of diffusion in biological tissues

Molecular-level computer simulations of restricted water diffusion can be used to develop models for relating diffusion tensor imaging measurements of anisotropic tissue to microstructural tissue characteristics. The diffusion tensors resulting from these simulations can then be analyzed in terms of their relationship to the structural anisotropy of the model used. As the translational motion of water molecules is essentially random, their dynamics can be effectively simulated using computers. In addition to modeling water dynamics and water-tissue interactions, the simulation software of the present study was developed to automatically generate collagen fiber networks from user-defined parameters. This flexibility provides the opportunity for further investigations of the relationship between the diffusion tensor of water and morphologically different models representing different anisotropic tissues.