On a difference scheme of exponential type for a nonlinear singular perturbation problem

Summary A difference scheme of exponential type for solving a nonlinear singular perturbation problem is analysed. Although this scheme is not of monotone type, aL ¹ convergence result is obtained. Relations between this scheme and Engquist-Osher scheme are also discussed.     

On an iterative method for saddle point problems

Convergence behavior of a nested iterative scheme presented in a paper by Bank, Welfert and Yserentant is studied. It is shown that this scheme converges under conditions weaker than that stated in their paper. Received November 18, 1996 / Revised version received March 18, 1997     

A uniformly convergent difference scheme for a semilinear singular perturbation problem

Summary We present a difference scheme for solving a semilinear singular perturbation problem with any number of turning points of arbitrary orders. It is shown that a solution of the scheme converges, uniformly in a perturbation parameter, to that of the continuous problem.     

Relative impatience determines preference between contract bargaining and repeated bargaining

In this paper we consider the effect of the `impatience ratio' I (of the worker discount factor to the firm discount factor) on the preferences of the players between two bargaining schemes in an asymmetric information wage bargaining context. The firm has private information about the worker's value and the worker makes wage demands. In the contact bargaining scheme, a wage demand which is accepted in one period is binding for all future periods (and hence the bargaining ends after acceptance of a wage demand). In the repeated bargaining scheme, the parties continue to bargain irrespective of whether the worker has been hired or not, and any accepted wage demand is only valid for the period in which it was accepted. We establish the following results under the assumption that the worker's value is uniformly distributed on an interval: When the firm is more patient than the worker (I<1) both parties prefer contract bargaining, and when the worker is more patient than the firm (I >1) both prefer repeated bargaining. For any value of I, the preferred type of bargaining gives the lower unemployment.  The work of Bae has already shown that when players are equally patient (I=1) the players are indifferent between the two schemes, regardless of the distribution of the worker's value. This paper shows that Bae's indifference result (Bae, 1991) cannot be extended to unequally patient players. Received: December 1996/Final version: October 1998     

The finite element approximation of a coupled reaction-diffusion problem with non-Lipschitz nonlinearities

Summary. A coupled semilinear elliptic problem modelling an irreversible, isothermal chemical reaction is introduced, and discretised using the usual piecewise linear Galerkin finite element approximation. An interesting feature of the problem is that a reaction order of less than one gives rise to a "dead core" region. Initially, one reactant is assumed to be acting as a catalyst and is kept constant. It is shown that error bounds previously obtained for a scheme involving numerical integration can be improved upon by considering a quadratic regularisation of the nonlinear term. This technique is then applied to the full coupled problem, and optimal and error bounds are proved in the absence of quadrature. For a scheme involving numerical integration, bounds similar to those obtained for the catalyst problem are shown to hold. Received May 25, 1993 / Revised version received July 5, 1994     

On the integration of stiff systems of O.D.E.s using extended backward differentiation formulae

Summary A class of extended backward differentiation formulae suitable for the approximate numerical integration of stiff systems of first order ordinary differential equations is derived. An algorithm is described whereby the required solution is predicted using a conventional backward differentiation scheme and then corrected using an extended backward differentiation scheme of higher order. This approach allows us to developL-stable schemes of order up to 4 andL()-stable schemes of order up to 9. An algorithm based on the integration formulae derived in this paper is illustrated by some numerical examples and it is shown that it is often superior to certain existing algorithms.     

Numerical solutions for some coupled systems of nonlinear boundary value problems

Summary In the well-known Volterra-Lotka model concerning two competing species with diffusion, the densities of the species are governed by a coupled system of reaction diffusion equations. The aim of this paper is to present an iterative scheme for the steady state solutions of a finite difference system which corresponds to the coupled nonlinear boundary value problems. This iterative scheme is based on the method of upper-lower solutions which leads to two monotone sequences from some uncoupled linear systems. It is shown that each of the two sequences converges to a nontrivial solution of the discrete equations. The model under consideration may have one, two or three nonzero solutions and each of these solutions can be computed by a suitable choice of initial iteration. Numerical results are given for these solutions under both the Dirichlet boundary condition and the mixed type boundary condition.     

A stable and conservative finite difference scheme for the Cahn-Hilliard equation

Summary. We propose a stable and conservative finite difference scheme to solve numerically the Cahn-Hilliard equation which describes a phase separation phenomenon. Numerical solutions to the equation is hard to obtain because it is a nonlinear and nearly ill-posed problem. We design a new difference scheme based on a general strategy proposed recently by Furihata and Mori. The new scheme inherits characteristic properties, the conservation of mass and the decrease of the total energy, from the equation. The decrease of the total energy implies boundedness of discretized Sobolev norm of the solution. This in turn implies, by discretized Sobolev's lemma, boundedness of max norm of the solution, and hence the stability of the solution. An error estimate for the solution is obtained and the order is . Numerical examples demonstrate the effectiveness of the proposed scheme. Received July 22, 1997 / Revised version received October 19, 1999 / Published online August 2, 2000     

Convergence of the multigrid full approximation scheme for a class of elliptic mildly nonlinear boundary value problems

Summary The multigrid full approximation scheme (FAS MG) is a well-known solver for nonlinear boundary value problems. In this paper we restrict ourselves to a class of second order elliptic mildly nonlinear problems and we give local conditions, e.g. a local Lipschitz condition on the derivative of the continuous operator, under which the FAS MG with suitably chosen parameters locally converges. We prove quantitative convergence statements and deduce explicit bounds for important quantities such as the radius of a ball of guaranteed convergence, the number of smoothings needed, the number of coarse grid corrections needed and the number of FAS MG iterations needed in a nested iteration. These bounds show well-known features of the FAS MG scheme.     

An asymptotic finite element method for improvement of solutions of boundary layer problems

Summary A scheme that uses singular perturbation theory to improve the performance of existing finite element methods is presented. The proposed scheme improves the error bounds of the standard Galerkin finite element scheme by a factor of O(ⁿ⁺¹) (where is the small parameter andn is the order of the asymptotic approximation). Numerical results for linear second order O.D.E.'s are given and are compared with several other schemes.     

On the numerical analysis of the Von Karman equations: Mixed finite element approximation and continuation techniques

Summary The purpose of this paper is to study the approximation of the Von Karman equations by the mixed finite element scheme of Miyoshi and to follow the solutions arcs at a neighbourhood of the first eigenvalue of the linearized problem. This last problem is solved by a continuation method.     

Analysis of finite difference schemes for unsteady Navier-Stokes equations in vorticity formulation

Summary. In this paper, we provide stability and convergence analysis for a class of finite difference schemes for unsteady incompressible Navier-Stokes equations in vorticity-stream function formulation. The no-slip boundary condition for the velocity is converted into local vorticity boundary conditions. Thom's formula, Wilkes' formula, or other local formulas in the earlier literature can be used in the second order method; while high order formulas, such as Briley's formula, can be used in the fourth order compact difference scheme proposed by E and Liu. The stability analysis of these long-stencil formulas cannot be directly derived from straightforward manipulations since more than one interior point is involved in the formula. The main idea of the stability analysis is to control local terms by global quantities via discrete elliptic regularity for stream function. We choose to analyze the second order scheme with Wilkes' formula in detail. In this case, we can avoid the complicated technique necessitated by the Strang-type high order expansions. As a consequence, our analysis results in almost optimal regularity assumption for the exact solution. The above methodology is very general. We also give a detailed analysis for the fourth order scheme using a 1-D Stokes model. Received December 10, 1999 / Revised version received November 5, 2000 / Published online August 17, 2001     

L 1 andL ∞ Uniform convergence of a difference scheme for a semilinear singular perturbation problem

Summary A nonlinear difference scheme is given for solving a semilinear singularly perturbed two-point boundary value problem. Without any restriction on turning points, the solution of the scheme is shown to be first order accurate in the discreteL ¹ norm, uniformly in the perturbation parameter. When turning points are excluded, the scheme is first order accurate in the discreteL norm, uniformly in the perturbation parameter.Partly supported by the Arts Faculty Research Fund of University College, Cork     

A new extension of boundary elements method for classical elliptic partial differential equations with constant coefficients

This article is devoted to an extension of boundary elements method (BEM) for solving elliptic partial differential equations of general type with constant coefficients. As the fundamental solution of these equations was not available in the literature, BEM was not able to handle them, directly. So the dual reciprocity method (DRM) has been applied to tackle these problems. In this work, a fundamental solution for these equations is obtained and a new formulation is derived to solve them. Besides, we show that the rate of convergence of the new scheme is quadratic when singular (boundary and domain) integrals are calculated, accurately. The new scheme is applicable on complex domains, without needing internal nodes, just same as conventional BEM. So the CPU time of the new scheme is much less than that of the DRM. Numerical examples presented in the article show ability and efficiency of the new scheme in solving two‐dimensional nonhomogenous elliptic boundary value problems, clearly. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2027–2042, 2015     

Thirteen ways to estimate global error

Summary Various techniques that have been proposed for estimating the accumulated discretization error in the numerical solution of differential equations, particularly ordinary differential equations, are classified, described, and compared. For most of the schemes either an outline of an error analysis is given which explains why the scheme works or a weakness of the scheme is illustrated.This research is partially supported by NSF Grant No. MCS-8107046     

Multigrid methods for boundary integral equations

Summary Multigrid methods are applied for solving algebraic systems of equations that occur to the numerical treatment of boundary integral equations of the first and second kind. These methods, originally formulated for partial differential equations of elliptic type, combine relaxation schemes and coarse grid corrections. The choice of the relaxation scheme is found to be essential to attain a fast convergent iterative process. Theoretical investigations show that the presented relaxation scheme provides a multigrid algorithm of which the rate of convergence increases with the dimension of the finest grid. This is illustrated for the calculation of potential flow around an aerofoil.     

Numerical solution of retarded initial value problems: Local and global error and stepsize control

Summary Retarded initial value problems are routinely replaced by an initial value problem of ordinary differential equations along with an appropriate interpolation scheme. Hence one can control the global error of the modified problem but not directly the actual global error of the original problem. In this paper we give an estimate for the actual global error in terms of controllable quantities. Further we show that the notion of local error as inherited from the theory of ordinary differential equations must be generalized for retarded problems. Along with the new definition we are led to developing a reliable basis for a step selection scheme.     

不确定混沌系统的混合投影同步

本文研究了一类不确定混沌(超混沌)系统的混合投影问题.利用自适应方法和Lyapunov稳定性理论,获得了两个恒同或不同混沌系统实现混沌投影同步的一般方法.最后,数值仿真的结果验证了方法的有效性和鲁棒性.     

A third-order semi-discrete genuinely multidimensional central scheme for hyperbolic conservation laws and related problems

Summary. We construct a new third-order semi-discrete genuinely multidimensional central scheme for systems of conservation laws and related convection-diffusion equations. This construction is based on a multidimensional extension of the idea, introduced in [17] – the use of more precise information about the local speeds of propagation, and integration over nonuniform control volumes, which contain Riemann fans. As in the one-dimensional case, the small numerical dissipation, which is independent of , allows us to pass to a limit as . This results in a particularly simple genuinely multidimensional semi-discrete scheme. The high resolution of the proposed scheme is ensured by the new two-dimensional piecewise quadratic non-oscillatory reconstruction. First, we introduce a less dissipative modification of the reconstruction, proposed in [29]. Then, we generalize it for the computation of the two-dimensional numerical fluxes. Our scheme enjoys the main advantage of the Godunov-type central schemes –simplicity, namely it does not employ Riemann solvers and characteristic decomposition. This makes it a universal method, which can be easily implemented to a wide variety of problems. In this paper, the developed scheme is applied to the Euler equations of gas dynamics, a convection-diffusion equation with strongly degenerate diffusion, the incompressible Euler and Navier-Stokes equations. These numerical experiments demonstrate the desired accuracy and high resolution of our scheme. Received February 7, 2000 / Published online December 19, 2000     

A constraint linearization method for nondifferentiable convex minimization

Summary This paper presents a readily implementable algorithm for solving constrained minimization problems involving (possibly nonsmooth) convex functions. The constraints are handled as in the successive quadratic approximations methods for smooth problems. An exact penalty function is employed for stepsize selection. A scheme for automatic limitation of penalty growth is given. Global convergence of the algorithm is established, as well as finite termination for piecewise linear problems. Numerical experience is reported.Sponsored by Program CPBP 02.15