A meshless method for solving three-dimensional time fractional diffusion equation with variable-order derivatives

In this study a new framework for solving three-dimensional (3D) time fractional diffusion equation with variable-order derivatives is presented. Firstly, a θ-weighted finite difference scheme with second-order accuracy is introduced to perform temporal discretization. Then a meshless generalized finite difference (GFD) scheme is employed for the solutions of remaining problems in the space domain. The proposed scheme is truly meshless and can be used to solve problems defined on an arbitrary domain in three dimensions. Preliminary numerical examples illustrate that the new method proposed here is accurate and efficient for time fractional diffusion equation in three dimensions, particularly when high accuracy is desired.     

The method of approximate particular solutions to simulate an anomalous mobile-immobile transport process

In the current work, a generalized mathematical model based on the Coimbra time fractional derivative of variable order, which describes an anomalous mobile-immobile transport process in complex systems is investigated numerically. A robust numerical technique based on the meshfree strong form method combined with an efficient time-stepping scheme is performed to compute the approximate solution of the problem with high accuracy. For this purpose, firstly, an effective implicit time discretization approach is used for discretizing the variable-order time fractional problem in the time direction. Then a global meshless technique based on the method of approximate particular solutions is performed to fully discretize the model in the spatial domain. The validity and performance of the procedure to numerically simulate the proposed generalized solute transport model on regular and irregular domains are demonstrated through some numerical examples.     

The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation

In this article we describe a numerical method to solve a nonhomogeneous diffusion equation with arbitrary geometry by combining the method of fundamental solutions (MFS), the method of particular solutions (MPS), and the eigenfunction expansion method (EEM). This forms a meshless numerical scheme of the MFS‐MPS‐EEM model to solve nonhomogeneous diffusion equations with time‐independent source terms and boundary conditions for any time and any shape. Nonhomogeneous diffusion equation with complex domain can be separated into a Poisson equation and a homogeneous diffusion equation using this model. The Poisson equation is solved by the MFS‐MPS model, in which the compactly supported radial basis functions are adopted for the MPS. On the other hand, utilizing the EEM the diffusion equation is first translated to a Helmholtz equation, which is then solved by the MFS together with the technique of the singular value decomposition (SVD). Since the present meshless method does not need mesh generation, nodal connectivity, or numerical integration, the computational effort and memory storage required are minimal as compared with other numerical schemes. Test results for two 2D diffusion problems show good comparability with the analytical solutions. The proposed algorithm is then extended to solve a problem with irregular domain and the results compare very well with solutions of a finite element scheme. Therefore, the present scheme has been proved to be very promising as a meshfree numerical method to solve nonhomogeneous diffusion equations with time‐independent source terms of any time frame, and for any arbitrary geometry. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Solitary wave solutions of the MRLW equation using radial basis functions

In this study, traveling wave solutions of the modified regularized long wave (MRLW) equation are simulated by using the meshless method based on collocation with well‐known radial basis functions. The method is tested for three test problems which are single solitary wave motion, interaction of two solitary waves and interaction of three solitary waves. Invariant values for all test problems are calculated, also L₂, L_∞ norms and values of the absolute error for single solitary wave motion are calculated. Numerical results by using the meshless method with different radial basis functions are presented. Figures of wave motions for all test problems are shown. Altogether, meshless methods with radial basis functions solve the MRLW equation very satisfactorily.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 235–247, 2012     

Numerical solution of two-dimensional time fractional cable equation with Mittag-Leffler kernel

