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.

A Novel High Order Space-Time Spectral Method for the Time Fractional Fokker--Planck Equation

The fractional Fokker-Planck equation is an important physical model for simulating anomalous diffusions with external forces. Because of the non-local property of the fractional derivative an interesting problem is to explore high accuracy numerical methods for fractional differential equations. In this paper, a space-time spectral method is presented for the numerical solution of the time fractional Fokker-Planck initial-boundary value problem. The proposed method employs the Jacobi polynomials for the temporal discretization and Fourier-like basis functions for the spatial discretization. Due to the diagonalizable trait of the Fourier-like basis functions, this leads to a reduced representation of the inner product in the Galerkin analysis. We prove that the time fractional Fokker-Planck equation attains the same approximation order as the time fractional diffusion equation developed in [23] by using the present method. That indicates an exponential decay may be achieved if the exact solution is sufficiently smooth. Finally, some numerical results are given to demonstrate the high order accuracy and efficiency of the new numerical scheme. The results show that the errors of the numerical solutions obtained by the space-time spectral method decay exponentially.

H2 purification by functionalized graphdiyne – role of nitrogen doping

A nitrogen modified graphdiyne is investigated concerning its performance for hydrogen purification from CH4 and CO by density functional theory with dispersion correction and transition state theory. After nitrogen doping, the porous N-graphdiyne nano-mesh shows a reduced H2 diffusion barrier and increased CH4/CO diffusion barriers, hence leading to an enhanced hydrogen purification capability.

Performance dispersion for evidence-based classification of stationary throwers

Background There is a need for better understanding of the dispersion of classification-related variable to develop an evidence-based classification of athletes with a disability participating in stationary throwing events. Objectives The purposes of this study are (A) to describe tools designed to comprehend and represent the dispersion of the performance between successive classes, and (B) to present this dispersion for the elite male and female stationary shot-putters who participated in Beijing 2008 Paralympic Games. Study design Retrospective study Methods This study analysed a total of 479 attempts performed by 114 male and female stationary shot-putters in three F30s (F32-F34) and six F50s (F52-F58) classes during the course of eight events during Beijing 2008 Paralympic Games. Results The average differences of best performance were 1.46±0.46 m for males between F54 and F58 classes as well as 1.06±1.18 m for females between F55 and F58 classes. The results demonstrated a linear relationship between best performance and classification while revealing two male Gold Medallists in F33 and F52 classes were outliers. Conclusions This study confirms the benefits of the comparative matrices, performance continuum and dispersion plots to comprehend classification-related variables. The work presented here represents a stepping stone into biomechanical analyses of stationary throwers, particularly on the eve of the London 2012 Paralympic Games where new evidences could be gathered.

Quantitative characterization of kaolinite dispersibility in styrene–butadiene rubber composites by fractal dimension

A fractal method was introduced to quantitatively characterize the dispersibility of modified kaolinite (MK) and precipitated silica (PS) in styrene–butadiene rubber (SBR) matrix based on the lower magnification transmission electron microscopic images. The fractal dimension (FD) is greater, and the dispersion is worse. The fractal results showed that the dispersibility of MK in the latex blending sample is better than that in the mill blending samples. With the increase of kaolinite content, the FD increases from 1.713 to 1.800, and the dispersibility of kaolinite gradually decreases. There is a negative correlation between the dispersibility and loading content. With the decrease of MK and increase of PS, the FD significantly decreases from 1.735 to 1.496 and the dipersibility of kaolinite remarkably increases. The hybridization can improve the dispersibility of fillers in polymer matrix. The FD can be used to quantitatively characterize the aggregation and dispersion of kaolinite sheets in rubber matrix.

Effect of carrier size on the dispersion of salmeterol xinafoate from interactive mixtures

