Joint additive Kullback–Leibler residual minimization and regularization for linear inverse problems

For the approximate solution of ill‐posed inverse problems, the formulation of a regularization functional involves two separate decisions: the choice of the residual minimizer and the choice of the regularizor. In this paper, the Kullback–Leibler functional is used for both. The resulting regularization method can solve problems for which the operator and the observational data are positive along with the solution, as occur in many inverse problem applications. Here, existence, uniqueness, convergence and stability for the regularization approximations are established under quite natural regularity conditions. Convergence rates are obtained by using an a priori strategy. Copyright © 2007 John Wiley & Sons, Ltd.     

The quasi‐reversibility method for a final value problem of the time‐fractional diffusion equation with inhomogeneous source

This paper is devoted to discuss a multidimensional backward heat conduction problem for time‐fractional diffusion equation with inhomogeneous source. This problem is ill‐posed. We use quasi‐reversibility regularization method to solve this inverse problem. Moreover, the convergence estimates between regularization solution and the exact solution are obtained under the a priori and the a posteriori choice rules. Finally, the numerical examples for one‐dimensional and two‐dimensional cases are presented to show that our method is feasible and effective.     

Identification of a time-dependent source term for a time fractional diffusion problem

In this paper, we study an inverse problem of identifying a time-dependent term of an unknown source for a time fractional diffusion equation using nonlocal measurement data. Firstly, we establish the conditional stability for this inverse problem. Then two regularization methods are proposed to for reconstructing the time-dependent source term from noisy measurements. The first method is an integral equation method which formulates the inverse source problem into an integral equation of the second kind; and a prior convergence rate of regularized solutions is derived with a suitable choice strategy of regularization parameters. The second method is a standard Tikhonov regularization method and formulates the inverse source problem as a minimizing problem of the Tikhonov functional. Based on the superposition principle and the technique of finite-element interpolation, a numerical scheme is proposed to implement the second regularization method. One- and two-dimensional examples are carried out to verify efficiency and stability of the second regularization method.     

A backward problem for the time‐fractional diffusion equation

In this paper, we are concerned with the backward problem of reconstructing the initial condition of a time‐fractional diffusion equation from interior measurements. We establish uniqueness results and provide stability analysis. Our method is based on the eigenfunction expansion of the forward solution and the Tikhonov regularization to tackle the ill‐posedness issue of the underlying inverse problem. Some numerical examples are included to illustrate the effectiveness of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

A modified Chebyshev ϑ‐weighted Crank–Nicolson method for analyzing fractional sub‐diffusion equations

This study presents a robust modification of Chebyshev ? ‐weighted Crank–Nicolson method for analyzing the sub‐diffusion equations in the Caputo fractional sense. In order to solve the problem, by discretization of the sub‐fractional diffusion equations using Taylor's expansion a linear system of algebraic equations that can be analyzed by numerical methods is presented. Furthermore, consistency, convergence, and stability analysis of the suggested method are discussed. In this framework, compact structures of sub‐diffusion equations are considered as prototype examples. The main advantage of the proposed method is that, it is more efficient in terms of CPU time, computational cost and accuracy in comparing with the existing ones in open literature.     

Polynomial spectral collocation method for space fractional advection–diffusion equation

This article discusses the spectral collocation method for numerically solving nonlocal problems: one‐dimensional space fractional advection–diffusion equation; and two‐dimensional linear/nonlinear space fractional advection–diffusion equation. The differentiation matrixes of the left and right Riemann–Liouville and Caputo fractional derivatives are derived for any collocation points within any given bounded interval. Several numerical examples with different boundary conditions are computed to verify the efficiency of the numerical schemes and confirm the exponential convergence; the physical simulations for Lévy–Feller advection–diffusion equation and space fractional Fokker–Planck equation with initial δ‐peak and reflecting boundary conditions are performed; and the eigenvalue distributions of the iterative matrix for a variety of systems are displayed to illustrate the stabilities of the numerical schemes in more general cases. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 514–535, 2014     

A linear,decoupled fractional time‐stepping method for the nonlinear fluid–fluid interaction

