Numerical methods for fourth‐order nonlinear elliptic boundary value problems

The aim of this article is to present several computational algorithms for numerical solutions of a nonlinear finite difference system that represents a finite difference approximation of a class of fourth‐order elliptic boundary value problems. The numerical algorithms are based on the method of upper and lower solutions and its associated monotone iterations. Three linear monotone iterative schemes are given, and each iterative scheme yields two sequences, which converge monotonically from above and below, respectively, to a maximal solution and a minimal solution of the finite difference system. This monotone convergence property leads to upper and lower bounds of the solution in each iteration as well as an existence‐comparison theorem for the finite difference system. Sufficient conditions for the uniqueness of the solution and some techniques for the construction of upper and lower solutions are obtained, and numerical results for a two‐point boundary‐value problem with known analytical solution are given. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:347–368, 2001     

Finite difference reaction diffusion equations with nonlinear boundary conditions

This article is concerned with numerical solutions of finite difference systems of reaction diffusion equations with nonlinear internal and boundary reaction functions. The nonlinear reaction functions are of general form and the finite difference systems are for both time-dependent and steady-state problems. For each problem a unified system of nonlinear equations is treated by the method of upper and lower solutions and its associated monotone iterations. This method leads to a monotone iterative scheme for the computation of numerical solutions as well as an existence-comparison theorem for the corresponding finite difference system. Special attention is given to the dynamical property of the time-dependent solution in relation to the steady-state solutions. Application is given to a heat-conduction problem where a nonlinear radiation boundary condition obeying the Boltzmann law of cooling is considered. This application demonstrates a bifurcation property of two steady-state solutions, and determines the dynamic behavior of the time-dependent solution. Numerical results for the heat-conduction problem, including a test problem with known analytical solution, are presented to illustrate the various theoretical conclusions. © 1995 John Wiley & Sons, Inc.     

Numerical methods for coupled systems of quasilinear elliptic equations with nonlinear boundary conditions

This paper is concerned with numerical solutions of a coupled system of arbitrary number of quasilinear elliptic equations under combined Dirichlet and nonlinear boundary conditions. A finite difference system for a transformed system of the quasilinear equations is formulated, and three monotone iterative schemes for the computation of numerical solutions are given using the method of upper and lower solutions. It is shown that each of the three monotone iterations converges to a minimal solution or a maximal solution depending on whether the initial iteration is a lower solution or an upper solution. A comparison result among the three iterative schemes is given. Also shown is the convergence of the minimal and maximal discrete solutions to the corresponding minimal and maximal solutions of the continuous system as the mesh size tends to zero. These results are applied to a heat transfer problem with temperature dependent thermal conductivity and a Lotka-Volterra cooperation system with degenerate diffusion. This degenerate property leads to some interesting distinct property of the system when compared with the non-degenerate semilinear systems. Numerical results are given to the above problems, and in each problem an explicit continuous solution is constructed and is used to compare with the computed solution     

Blow‐up behavior of the Solution for the problem in a subdiffusive mediums

The diffusion problem in a subdiffusive medium is formulated by using the fractional differential operator. In this paper, we consider a fractional differential equation with concentrated source. The existence of the solution in a finite time is given. The finite time blow‐up criteria for the solution of the problem is established, and the location of the blow‐up point is investigated.     

Accelerated monotone iterative methods for finite difference equations of reaction-diffusion

This paper is concerned with numerical methods for a finite difference system of reaction-diffusion-convection equation under nonlinear boundary condition. Various monotone iterative methods are presented, and each of these methods leads to an existence-comparison theorem as well as a computational algorithm for numerical solutions. The monotone property of the iterations gives improved upper and lower bounds of the solution in each iteration, and the rate of convergence of the iterations is either quadratic or nearly quadratic depending on the property of the nonlinear function. Application is given to a model problem from chemical engineering, and some numerical results, including a test problem with known analytical solution, are presented to illustrate the various rates of convergence of the iterations. Received November 2, 1995 / Revised version received February 10, 1997     

Nonstationary monotone iterative methods for nonlinear partial differential equations

A simple technique is given in this paper for the construction and analysis of monotone iterative methods for a class of nonlinear partial differential equations. With the help of the special nonlinear property we can construct nonstationary parameters which can speed up the iterative process in solving the nonlinear system. Picard, Gauss–Seidel, and Jacobi monotone iterative methods are presented and analyzed for the adaptive solutions. The adaptive meshes are generated by the 1-irregular mesh refinement scheme which together with the M-matrix of the finite element stiffness matrix lead to existence–uniqueness–comparison theorems with simple upper and lower solutions as initial iterates. Some numerical examples, including a test problem with known analytical solution, are presented to demonstrate the accuracy and efficiency of the adaptive and monotone properties. Numerical results of simulations on a MOSFET with the gate length down to 34 nm are also given.     