The objective of this study was to determine the influence of lactose carrier size on drug dispersion of salmeterol xinafoate (SX) from interactive mixtures. SX dispersion was measured by using the fine particle fractions determined by a twin stage impinger attached to a Rotahaler1. The particle size of the lactose carrier in the SX interactive mixtures was varied using a range of commercial inhalation-grade lactoses. In addition, differing size fractions of individual lactose samples were achieved by dry sieving. The dispersion ofSXappeared to increase as the particle size of the lactose carrier decreased for the mixtures prepared from different particle size commercial samples of lactose and from different sieve fractions of the same lactose. Fine particles of lactose (<5 mm) associated with the lactose carrier were removed from the carrier surface by a wet decantation process to produce lactose samples with low but similar concentrations of fine lactose particles. The fine particle fractions of SX in mixtures prepared with the decanted lactose decreased significantly (analysis of variance, p<0.001) and the degree of dispersion became independent of the volume mean diameter of the carriers (analysis of variance, p<0.05). The dispersion behavior is therefore associated with the presence of fine adhered particles associated with the carriers and the inherent size of the carrier itself has little influence on dispersion.

In-situ carbon control in the preparation of precursors to boron carbide by a non-aqueous solution technique

Synthesis of high quality boron carbide (B4C) powders is achieved by carbothermal reduction of boron oxide (B2O3) from a condensed boric acid (H3BO3)/polyvinyl acetate (PVAc) product. Precursor solutions are prepared via free radical polymerisation of vinyl acetate (VA) monomer in methanol in the presence of dissolved H3BO3. A condensed product is then formed by flash evaporation under vacuum. As excess VA monomer is removed at the evaporation step, the polymerisation time is used to manage availability of carbon for reaction. This control of carbon facilitates dispersion of H3BO3 in solution due to the presence of residual VA monomer. B4C powders with very low residual carbon are formed at temperatures as low as 1,250 °C with a 4 hour residence time.

Estimate of Lagrangian integral scales in shallow tidal water using high resolution GPS-tracked drifters

In estuaries and natural water channels, the estimate of velocity and dispersion coefficients is critical to the knowledge of scalar transport and mixing. This estimate is rarely available experimentally at sub-tidal time scale in shallow water channels where high frequency is required to capture its spatio-temporal variation. This study estimates Lagrangian integral scales and autocorrelation curves, which are key parameters for obtaining velocity fluctuations and dispersion coefficients, and their spatio-temporal variability from deployments of Lagrangian drifters sampled at 10 Hz for a 4-hour period. The power spectral densities of the velocities between 0.0001 and 0.8 Hz were well fitted with a slope of 5/3 predicted by Kolmogorov’s similarity hypothesis within the inertial subrange, and were similar to the Eulerian power spectral previously observed within the estuary. The result showed that large velocity fluctuations determine the magnitude of the integral time scale, TL. Overlapping of short segments improved the stability of the estimate of TL by taking advantage of the redundant data included in the autocorrelation function. The integral time scales were about 20 s and varied by up to a factor of 8. These results are essential inputs for spatial binning of velocities, Lagrangian stochastic modelling and single particle analysis of the tidal estuary.

Characterization of alluvial formation by stochastic modelling of paleo-fluvial processes: The concept and method

Modelling fluvial processes is an effective way to reproduce basin evolution and to recreate riverbed morphology. However, due to the complexity of alluvial environments, deterministic modelling of fluvial processes is often impossible. To address the related uncertainties, we derive a stochastic fluvial process model on the basis of the convective Exner equation that uses the statistics (mean and variance) of river velocity as input parameters. These statistics allow for quantifying the uncertainty in riverbed topography, river discharge and position of the river channel. In order to couple the velocity statistics and the fluvial process model, the perturbation method is employed with a non-stationary spectral approach to develop the Exner equation as two separate equations: the first one is the mean equation, which yields the mean sediment thickness, and the second one is the perturbation equation, which yields the variance of sediment thickness. The resulting solutions offer an effective tool to characterize alluvial aquifers resulting from fluvial processes, which allows incorporating the stochasticity of the paleoflow velocity.

Transition to an unsteady flow induced by a fin on the sidewall of a differentially heated air-filled square cavity and heat transfer