In this paper, a linear decoupled fractional time stepping method is proposed and developed for the nonlinear fluid–fluid interaction governed by the two Navier–Stokes equations. Partitioned time stepping method is applied to two‐physics problems with stiffness of the coupling terms being treated explicitly and is also unconditionally stable. As for each fluid, the velocity and pressure are respectively determined by just solving one vector‐valued quasi‐elliptic equation and the Possion equation with homogeneous Neumann boundary condition per time step. Therefore, the cost of the fluid–fluid interaction is dominant to solve four simple linear equations, which greatly reduces the computational cost of the whole system. The method exploits properties of the fluid–fluid system to establish its stability and convergence with the same results as the standard scheme. Finally, numerical experiments are presented to show the performance of the proposed method.     

A second‐order finite difference method for fractional diffusion equation with Dirichlet and fractional boundary conditions

In this paper, we develop a practical numerical method to approximate a fractional diffusion equation with Dirichlet and fractional boundary conditions. An approach based on the classical Crank–Nicolson method combined with spatial extrapolation is used to obtain temporally and spatially second‐order accurate numerical estimates. The solvability, stability, and convergence of the proposed numerical scheme are proved via the Gershgorin theorem. Numerical experiments are performed to confirm the accuracy and efficiency of our scheme.     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

An inverse problem for space‐fractional backward diffusion problem

In this paper, an inverse problem for space‐fractional backward diffusion equation, which is highly ill‐posed, is considered. This problem is obtained from the classical diffusion equation by replacing the second‐order space derivative with a Riesz–Feller derivative of order α ∈ (0,2]. We show that such a problem is severely ill‐posed, and further present a simplified Tikhonov regularization method to deal with this problem. Convergence estimate is presented under a priori choice of regularization parameter. Numerical experiments are given to illustrate the accuracy and efficiency of the proposed method. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical study of the SVD–MFS solution of inverse boundary value problems in two‐dimensional steady‐state linear thermoelasticity

We study the reconstruction of the missing thermal and mechanical data on an inaccessible part of the boundary in the case of two‐dimensional linear isotropic thermoelastic materials from overprescribed noisy measurements taken on the remaining accessible boundary part. This inverse problem is solved by using the method of fundamental solutions together with the method of particular solutions. The stabilization of this inverse problem is achieved using several singular value decomposition (SVD)‐based regularization methods, such as the Tikhonov regularization method (Tikhonov and Arsenin, Methods for solving ill‐posed problems, Nauka, Moscow, 1986), the damped SVD and the truncated SVD (Hansen, Rank‐deficient and discrete ill‐posed problems: numerical aspects of linear inversion, SIAM, Philadelphia, 1998), whilst the optimal regularization parameter is selected according to the discrepancy principle (Morozov, Sov Math Doklady 7 (1966), 414–417), generalized cross‐validation criterion (Golub et al. Technometrics 22 (1979), 1–35) and Hansen's L‐curve method (Hansen and O'Leary, SIAM J Sci Comput 14 (1993), 1487–503). © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 168–201, 2015     

A backward euler orthogonal spline collocation method for the time‐fractional Fokker–Planck equation

We formulate and analyze a novel numerical method for solving a time‐fractional Fokker–Planck equation which models an anomalous subdiffusion process. In this method, orthogonal spline collocation is used for the spatial discretization and the time‐stepping is done using a backward Euler method based on the L1 approximation to the Caputo derivative. The stability and convergence of the method are considered, and the theoretical results are supported by numerical examples, which also exhibit superconvergence. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1534–1550, 2015     

Vieta–Fibonacci wavelets: Application in solving fractional pantograph equations

In this paper, the Vieta–Fibonacci wavelets as a new family of orthonormal wavelets are generated. An operational matrix concerning fractional integration of these wavelets is extracted. A numerical scheme is established based on these wavelets and their fractional integral matrix together with the collocation technique to solve fractional pantograph equations. The presented method reduces solving the problem under study into solving a system of algebraic equations. Several examples are provided to show the accuracy of the method.     

Numerical solution of the Rosenau–KdV–RLW equation by operator splitting techniques based on B‐spline collocation method

In the present study, the operator splitting techniques based on the quintic B‐spline collocation finite element method are presented for calculating the numerical solutions of the Rosenau–KdV–RLW equation. Two test problems having exact solutions have been considered. To demonstrate the efficiency and accuracy of the present methods, the error norms L₂ and L_∞ with the discrete mass Q and energy E conservative properties have been calculated. The results obtained by the method have been compared with the exact solution of each problem and other numerical results in the literature, and also found to be in good agreement with each other. A Fourier stability analysis of each presented method is also investigated.     

