Booster Method for Singularly-Perturbed One-Dimensional Reaction-Diffusion Neumann Problems

A numerical method for singularly-perturbed self-adjoint boundary-value problems for second-order ordinary differential equations subject to Neumann boundary conditions is proposed. In this method (booster method), an asymptotic approximation is incorporated into a finite-difference scheme to improve the numerical solution. Uniform error estimates are derived for this method when implemented in known difference schemes. Numerical examples are presented to illustrate the present method.     

Numerical Analysis of Boundary-Value Problems for Singularly-Perturbed Differential-Difference Equations with Small Shifts of Mixed Type

In this paper, we use a numerical method to solve boundary-value problems for a singularly-perturbed differential-difference equation of mixed type, i.e., containing both terms having a negative shift and terms having a positive shift. Similar boundary-value problems are associated with expected first exit time problems of the membrane potential in models for the neuron. The stability and convergence analysis of the method is given. The effect of a small shift on the boundary-layer solution is shown via numerical experiments. The numerical results for several test examples demonstrate the efficiency of the method.     

Solution of Two-Point Boundary-Value Problems Using the Differential Transformation Method

This paper considers two-point boundary-value problems using the differential transformation method. An iterative procedure is proposed for both the linear and nonlinear cases. Using the proposed approach, an analytic solution of the two-point boundary-value problem, represented by an mth-order Taylor series expansion, can be obtained throughout the prescribed range.     

Spline Techniques for Solving Singularly-Perturbed Nonlinear Problems on Nonuniform Grids

A numerical method based on cubic splines with nonuniform grid is given for singularly-perturbed nonlinear two-point boundary-value problems. The original nonlinear equation is linearized using quasilinearization. Difference schemes are derived for the linear case using a variable-mesh cubic spline and are used to solve each linear equation obtained via quasilinearization. Second-order uniform convergence is achieved. Numerical examples are given in support of the theoretical results.     

Asymptotic Initial-Value Method for Singularly-Perturbed Boundary-Value Problems for Second-Order Ordinary Differential Equations

A computational method is presented to solve a class of nonturning-point singularly-perturbed two-point boundary-value problems for second-order ordinary differential equations with a small parameter multiplying the highest derivative, subject to Dirichlet-type boundary conditions. In this method, first we construct a zeroth order asymptotic expansion for the solution of the given boundary-value problem. Then, this problem is integrated to get an equivalent initial-value problem for first-order ordinary differential equations. This initial-value problem is solved by either a classical method or a fitted operator method after approximating some of the terms in the differential equations by using the zeroth order asymptotic expansion. This method is effective and easy to implement. An error estimate is derived for the numerical solution. Examples are given to illustrate the method.     

二阶常微分方程初值问题的Laguerre-Gauss配置法

研究二阶常微分方程初值问题的数值解法．该文中基于Laguerre-Gauss插值设计了一类新的配置法, 它易于计算,且特别适用于非线性问题．该文中分析了二种不同情况时的收敛性,并应用Laguerre-Gauss插值的最新结果,证明了它的谱精度．该文还提供了一种多步配置法,它既简化了计算,又保持同样的谱精度．数值结果显示了这些算法的高精度．     

Quasilinearization Method for Nonlinear Elliptic Boundary-Value Problems

The quasilinearization method is developed for strong solutions of semilinear and nonlinear elliptic boundary-value problems. We obtain two monotone, L^p-convergent sequences of approximate solutions. The order of convergence is two. The tools are some results on the abstract quasilinearization method and from weakly–near operators theory.     

非线性Klein-Gordon方程的广义Jacobi谱配置方法

构造非线性Klein-Gordon方程的广义Jacobi谱配置格式,并给出相应收敛性分析.文中的方法和技巧为设计和分析各类线性与非线性偏微分方程的谱配置格式提供了有效的框架.     

Two-Point Boundary Value Problems for Duffing Equations across Resonance

In this paper, we consider the equation y″+f(x,y)=0 with a nonresonance condition of the form A≤f _y (x,y)≤B, where (k−1)²<A≤k ²<⋅⋅⋅<m ²≤B<(m+1)², k,m∈ℤ⁺. With optimal control theory and the Schauder fixed-point theorem, by introducing a new cost functional, we obtain a new existence and uniqueness result for the above equation with two-point boundary-value conditions. This work was supported by NSFC Grant 10501017 and 985 Project of Jilin University.     

Asymptotic Solution of a Boundary-Value Problem for Linear Singularly-Perturbed Functional Differential Equations Arising in Optimal Control Theory

The Hamiltonian boundary-value problem, associated with a singularly-perturbed linear-quadratic optimal control problem with delay in the state variables, is considered. A formal asymptotic solution of this boundary-value problem is constructed by application of the boundary function method. The justification of this asymptotic solution is done. The asymptotic solution of the Hamiltonian boundary-value problem is constructed and justified assuming boundary-layer stabilizability and detectability.     

