<i>L</i>-Stable Block Hybrid Second Derivative Algorithm for Parabolic Partial Differential Equations

An L-stable block method based on hybrid second derivative algorithm (BHSDA) is provided by a continuous second derivative method that is defined for all values of the independent variable and applied to parabolic partial differential equations (PDEs). The use of the BHSDA to solve PDEs is facilitated by the method of lines which involves making an approximation to the space derivatives, and hence reducing the problem to that of solving a time-dependent system of first order initial value ordinary differential equations. The stability properties of the method is examined and some numerical results presented.

Data-driven and physical-based identification of partial differential equations for multivariable system

Data-driven partial differential equation identification is a potential breakthrough to solve the lack of physical equations in complex dynamic systems.However,existing equation identification methods still cannot effectively identify equations from multivariable complex systems.In this work,we combine physical constraints such as dimension and direction of equation with data-driven method,and successfully identify the Navier-Stocks equations from the flow field data of Karman vortex street.This method provides an effective approach to identify partial differential equations of multivariable complex systems.

THE HIGH ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING PARABOLIC EQUATIONS 3-DIMENSION

In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.

A SINGULAR PERTURBATION PROBLEM FOR PERIODIC BOUNDARY PARTIAL DIFFERENTIAL EQUATION

In this paper, we consider a singular perturbation elliptic-parabolic partial differential equation for periodic boundary value problem, and construct a difference scheme. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, we prove that the scheme constructed by this paper converges uniformly to the solution of its original problem with O(r+h2).

Optimality Conditions and Algorithms for Direct Optimizing the Partial Differential Equations

New form of necessary conditions for optimality (NCO) is considered. They can be useful for design the direct infinite- dimensional optimization algorithms for systems described by partial differential equations (PDE). Appropriate algo-rithms for unconstrained minimizing a functional are considered and tested. To construct the algorithms, new form of NCO is used. Such approach demonstrates fast uniform convergence at optimal solution in infinite-dimensional space.

Classification of All Single Traveling Wave Solutions to (3 + 1)-Dimensional Breaking Soliton Equation

In order to get the exact traveling wave solutions to nonlinear partial differential equation, the complete discrimination system for polynomial and direct integral method are applied to the considered equation. All single traveling wave solutions to the equation can be obtained. As an example, we give the solutions to (3 + 1)-dimensional breaking soliton equation.

EFFICIENT NUMERICAL ALGORITHMS FOR THREE-DIMENSIONAL FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS

This paper detailedly discusses the locally one-dimensional numerical methods for ef- ficiently solving the three-dimensional fractional partial differential equations, including fractional advection diffusion equation and Riesz fractional diffusion equation. The second order finite difference scheme is used to discretize the space fractional derivative and the Crank-Nicolson procedure to the time derivative. We theoretically prove and numerically verify that the presented numerical methods are unconditionally stable and second order convergent in both space and time directions. In particular, for the Riesz fractional dif- fusion equation, the idea of reducing the splitting error is used to further improve the algorithm, and the unconditional stability and convergency are also strictly proved and numerically verified for the improved scheme.

A POSTERIORI ERROR ESTIMATE FOR BOUNDARY CONTROL PROBLEMS GOVERNED BY THE PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

In this paper, we discuss the a posteriori error estimate of the finite element approximation for the boundary control problems governed by the parabolic partial differential equations. Three different a posteriori error estimators are provided for the parabolic boundary control problems with the observations of the distributed state, the boundary state and the final state. It is proven that these estimators are reliable bounds of the finite element approximation errors, which can be used as the indicators of the mesh refinement in adaptive finite element methods.

OSCILLATION THEOREM TO SYSTEMS OF IMPULSIVE NEUTRAL DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

In this paper, we study the oscillation of solutions to the systems of impulsive neutral delay parabolic partial differential equations. Under two different boundary conditions, we obtain some sufficient conditions for oscillation of all solutions to the systems.

A Fitted Numerov Method for Singularly Perturbed Parabolic Partial Differential Equation with a Small Negative Shift Arising in Control Theory

In this paper,a fitted Numerov method is constructed for a class of singularly perturbed one-dimensional parabolic partial differential equations with a small negative shift in the temporal variable.Similar boundary value problems are associated with a furnace used to process a metal sheet in control theory.Here,the study focuses on the effect of shift on the boundary layer behavior of the solution via finite difference approach.When the shift parameter is smaller than the perturbation parameter,the shifted term is expanded in Taylor series and an exponentially fitted tridiagonal finite difference scheme is developed.The proposed finite difference scheme is unconditionally stable.When the shift parameter is larger than the perturbation parameter,a special type of mesh is used for the temporal variable so that the shift lies on the nodal points and an exponentially fitted scheme is developed.This scheme is also unconditionally stable.The applicability of the proposed methods is demonstrated by means of two examples.

ODE-Based Multistep Schemes for Backward Stochastic Differential Equations