Numerical and theoretical discussions for solving nonlinear generalized Benjamin–Bona–Mahony–Burgers equation based on the Legendre spectral element method

This article is devoted to solving numerically the nonlinear generalized Benjamin–Bona–Mahony–Burgers (GBBMB) equation that has several applications in physics and applied sciences. First, the time derivative is approximated by using a finite difference formula. Afterward, the stability and convergence analyses of the obtained time semi‐discrete are proven by applying the energy method. Also, it has been demonstrated that the convergence order in the temporal direction is O(dt) . Second, a fully discrete formula is acquired by approximating the spatial derivatives via Legendre spectral element method. This method uses Lagrange polynomial based on Gauss–Legendre–Lobatto points. An error estimation is also given in detail for full discretization scheme. Ultimately, the GBBMB equation in the one‐ and two‐dimension is solved by using the proposed method. Also, the calculated solutions are compared with theoretical solutions and results obtained from other techniques in the literature. The accuracy and efficiency of the mentioned procedure are revealed by numerical samples.     

Variational formulations for Vlasov–Poisson–Fokker–Planck systems

Time‐discrete variational schemes are introduced for both the Vlasov–Poisson–Fokker–Planck (VPFP) system and a natural regularization of the VPFP system. The time step in these variational schemes is governed by a certain Kantorovich functional (or scaled Wasserstein metric). The discrete variational schemes may be regarded as discretized versions of a gradient flow, or steepest descent, of the underlying free energy functionals for these systems. For the regularized VPFP system, convergence of the variational scheme is rigorously established. Copyright © 2000 John Wiley & Sons, Ltd.     

Poisson方程热源识别反问题

探讨了半带状区域上二维Poisson方程只含有一个空间变量的热源识别反问题.这类问题是不适定的,即问题的解(如果存在的话)不连续依赖于测量数据.利用Carasso-Tikhonov正则化方法,得到了问题的一个正则近似解,并且给出了正则解和精确解之间具有Holder型误差估计.数值实验表明Carasso-Tikhonov正则化方法对于这种热源识别是非常有效的.     

Multi‐switching combination–combination synchronization of non‐identical fractional‐order chaotic systems

In this paper, multi‐switching combination–combination synchronization scheme has been investigated between a class of four non‐identical fractional‐order chaotic systems. The fractional‐order Lorenz and Chen's systems are taken as drive systems. The combination–combination of multi drive systems is then synchronized with the combination of fractional‐order Lü and Rössler chaotic systems. In multi‐switching combination–combination synchronization, the state variables of two drive systems synchronize with different state variables of two response systems simultaneously. Based on the stability of fractional‐order chaotic systems, the multi‐switching combination–combination synchronization of four fractional‐order non‐identical systems has been investigated. For the synchronization of four non‐identical fractional‐order chaotic systems, suitable controllers have been designed. Theoretical analysis and numerical results are presented to demonstrate the validity and feasibility of the applied method. Copyright © 2017 John Wiley & Sons, Ltd.     

Numerical solutions of distributed order fractional differential equations in the time domain using the Müntz–Legendre wavelets approach

In this paper, a numerical method is presented to obtain and analyze the behavior of numerical solutions of distributed order fractional differential equations of the general form in the time domain with the Caputo fractional derivative. The suggested method is based on the Müntz–Legendre wavelet approximation. We derive a new operational vector for the Riemann–Liouville fractional integral of the Müntz–Legendre wavelets by using the Laplace transform method. Applying this operational vector and collocation method in our approach, the problem can be reduced to a system of linear and nonlinear algebraic equations. The arising system can be solved by the Newton method. Discussion on the error bound and convergence analysis for the proposed method is presented. Finally, seven test problems are considered to compare our results with other well‐known methods used for solving these problems. The results in the tabulated tables highlighted that the proposed method is an efficient mathematical tool for analyzing distributed order fractional differential equations of the general form.     

A three‐level linearized finite difference scheme for the camassa–holm equation

The Camassa–Holm (CH) system is a strong nonlinear third‐order evolution equation. So far, the numerical methods for solving this problem are only a few. This article deals with the finite difference solution to the CH equation. A three‐level linearized finite difference scheme is derived. The scheme is proved to be conservative, uniquely solvable, and conditionally second‐order convergent in both time and space in the discrete L_∞ norm. Several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 451–471, 2014