半线性抛物型方程奇异摄动问题的数值解

本文讨论具有抛物边界层的半线性抛物型方程奇异摄动问题的数值解法,在非均匀网格上构造了两层非线性差分格式,证明了差分格式是一致收敛的,给出了一些数值例子.     

Stability of Rothe difference scheme for the reverse parabolic problem with integral boundary condition

In this paper, we study the approximation of reverse parabolic problem with integral boundary condition. The Rothe difference scheme for an approximate solution of reverse problem is discussed. We establish stability and coercive stability estimates for the solution of the Rothe difference scheme. In sequel, we investigate the first order of accuracy difference scheme for approximation of boundary value problem for multidimensional reverse parabolic equation and obtain stability estimates for its solution. Finally, we give numerical results together with an explanation on the realization in one- and two-dimensional test examples.     

CONVERGENCE OF THE CRANK-NICOLSON/NEWTON SCHEME FOR NONLINEAR PARABOLIC PROBLEM

In this paper, the Crank-Nicolson/Newton scheme for solving numerically secondorder nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the spatial discretization of the problem and the Crank-Nicolson/Newton scheme is applied to the time discretization of the resulted finite element equations. Moreover, assuming the appropriate regularity of the exact solution and the finite element solution, we obtain optimal error estimates of the fully discrete CrankNicolson/Newton scheme of nonlinear parabolic problem. Finally, numerical experiments are presented to show the efficient performance of the proposed scheme.     

Convergence of Iterative Difference Schemes for Two and Three Dimensional Nonlinear Parabolic Systems

In this paper, we study the general difference schemes of the boundary value problem for the nonlinear parabolic systems with two and three space dimensions. To solve the nonlinear difference schemes, we construct an iterative sequence from the solutions or the linearized difference schemes. We shall prove the convergence of the difference solutions for the iterative difference schemes to the solution of the original boundary value problem or the nonlinear parabolic systems.     

Numerical solution of a quasilinear parabolic problem

A combined approach of linearisation techniques and finite difference method is presented for obtaining the numerical solution of a quasilinear parabolic problem. The given problem is reduced to a sequence of linear problems by using the Picard or Newton methods. Each problem of this sequence is approximated by Crank-Nicolson difference scheme. The solutions of the resulting system of algebraic equations are obtained by using Block-Gaussian elimination method. Two numerical examples are solved by using both linearisation procedures to illustrate the method. For these examples, the Newton method is found to be more effective, especially when the given nonlinear problem has steep gradients.     

Stability with respect to the input data and monotonicity of an implicit finite-difference scheme for a quasilinear parabolic equation

We analyze the stability and monotonicity of a conservative difference scheme approximating an initial-boundary value problem for a quasilinear parabolic equation under specific conditions imposed solely on the problem input data. We prove some kinds of the maximum principle for the nonlinear equations that are used in the derivation of a priori estimates for the solution; we also prove estimates for some kinds of recursion inequalities that are used in the derivation of a priori estimates for higher-order derivatives, these estimates being necessary for proving the continuous dependence of the solution on small perturbations of the input data and for analyzing monotonicity in the nonlinear case. We show that, depending on the properties of the input data, higher derivatives can become infinite in finite critical time. We obtain conditions on the input data guaranteeing the stability of the difference scheme on the entire time interval.     

Compact difference methods applied to initial-boundary value problems for mixed systems

Summary. An initial--boundary value problem to a system of nonlinear partial differential equations, which consists of a hyperbolic and a parabolic part, is taken into consideration. The problem is discretised by a compact finite difference method. An approximation of the numerical solution is constructed, at which the difference scheme is linearised. Nonlinear convergence is proved using the stability of the linearised scheme. Finally, a computational experiment for a noncompact scheme is presented. Received May 20, 1995     

Uniform quadratic convergence of monotone iterates for nonlinear singularly perturbed parabolic problems

