Simplest equation method for some time-fractional partial differential equations with conformable derivative

The conformable fractional derivative was proposed by R. Khalil et al. in 2014, which is natural and obeys the Leibniz rule and chain rule. Based on the properties, a class of time-fractional partial differential equations can be reduced into ODEs using traveling wave transformation. Then the simplest equation method is applied to find exact solutions of some time-fractional partial differential equations. The exact solutions (solitary wave solutions, periodic function solutions, rational function solutions) of time-fractional generalized Burgers equation, time-fractional generalized KdV equation, time-fractional generalized Sharma–Tasso–Olver (FSTO) equation and time-fractional fifth-order KdV equation, $(3+1)$-dimensional time-fractional KdV–Zakharov–Kuznetsov (KdV–ZK) equation are constructed. This method presents a wide applicability to solve some nonlinear time-fractional differential equations with conformable derivative.     

Homotopy perturbation method for fractional Fornberg–Whitham equation

This article presents the approximate analytical solutions to solve the nonlinear Fornberg–Whitham equation with fractional time derivative. By using initial values, the explicit solutions of the equations are solved by using a reliable algorithm like homotopy perturbation method. The fractional derivatives are taken in the Caputo sense. Numerical results show that the HPM is easy to implement and accurate when applied to time-fractional PDEs.     

Variable-coefficient auxiliary equation method for exact solutions of non-linear evolution equations

In this paper, a variable-coefficient auxiliary equation method is proposed to seek more general exact solutions of non-linear evolution equations. Being concise and straightforward, this method is applied to the Kawahara equation, Sawada–Kotera equation and (2+1)-dimensional Korteweg–de Vries equations. As a result, many new and more general exact solutions are obtained including Jacobi elliptic, hyperbolic and trigonometric function solutions. It is shown that the proposed method provides a straightforward and effective method for non-linear evolution equations in mathematical physics.     

分数阶非线性方程精确解的新方法

将分数阶复变换方法和[(G/G)]方法相结合得到了一种辅助方程方法,用来求解分数阶非线性微分方程。利用该方法并借助于软件Mathematica的符号计算功能求解了分数阶Calogero KDV方程,得到了该方程新的精确解。     

Particle tracking for fractional diffusion with two time scales

Previous work [Y. Zhang, M.M. Meerschaert, B. Baeumer, Particle tracking for time-fractional diffusion, Phys. Rev. E 78 (2008) 036705] showed how to solve time-fractional diffusion equations by particle tracking. This paper extends the method to the case where the order of the fractional time derivative is greater than one. A subordination approach treats the fractional time derivative as a random time change of the corresponding Cauchy problem, with a first derivative in time. One novel feature of the time-fractional case of order greater than one is the appearance of clustering in the operational time subordinator, which is non-Markovian. Solutions to the time-fractional equation are probability densities of the underlying stochastic process. The process models movement of individual particles. The evolution of an individual particle in both space and time is captured in a pair of stochastic differential equations, or Langevin equations. Monte Carlo simulation yields particle location, and the ensemble density approximates the solution to the variable coefficient time-fractional diffusion equation in one or several spatial dimensions. The particle tracking code is validated against inverse transform solutions in the simplest cases. Further applications solve model equations for fracture flow, and upscaling flow in complex heterogeneous porous media. These variable coefficient time-fractional partial differential equations in several dimensions are not amenable to solution by any alternative method, so that the grid-free particle tracking approach presented here is uniquely appropriate.     

Burgers方程的精确解

引入一个变换,将二阶非线性偏微分方程—Burgers方程降阶为一阶的非线性方程,再直接求解该方程,得出了Burgers方程精确解的新形式,并与已有结果完全吻合.这种方法也适合于求解其他非线性偏微分方程.     

The Unstructured Mesh Finite Element Method for the Two-Dimensional Multi-term Time–Space Fractional Diffusion-Wave Equation on an Irregular Convex Domain

In this paper, the two-dimensional multi-term time-space fractional diffusion-wave equation on an irregular convex domain is considered as a much more general case for wider applications in fluid mechanics. A novel unstructured mesh finite element method is proposed for the considered equation. In most existing works, the finite element method is applied on regular domains using uniform meshes. The case of irregular convex domains, which would require subdivision using unstructured meshes, is mostly still open. Furthermore, the orders of the multi-term time-fractional derivatives have been considered to belong to (0, 1] or (1, 2] separately in existing models. In this paper, we consider two-dimensional multi-term time-space fractional diffusion-wave equations with the time fractional orders belonging to the whole interval (0, 2) on an irregular convex domain. We propose to use a mixed difference scheme in time and an unstructured mesh finite element method in space. Detailed implementation and the stability and convergence analyses of the proposed numerical scheme are given. Numerical examples are conducted to evaluate the theoretical analysis.     

New analytical method for gas dynamics equation arising in shock fronts