Block monotone iterative methods for numerical solutions of nonlinear elliptic equations

Summary. Two block monotone iterative schemes for a nonlinear algebraic system, which is a finite difference approximation of a nonlinear elliptic boundary-value problem, are presented and are shown to converge monotonically either from above or from below to a solution of the system. This monotone convergence result yields a computational algorithm for numerical solutions as well as an existence-comparison theorem of the system, including a sufficient condition for the uniqueness of the solution. An advantage of the block iterative schemes is that the Thomas algorithm can be used to compute numerical solutions of the sequence of iterations in the same fashion as for one-dimensional problems. The block iterative schemes are compared with the point monotone iterative schemes of Picard, Jacobi and Gauss-Seidel, and various theoretical comparison results among these monotone iterative schemes are given. These comparison results demonstrate that the sequence of iterations from the block iterative schemes converges faster than the corresponding sequence given by the point iterative schemes. Application of the iterative schemes is given to a logistic model problem in ecology and numerical ressults for a test problem with known analytical solution are given. Received August 1, 1993 / Revised version received November 7, 1994     

分数阶微分方程的一些新结果

第一部分,介绍分数阶导数的定义和著名的Mittag—Leffler函数的性质．第二部分,利用单调迭代方法给出了具有2序列Riemann—Liouville分数阶导数微分方程初值问题解的存在性和唯一性．第三部分,利用上下解方法和Schauder不动点定理给出了具有2序列Riemann—Liouville分数阶导数微分方程周期边值问题解的存在性．第四部分,利用Leray—Schauder不动点定理和Banach压缩映像原理建立了具有n序列Riemann—Liouville分数阶导数微分方程初值问题解的存在性、唯一性和解对初值的连续依赖性．第五部分,利用锥上的不动点定理给出了具有Caputo分数阶导数微分方程边值问题,在超线性（次线性）条件下C310,11正解存在的充分必要条件．最后一部分,通过建立比较定理和利用单调迭代方法给出了具有Caputo分数阶导数脉冲微分方程周期边值问题最大解和最小解的存在性．     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

A modified accelerated monotone iterative method for finite difference reaction-diffusion-convection equations

This paper is concerned with monotone algorithms for the finite difference solutions of a class of nonlinear reaction-diffusion-convection equations with nonlinear boundary conditions. A modified accelerated monotone iterative method is presented to solve the finite difference systems for both the time-dependent problem and its corresponding steady-state problem. This method leads to a simple and yet efficient linear iterative algorithm. It yields two sequences of iterations that converge monotonically from above and below, respectively, to a unique solution of the system. The monotone property of the iterations gives concurrently improving upper and lower bounds for the solution. It is shown that the rate of convergence for the sum of the two sequences is quadratic. Under an additional requirement, quadratic convergence is attained for one of these two sequences. In contrast with the existing accelerated monotone iterative methods, our new method avoids computing local maxima in the construction of these sequences. An application using a model problem gives numerical results that illustrate the effectiveness of the proposed method.     

On the blow‐up phenomena for a 1‐dimensional equation of ion sound waves in a plasma: Analytical and numerical investigation

The initial‐boundary value problem for an equation of ion sound waves in plasma is considered. A theorem on nonextendable solution is proved. The blow‐up phenomena are studied. The sufficient blow‐up conditions and the blow‐up time are analysed by the method of the test functions. This analytical a priori information is used in the numerical experiments, which are able to determine the process of the solution's blow‐up more accurately.     

Nonlinear fourth-order elliptic equations with nonlocal boundary conditions

This paper is concerned with a class of fourth-order nonlinear elliptic equations with nonlocal boundary conditions, including a multi-point boundary condition in a bounded domain of Rn. Also considered is a second-order elliptic equation with nonlocal boundary condition, and the usual multi-point boundary problem in ordinary differential equations. The aim of the paper is to show the existence of maximal and minimal solutions, the uniqueness of a positive solution, and the method of construction for these solutions. Our approach to the above problems is by the method of upper and lower solutions and its associated monotone iterations. The monotone iterative schemes can be developed into computational algorithms for numerical solutions of the problem by either the finite difference method or the finite element method.     

Finite difference reaction–diffusion equations with nonlinear diffusion coefficients

Summary. A monotone iterative method for numerical solutions of a class of finite difference reaction-diffusion equations with nonlinear diffusion coefficient is presented. It is shown that by using an upper solution or a lower solution as the initial iteration the corresponding sequence converges monotonically to a unique solution of the finite difference system. It is also shown that the solution of the finite difference system converges to the solution of the continuous equation as the mesh size decreases to zero. Received February 18, 1998 / Revised version received April 21, 1999 / Published online February 17, 2000     

Series solutions of a semilinear wave equation