This paper deals with a monotone iterative method for solving nonlinear singularly perturbed parabolic problems. Monotone sequences, based on the method of upper and lower solutions, are constructed for a nonlinear difference scheme which approximates the nonlinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. The monotone sequences possess quadratic convergence rate. An analysis of uniform convergence of the monotone iterative method to the solutions of the nonlinear difference scheme and to the continuous problem is given. Numerical experiments are presented.     

Uniform convergent monotone iterates for semilinear singularly perturbed parabolic problems

This paper deals with discrete monotone iterative methods for solving semilinear singularly perturbed parabolic problems. Monotone sequences, based on the accelerated monotone iterative method, are constructed for a nonlinear difference scheme which approximates the semilinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. An analysis of uniform convergence of the monotone iterative method to the solutions of the nonlinear difference scheme and continuous problem is given. Numerical experiments are presented.     

Exact three-point difference scheme for a nonlinear boundary-value problem on the semiaxis

For the numerical solution of boundary-value problems on the semiaxis for second-order nonlinear ordinary differential equations, an exact three-point difference scheme is constructed and substantiated. Under the conditions of existence and uniqueness of solution of a boundary-value problem, we prove the existence and uniqueness of solution of the exact three-point difference scheme and convergence of the method of successive approximations for its solution.     

An Improved Two-Grid Technique for the Nonlinear Time-Fractional Parabolic Equation Based on the Block-Centered Finite Difference Method

A combined scheme of the improved two-grid technique with the block-centered finite difference method is constructed and analyzed to solve the nonlinear time-fractional parabolic equation. This method is considered where the nonlinear problem is solved only on a coarse grid of size $H$ and two linear problems based on the coarse-grid solutions and one Newton iteration is considered on a fine grid of size $h$. We provide the rigorous error estimate, which demonstrates that our scheme converges with order $\mathcal{O}(\Delta t^{2-\alpha}+h^2+H^4)$ on non-uniform rectangular grid. This result indicates that the improved two-grid method can obtain asymptotically optimal approximation as long as the mesh sizes satisfy $h=\mathcal{O}(H^2).$ Finally, numerical tests confirm the theoretical results of the presented method.     

The Nonlinear Initial-Boundary Value Problem and the Existence of Multi-Dimension al Shook Wave for Quasilinear Hyperbolic-Parabolic Coupled Systems

For the quasilinear hyperbolie-parabolio coupled system, the nonlinear initial- boundary value problem and the shook wave free boundary problem are considered. By linear iteration, the existence and uniqueness of the local H^m (m\geq [N+1/2]+4) solution are obtained under the assumption that for the fixed boundary problem, the boundary conditions are uniformly Lopatinski well-posed with respect to the hyperbolic and parabolic part, and for the free boundary problem, there exists a linear stable shock front structure. In particular, the local existence of the isothermal shock wave solution for radiative hydrodynamic eqations is proved.     

A FINITE DIFFERENCE METHOD FOR THE HEAT EQUATION WITH A NONLINEAR BOUNDARY CONDITION

This article is concerned with the numerical solution to a parabolic equation with a kind of nonlinear boundary conditions. A difference scheme is constructed by the method of reduction of order on uniform mesh to solve the problem. It is proved that the difference scheme is uniquely solvable and uncon-ditionaUy convergent with the convergence order 2 in both space and time in an energy norm. An effective iterative algorithm is given and a numerical example is presented to demonstrate the theoretical results.     

Monotone iterates with quadratic convergence rate for solving weighted average approximations to semilinear parabolic problems

This paper deals with discrete monotone iterative methods for solving semilinear singularly perturbed parabolic problems. Monotone sequences, based on the accelerated monotone iterative method, are constructed for a nonlinear difference scheme which approximates the semilinear parabolic problem. This monotone convergence leads to the existence-uniqueness theorem. An analysis of convergence of the monotone iterative method to the solutions of the nonlinear difference scheme is given. Numerical experiments are presented.     

热方程非线性边界条件问题的有限差分方法

This article is concerned with the numerical solution to a parabolic equation with a kind of nonlinear boundary conditions. A difference scheme is constructed by the method of reduction of order on uniform mesh to solve the problem. It is proved that the difference scheme is uniquely solvable and uncon-ditionaUy convergent with the convergence order 2 in both space and time in an energy norm. An effective iterative algorithm is given and a numerical example is presented to demonstrate the theoretical results.     

