Efficient preconditioning of the method of lines for solving nonlinear two-sided space-fractional diffusion equations

A standard method for the numerical solution of partial differential equations (PDEs) is the method of lines. In this approach the PDE is discretised in space using �finite di�fferences or similar techniques, and the resulting semidiscrete problem in time is integrated using an initial value problem solver. A significant challenge when applying the method of lines to fractional PDEs is that the non-local nature of the fractional derivatives results in a discretised system where each equation involves contributions from many (possibly every) spatial node(s). This has important consequences for the effi�ciency of the numerical solver. First, 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. Second, since the Jacobian matrix of the system is dense (partially or fully), methods that avoid the need to form and factorise this matrix are preferred. In this paper, we consider a nonlinear two-sided space-fractional di�ffusion equation in one spatial dimension. A key contribution of this paper is to demonstrate how an eff�ective preconditioner is crucial for improving the effi�ciency of the method of lines for solving this equation. In particular, we 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.

Alternating direction implicit numerical method for the space and time fractional Bloch-Torrey equation in 3-D

Recently, some authors have considered a new diffusion model–space and time fractional Bloch-Torrey equation (ST-FBTE). Magin et al. (2008) have derived analytical solutions with fractional order dynamics in space (i.e., _ = 1, β an arbitrary real number, 1 < β ≤ 2) and time (i.e., 0 < α < 1, and β = 2), respectively. Yu et al. (2011) have derived an analytical solution and an effective implicit numerical method for solving ST-FBTEs, and also discussed the stability and convergence of the implicit numerical method. However, due to the computational overheads necessary to perform the simulations for nuclear magnetic resonance (NMR) in three dimensions, they present a study based on a two-dimensional example to confirm their theoretical analysis. Alternating direction implicit (ADI) schemes have been proposed for the numerical simulations of classic differential equations. The ADI schemes will reduce a multidimensional problem to a series of independent one-dimensional problems and are thus computationally efficient. In this paper, we consider the numerical solution of a ST-FBTE on a finite domain. The time and space derivatives in the ST-FBTE are replaced by the Caputo and the sequential Riesz fractional derivatives, respectively. A fractional alternating direction implicit scheme (FADIS) for the ST-FBTE in 3-D is proposed. Stability and convergence properties of the FADIS are discussed. Finally, some numerical results for ST-FBTE are given.

Numerical simulation of a new two-dimensional variable-order fractional percolation equation in non-homogeneous porous media

Percolation flow problems are discussed in many research fields, such as seepage hydraulics, groundwater hydraulics, groundwater dynamics and fluid dynamics in porous media. Many physical processes appear to exhibit fractional-order behavior that may vary with time, or space, or space and time. The theory of pseudodifferential operators and equations has been used to deal with this situation. In this paper we use a fractional Darcys law with variable order Riemann-Liouville fractional derivatives, this leads to a new variable-order fractional percolation equation. In this paper, a new two-dimensional variable-order fractional percolation equation is considered. A new implicit numerical method and an alternating direct method for the two-dimensional variable-order fractional model is proposed. Consistency, stability and convergence of the implicit finite difference method are established. Finally, some numerical examples are given. The numerical results demonstrate the effectiveness of the methods. This technique can be used to simulate a three-dimensional variable-order fractional percolation equation.

A finite volume method for solving the two-sided time-space fractional advection-dispersion equation

The field of fractional differential equations provides a means for modelling transport processes within complex media which are governed by anomalous transport. Indeed, the application to anomalous transport has been a significant driving force behind the rapid growth and expansion of the literature in the field of fractional calculus. In this paper, we present a finite volume method to solve the time-space two-sided fractional advection dispersion equation on a one-dimensional domain. Such an equation allows modelling different flow regime impacts from either side. The finite volume formulation provides a natural way to handle fractional advection-dispersion equations written in conservative form. The novel spatial discretisation employs fractionally-shifted Gr¨unwald formulas to discretise the Riemann-Liouville fractional derivatives at control volume faces in terms of function values at the nodes, while the L1-algorithm is used to discretise the Caputo time fractional derivative. Results of numerical experiments are presented to demonstrate the effectiveness of the approach.

Their space, my space, our space? Report of the 2009 Redland City Council Youth and Public Space Survey

This report analyses data collected through the Redland City Council’s Young People and Public Space Survey of 2148 high school students aged 12-19. The survey conducted in 2009 explored their sense of safety and experiences in public spaces across the City, and views on what Council could do to improve these. It is apparent they base their assessment of a space as ‘public’ on their ‘use’ of a space alone or with friends, and where strangers may be present, rather than on a type of ownership of a space (public/ private). The findings of the survey are summarised according to the themes of safety, community attitudes towards young people being in public spaces, young people and authorities, young people’s views of what is needed, and understanding different young people’s experiences of public space.