The transition from a steady to an unsteady flow induced by an adiabatic fin on the sidewall of a differentially heated air-filled cavity is numerically investigated. Numerical simulations have been performed over the range of Rayleigh numbers from Ra = 105–109. The temporal development and spatial structures of natural convection flows in the cavity with a fin are described. It has been demonstrated that the fin may induce the transition to an unsteady flow and the critical Rayleigh number for the occurrence of the transition is between 3.72 × 106 and 3.73 × 106. Furthermore, the peak frequencies of the oscillations triggered by different mechanisms are obtained through spectral analysis. It has been found that the flow rate through the cavity with a fin is larger than that without a fin under the unsteady flow, indicating that the fin may improve the unsteady flow in the cavity.

Fourier spectral methods for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reaction-diffusion equations described by the fractional Laplacian in bounded rectangular domains ofRn. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is illustrated by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models, together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator.

Fractional diffusion models of cardiac electrical propagation: Role of structural heterogeneity in dispersion of repolarization

Impulse propagation in biological tissues is known to be modulated by structural heterogeneity. In cardiac muscle, improved understanding on how this heterogeneity influences electrical spread is key to advancing our interpretation of dispersion of repolarization. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a means of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, analysed against in vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies relevant characteristics of cardiac electrical propagation at tissue level. These include conduction effects on action potential (AP) morphology, the shortening of AP duration along the activation pathway and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media.

Numerical analysis of a new space-time variable fractional order advection-dispersion equation

Many physical processes appear to exhibit fractional order behavior that may vary with time and/or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. In this paper, we consider a new space–time variable fractional order advection–dispersion equation on a finite domain. The equation is obtained from the standard advection–dispersion equation by replacing the first-order time derivative by Coimbra’s variable fractional derivative of order α(x)∈(0,1]α(x)∈(0,1], and the first-order and second-order space derivatives by the Riemann–Liouville derivatives of order γ(x,t)∈(0,1]γ(x,t)∈(0,1] and β(x,t)∈(1,2]β(x,t)∈(1,2], respectively. We propose an implicit Euler approximation for the equation and investigate the stability and convergence of the approximation. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.

Stress analysis on human arterial plaques by fluid structure interactions - Multi-case study

Atherosclerotic plaque rupture has been extensively considered as the leading cause of death in the world. It is believed that high stress within plaque can be an important factor which can trigger the rupture of the plaque. High resolution multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components (arterial wall, lipids, and fibrous cap) to be visualized in vivo [1]. The patient specific finite element model can be generated from the image data to perform stress analysis and provide critical information on understanding plaque rupture mechanisms [2]. The present work is to apply the procedure to a total of 14 patients (S1 ∼ S14), to study the stress distributions on carotid artery plaque reconstructed from multi-spectral magnetic resonance images, and the possible relationships between stress and plaque burdens.

On the order of the fractional Laplacian in determining the spatio-temporal evolution of a space-fractional model of cardiac electrophysiology

Space-fractional operators have been used with success in a variety of practical applications to describe transport processes in media characterised by spatial connectivity properties and high structural heterogeneity altering the classical laws of diffusion. This study provides a systematic investigation of the spatio-temporal effects of a space-fractional model in cardiac electrophysiology. We consider a simplified model of electrical pulse propagation through cardiac tissue, namely the monodomain formulation of the Beeler-Reuter cell model on insulated tissue fibres, and obtain a space-fractional modification of the model by using the spectral definition of the one-dimensional continuous fractional Laplacian. The spectral decomposition of the fractional operator allows us to develop an efficient numerical method for the space-fractional problem. Particular attention is paid to the role played by the fractional operator in determining the solution behaviour and to the identification of crucial differences between the non-fractional and the fractional cases. We find a positive linear dependence of the depolarization peak height and a power law decay of notch and dome peak amplitudes for decreasing orders of the fractional operator. Furthermore, we establish a quadratic relationship in conduction velocity, and quantify the increasingly wider action potential foot and more pronounced dispersion of action potential duration, as the fractional order is decreased. A discussion of the physiological interpretation of the presented findings is made.