Further development of discrete computational techniques for calculation of restricted diffusion propagators in porous media

Magnetic resonance is a well-established tool for structural characterisation of porous media. Features of pore-space morphology can be inferred from NMR diffusion-diffraction plots or the time-dependence of the apparent diffusion coefficient. Diffusion NMR signal attenuation can be computed from the restricted diffusion propagator, which describes the distribution of diffusing particles for a given starting position and diffusion time. We present two techniques for efficient evaluation of restricted diffusion propagators for use in NMR porous-media characterisation. The first is the Lattice Path Count (LPC). Its physical essence is that the restricted diffusion propagator connecting points A and B in time t is proportional to the number of distinct length-t paths from A to B. By using a discrete lattice, the number of such paths can be counted exactly. The second technique is the Markov transition matrix (MTM). The matrix represents the probabilities of jumps between every pair of lattice nodes within a single timestep. The propagator for an arbitrary diffusion time can be calculated as the appropriate matrix power. For periodic geometries, the transition matrix needs to be defined only for a single unit cell. This makes MTM ideally suited for periodic systems. Both LPC and MTM are closely related to existing computational techniques: LPC, to combinatorial techniques; and MTM, to the Fokker-Planck master equation. The relationship between LPC, MTM and other computational techniques is briefly discussed in the paper. Both LPC and MTM perform favourably compared to Monte Carlo sampling, yielding highly accurate and almost noiseless restricted diffusion propagators. Initial tests indicate that their computational performance is comparable to that of finite element methods. Both LPC and MTM can be applied to complicated pore-space geometries with no analytic solution. We discuss the new methods in the context of diffusion propagator calculation in porous materials and model biological tissues.

A meshless method based on Point Interpolation Method (PIM) for the space fractional diffusion equation

This paper aims to develop a meshless approach based on the Point Interpolation Method (PIM) for numerical simulation of a space fractional diffusion equation. Two fully-discrete schemes for the one-dimensional space fractional diffusion equation are obtained by using the PIM and the strong-forms of the space diffusion equation. Numerical examples with different nodal distributions are studied to validate and investigate the accuracy and efficiency of the newly developed meshless approach.

A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients

In this paper, we derive a new nonlinear two-sided space-fractional diffusion equation with variable coefficients from the fractional Fick’s law. A semi-implicit difference method (SIDM) for this equation is proposed. The stability and convergence of the SIDM are discussed. For the implementation, we develop a fast accurate iterative method for the SIDM by decomposing the dense coefficient matrix into a combination of Toeplitz-like matrices. This fast iterative method significantly reduces the storage requirement of O(n2)O(n2) and computational cost of O(n3)O(n3) down to n and O(nlogn)O(nlogn), where n is the number of grid points. The method retains the same accuracy as the underlying SIDM solved with Gaussian elimination. Finally, some numerical results are shown to verify the accuracy and efficiency of the new method.

A Crank--Nicolson ADI spectral method for a two-dimensional riesz space fractional nonlinear reaction-diffusion equation

In this paper, a new alternating direction implicit Galerkin--Legendre spectral method for the two-dimensional Riesz space fractional nonlinear reaction-diffusion equation is developed. The temporal component is discretized by the Crank--Nicolson method. The detailed implementation of the method is presented. The stability and convergence analysis is strictly proven, which shows that the derived method is stable and convergent of order $2$ in time. An optimal error estimate in space is also obtained by introducing a new orthogonal projector. The present method is extended to solve the fractional FitzHugh--Nagumo model. Numerical results are provided to verify the theoretical analysis.

Impact of mobility and topology on information diffusion in MANETs

In some delay-tolerant communication systems such as vehicular ad-hoc networks, information flow can be represented as an infectious process, where each entity having already received the information will try to share it with its neighbours. The random walk and random waypoint models are popular analysis tools for these epidemic broadcasts, and represent two types of random mobility. In this paper, we introduce a simulation framework investigating the impact of a gradual increase of bias in path selection (i.e. reduction of randomness), when moving from the former to the latter. Randomness in path selection can significantly alter the system performances, in both regular and irregular network structures. The implications of these results for real systems are discussed in details.

Stability and convergence of a new finite volume method for a two-sided space-fractional diffusion equation

In this paper, we consider a two-sided space-fractional diffusion equation with variable coefficients on a finite domain. Firstly, based on the nodal basis functions, we present a new fractional finite volume method for the two-sided space-fractional diffusion equation and derive the implicit scheme and solve it in matrix form. Secondly, we prove the stability and convergence of the implicit fractional finite volume method and conclude that the method is unconditionally stable and convergent. Finally, some numerical examples are given to show the effectiveness of the new numerical method, and the results are in excellent agreement with theoretical analysis.