Three-dimensional medical ultrasound

A number of groups around the world are working in the field of three dimensional(3D) ultrasound (US) in order to obtain higher quality diagnostic information. 3D US, in general, involves collecting a sequence of conventional 2D US images along with information on the position and orientation of each image plane. A transformation matrix is calculated relating image space to real world space. This allows image pixels and region of interest (ROI) points drawn on the image to be displayed in 3D. The 3D data can be used for the production of volume or surface rendered images, or for the direct calculation of ROI volumes.

Serious playground : using Second Life to engage high school students in urban planning

Virtual world platforms such as Second Life have been successfully used in educational contexts to motivate and engage learners. This article reports on an exploratory workshop involving a group of high school students using Second Life for an urban planning project. Young people are traditionally an under-represented demographic when it comes to participating in urban planning and decision making processes. The research team developed activities that combined technology with a constructivist approach to learning. Real world experiences and purposes ensured that the workshop enabled students to see the relevance of their learning. Our design also ensured that play remained an important part of the learning. By conceiving of the workshop as a ‘serious playground’ we investigated the ludic potential of learning in a virtual world.

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.

Experiences of neighbourhood walkability among older Australians living in high density inner-city areas

Walking as an out-of-home mobility activity is recognised for its contribution to healthy and active ageing. The environment can have a powerful effect on the amount of walking activity undertaken by older people, thereby influencing their capacity to maintain their wellbeing and independence. This paper reports the findings from research examining the experiences of neighbourhood walking for 12 older people from six different inner-city high density suburbs, through analysis of data derived from travel diaries, individual time/space activity maps (created via GPS tracking over a seven-day period and GIS technology), and in-depth interviews. Reliance on motor vehicles, the competing interests of pedestrians and cyclists on shared pathways and problems associated with transit systems, public transport, and pedestrian infrastructure emerged as key barriers to older people venturing out of home on foot. GPS and GIS technology provide new opportunities for furthering understanding of the out-of-home mobility of older populations.

Orthogonality defect and reduced search-space size for solving integer least-squares problems

In the context of ambiguity resolution (AR) of Global Navigation Satellite Systems (GNSS), decorrelation among entries of an ambiguity vector, integer ambiguity search and ambiguity validations are three standard procedures for solving integer least-squares problems. This paper contributes to AR issues from three aspects. Firstly, the orthogonality defect is introduced as a new measure of the performance of ambiguity decorrelation methods, and compared with the decorrelation number and with the condition number which are currently used as the judging criterion to measure the correlation of ambiguity variance-covariance matrix. Numerically, the orthogonality defect demonstrates slightly better performance as a measure of the correlation between decorrelation impact and computational efficiency than the condition number measure. Secondly, the paper examines the relationship of the decorrelation number, the condition number, the orthogonality defect and the size of the ambiguity search space with the ambiguity search candidates and search nodes. The size of the ambiguity search space can be properly estimated if the ambiguity matrix is decorrelated well, which is shown to be a significant parameter in the ambiguity search progress. Thirdly, a new ambiguity resolution scheme is proposed to improve ambiguity search efficiency through the control of the size of the ambiguity search space. The new AR scheme combines the LAMBDA search and validation procedures together, which results in a much smaller size of the search space and higher computational efficiency while retaining the same AR validation outcomes. In fact, the new scheme can deal with the case there are only one candidate, while the existing search methods require at least two candidates. If there are more than one candidate, the new scheme turns to the usual ratio-test procedure. Experimental results indicate that this combined method can indeed improve ambiguity search efficiency for both the single constellation and dual constellations respectively, showing the potential for processing high dimension integer parameters in multi-GNSS environment.

Micromarrows—Three-dimensional coculture of hematopoietic stem cells and mesenchymal stromal cells

Hematopoietic stem cell (HSC) transplant is a well established curative therapy for some hematological malignancies. However, achieving adequate supply of HSC from some donor tissues can limit both its application and ultimate efficacy. The theory that this limitation could be overcome by expanding the HSC population before transplantation has motivated numerous laboratories to develop ex vivo expansion processes. Pioneering work in this field utilized stromal cells as support cells in cocultures with HSC to mimic the HSC niche. We hypothesized that through translation of this classic coculture system to a three-dimensional (3D) structure we could better replicate the niche environment and in turn enhance HSC expansion. Herein we describe a novel high-throughput 3D coculture system where murine-derived HSC can be cocultured with mesenchymal stem/stromal cells (MSC) in 3D microaggregates—which we term “micromarrows.” Micromarrows were formed using surface modified microwells and their ability to support HSC expansion was compared to classic two-dimensional (2D) cocultures. While both 2D and 3D systems provide only a modest total cell expansion in the minimally supplemented medium, the micromarrow system supported the expansion of approximately twice as many HSC candidates as the 2D controls. Histology revealed that at day 7, the majority of bound hematopoietic cells reside in the outer layers of the aggregate. Quantitative polymerase chain reaction demonstrates that MSC maintained in 3D aggregates express significantly higher levels of key hematopoietic niche factors relative to their 2D equivalents. Thus, we propose that the micromarrow platform represents a promising first step toward a high-throughput HSC 3D coculture system that may enable in vitro HSC niche recapitulation and subsequent extensive in vitro HSC self-renewal.