This work suggests a new analytical technique called the fractional homotopy analysis transform method (FHATM) for solving nonlinear homogeneous and nonhomogeneous time-fractional gas dynamics equations. The FHATM is an innovative adjustment in Laplace transform algorithm (LTA) and makes the calculation much simpler. The proposed technique solves the nonlinear problems without using Adomian polynomials and He’s polynomials which can be considered as a clear advantage of this new algorithm over decomposition and the homotopy perturbation transform method. In this paper, it can be observed that the auxiliary parameter ?$?$, which controls the convergence of the HATM approximate series solutions, also can be used in predicting and calculating multiple solutions. This is a basic and more qualitative difference in analysis between HATM and other methods. The solutions obtained by the proposed method indicate that the approach is easy to implement and computationally very attractive. The proposed method is illustrated by solving some numerical examples.     

Well-posedness,optimal control and discretization for time-fractional parabolic equations with time-dependent coefficients on metric graphs

We study distributed optimal control problems governed by time-fractional parabolic equations with time dependent coefficients on metric graphs, where the fractional derivative is considered in the Caputo sense. Using the Galerkin method and compactness results, for the spatial part, and approximating the kernel of the time-fractional Caputo derivative by a sequence of more regular kernel functions, we first prove the well-posedness of the system. We then turn to the existence and uniqueness of solutions to the distributed optimal control problem. By means of the Lagrange multiplier method, we develop an adjoint calculus for the right Caputo derivative and derive the corresponding first order optimality system. Moreover, we propose a finite difference scheme to find the approximate solution of the state equation and the resulting optimality system on metric graphs. Finally, examples are provided on two different graphs to illustrate the performance of the proposed difference scheme.     

Exact dynamic and static element stiffness matrices of nonsymmetric thin-walled beam-columns

An improved numerical method to exactly evaluate 14 × 14 dynamic and static element stiffness matrices is proposed for the spatial free vibration and stability analysis of nonsymmetric thin-walled straight beams subjected to eccentrically axial loads. Firstly equations of motion and force-deformation relations are rigorously derived from the total potential energy for a uniform beam element with nonsymmetric thin-walled cross-section. Next a system of linear algebraic equations with nonsymmetric matrices is constructed by introducing 14 displacement parameters and transforming the higher order simultaneous differential equation into the first order simultaneous equation. And then explicit expressions for displacement parameters are exactly evaluated by solving a generalized eigenproblem with complex eigenvalues. Finally exact element stiffness matrices are determined using force-deformation relations. Particularly straightforward application of the present method may not give the exact static stiffness because of existence of multiple zero eigenvalues in case of static buckling problems. Accordingly, a modified numerical method to resolve this difficulty is developed for two cases depending on the initial state of stress resultants. In order to demonstrate the validity and the accuracy of this method, the natural frequencies and buckling loads of nonsymmetric thin-walled beam-columns having bending-torsional deformation modes are evaluated and compared with analytical and F.E. solutions or results analyzed by ABAQUS’s shell element.     

求一类非线性偏微分方程精确解的简化试探函数法

利用试探函数法,将一个难于求解的非线性偏微分方程化为一个易于求解的代数方程,然后用待定系数法确定相应的常数,简洁地求得了一类非线性偏微分方程的精确解．将此方法应用到Burgers方程、KdV方程和KdV—Burgers方程,所得结果与已有结果完全吻合．本方法可望进一步推广用于求解其它非线性偏微分方程．     

He's homotopy perturbation method for solving the space- and time-fractional telegraph equations

In this study, homotopy perturbation method (HPM) is used to obtain analytic and approximate solutions of the space- and time-fractional telegraph equations. The space- and time-fractional derivatives are considered in the Caputo sense. The analytic solutions are calculated in the form of a series with easily computable terms. Some examples are given. The results reveal that HPM is very effective and convenient.     

The non-standard finite difference scheme for linear fractional PDEs in fluid mechanics

A non-standard finite difference scheme is developed to solve the linear partial differential equations with time- and space-fractional derivatives. The Grunwald–Letnikov method is used to approximate the fractional derivatives. Numerical illustrations that include the linear inhomogeneous time-fractional equation, linear space-fractional telegraph equation, linear inhomogeneous fractional Burgers equation and the fractional wave equation are investigated to show the pertinent features of the technique. Numerical results are presented graphically and reveal that the non-standard finite difference scheme is very effective and convenient for solving linear partial differential equations of fractional order.     

On the numerical treatment of an ordinary differential equation arising in one-dimensional non-Fickian diffusion problems