The main motive of this article is to study the recently developed Atangana-Baleanu Caputo (ABC) fractional operator that is obtained by replacing the classical singular kernel by Mittag-Leffler kernel in the definition of the fractional differential operator. We investigate a novel numerical method for the nonlinear two-dimensional cable equation in which time-fractional derivative is of Mittag-Leffler kernel type. First, we derive an approximation formula of the fractional-order ABC derivative of a function t^k using a numerical integration scheme. Using this approximation formula and some properties of shifted Legendre polynomials, we derived the operational matrix of ABC derivative. In the author of knowledge, this operational matrix of ABC derivative is derived the first time. We have shown the efficiency of this newly derived operational matrix by taking one example. Then we solved a new class of fractional partial differential equations (FPDEs) by the implementation of this ABC operational matrix. The two-dimensional model of the time-fractional model of the cable equation is solved and investigated by this method. We have shown the effectiveness and validity of our proposed method by giving the solution of some numerical examples of the two-dimensional fractional cable equation. We compare our obtained numerical results with the analytical results, and we conclude that our proposed numerical method is feasible and the accuracy can be seen by error tables. We see that the accuracy is so good. This method will be very useful to investigate a different type of model that have Mittag-Leffler fractional derivative.     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method together with the fractional integral operator, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved, and the solution of fractional nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

Analysis of a meshless method for the time fractional diffusion-wave equation

In this paper a numerical technique is proposed for solving the time fractional diffusion-wave equation. We obtain a time discrete scheme based on finite difference formula. Then, we prove that the time discrete scheme is unconditionally stable and convergent using the energy method and the convergence order of the time discrete scheme is \(\mathcal {O}(\tau ^{3-\alpha })\). Firstly, we change the main problem based on Dirichlet boundary condition to a new problem based on Robin boundary condition and then, we consider a semi-discrete scheme with Robin boundary condition and show when \(\beta \rightarrow +\infty \) solution of the main semi-discrete problem with Dirichlet boundary condition is convergent to the solution of the new semi-discrete problem with Robin boundary condition. We consider the new semi-discrete problem with Robin boundary condition and use the meshless Galerkin method to approximate the spatial derivatives. Finally, we obtain an error bound for the new problem. We prove that convergence order of the numerical scheme based on Galekin meshless is \(\mathcal {O}(h)\). In the considered method the appeared integrals are approximated using Gauss Legendre quadrature formula. The main aim of the current paper is to obtain an error estimate for the meshless Galerkin method based on the radial basis functions. Numerical examples confirm the efficiency and accuracy of the proposed scheme.     

Analytical approximate solutions of Riesz fractional diffusion equation and Riesz fractional advection–dispersion equation involving nonlocal space fractional derivatives

In this paper, we consider the analytical solutions of fractional partial differential equations (PDEs) with Riesz space fractional derivatives on a finite domain. Here we considered two types of fractional PDEs with Riesz space fractional derivatives such as Riesz fractional diffusion equation (RFDE) and Riesz fractional advection–dispersion equation (RFADE). The RFDE is obtained from the standard diffusion equation by replacing the second‐order space derivative with the Riesz fractional derivative of order α∈(1,2]. The RFADE is obtained from the standard advection–dispersion equation by replacing the first‐order and second‐order space derivatives with the Riesz fractional derivatives of order β∈(0,1] and of order α∈(1,2] respectively. Here the analytic solutions of both the RFDE and RFADE are derived by using modified homotopy analysis method with Fourier transform. Then, we analyze the results by numerical simulations, which demonstrate the simplicity and effectiveness of the present method. Here the space fractional derivatives are defined as Riesz fractional derivatives. Copyright © 2015 John Wiley & Sons, Ltd.     

Spectral meshless radial point interpolation (SMRPI) method to two‐dimensional fractional telegraph equation

