Implicit difference methods for Hamilton-Jacobi functional differential equations

Classical solutions of initial boundary value problems are approximated by solutions of associated implicit difference functional equations. A stability result is proved by using a comparison technique with nonlinear estimates of the Perron type for given functions. The Newton method is used to numerically solve nonlinear equations generated by implicit difference schemes. It is shown that there are implicit difference schemes which are convergent whereas the corresponding explicit difference methods are not. The results obtained can be applied to differential integral problems and differential equations with deviated variables.     

Numerical methods for Hamilton Jacobi functional differential equations

Initial and initial boundary value problems for first order partial functional differential equations are considered. Explicit difference schemes of the Euler type and implicit difference methods are investigated. The following theoretical aspects of the methods are presented. Sufficient conditions for the convergence of approximate solutions are given and comparisons of the methods are presented. It is proved that assumptions on the regularity of given functions are the same for both the methods. It is shown that conditions on the mesh for explicit difference schemes are more restrictive than suitable assumptions for implicit methods. There are implicit difference schemes which are convergent and corresponding explicit difference methods are not convergent. Error estimates for both the methods are construted.     

Exact finite-difference methods based on nonlinear means

The necessary and sufficient condition for an autonomous, ordinary differential equation to be exactly solved by a given Evans-Sanugi, nonlinear, one-step, finite difference method based on ‘classical means’ are derived. A necessary, but not sufficient, condition yields the most general nonlinearity for the differential equation which is independent of the step size. Examples of differential equations for which either nonlinear trapezoidal or nonlinear implicit midpoint methods based on arithmetic, harmonic, contraharmonic, quadratic, geometric, Heronian, centroidal and logarithmic means are exact, are presented. These new exact difference schemes may be useful in future developments of new Denk-Bulirsch, Le-Roux or Kojouhavov-Chen schemes for nonlinear evolution equations with or without blow-up.     

Nonlinear Implicit One-Step Schemes for Solving Initial Value Problems for Ordinary Differential Equation with Steep Gradients

A general theory for nonlinear implicit one-step schemes for solving initial value problems for ordinary differential equations is presented in this paper. The general expansion of "symmetric" implicit one-step schemes having second-order is derived and stability and convergence are studied. As examples, some geometric schemes are given. Based on previous work of the first author on a generalization of means, a fourth-order nonlinear implicit one-step scheme is presented for solving equations with steep gradients. Also, a hybrid method based on the GMS and a fourth-order linear scheme is discussed. Some numerical results are given.     

On nonlinear difference and differential equations

A generalization of the logarithmic norm to nonlinear operators, the Dahlquist constant is introduced as a useful tool for the estimation and analysis of error propagation in general nonlinear first-order ODE's. It is a counterpart to the Lipschitz constant which has similar applications to difference equations. While Lipschitz constants can also be used for ODE's, estimates based on the Dahlquist constant always give sharper results.The analogy between difference and differential equations is investigated, and some existence and uniqueness results for nonlinear (algebraic) equations are given. We finally apply the formalism to the implicit Euler method, deriving a rigorous global error bound for stiff nonlinear problems.Dedicated to my teacher and friend, Professor Germund Dahlquist, on the occasion of his 60th birthday.     

Generalized implicit Euler method for hyperbolic functional differential equations