What makes inner city high density living attractive? Insight from residents in Brisbane, Australia

Australia’s urban form and planning has shifted from traditional individual dwellings on spacious suburban blocks towards higher density urban consolidation. Despite relatively strong market demand for inner city high density (ICHD) living, there is ongoing need to explore and understand the aspects that make this urban form liveable and sustainable. The purpose of this research is to explore the viewpoints of current ICHD residents to better understand the liveability and sustainability matters that affect their everyday experiences and perceptions of this urban form. Qualitative interviews with 24 ICHD Brisbane (Australia) residents illustrates their perceptions and experiences of liveability and the ways in which it is broadly understood within three main domains and nine key sub-concepts, including: individual dwelling (thermal comfort, natural light and balconies, noise mitigation), building complex (shared space, good neighbour protocols, environmental sustainability) and the community (transport, amenities, sense of community). Focussing on the experience of ICHD residents, this research highlights the ways in which multiple aspects of the immediate living environment, the dwelling, building complex and the community intertwine to provide residents with a liveable space. The results show that urban features that reflect current societal pressure for greater sustainability such as lower energy use are the exact same features sought by ICHD residents in determining their liveability. By highlighting the aspects current ICHD residents value most about their dwellings, buildings and communities, these findings will help inform policy-makers, planners, developers and designers as they create urban spaces and dwellings that are more liveable and sustainable.

An evolutionary computation approach to three-dimensional path planning for unmanned aerial vehicles with tactical and kinematic constraints

This paper presents a novel evolutionary computation approach to three-dimensional path planning for unmanned aerial vehicles (UAVs) with tactical and kinematic constraints. A genetic algorithm (GA) is modified and extended for path planning. Two GAs are seeded at the initial and final positions with a common objective to minimise their distance apart under given UAV constraints. This is accomplished by the synchronous optimisation of subsequent control vectors. The proposed evolutionary computation approach is called synchronous genetic algorithm (SGA). The sequence of control vectors generated by the SGA constitutes to a near-optimal path plan. The resulting path plan exhibits no discontinuity when transitioning from curve to straight trajectories. Experiments and results show that the paths generated by the SGA are within 2% of the optimal solution. Such a path planner when implemented on a hardware accelerator, such as field programmable gate array chips, can be used in the UAV as on-board replanner, as well as in ground station systems for assisting in high precision planning and modelling of mission scenarios.

Mechanism of dehydroxylation temperature decrease and high temperature phase transition of coal-bearing strata kaolinite intercalated by potassium acetate

The thermal decomposition and dehydroxylation process of coal-bearing strata kaolinite–potassium acetate intercalation complex (CSKK) has been studied using X-ray diffraction (XRD), infrared spectroscopy (IR), thermal analysis, mass spectrometric analysis and infrared emission spectroscopy. The XRD results showed that the potassium acetate (KAc) have been successfully intercalated into coal-bearing strata kaolinite with an obvious basal distance increase of the first basal peak, and the positive correlation was found between the concentration of intercalation regent KAc and the degree of intercalation. As the temperature of the system is raised, the formation of KHCO3, KCO3 and KAlSiO4, which is derived from the thermal decomposition or phase transition of CSKK, is observed in sequence. The IR results showed that new bands appeared, the position and intensities shift can also be found when the concentration of intercalation agent is raised. The thermal analysis and mass spectrometric analysis results revealed that CSKK is stable below 300 °C, and the thermal decomposition products (H2O and CO2) were further proved by the mass spectrometric analysis. A comparison of thermal analysis results of original coal-bearing strata kaolinite and its intercalation complex gives new discovery that not only a new mass loss peak is observed at 285 °C, but also the temperature of dehydroxylation and dehydration of coal bearing strata kaolinite is decreased about 100 °C. This is explained on the basis of the interlayer space of the kaolinite increased obviously after being intercalated by KAc, which led to the interlayer hydrogen bonds weakened, enables the dehydroxylation from kaolinite surface more easily. Furthermore, the possible structural model for CSKK has been proposed, with further analysis required in order to prove the most possible structures.

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.