H. Ammari In this article, an innovative technique so‐called spectral meshless radial point interpolation (SMRPI) method is proposed and, as a test problem, is applied to a classical type of two‐dimensional time‐fractional telegraph equation defined by Caputo sense for (1 < α≤2). This new methods is based on meshless methods and benefits from spectral collocation ideas, but it does not belong to traditional meshless collocation methods. The point interpolation method with the help of radial basis functions is used to construct shape functions, which play as basis functions in the frame of SMRPI method. These basis functions have Kronecker delta function property. Evaluation of high‐order derivatives is not difficult by constructing operational matrices. In SMRPI method, it does not require any kind of integration locally or globally over small quadrature domains, which is essential of the finite element method (FEM) and those meshless methods based on Galerkin weak form. Also, it is not needed to determine strict value for the shape parameter, which plays an important role in collocation method based on the radial basis functions (Kansa's method). Therefore, computational costs of SMRPI method are less expensive. Two numerical examples are presented to show that SMRPI method has reliable rates of convergence. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of oblique and multidirectional wave propagation and breaking on steep slope based on FEM model of Boussinesq equations

Creating a representative numerical simulation of the propagation and breaking of waves along slopes is an important problem in engineering design. Most studies on wave breaking have focused on the propagation of normal incident waves on gentle slopes. In practice, however, waves on steep slopes are obliquely incident or multidirectional irregular waves. In this paper, the eddy viscosity term is introduced to the momentum equation of the improved Boussinesq equations to model wave dissipation caused by breaking and friction, and a numerical model based on an unstructured finite element method (FEM) is established based on the governing equations. It is applied to simulate wave propagation on a steep slope of 1:5. Parallel physical experiments are conducted for comparative analysis that considered a large number of cases, including those featuring of normal and oblique incident regular and irregular waves, and multidirectional waves. The heights of the incident wave increase for different periods to represent different kinds of waves breaking. Based on examination, the effectiveness and accuracy of the numerical model is verified through a comprehensive comparison between the numerical and the experimental results, including in terms of variation in wave height, wave spectrum, and nonlinear parameters. Satisfactory agreement between the numerical and experimental values shows that the proposed model is effective in representing the breaking of oblique incident regular waves, irregular waves, and multidirectional incident irregular waves. However, the initial threshold of the breaking parameter η_t^(I) takes different values for oblique and multidirectional waves. This needs to be paid attention when the breaking of waves is simulated using the Boussinesq equations.     

基于耦合纯无网格方法时间分数阶下孤立子波碰撞过程的数值模拟研究北大核心CSCD

为数值预测时间分数阶耦合非线性Schrödinger(TF-CNLS)方程描述的孤立子波非弹性碰撞过程,首次发展了一种耦合纯无网格有限点集法(coupled finite pointset method,CFPM).其构造过程为:1)对时间分数阶Caputo导数项采用一种高精度的差分格式;2)对空间导数采用基于Taylor展开和加权最小二乘法的有限粒子法(FPM)离散格式;3)对区域进行局部加密和采用稳定性好的双曲余弦核函数以提高数值精度.数值研究中,首先,运用CFPM对有解析解的一维TF-CNLS方程进行求解,分析了节点均匀分布或局部加密情况下的误差和收敛阶,表明给出的耦合无网格法具有近似二阶精度和易局部加密求解的灵活性;其次,运用CFPM对无解析解一维TF-CNLS方程描述的孤立子波非弹性碰撞过程进行了数值预测,其出现的波塌缩现象与整数阶下出现的多波现象截然不同;最后,与有限差分结果作对比,表明CFPM数值预测时间分数阶下孤立子波非弹性碰撞过程的复杂传播现象是可靠的.     

Numerical Solution of the Bagley-Torvik Equation

We consider the numerical solution of the Bagley-Torvik equation Ay(t) + BD _* ^3/2 y(t) + Cy(t) = f(t), as a prototype fractional differential equation with two derivatives. Approximate solutions have recently been proposed in the book and papers of Podlubny in which the solution obtained with approximate methods is compared to the exact solution. In this paper we consider the reformulation of the Bagley-Torvik equation as a system of fractional differential equations of order 1/2. This allows us to propose numerical methods for its solution which are consistent and stable and have arbitrarily high order. In this context we specifically look at fractional linear multistep methods and a predictor-corrector method of Adams type.     

Meshless local Petrov-Galerkin (MLPG) approximation to the two dimensional sine-Gordon equation

During the past few years, the idea of using meshless methods for numerical solution of partial differential equations (PDEs) has received much attention throughout the scientific community, and remarkable progress has been achieved on meshless methods. The meshless local Petrov-Galerkin (MLPG) method is one of the “truly meshless” methods since it does not require any background integration cells. The integrations are carried out locally over small sub-domains of regular shapes, such as circles or squares in two dimensions and spheres or cubes in three dimensions. In this paper the MLPG method for numerically solving the non-linear two-dimensional sine-Gordon (SG) equation is developed. A time-stepping method is employed to deal with the time derivative and a simple predictor-corrector scheme is performed to eliminate the non-linearity. A brief discussion is outlined for numerical integrations in the proposed algorithm. Some examples involving line and ring solitons are demonstrated and the conservation of energy in undamped SG equation is investigated. The final numerical results confirm the ability of proposed method to deal with the unsteady non-linear problems in large domains.     