We are concerned with the reconstruction of series solutions of a semilinear wave equation with a quadratic nonlinearity. The solution which may blow up in finite time is sought as a sum of exponential functions and is shown to be a classical one. The constructed solutions can be used to benchmark numerical methods used to approximate solutions of nonlinear equations.     

求解带非线性边界条件的反应扩散方程数值解的组合迭代方法

1　引　　言我们首先考虑如下抛物型方程ｕｔ-ＤΔｕ =ｆ(ｘ ,ｔ ,ｕ) (ｔ∈ ( 0 ,Ｔ],ｘ∈Ω ) ｕ/ ν+ βｕ =ｇ(ｘ ,ｔ ,ｕ) (ｔ∈ ( 0 ,Ｔ],ｘ∈ Ω )ｕ(ｘ ,0 ) =ψ(ｘ) (ｘ∈Ω )( 1 .1 )其中Ｔ为正常数 ,Ω 是ＲＰ 空间的有界区域 记ＱＴ=Ω × ( 0 ,Ｔ],ＳＴ= Ω × ( 0 ,Ｔ],假设在ＱＴ上Ｄ≡ｄ(ｘ ,ｔ) >0 ,在ＳＴ 上β≡β(ｘ ,ｔ)≥ 0 又设 ｆ(ｘ ,ｔ,ｕ) ,ｇ(ｘ ,ｔ,ｕ)为关于ｕ的非线性函数 ,且对ｘ ,ｔ各参数满足Ｈ¨ｏｌｄｅｒ连续条件 将 ( 1 .1 )离散化之后我们得到相应的有限差分系统 ,当 ｇ(ｘ ,ｔ,ｕ)为ｕ的线性…     

一类非线性抛物型方程组解的整体存在及爆破

研究了一类带有非线性边界条件的非线性抛物型方程组解的整体存在及解在有限时刻爆破问题.通过构造方程组的上、下解.得到了解整体存在及解在有限时刻爆破的充分条件.对指数型反应项和边界流采用了常微分方程方法构造其上下解,而其它例如第一特征值等方法运用于该方程就比较困难.     

Higher-order compact finite difference method for systems of reaction–diffusion equations

This paper is concerned with a compact finite difference method for solving systems of two-dimensional reaction–diffusion equations. This method has the accuracy of fourth-order in both space and time. The existence and uniqueness of the finite difference solution are investigated by the method of upper and lower solutions, without any monotone requirement on the nonlinear term. Three monotone iterative algorithms are provided for solving the resulting discrete system efficiently, and the sequences of iterations converge monotonically to a unique solution of the system. A theoretical comparison result for the various monotone sequences is given. The convergence of the finite difference solution to the continuous solution is proved, and Richardson extrapolation is used to achieve fourth-order accuracy in time. An application is given to an enzyme–substrate reaction–diffusion problem, and some numerical results are presented to demonstrate the high efficiency and advantages of this new approach.     

Finding a Zero of The Sum of Two Maximal Monotone Operators

In this paper, the equivalence between variational inclusions and a generalized type of Weiner–Hopf equation is established. This equivalence is then used to suggest and analyze iterative methods in order to find a zero of the sum of two maximal monotone operators. Special attention is given to the case where one of the operators is Lipschitz continuous and either is strongly monotone or satisfies the Dunn property. Moreover, when the problem has a nonempty solution set, a fixed-point procedure is proposed and its convergence is established provided that the Brézis–Crandall–Pazy condition holds true. More precisely, it is shown that this allows reaching the element of minimal norm of the solution set.     

Monotone iterative methods for finite difference system of reaction-diffusion equations

Summary This paper presents an existence-comparison theorem and an iterative method for a nonlinear finite difference system which corresponds to a class of semilinear parabolic and elliptic boundary-value problems. The basic idea of the iterative method for the computation of numerical solutions is the monotone approach which involves the notion of upper and lower solutions and the construction of monotone sequences from a suitable linear discrete system. Using upper and lower solutions as two distinct initial iterations, two monotone sequences from a suitable linear system are constructed. It is shown that these two sequences converge monotonically from above and below, respectively, to a unique solution of the nonlinear discrete equations. This formulation leads to a well-posed problem for the nonlinear discrete system. Applications are given to several models arising from physical, chemical and biological systems. Numerical results are given to some of these models including a discussion on the rate of convergence of the monotone sequences.     

Numerical blow‐up for the porous medium equation with a source

We study numerical approximations of positive solutions of the porous medium equation with a nonlinear source, where m > 1, p > 0 and L > 0 are parameters. We describe in terms of p, m, and L when solutions of a semidiscretization in space exist globally in time and when they blow up in a finite time. We also find the blow‐up rates and the blow‐up sets, proving that there is no regional blow‐up for the numerical scheme. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004