Initial-Value Technique for Singularly Perturbed Boundary-Value Problems for Second-Order Ordinary Differential Equations Arising in Chemical Reactor Theory

An initial-value technique is presented for solving singularly perturbed two-point boundary-value problems for linear and semilinear second-order ordinary differential equations arising in chemical reactor theory. In this technique, the required approximate solution is obtained by combining solutions of two terminal-value problems and one initial-value problem which are obtained from the original boundary-value problem through asymptotic expansion procedures. Error estimates for approximate solutions are obtained. Numerical examples are presented to illustrate the present technique.     

A Collocation Code for Singular Boundary Value Problems in Ordinary Differential Equations

We present a MATLAB package for boundary value problems in ordinary differential equations. Our aim is the efficient numerical solution of systems of ODEs with a singularity of the first kind, but the solver can also be used for regular problems. The basic solution is computed using collocation methods and a new, efficient estimate of the global error is used for adaptive mesh selection. Here, we analyze some of the numerical aspects relevant for the implementation, describe measures to increase the efficiency of the code and compare its performance with the performance of established standard codes for boundary value problems.     

Booster Method for Singularly-Perturbed One-Dimensional Convection-Diffusion Neumann Problems

An improved numerical method for singularly-perturbed two-point boundary-value problems for second-order ordinary differential equations subject to Neumann-type boundary conditions is proposed. In this method, an asymptotic approximation is incorporated into a finite-difference scheme to improve the numerical solution. Uniform error estimates are derived when implemented in known difference schemes. Numerical results are presented in support of the proposed method.     

Efficient Collocation Schemes for Singular Boundary Value Problems

We discuss an error estimation procedure for the global error of collocation schemes applied to solve singular boundary value problems with a singularity of the first kind. This a posteriori estimate of the global error was proposed by Stetter in 1978 and is based on the idea of Defect Correction, originally due to Zadunaisky. Here, we present a new, carefully designed modification of this error estimate which not only results in less computational work but also appears to perform satisfactorily for singular problems. We give a full analytical justification for the asymptotical correctness of the error estimate when it is applied to a general nonlinear regular problem. For the singular case, we are presently only able to provide computational evidence for the full convergence order, the related analysis is still work in progress. This global estimate is the basis for a grid selection routine in which the grid is modified with the aim to equidistribute the global error. This procedure yields meshes suitable for an efficient numerical solution. Most importantly, we observe that the grid is refined in a way reflecting only the behavior of the solution and remains unaffected by the unsmooth direction field close to the singular point.     

Efficient Implementation of a Spectral-Tau Method for a Class of Two-Point Boundary-Value and Initial-Boundary-Value Problems

A simple and efficient spectral method for solving the second, third order and fourth order elliptic equations with variable coefficients and nonlinear differential equations is presented. It is different from spectral-collocation method which leads to dense, ill-conditioned matrices. The spectral method in this paper solves for the coefficients of the solution in a Chebyshev series, leads to discrete systems with special structured matrices which can be factorized and solved efficiently. We also extend the method to boundary value problems in two space dimensions and solve 2-D separable equation with variable coefficients. As an application, we solve Cahn-Hilliard equation iteratively via first-order implicit time discretization scheme. Ample numerical results indicate that the proposed method is extremely accurate and efficient.     

n阶非线性常微分方程的非线性两点边值问题

本文利用打靶法,给出了ｎ阶非线性常微分方程具有非线性两点边界条件的边值问题存在解与存在唯一解的一般性结果,并将所得结果应用于Ｌｉｐｓｃｈｉｔｚ方程的两点过值问题,给出了存在解与存在唯一解的具体的充分条件。     

Legendre-Gauss Spectral Collocation Method for Second Order Nonlinear Delay Differential Equations

In this paper, we present and analyze a single interval Legendre-Gauss spectral collocation method for solving the second order nonlinear delay differential equations with variable delays. We also propose a novel algorithm for the single interval scheme and apply it to the multiple interval scheme for more efficient implementation. Numerical examples are provided to illustrate the high accuracy of the proposed methods.     

Shooting Methods for Numerical Solution of Nonlinear Stochastic Boundary-Value Problems

Abstract

In the present investigation, shooting methods are described for numerically solving nonlinear stochastic boundary-value problems. These stochastic shooting methods are analogous to standard shooting methods for numerical solution of ordinary deterministic boundary-value problems. It is shown that the shooting methods provide accurate approximations. An error analysis is performed and computational simulations are described.     

A Chebyshev-Gauss Spectral Collocation Method for Ordinary Differential Equations

In this paper, we introduce an efficient Chebyshev-Gauss spectral collocation method for initial value problems of ordinary differential equations. We first propose a single interval method and analyze its convergence. We then develop a multi-interval method. The suggested algorithms enjoy spectral accuracy and can be implemented in stable and efficient manners. Some numerical comparisons with some popular methods are given to demonstrate the effectiveness of this approach.