Convergence analysis of the Jacobi spectral-collocation method for fractional integro-differential equations

We propose and analyze a spectral Jacobi-collocation approximation for fractional order integro-differential equations of Volterra type. The fractional derivative is described in the Caputo sense. We provide a rigorous error analysis for the collection method, which shows that the errors of the approximate solution decay exponentially in L^∞ norm and weighted L²-norm. The numerical examples are given to illustrate the theoretical results.     

MIXED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR ONE-DIMENSIONAL FULLY NONLINEAR SECOND ORDER ELLIPTIC AND PARABOLIC EQUATIONS

This paper is concerned with developing accurate and efficient numerical methods for one-dimensional fully nonlinear second order elliptic and parabolic partial differential equations （PDEs）. In the paper we present a general framework for constructing high order interior penalty discontinuous Galerkin （IP-DG） methods for approximating viscosity solutions of these fully nonlinear PDEs. In order to capture discontinuities of the second order derivative uxx of the solution u, three independent functions p1,p2 and p3 are introduced to represent numerical derivatives using various one-sided limits. The proposed DG frame- work, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a nonlinear problem into a mostly linear system of equations where the nonlinearity has been modified to include multiple values of the second order derivative uxz. The proposed framework extends a companion finite difference framework developed by the authors in [9] and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes. In addition to the nonstandard mixed formulation setting, another main idea is to replace the fully nonlinear differential operator by a numerical operator which is consistent with the differential operator and satisfies certain monotonicity （called g-monotonicity） properties. To ensure such a g-monotonicity, the crux of the construction is to introduce the numerical moment, which plays a critical role in the proposed DG frame- work. The g-monotonicity gives the DG methods the ability to select the mathematically ＂correct＂ solution （i.e., the viscosity solution） among all possible solutions. Moreover, the g-monotonicity allows for the possible development of more efficient nonlinear solvers as the special nonlinearity of the algebraic systems can be explored to decouple the equations. This paper also presents and analyzes numerical results for several numerical test problems which are used to guage the accuracy and efficiency of the proposed DG methods.     

OPTIMAL POINT-WISE ERROR ESTIMATE OF A COMPACT FINITE DIFFERENCE SCHEME FOR THE COUPLED NONLINEAR SCHRODINGER EQUATIONS

In this paper, we analyze a compact finite difference scheme for computing a coupled nonlinear SchrSdinger equation. The proposed scheme not only conserves the totM mass and energy in the discrete level but also is decoupled and linearized in practical computa- tion. Due to the difficulty caused by compact difference on the nonlinear term, it is very hard to obtain the optimal error estimate without any restriction on the grid ratio. In order to overcome the difficulty, we transform the compact difference scheme into a special and equivalent vector form, then use the energy method and some important lemmas to obtain the optimal convergent rate, without any restriction on the grid ratio, at the order of O（h4 ＋r2） in the discrete L∞ -norm with time step - and mesh size h. Finally, numerical results are reported to test our theoretical results of the proposed scheme.     

A new fractional finite volume method for solving the fractional diffusion equation

The inherent heterogeneities of many geophysical systems often gives rise to fast and slow pathways to water and chemical movement. One approach to model solute transport through such media is by fractional diffusion equations with a space–time dependent variable coefficient. In this paper, a two-sided space fractional diffusion model with a space–time dependent variable coefficient and a nonlinear source term subject to zero Dirichlet boundary conditions is considered.Some finite volume methods to solve a fractional differential equation with a constant dispersion coefficient have been proposed. The spatial discretisation employs fractionally-shifted Grünwald formulas to discretise the Riemann–Liouville fractional derivatives at control volume faces in terms of function values at the nodes. However, these finite volume methods have not been extended to two-dimensional and three-dimensional problems in a natural manner. In this paper, a new weighted fractional finite volume method with a nonlocal operator (using nodal basis functions) for solving this two-sided space fractional diffusion equation is proposed. Some numerical results for the Crank–Nicholson fractional finite volume method are given to show the stability, consistency and convergence of our computational approach. This novel simulation technique provides excellent tools for practical problems even when a complex transition zone is involved. This technique can be extend to two-dimensional and three-dimensional problems with complex regions.     

Global stability and stabilization of more general stochastic nonlinear systems

This paper considers the global stability and stabilization of more general stochastic nonlinear systems. Due to the absence of the conventional assumptions (e.g., Lipschitz condition), the stochastic nonlinear systems under investigation may have more than one weak solution. However, the most associated results are only applicable to the stochastic systems having a unique strong solution, and therefore, it is meaningful to refine and extend the relevant concepts and methods to the more general case. In this paper, the concepts of stochastic stability in the more general sense are first introduced to cover the stochastic nonlinear systems having more than one weak solution. Then, the generalized stochastic Barbashin–Krasovskii theorem and LaSalle theorem are established, which present the criterions of stochastic stability for more general stochastic nonlinear systems. As one of the main contributions in this paper, we rigorously prove the generalized stochastic Barbashin–Krasovskii theorem. Moreover, based on the generalized theorems, the output-feedback and state-feedback stabilization are accomplished respectively for two classes of high-order stochastic nonlinear systems under rather weaker assumptions comparing to the existing literature.     

On super-linear Emden–Fowler type differential equations

We study the second order Emden–Fowler type differential equation

(a(t)|x^′|^αsgnx^′)^′+b(t)|x|^βsgnx=0${(a(t){\left|x^{\prime }\right|}^{\alpha }\mathrm{sgn}{x}^{\prime })}^{\prime }+b(t){|x|}^{\beta }\mathrm{sgn}x=0$

in the super-linear case α<β$\alpha <\beta$. Using a Hölder-type inequality, we resolve the open problem on the possible coexistence on three possible types of nononscillatory solutions (subdominant, intermediate, and dominant solutions). Jointly with this, sufficient conditions for the existence of globally positive intermediate solutions are established. Some of our results are new also for the Emden–Fowler equation.     

EXISTENCE OF POSITIVE MULTIPLE SOLUTIONS TO FRACTIONAL DIFFERENTIAL EQUATIONS WITH p-LAPLACIAN OPERATOR

In this paper, we consider a two-point fractional boundary value problem. We provide sufficient conditions for the existence of multiple positive solutions to the boundary value problem by Krasnosel＇skii fixed point theorem on the cone.     

THE FRACTIONAL QUADRATIC-FORM IDENTITY AND BI-HAMILTONIAN STRUCTURES OF AN INTEGRABLE COUPLING OF THE FRACTIONAL COUPLED BURGERS HIERARCHY

Based on a general isospectral problem of fractional order and the fractional quadratic-form identity by Yue and Xia, the new integrable coupling of fractional coupled Burgers hierarchy and its fractional bi-Hamiltonian structures are obtained.     

带参数共轭梯度法簇的全局收敛性

共轭梯度法是最优化中最常用的方法之一,广泛地应用于求解大规模优化问题,其中参数β_k的不同选取可以构成不同的共轭梯度法.给出了一类含有三个参数的共轭梯度算法,这种算法能够在给定的条件下证明选定的β_k在每一步都能产生一个下降方向,同时在强Wolfe线搜索下,这种算法具有全局收敛性.