Numerical simulation of a class of fractional subdiffusion equations via the alternating direction implicit method

In this article, a new numerical technique is proposed for solving the two‐dimensional time fractional subdiffusion equation with nonhomogeneous terms. After a transformation of the original problem, standard central difference approximation is used for the spatial discretization. For the time step, a new fractional alternating direction implicit (FADI) scheme based on the L₁ approximation is considered. This FADI scheme is constructed by adding a small term, so it is different from standard FADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Finally, numerical examples show the effectiveness and accuracy of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 531–547, 2016     

Local meshless method for PDEs arising from models of wound healing

A meshless collocation procedure is proposed for one- and two-dimensional partial differential equations arising from modeling of wound healing processes (Sherratt and Murray, 1991). Main motivation of this choice is its straightforward application in higher dimensions for both regular and irregular domains on various nodal points distributions. In the case of numerical solution of convection-dominated wound healing PDE models, a stencil based upwind stabilization technique is coupled with the local meshless method to counter instabilities of the computed solution. To assess efficacy, efficiency and accuracy of the proposed method on regular and irregular domains, numerical approximations of different wound healing models are obtained and validated against the exact solution and medically tested healing time duration.     

A hybrid radial basis functions collocation technique to numerically solve fractional advection–diffusion models

In this work, we propose a hybrid radial basis functions (RBFs) collocation technique for the numerical solution of fractional advection–diffusion models. In the formulation of hybrid RBFs (HRBFs), there exist shape parameter (c* ) and weight parameter (ϵ) that control numerical accuracy and stability. For these parameters, an adaptive algorithm is developed and validated. The proposed HRBFs method is tested for numerical solutions of some fractional Black–Sholes and diffusion models. Numerical simulations performed for several benchmark problems verified the proposed method accuracy and efficiency. The quantitative analysis is made in terms of L_∞, L₂, L_rms , and L_rel error norms as well as number of nodes N over space domain and time-step δt. Numerical convergence in space and time is also studied for the proposed method. The unconditional stability of the proposed HRBFs scheme is obtained using the von Neumann methodology. It is observed that the HRBFs method circumvented the ill-conditioning problem greatly, a major issue in the Kansa method.     

A novel attempt for finding comparatively accurate solution for sine‐Gordon equation comprising Riesz space fractional derivative

In this paper, a numerical procedure involving Chebyshev wavelet method has been implemented for computing the approximate solution of Riesz space fractional sine‐Gordon equation (SGE). Two‐dimensional Chebyshev wavelet method is implemented to calculate the numerical solution of space fractional SGE. The fractional SGE is considered as an interpolation between the classical SGE (corresponding to α = 2) and nonlocal SGE (corresponding to α = 1). As a consequence, the approximate solutions of fractional SGE obtained by using Chebyshev wavelet approach were compared with those derived by using modified homotopy analysis method with Fourier transform. Copyright © 2015 John Wiley & Sons, Ltd.     

Existence and uniqueness of minimizers for L2-constrained problems related to fractional Kirchhoff equation

In this paper, we discuss the existence and uniqueness of solutions of the constrained variational problem with respect to the fractional Kirchhoff equation. For the exponent p<p_*(s,N), a complete classification with respect to p for the existence of solutions of the fractional Kirchhoff functional on the L²-normalized manifold was given. Furthermore, all these solutions are unique up to translations, and our methods depend only on some simple energy estimates.     

关于薄板的无网格局部边界积分方程方法中的友解

无网格局部边界积分方程方法是最近发展起来的一种新的数值方法,这种方法综合了伽辽金有限元、边界元和无单元伽辽金法的优点,是一种具有广阔应用前景的、真正的无网格方法．把无网格局部边界积分方程方法应用于求解薄板问题,给出了薄板无网格局部边界积分方程方法所需要的友解及其全部公式．