Stability and convergence of implicit numerical methods for a class of fractional advection-dispersion models

In this paper, a class of fractional advection-dispersion models (FADM) is investigated. These models include five fractional advection-dispersion models: the immobile, mobile/immobile time FADM with a temporal fractional derivative 0 < γ < 1, the space FADM with skewness, both the time and space FADM and the time fractional advection-diffusion-wave model with damping with index 1 < γ < 2. They describe nonlocal dependence on either time or space, or both, to explain the development of anomalous dispersion. These equations can be used to simulate regional-scale anomalous dispersion with heavy tails, for example, the solute transport in watershed catchments and rivers. We propose computationally effective implicit numerical methods for these FADM. The stability and convergence of the implicit numerical methods are analyzed and compared systematically. Finally, some results are given to demonstrate the effectiveness of our theoretical analysis.

Dry powder inhaler formulations : effect of carrier size on the drug dispersion

Background: The size of the carrier influences drug aerosolization from a dry powder inhaler (DPI) formulation. Lactose particles with irregular shape and rough surface in a variety of sizes are additionally used as carriers; however, contradictory reports exist regarding the effect of carrier size on the dispersion of drug. We examined the influence of the spherical particle size of the biodegradable polylactide-co-glycolide (PLGA) carrier on the aerosolization of a model drug, salbutamol sulphate (SS). Methods: Four different sizes (20-150 µm) of polymer carriers were fabricated using solvent evaporation technique and the dispersion of SS from these carriers was measured by a Twin Stage Impinger (TSI). The size and morphological properties of polymer carriers were determined by laser diffraction and SEM, respectively. Results: The FPF was found to increase from 5.6% to 21.3% with increasing carrier sizeup to150 µm. Conclusions: The aerosolization of drug increased linearly with the size of polymer carriers. For a fixed mass of drug particles in a formulation, the mass of drug particles per unit area of carriers is higher in formulations containing the larger carriers, which leads to an increase in the dispersion of drug due to the increased mechanical forces occurred between the carriers and the device walls.

Numerical techniques for simulating a fractional mathematical model of epidermal wound healing

A number of mathematical models investigating certain aspects of the complicated process of wound healing are reported in the literature in recent years. However, effective numerical methods and supporting error analysis for the fractional equations which describe the process of wound healing are still limited. In this paper, we consider the numerical simulation of a fractional mathematical model of epidermal wound healing (FMM-EWH), which is based on the coupled advection-diffusion equations for cell and chemical concentration in a polar coordinate system. The space fractional derivatives are defined in the Left and Right Riemann-Liouville sense. Fractional orders in the advection and diffusion terms belong to the intervals (0, 1) or (1, 2], respectively. Some numerical techniques will be used. Firstly, the coupled advection-diffusion equations are decoupled to a single space fractional advection-diffusion equation in a polar coordinate system. Secondly, we propose a new implicit difference method for simulating this equation by using the equivalent of Riemann-Liouville and Grünwald-Letnikov fractional derivative definitions. Thirdly, its stability and convergence are discussed, respectively. Finally, some numerical results are given to demonstrate the theoretical analysis.

Dry powder inhaler formulations-effect of polymer carrier size on the drug dispersion

Background The size of the carrier influences the aerosolization of drug from a dry powder inhaler (DPI) formulation. Currently, lactose monohydrate particles in a variety of sizes are preferably used in carrier based DPI formulations of various drugs; however, contradictory reports exist regarding the effect of the size of the carrier on the dispersion of drug. In this study we examined the influence of the intrinsic particle size of the polymeric carrier on the aerosolization of a model drug salbutamol sulphate (SS). Methods Four different sizes (20–150 lm) of polymer carriers were fabricated using solvent evaporation technique and the dispersion of SS particles from these carriers was measured by a Twin Stage Impinger (TSI). The size and morphological properties of polymer carriers were by laser diffraction and SEM, respectively. Results The FPF from these carriers was found to be increasing from 5.6% to 21.3% with increasing the carrier size. The FPF was found to be greater (21%) with the highest particle size of the carrier (150 lm). Conclusions The aerosolization of drug was dependent on the size of polymer carriers. The smaller size of the carrier resulted in lower FPF which was increased with increasing the carrier size. For a fixed mass of drug particles in a formulation, the mass of drug particles per unit area of carriers is higher in formulations containing the larger carriers, which leads to an increase in the dispersion of drug due to the increased mechanical forces occurred between the carriers and the device walls.

Dispersion properties of surface waves in a two-layer plasma structure bounded by a metal

The effect of the nonuniformity of the electron density on the dispersion properties of surface waves propagating in a direction transverse to an external magnetic field is studied for the model of a two-layer plasma structure bounded by a metal. It is shown that the spectra of the waves can be effectively controlled by varying the degree of nonuniformity of the density and the dimensions of the layers.

An investigation of Nano-particle deposition in cylindrical tubes under Laminar condition using Lagrangian transport model

Aerosol deposition in cylindrical tubes is a subject of interest to researchers and engineers in many applications of aerosol physics and metrology. Investigation of nano-particles in different aspects such as lungs, upper airways, batteries and vehicle exhaust gases is vital due the smaller size, adverse health effect and higher trouble for trapping than the micro-particles. The Lagrangian particle tracking provides an effective method for simulating the deposition of nano-particles as well as micro-particles as it accounts for the particle inertia effect as well as the Brownian excitation. However, using the Lagrangian approach for simulating ultrafine particles has been limited due to computational cost and numerical difficulties. In this paper, the deposition of nano-particles in cylindrical tubes under laminar condition is studied using the Lagrangian particle tracking method. The commercial Fluent software is used to simulate the fluid flow in the pipes and to study the deposition and dispersion of nano-particles. Different particle diameters as well as different flow rates are examined. The point analysis in a uniform flow is performed for validating the Brownian motion. The results show good agreement between the calculated deposition efficiency and the analytic correlations in the literature. Furthermore, for the nano-particles with the diameter more than 40 nm, the calculated deposition efficiency by the Lagrangian method is less than the analytic correlations based on Eulerian method due to statistical error or the inertia effect.

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.

An approximate method for the solution of an influence function foil rolling model

Fleck and Johnson (Int. J. Mech. Sci. 29 (1987) 507) and Fleck et al. (Proc. Inst. Mech. Eng. 206 (1992) 119) have developed foil rolling models which allow for large deformations in the roll profile, including the possibility that the rolls flatten completely. However, these models require computationally expensive iterative solution techniques. A new approach to the approximate solution of the Fleck et al. (1992) Influence Function Model has been developed using both analytic and approximation techniques. The numerical difficulties arising from solving an integral equation in the flattened region have been reduced by applying an Inverse Hilbert Transform to get an analytic expression for the pressure. The method described in this paper is applicable to cases where there is or there is not a flat region.