In this paper,we explore a new approach to design and analyze numerical schemes for backward stochastic differential equations(BSDEs).By the nonlinear Feynman-Kac formula,we reformulate the BSDE into a pair of reference ordinary differential equations(ODEs),which can be directly discretized by many standard ODE solvers,yielding the corresponding numerical schemes for BSDEs.In particular,by applying strong stability preserving(SSP)time discretizations to the reference ODEs,we can propose new SSP multistep schemes for BSDEs.Theoretical analyses are rigorously performed to prove the consistency,stability and convergency of the proposed SSP multistep schemes.Numerical experiments are further carried out to verify our theoretical results and the capacity of the proposed SSP multistep schemes for solving complex associated problems.

Zakharov-Rubenchik方程组的格子Boltzmann方法

Zakharov-Rubenchik方程组常用于描述非线性介质中高、低频波间相互作用的波耦合现象。本文针对该方程组的数值求解问题,构建了一种格子Boltzmann方法的D1Q3演化模型,并利用Chapman-Enskog展开和多尺度分析技术,推导出了各个方向的平衡态分布函数和修正函数的具体表达式,从而将所建的演化模型准确恢复到宏观方程组。最后,通过数值算例证明了该方法的有效性。

Adaptive Sparse Grid Discontinuous Galerkin Method:Review and Software Implementation

This paper reviews the adaptive sparse grid discontinuous Galerkin(aSG-DG)method for computing high dimensional partial differential equations(PDEs)and its software implementation.The C++software package called AdaM-DG,implementing the aSG-DG method,is available on GitHub at https://github.com/JuntaoHuang/adaptive-multiresolution-DG.The package is capable of treating a large class of high dimensional linear and nonlinear PDEs.We review the essential components of the algorithm and the functionality of the software,including the multiwavelets used,assembling of bilinear operators,fast matrix-vector product for data with hierarchical structures.We further demonstrate the performance of the package by reporting the numerical error and the CPU cost for several benchmark tests,including linear transport equations,wave equations,and Hamilton-Jacobi(HJ)equations.

OSCILLATION OF SYSTEMS OF IMPULSIVE DELAY PARABOLIC EQUATIONS ABOUT BOUNDARY VALUE PROBLEMS

In this paper, oscillation of solutions to a class of impulsive delay parabolic partial differential equations system with higher order Laplace operator is studied. Under two different boundary value conditions, we establish some sufficient criteria with respect to the oscillations of such systems, employing first-order impulsive delay differential inequalities. The results fully reflect the influence action of impulsive and delay in oscillation.

Difference Approximation of Stochastic Elastic Equation Driven by Infinite Dimensional Noise

An explicit differencescheme is described,analyzed and tested for numer-ically approximating stochastic elastic equation driven by infinite dimensional noise.The noise processes are approximated by piecewise constant random processes and the integral formula of the stochastic elastic equation is approximated by a truncated series.Error analysis of the numerical method yields estimate of convergence rate.The rate of convergence is demonstrated with numerical experiments.

EXISTENCE AND UNIQUENESS OF WEAK SOLUTIONS FOR A NONLINEAR PARABOLIC EQUATION RELATED TO IMAGE ANALYSIS

In this paper we establish the existence and uniqueness of weak solutions for the initial-boundary value problem of a nonlinear parabolic partial differential equation, which is related to the Malik-Perona model in image analysis.

A QUASI-NEWTON METHOD IN INFINITE-DIMENSIONAL SPACES AND ITS APPLICATION FOR SOLVING A PARABOLIC INVERSE PROBLEM

A Quasi-Newton method in Infinite-dimensional Spaces (QNIS) for solving operator equations is presellted and the convergence of a sequence generated by QNIS is also proved in the paper. Next, we suggest a finite-dimensional implementation of QNIS and prove that the sequence defined by the finite-dimensional algorithm converges to the root of the original operator equation providing that the later exists and that the Frechet derivative of the governing operator is invertible. Finally, we apply QNIS to an inverse problem for a parabolic differential equation to illustrate the efficiency of the finite-dimensional algorithm.

On the speed of convergence of Picard iterations of backward stochastic differential equations

It is a well-established fact in the scientific literature that Picard iterations of backward stochastic differential equations with globally Lipschitz continuous nonlinearities converge at least exponentially fast to the solution.In this paper we prove that this convergence is in fact at least square-root factorially fast.We show for one example that no higher convergence speed is possible in general.Moreover,if the nonlinearity is zindependent,then the convergence is even factorially fast.Thus we reveal a phase transition in the speed of convergence of Picard iterations of backward stochastic differential equations.

An efficient technique for solving fractional-order diffusion equations arising in oil pollution

In this article,non-linear time-fractional diffusion equations are considered to describe oil pollution in the water.The latest technique,fractional reduced differential transform method(FRDTM),is used to ac-quire approximate solutions of the time fractional-order diffusion equation and two cases of Allen-Cahn equations.The acquired results are collated with the exact solutions and other results from literature for integer-orderα,which reveal that the proposed method is effective.Hence,FRDTM can be employed to obtain solutions for different types of nonlinear fractional-order IVPs arising in engineering and science.

具连续分布滞量的非线性中立型抛物偏泛函微分方程的振动性

考虑一类具连续分布滞量的非线性中立型抛物偏泛函微分方程解的振动性,借助Green定理和时滞微分不等式获得了这类方程在Robin,Dirichlet边值条件下所有解振动的若干充分条件.