Nonlinear hyperbolic functional differential equations with initial boundary conditions are considered. Theorems on the convergence of difference schemes and error estimates of approximate solutions are presented. The proof of the stability of the difference functional problem is based on a comparison technique. Nonlinear estimates of the Perron type with respect to the functional variable for given functions are used. Numerical examples are given (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

非线性波动方程的弱隐式与显式差分方法

广泛出现于物理、化学、机械动力学、生物、几何学等领域的非线性波动方程已经有很多的研究工作,Sine-Gordon方程和非线性受迫振动方程就是典型的例子.周毓麟教授在[1]中研究了非线性波动方程组     

Estimate of solutions for differential and difference functional equations with applications to difference methods

The theorems on the estimate of solutions for nonlinear second-order partial differential functional equations mainly of parabolic type with Dirichlet’s condition and for the suitable explicit finite difference functional schemes are proved. The proofs are based on the comparison technique. The convergent difference method given is considered without an assumption of the global generalized Perron condition on the functional variable but with local one in some sense only. It is a consequence of our estimate theorems. The functional dependence is of the Volterra type.     

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.     

On compact approximations of divergence differential equations

A method for the construction of compact difference schemes approximating divergence differential equations is proposed. The schemes have an arbitrarily prescribed order of approximation on general stencils. It is shown that the construction of such schemes for partial differential equations is based on special compact schemes approximating ordinary differential equations in several independent functions. Necessary and sufficient conditions on the coefficients of these schemes with high order of approximation are obtained. Examples of reconstruction of compact difference schemes for partial differential equations with these schemes are given. It is shown that such compact difference schemes have the same order of accuracy both for classical approximations on smooth solutions and weak approximations on discontinuous solutions.     

Numerical algorithm for solving diffusion-type equations on the basis of multigrid methods

A new numerical algorithm based on multigrid methods is proposed for solving equations of the parabolic type. Theoretical error estimates are obtained for the algorithm as applied to a two-dimensional initial-boundary value model problem for the heat equation. The good accuracy of the algorithm is demonstrated using model problems including ones with discontinuous coefficients. As applied to initial-boundary value problems for diffusion equations, the algorithm yields considerable savings in computational work compared to implicit schemes on fine grids or explicit schemes with a small time step on fine grids. A parallelization scheme is given for the algorithm.     

The method of contractor directions for partial differential equations

A further development of the concept of contractor directions is presented. An application of the method of specialized contractor directions to local existence theorems for general nonlinear equations and to solutions of some partial differential equations is discussed. An implicit function theorem of the type investigated by Nash, Schwartz and Moser follows.     

Second-order characteristic schemes in time and age for a nonlinear age-structured population model

In this paper, we analyze two new second-order characteristic schemes in time and age for an age-structured population model with nonlinear diffusion and reaction. By using the characteristic difference to approximate the transport term and the average along the characteristics to treat the nonlinear spatial diffusion and reaction terms, an implicit second-order characteristic scheme is proposed. To compute the nonlinear approximation system, an explicit second-order characteristic scheme in time and age is further proposed by using the extrapolation technique. The global existence and uniqueness of the solution of the nonlinear approximation scheme are established by using the theory of variation methods, Schauder’s fixed point theorem, and the technique of prior estimates. The optimal error estimates of second order in time and age are strictly proved for both the implicit and the explicit characteristic schemes. Numerical examples are given to illustrate the performance of the methods.     

Linearly implicit methods for nonlinear parabolic equations

We construct and analyze combinations of rational implicit and explicit multistep methods for nonlinear parabolic equations. The resulting schemes are linearly implicit and include as particular cases implicit-explicit multistep schemes as well as the combination of implicit Runge-Kutta schemes and extrapolation. An optimal condition for the stability constant is derived under which the schemes are locally stable. We establish optimal order error estimates.

     

Conservative finite difference methods for fractional Schrödinger–Boussinesq equations and convergence analysis

In this paper, two conservative finite difference schemes for fractional Schrödinger–Boussinesq equations are formulated and investigated. The convergence of the nonlinear fully implicit scheme is established via discrete energy method, while the linear semi‐implicit scheme is analyzed by means of mathematical induction method. Our schemes are proved to preserve the total mass and energy in discrete level. The numerical results are given to confirm the theoretical analysis.     

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.     

Semi‐implicit schemes for transient Navier–Stokes equations and eddy viscosity models

This study presents two computational schemes for the numerical approximation of solutions to eddy viscosity models as well as transient Navier–Stokes equations. The eddy viscosity model is one example of a class of Large Eddy Simulation models, which are used to simulate turbulent flow. The first approximation scheme is a first order single step method that treats the nonlinear term using a semi‐implicit discretization. The second scheme employs a two step approach that applies a Crank–Nicolson method for the nonlinear term while also retaining the semi‐implicit treatment used in the first scheme. A finite element approximation is used in the spatial discretization of the partial differential equations. The convergence analysis for both schemes is discussed in detail, and numerical results are given for two test problems one of which is the two dimensional flow around a cylinder. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Variational Principle for General Diffusion Problems

We employ the Monge–Kantorovich mass transfer theory to study the existence of solutions for a large class of parabolic partial differential equations. We deal with nonhomogeneous nonlinear diffusion problems (of Fokker–Planck type) with time-dependent coefficients. This work greatly extends the applicability of known techniques based on constructing weak solutions by approximation with time-interpolants of minimizers arising from Wasserstein-type implicit schemes. It also generalizes previous results of the authors, where proofs of convergence in the case of a right-hand side in the equation is given by these methods. To prove the existence of weak solutions we establish an interesting maximum principle for such equations. This involves comparison with the solution for the corresponding homogeneous, time-independent equation.     

Numerical Solutions to the Darboux Problem with Functional Dependence

The paper deals with the Darboux problem for the equation D _{xy z} (x,y) = f(x,y,z( _x,y ) where z( _x,y ) is a function defined by . We construct a general class of difference methods for this problem. We prove the existence and uniqueness of solutions to implicit functional difference equations by means of a comparison method; moreover we give an error estimate. The convergence of explicit difference schemes is proved under a general assumption that given functions satisfy nonlinear estimates of the Perron type. Our results are illustrated by a numerical example.     

Conservative finite difference schemes for the generalized Zakharov–Kuznetsov equations

This paper is concerned with the construction of conservative finite difference schemes by means of discrete variational method for the generalized Zakharov–Kuznetsov equations and the numerical solvability of the two-dimensional nonlinear wave equations. A finite difference scheme is proposed such that mass and energy conservation laws associated with the generalized Zakharov–Kuznetsov equations hold. Our arguments are based on the procedure that D. Furihata has recently developed for real-valued nonlinear partial differential equations. Numerical results are given to confirm the accuracy as well as validity of the numerical solutions and then exhibit remarkable nonlinear phenomena of the interaction and behavior of pulse wave solutions.