Proportionality : a valid alternative to correlation for relative data

Relative abundance data is common in the life sciences, but appreciation that it needs special analysis and interpretation is scarce. Correlation is popular as a statistical measure of pairwise association but should not be used on data that carry only relative information. Using timecourse yeast gene expression data, we show how correlation of relative abundances can lead to conclusions opposite to those drawn from absolute abundances, and that its value changes when different components are included in the analysis. Once all absolute information has been removed, only a subset of those associations will reliably endure in the remaining relative data, specifically, associations where pairs of values behave proportionally across observations. We propose a new statistic φ to describe the strength of proportionality between two variables and demonstrate how it can be straightforwardly used instead of correlation as the basis of familiar analyses and visualization methods.

Different correlation of Sphingosine-1-Phosphate receptor 1 with receptor activator of nuclear factor kappa B ligand and regulatory T cells in rat periapical lesions

Introduction Sphingosine-1-phosphate receptor 1 (S1P1) is crucial for regulation of immunity and bone metabolism. This study aimed to investigate the expression of S1P1 in rat periapical lesions and its relationship with receptor activator of nuclear factor kappa B ligand (RANKL) and regulatory T (Treg) cells. Methods Periapical lesions were induced by pulp exposure in the first lower molars of 55 Wistar rats. Thirty rats were killed on days 0, 7, 14, 21, 28, and 35, and their mandibles were harvested for x-ray imaging, micro–computed tomography scanning, histologic observation, immunohistochemistry, enzyme histochemistry, and double immunofluorescence analysis. The remaining 25 rats were killed on days 0, 14, 21, 28, and 35, and mandibles were harvested for flow cytometry. Results The volume and area of the periapical lesions increased from day 0 to day 21 and then remained comparably stable after day 28. S1P1-positive cells were observed in the inflammatory periapical regions; the number of S1P1-positive cells peaked at day 14 and then decreased from day 21 to day 35. The distribution of S1P1-positive cells was positively correlated with the dynamics of RANKL-positive cells but was negatively correlated with that of Treg cells. Conclusions S1P1 expression was differentially correlated with RANKL and Treg cell infiltration in the periapical lesions and is therefore a contributing factor to the pathogenesis of such lesions.

Performance evaluation of non-thermal plasma on particulate matter, ozone and CO2 correlation for diesel exhaust emission reduction

This study is seeking to investigate the effect of non-thermal plasma technology in the abatement of particulate matter (PM) from the actual diesel exhaust. Ozone (O3) strongly promotes PM oxidation, the main product of which is carbon dioxide (CO2). PM oxidation into the less harmful product (CO2) is the main objective whiles the correlation between PM, O3 and CO2 is considered. A dielectric barrier discharge reactor has been designed with pulsed power technology to produce plasma inside the diesel exhaust. To characterise the system under varied conditions, a range of applied voltages from 11 kVPP to 21kVPP at repetition rates of 2.5, 5, 7.5 and 10 kHz, have been experimentally investigated. The results show that by increasing the applied voltage and repetition rate, higher discharge power and CO2 dissociation can be achieved. The PM removal efficiency of more than 50% has been achieved during the experiments and high concentrations of ozone on the order of a few hundreds of ppm have been observed at high discharge powers. Furthermore, O3, CO2 and PM concentrations at different plasma states have been analysed for time dependence. Based on this analysis, an inverse relationship between ozone concentration and PM removal has been found and the role of ozone in PM removal in plasma treatment of diesel exhaust has been highlighted.

Numerical treatment of a two-dimensional variable-order fractional nonlinear reaction-diffusion model

A two-dimensional variable-order fractional nonlinear reaction-diffusion model is considered. A second-order spatial accurate semi-implicit alternating direction method for a two-dimensional variable-order fractional nonlinear reaction-diffusion model is proposed. Stability and convergence of the semi-implicit alternating direct method are established. Finally, some numerical examples are given to support our theoretical analysis. These numerical techniques can be used to simulate a two-dimensional variable order fractional FitzHugh-Nagumo model in a rectangular domain. This type of model can be used to describe how electrical currents flow through the heart, controlling its contractions, and are used to ascertain the effects of certain drugs designed to treat arrhythmia.

GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

Lower mantle diamonds from Rio Soriso (Juina area, Mato Grosso, Brazil)