A linearized difference scheme for semilinear parabolic equations with nonlinear absorbing boundary conditions

A novel three level linearized difference scheme is proposed for the semilinear parabolic equation with nonlinear absorbing boundary conditions. The solution of this problem will blow up in finite time. Hence this difference scheme is coupled with an adaptive time step size, i.e., when the solution tends to infinity, the time step size will be smaller and smaller. Furthermore, the solvability, stability and convergence of the difference scheme are proved by the energy method. Numerical experiments are also given to demonstrate the theoretical second order convergence both in time and in space in L_∞-norm.     

On the role of conservation laws in the problem on the occurrence of unstable solutions for quasilinear parabolic equations and their approximations

For a broad class of initial-boundary value problems for quasilinear parabolic equations with nonlinear source and their approximations, we show that if the initial energy is negative, then the solution always blows up in finite time. This is especially important for finding sufficiently simple and easy-to-verify conditions guaranteeing the presence of physical effects such as heat localization in peaking modes or thermal explosions and for deriving two-sided estimates of the solution lifespan. We construct the corresponding new classes of difference schemes for which grid analogs of integral conservation laws hold. We show that, to obtain efficient two-sided estimates for the blow-up time of the solution of the differential problem, in practice, one should use difference schemes with explicit as well as implicit approximation to the source.     

A Trilayer Difference Scheme for One-Dimensional Parabolic Systems

In order to obtain the numerical solution for a one-dimensional parabolic system, an unconditionally stable difference method is investigated in [1]. If the number of unknown functions is M, for each time step only M times of calculation are needed. The rate of convergence is $O(\tau+h^2)$. On the basis of [1], an alternating calculation difference scheme is presented in [2]; the rate of the convergence is $O(\tau^2+h^2)$. The difference schemes in [1] and [2] are economic ones. For the $\alpha$-$th$ equation, only $U_{\alpha}$ is an unknown function; the others $U_{\beta}$ are given evaluated either in the last step or in the present step. So the practical calculation is quite convenient. The purpose of this paper is to derive a trilayer difference scheme for one-dimensional parabolic systems. It is known that the scheme is also unconditionally stable and the rate of convergence is $O(\tau^2+h^2)$.     

Numerical method of lines for parabolic functional differential equations

Initial boundary value problems for nonlinear parabolic functional differential equations are transformed by discretization in space variables into systems of ordinary functional differential equations. A comparison theorem for differential difference inequalities is proved. Sufficient conditions for the convergence of the numerical method of lines are given. An explicit Euler method is proposed for the numerical solution of systems thus obtained. This leads to difference scheme for the original problem. A complete convergence analysis for the method is given.     

Discrete Approximations for Singularly Perturbed Boundary Value Problems with Parabolic Layers,I

In this series of three papers we study singularly perturbed (SP) boundary value problems for equations of elliptic and parabolic type. For small values of the perturbation parameter parabolic boundary and interior layers appear in these problems. If classical discretisation methods are used, the solution of the finite difference scheme and the approximation of the diffusive flux do not converge uniformly with respect to this parameter. Using the method of special, adapted grids, we can construct difference schemes that allow approximation of the solution and the normalised diffusive flux uniformly with respect to the small parameter. We also consider singularly perturbed boundary value problems for convection-diffusion equations. Also for these problems we construct special finite difference schemes, the solution of which converges $ε$-uniformly. We study what problems appear, when classical schemes are used for the approximation of the spatial derivatives. We compare the results with those obtained by the adapted approach. Results of numerical experiments are discussed. In the three papers we first give an introduction on the general problem, and then we consider respectively (i) Problems for SP parabolic equations, for which the solution and the normalised diffusive fluxes are required; (ii) Problems for SP elliptic equations with boundary conditions of Dirichlet, Neumann and Robin type; (iii) Problems for SP parabolic equation with discontinuous boundary conditions.