Numerical study of oxygen diffusion from capillary to tissues during hypoxia with external force effects

A mathematical model to simulate oxygen delivery through a capillary to tissues under the influence of an external force field is presented. The multi-term general fractional diffusion equation containing force terms and a time dependent absorbent term is taken into account. Fractional calculus is applied to describe the phenomenon of sub-diffusion of oxygen in both transverse and longitudinal directions. A new computational algorithm, i.e., the new iterative method (NIM) is employed to solve the spatio-temporal fractional partial differential equation subject to appropriate physical boundary conditions. Validation of NIM solutions is achieved with a modified Adomian decomposition method (MADM). A parametric study is conducted for three loading scenarios on the capillary-radial force alone, axial force alone and the combined case of both forces. The results demonstrate that the force terms markedly influence the oxygen diffusion process. For example, the radial force exerts a more profound effect than axial force on sub-diffusion of oxygen indicating that careful manipulation of these forces on capillary tissues may assist in the effective reduction of hypoxia or other oxygen depletion phenomena

Solutions of System of Fractional Partial Differential Equations

In this paper, system of fractional partial differential equation which has numerous applications in many fields of science is considered. Adomian decomposition method, a novel method is used to solve these type of equations. The solutions are derived in convergent series form which shows the effectiveness of the method for solving wide variety of fractional differential equations

An alternating direction implicit spectral method for solving two dimensional multi-term time fractional mixed diffusion and diffusion-wave equations

In this paper, we consider the initial boundary value problem of the two dimensional multi-term time fractional mixed diffusion and diffusion-wave equations. An alternating direction implicit (ADI) spectral method is developed based on Legendre spectral approximation in space and finite difference discretization in time. Numerical stability and convergence of the schemes are proved, the optimal error is $O(N^{-r}+\tau^2)$, where $N, \tau, r$ are the polynomial degree, time step size and the regularity of the exact solution, respectively. We also consider the non-smooth solution case by adding some correction terms. Numerical experiments are presented to confirm our theoretical analysis. These techniques can be used to model diffusion and transport of viscoelastic non-Newtonian fluids

On adomian based numerical schemes for euler and navier-stokes equations, and application to aeroacoustic propagation

140 p.En esta tesis se ha desarrollado un nuevo método de integración en tiempo de tipo derivadas sucesivas (multiderivative), llamado ABS y basado en el algoritmo de Adomian. Su motivación radica en la reducción del coste de simulación para problemas en aeroacústica, muy costosos por su naturaleza transitoria y requisitos de alta precisión. El método ha sido satisfactoriamente empleado en ambas partes de un sistema híbrido, donde se distinguen la parte aerodinámica y la acústica.En la parte aerodinámica las ecuaciones de Navier-Stokes incompresibles son resueltas con unmétodo de proyección clásico. Sin embargo, la fase de predicción de velocidad ha sido modificadapara incluir el método ABS en combinación con dos métodos: una discretización espacial MAC devolúmenes finitos, y también con un método de alto orden basado en ADER. El método se ha validado respecto a los problemas (en 2D) de vórtices de Taylor-Green, y el desarrollo de vórticesde Karman en un cilindro cuadrado. La parte acústica resuelve la propagación de ondas descritaspor las ecuaciones linearizadas de Euler, empleando una discretización de Galerkin discontinua. El método se ha validado respecto a la ecuación de Burgers.El método ABS es sencillo de programar con una formulación recursiva. Los resultados demuestran que su sencillez junto con sus altas capacidades de adaptación lo convierten en un método fácilmente extensible a órdenes altos, a la vez que reduce el coste comparado con otros métodos clásicos

On Adomian Based Numerical Schemes for Euler and Navier-Stokes Equations, and Application to Aeroacoustic Propagation

In this thesis, an Adomian Based Scheme (ABS) for the compressible Navier-Stokes equations is constructed, resulting in a new multiderivative type scheme not found in the context of fluid dynamics. Moreover, this scheme is developed as a means to reduce the computational cost associated with aeroacoustic simulations, which are unsteady in nature with high-order requirements for the acoustic wave propagation. We start by constructing a set of governing equations for the hybrid computational aeroacoustics method, splitting the problem into two steps: acoustic source computation and wave propagation. The first step solves the incompressible Navier-Stokes equation using Chorin's projection method, which can be understood as a prediction-correction method. First, the velocity prediction is obtained solving the viscous Burgers' equation. Then, its divergence-free correction is performed using a pressure Poisson type projection. In the velocity prediction substep, Burgers' equation is solved using two ABS variants: a MAC type implementation, and a ``modern'' ADER method. The second step in the hybrid method, related to wave propagation, is solved combining ABS with the discontinuous Galerkin high-order approach. Described solvers are validated against several test cases: vortex shedding and Taylor-Green vortex problems for the first step, and a Gaussian wave propagation in the second case. Although ABS is a multiderivative type scheme, it is easily programmed with an elegant recursive formulation, even for the general Navier-Stokes equations. Results show that its simplicity combined with excellent adaptivity capabilities allows for a successful extension to very high-order accuracy at relatively low cost, obtaining considerable time savings in all test cases considered.Predoc Gobierno Vasc

Collocation Method using Compactly Supported Radial Basis Function for Solving Volterra's Population Model

In this paper, indirect collocation approach based on compactly supported radial basis function is applied for solving Volterras population model. The method reduces the solution of this problem to the solution of a system of algebraic equations. Volterras model is a non-linear integro-differential equation where the integral term represents the effect of toxin. To solve the problem, we use the well-known CSRBF: Wendland3,5. Numerical results and residual norm 2 show good accuracy and rate of convergence.Comment: 8 pages , 1 figure. arXiv admin note: text overlap with arXiv:1008.233

Numerical solution of fractional partial differential equations by spectral methods

Fractional partial differential equations (FPDEs) have become essential tool for the modeling of physical models by using spectral methods. In the last few decades, spectral methods have been developed for the solution of time and space dimensional FPDEs. There are different types of spectral methods such as collocation methods, Tau methods and Galerkin methods. This research work focuses on the collocation and Tau methods to propose an efficient operational matrix methods via Genocchi polynomials and Legendre polynomials for the solution of two and three dimensional FPDEs. Moreover, in this study, Genocchi wavelet-like basis method and Genocchi polynomials based Ritz- Galerkin method have been derived to deal with FPDEs and variable- order FPDEs. The reason behind using the Genocchi polynomials is that, it helps to generate functional expansions with less degree and small coefficients values to derive the operational matrix of derivative with less computational complexity as compared to Chebyshev and Legendre Polynomials. The results have been compared with the existing methods such as Chebyshev wavelets method, Legendre wavelets method, Adomian decomposition method, Variational iteration method, Finite difference method and Finite element method. The numerical results have revealed that the proposed methods have provided the better results as compared to existing methods due to minimum computational complexity of derived operational matrices via Genocchi polynomials. Additionally, the significance of the proposed methods has been verified by finding the error bound, which shows that the proposed methods have provided better approximation values for under consideration FPDEs

Homotopy analysis method for solving multi-term linear and nonlinear diffusion–wave equations of fractional order

AbstractIn this paper we have used the homotopy analysis method (HAM) to obtain solutions of multi-term linear and nonlinear diffusion–wave equations of fractional order. The fractional derivative is described in the Caputo sense. Some illustrative examples have been presented