The morphology, colour, fluorescence, cathodoluminescence, nitrogen content and aggregation state, internal structure and mineral inclusions have been studied for 69 alluvial diamonds from the Rio Soriso (Juina area, Mato Grosso State, Brazil). Nitrogen in most diamonds (53%) is fully aggregated as B centres, but there is also a large proportion of N-free stones (38%). A strong positive correlation between nitrogen and IR-active hydrogen concentrations is observed. The diamonds contain (in order of decreasing abundance) ferropericlase, CaSi-perovskite, magnetite, MgSi-perovskite, pyrrhotite, 'olivine', SiO2, perovskite, tetragonal almandine-pyrope phase and some other minerals represented by single grains. The Rio Soriso diamond suite is subdivided into several subpopulations that originated in upper and lower mantle of ultramafic and mafic compositions, with the largest subgroup forming in the ultramafic lower mantle. Analysed ferropericlase grains are enriched in Fe (Mg#=0.43-0.89), which is ascribed to their origin in the lowermost mantle. The Juina kimberlites may be unique in sampling the material from depths below 1,700 km that ascended in a plume formed at the core-mantle boundary.

Correlation between architectural variables and torque in the erector spinae muscle during maximal isometric contraction

This study analysed whether a significant relationship exists between the torque and muscle thickness and pennation angle of the erector spinae muscle during a maximal isometric lumbar extension with the lumbar spine in neutral position. This was a cross-sectional study in which 46 healthy adults performed three repetitions for 5 s of maximal isometric lumbar extension with rests of 90 s. During the lumbar extensions, bilateral ultrasound images of the erector spinae muscle (to measure pennation angle and muscle thickness) and torque were acquired. Reliability test analysis calculating the internal consistency (Cronbach's alpha) of the measure, correlation between pennation angle, muscle thickness and torque extensions were examined. Through a linear regression the contribution of each independent variable (muscle thickness and pennation angle) to the variation of the dependent variable (torque) was calculated. The results of the reliability test were: 0.976–0.979 (pennation angle), 0.980–0.980 (muscle thickness) and 0.994 (torque). The results show that pennation angle and muscle thickness were significantly related to each other with a range between 0.295 and 0.762. In addition, multiple regression analysis showed that the two variables considered in this study explained 68% of the variance in the torque. Pennation angle and muscle thickness have a moderate impact on the variance exerted on the torque during a maximal isometric lumbar extension with the lumbar spine in neutral position.

Extending genetic linkage analysis to diffusion tensor images to map single gene effects on brain fiber architecture

We extended genetic linkage analysis - an analysis widely used in quantitative genetics - to 3D images to analyze single gene effects on brain fiber architecture. We collected 4 Tesla diffusion tensor images (DTI) and genotype data from 258 healthy adult twins and their non-twin siblings. After high-dimensional fluid registration, at each voxel we estimated the genetic linkage between the single nucleotide polymorphism (SNP), Val66Met (dbSNP number rs6265), of the BDNF gene (brain-derived neurotrophic factor) with fractional anisotropy (FA) derived from each subject's DTI scan, by fitting structural equation models (SEM) from quantitative genetics. We also examined how image filtering affects the effect sizes for genetic linkage by examining how the overall significance of voxelwise effects varied with respect to full width at half maximum (FWHM) of the Gaussian smoothing applied to the FA images. Raw FA maps with no smoothing yielded the greatest sensitivity to detect gene effects, when corrected for multiple comparisons using the false discovery rate (FDR) procedure. The BDNF polymorphism significantly contributed to the variation in FA in the posterior cingulate gyrus, where it accounted for around 90-95% of the total variance in FA. Our study generated the first maps to visualize the effect of the BDNF gene on brain fiber integrity, suggesting that common genetic variants may strongly determine white matter integrity.

Brain fiber architecture, genetics, and intelligence: a high angular resolution diffusion imaging (HARDI) study

We developed an analysis pipeline enabling population studies of HARDI data, and applied it to map genetic influences on fiber architecture in 90 twin subjects. We applied tensor-driven 3D fluid registration to HARDI, resampling the spherical fiber orientation distribution functions (ODFs) in appropriate Riemannian manifolds, after ODF regularization and sharpening. Fitting structural equation models (SEM) from quantitative genetics, we evaluated genetic influences on the Jensen-Shannon divergence (JSD), a novel measure of fiber spatial coherence, and on the generalized fiber anisotropy (GFA) a measure of fiber integrity. With random-effects regression, we mapped regions where diffusion profiles were highly correlated with subjects' intelligence quotient (IQ). Fiber complexity was predominantly under genetic control, and higher in more highly anisotropic regions; the proportion of genetic versus environmental control varied spatially. Our methods show promise for discovering genes affecting fiber connectivity in the brain.