In a recent study, Chen and Liu [Comput. Phys. Comm. 150 (2003) 31] considered a one-dimensional, linear non-Fickian diffusion problem with a potential field, which, upon application of the Laplace transform, resulted in a second-order linear ordinary differential equation which was solved by means of a control-volume finite difference method that employs exponential shape functions. It is first shown that this formulation does not properly account for the spatial dependence of the drift forces and results in oscillatory solutions near the left boundary when these forces are large. A piecewise linearized method that provides piecewise analytical solutions, is exact in exact arithmetic for constant coefficients, homogeneous, second-order linear ordinary differential equations and results in three-point finite difference equations is then proposed. Numerical simulations indicate that the piecewise linearized method is free from unphysical oscillations and more accurate than that of Chen and Liu, especially for large drift forces. The method is then applied to non-Fickian diffusion problems with non-constant drift forces in order to determine the effects of the potential field on the concentration distribution.     

Numerical Solution of the Distributed-Order Fractional Bagley-Torvik Equation

In this paper, two numerical methods are proposed for solving distributed-order fractional Bagley-Torvik equation. This equation is used in modeling the motion of a rigid plate immersed in a Newtonian fluid with respect to the nonnegative density function. Using the composite Boole's rule the distributed-order Bagley-Torvik equation is approximated by a multi-term time-fractional equation, which is then solved by the Grunwald-Letnikov method (GLM) and the fractional differential transform method (FDTM). Finally, we compared our results with the exact results of some cases and show the excellent agreement between the approximate result and the exact solution.      

Analytical solution of gas lubricated slider microbearing

The analytical solution of a two-dimensional, isothermal, compressible gas flow in a slider microbearing is presented. A higher order accuracy of the solution is achieved by applying the boundary condition of Kn ² order for the velocity slip on the wall, together with the momentum equation of the same order (known as the Burnett equation). The analytical solution is obtained by the perturbation analysis. The order of all terms in continuum and momentum equations and in boundary conditions is evaluated by incorporating the exact relation between the Mach, Reynolds and Knudsen numbers in the modelling procedure. Low Mach number flows in microbearing with slowly varying cross-sections are considered, and it is shown that under these conditions the Burnett equation has the same form as the Navier–Stokes equation. Obtained analytical results for pressure distribution, load capacity and velocity field are compared with numerical solutions of the Boltzmann equation and some semi-analytical results, and excellent agreement is achieved. The model presented in this paper is a useful tool for the prediction of flow conditions in the microbearings. Also, its results are the benchmark test for the verifications of various numerical procedures.     

Piecewise quasilinearization techniques for singular boundary-value problems

Piecewise quasilinearization methods for singular boundary-value problems in second-order ordinary differential equations are presented. These methods result in linear constant-coefficients ordinary differential equations which can be integrated analytically, thus yielding piecewise analytical solutions. The accuracy of the globally smooth piecewise quasilinear method is assessed by comparisons with exact solutions of several Lane-Emden equations, a singular problem of non-Newtonian fluid dynamics and the Thomas-Fermi equation. It is shown that the smooth piecewise quasilinearization method provides accurate solutions even near the singularity and is more precise than (iterative) second-order accurate finite difference discretizations. It is also shown that the accuracy of the smooth piecewise quasilinear method depends on the kind of singularity, nonlinearity and inhomogeneities of singular ordinary differential equations. For the Thomas-Fermi equation, it is shown that the piecewise quasilinearization method that provides globally smooth solutions is more accurate than that which only insures global continuity, and more accurate than global quasilinearization techniques which do not employ local linearization.     

The Riccati equation with variable coefficients expansion algorithm to find more exact solutions of nonlinear differential equations

In this paper based on a system of Riccati equations with variable coefficients, we present a new Riccati equation with variable coefficients expansion method and its algorithm, which are direct and more powerful than the tanh-function method, sine-cosine method, the generalized hyperbolic-function method and the generalized Riccati equation with constant coefficient expansion method to construct more new exact solutions of nonlinear differential equations in mathematical physics. A pair of generalized Hamiltonian equations is chosen to illustrate our algorithm such that more families of new exact solutions are obtained which contain soliton-like solution and periodic solutions. This algorithm can also be applied to other nonlinear differential equations.     

Fully Petrov–Galerkin spectral method for the distributed-order time-fractional fourth-order partial differential equation

Distributed fractional derivative operators can be used for modeling of complex multiscaling anomalous transport, where derivative orders are distributed over a range of values rather than being just a fixed integer number. In this paper, we consider the space-time Petrov–Galerkin spectral method for a two-dimensional distributed-order time-fractional fourth-order partial differential equation. By applying a proper Gauss-quadrature rule to discretize the distributed integral operator, the problem is converted to a multi-term time-fractional equation. Then, the proposed method for solving the obtained equation is based on using Jacobi polyfractonomial, which are eigenfunctions of the first kind fractional Sturm–Liouville problem (FSLP), as temporal basis and Legendre polynomials for the spatial discretization. The eigenfunctions of the second kind FSLP are used as temporal basis in test space. This approach leads to finding the numerical solution of the problem through solving a system of linear algebraic equations. Finally, we provide some examples with smooth solutions and finite regular solutions to numerically demonstrate the efficiency, accuracy, and exponential convergence of the proposed method.