A NUMERICAL PERTURBATION EXPANSION METHOD FOR THE SOLUTION OF A LIENARD-TYPE INITIAL-VALUE PROBLEM WITH PERIODIC DECAYING FORCING TERM

A numerical perturbation expansion method is developed, analysed and implemented for the numerical solution of a second-order initial-value problem. The differential equation in this problem exhibits cubic damping, a cubic restoring force and a decaying forcing-term which is periodic with constant frequency. The method is compared with the numerical method by Twizell [1]. In fact, the later is first perturbation approximate solution in the present paper.     

Numerical solution of a nonlinear reaction diffusion equation

In this paper, the authors propose a numerical method to compute the solution of a Cauchy problem with blow-up of the solution. The problem is split in two parts: a hyperbolic problem which is solved by using Hopf and Lax formula and a parabolic problem solved by a backward linearized Euler method in time and a finite element method in space. It is proved that the numerical solution blows up in a finite time as the exact solution and the support of the approximation of a self-similar solution remains bounded. The convergence of the scheme is obtained.     

A recurrence method for a special class of continuous time linear programming problems

This article studies a numerical solution method for a special class of continuous time linear programming problems denoted by (SP). We will present an efficient method for finding numerical solutions of (SP). The presented method is a discrete approximation algorithm, however, the main work of computing a numerical solution in our method is only to solve finite linear programming problems by using recurrence relations. By our constructive manner, we provide a computational procedure which would yield an error bound introduced by the numerical approximation. We also demonstrate that the searched approximate solutions weakly converge to an optimal solution. Some numerical examples are given to illustrate the provided procedure.     

Finite-difference methods for solving the reaction-diffusion equations of a simple isothermal chemical system

Numerical methods are proposed for the numerical solution of a system of reaction-diffusion equations, which model chemical wave propagation. The reaction terms in this system of partial differential equations contain nonlinear expressions. Nevertheless, it is seen that the numerical solution is obtained by solving a linear algebraic system at each time step, as opposed to solving a nonlinear algebraic system, which is often required when integrating nonlinear partial differential equations. The development of each numerical method is made in the light of experience gained in solving the system of ordinary differential equations, which model the well-stirred analogue of the chemical system. The first-order numerical methods proposed for the solution of this initialvalue problem are characterized to be implicit. However, in each case it is seen that the numerical solution is obtained explicitly. In a series of numerical experiments, in which the ordinary differential equations are solved first of all, it is seen that the proposed methods have superior stability properties to those of the well-known, first-order, Euler method to which they are compared. Incorporating the proposed methods into the numerical solution of the partial differential equations is seen to lead to two economical and reliable methods, one sequential and one parallel, for solving the travelling-wave problem. © 1994 John Wiley & Sons, Inc.     

基于GA-Chebyshev神经网络的分数阶Bagley-Torvik方程数值解法

本文基于现有的切比雪夫神经网络,提出了一种利用遗传算法优化切比雪夫神经网络求解分数阶Bagley-Torvik方程数值解的新方法,结合多点处的泰勒公式原理,给出数值解的一般形式,将原问题转化为求解无约束最小化问题.与现有数值方法的数值结果进行比较表明了本文方法的可行性和有效性,为分数阶微分方程中类似问题的求解提供了新的思路.     

Two-dimensional integral-algebraic systems: Analysis and computational methods

In this article we formulate sufficient conditions for the existence and uniqueness of solution to systems of two-dimensional Volterra integral equations, in which the coefficient of the main term is a singular matrix. A numerical method is introduced which can be applied to approximate the solution when the given conditions are satisfied. The convergence of this method is proved and illustrated by numerical examples.     

Two-dimensional integral–algebraic systems: Analysis and computational methods

In this article we formulate sufficient conditions for the existence and uniqueness of solution to systems of two-dimensional Volterra integral equations, in which the coefficient of the main term is a singular matrix. A numerical method is introduced which can be applied to approximate the solution when the given conditions are satisfied. The convergence of this method is proved and illustrated by numerical examples.     

A starting method for the numerical solution of Volterra's integral equation of the second kind

The acquisition of starting values is one of the chief difficulties encountered in computing a numerical solution of Volterra's integral equation of the second kind by a multi-step method. The object of this note is to present a procedure which is derived from certain quadrature formulas and which provides these starting values, to provide a sufficient condition for the approximate solution to be unique, to bound the approximate solution and the error, and to give a numerical example.     

Analysis of traffic flow by the finite element method

A finite elernent methodology is developed for the numerical solution of traffic flow problems encountered in arterial streets. The simple continuum traffic flow model consisting of the equation of continuity and an equilibrium flow-density relationship is adopted. A Galerkin type finite element method is used to formulate the problem in discrete form and the solution is obtained by a step-by-step time integration in conjunction with the Newton-Raphson method. The proposed finite element methodology, which is of the shock capturing type, is applied to flow traffic problems. Two numerical examples illustrate the method and demonstrate its advantages over other analytical or numerical techniques.     

A research into the numerical method of Dirichlet's problem of complex Monge-Ampère equation on Cartan-Hartogs domain of the third type

Monge-Ampère equation is a nonlinear equation with high degree, therefore its numerical solution is very important and very difficult. In present paper the numerical method of Dirichlet's problem of Monge-Ampère equation on Cartan-Hartogs domain of the third type is discussed by using the analytic method. Firstly, the Monge-Ampère equation is reduced to the nonlinear ordinary differential equation, then the numerical method of the Dirichlet problem of Monge-Ampère equation becomes the numerical method of two point boundary value problem of the nonlinear ordinary differential equation. Secondly, the solution of the Dirichlet problem is given in explicit formula under the special case, which can be used to check the numerical solution of the Dirichlet problem.     

The Numerical Solution of the Fourth Boundary Value Problem for Parabolic Partial Differential Equations

A direct method for the numerical solution of the implicit finitedifference equations derived from a parabolic differential equationwith periodic spatial boundary conditions is presented in algorithmicfrom. Consideration is given to the stability of the roundingerrors involved in the solution process and numerical resultsare derived which compare favourably with those obtained fromthe analytical solution and a matrix spectral resolution methodwhich is closely allied to the method of lines.     

Convergence and stability of a numerical method for micromagnetics

The convergence and stability of a numerical method, which applies a nonconforming finite element method and an artificial boundary method to a multi-atomic Young measure relaxation model, for micromagnetics are analyzed. By revealing some key properties of the solution sets of both the continuous and discrete problems, we show that our numerical method is stable, and the solution set of the continuous problem is well approximated by those of the discrete problems. The performance of our method is also illustrated by some numerical examples. The research was supported in part by the Major State Basic Research Projects (2005CB321701), NSFC projects (10431050, 10571006, 10528102 and 10871011) and RFDP of China.     

The numerical solution of linear time-dependent partial differential equations by the Laplace transform and finite difference approximations

A feasible method is presented for the numerical solution of a large class of linear partial differential equations which may have source terms and boundary conditions which are time-varying. The Laplace transform is used to eliminate the time-dependency and to produce a subsidiary equation which is then solved in complex arithmetic by finite difference methods. An effective numerical Laplace transform inversion algorithm gives the final solution at each spatial mesh point for any specified set of values of t. The single-step property of the method obviates the need to evaluate the solution at a large number of unwanted intermediate time points. The method has been successfully applied to a variety of test problems and, with two alternative numerical Laplace transform inversion algorithms, has been found to give results of good to excellent accuracy. It is as accurate as other established finite difference methods using the same spatial grid. The algorithm is easily programmed and the same program handles equations of parabolic and hyperbolic type.     

A numerical method to compute the dissolution of second phases in ternary alloys

Dissolution of stoichiometric multi-component particles in ternary alloys is an important process occurring during the heat treatment of as-cast aluminium alloys prior to hot extrusion. A mathematical model is proposed to describe such a process. In this model an equation is given to determine the position of the particle interface in time, using two diffusion equations which are coupled by nonlinear boundary conditions at the interface. Some results concerning existence, uniqueness, and monotonicity are given. Furthermore, for an unbounded domain an analytical approximation is derived. The main part of this work is the development of a numerical solution method. Finite differences are used on a grid which changes in time. The discretization of the boundary conditions is important to obtain an accurate solution. The resulting nonlinear algebraic system is solved by the Newton-Raphson method. Numerical experiments illustrate the accuracy of the numerical method. The numerical solution is compared with the analytical approximation.     

带有不连续系数的线性输运方程差分格式的收敛性

提出了一个基于三角形网格的显式差分格式逼近带有不连续系数的线性输运方程. 通过对数值解的有界性、TVD(total variation decreasing)和空间、时间方向的平移估计, 利用Kolmogorov紧性原理证明了数值解在L¹_loc模下收敛于初值问题的唯一弱解.从而得到了初值问题解的存在唯一性和关于初值的稳定性. 数值算例表明本文提出的格式计算方便而且比 Lax-Friedrichs格式更有效.       

Numerical Method for Singularly Perturbed Third Order Ordinary Differential Equations of Convection-Diffusion Type

In this paper, we have proposed a numerical method for Singularly Perturbed Boundary Value Problems (SPBVPs) of convection-diffusion type of third order Ordinary Differential Equations (ODEs) in which the SPBVP is reduced into a weakly coupled system of two ODEs subject to suitable initial and boundary conditions. The numerical method combines boundary value technique, asymptotic expansion approximation, shooting method and finite difference scheme. In order to get a numerical solution for the derivative of the solution, the domain is divided into two regions namely inner region and outer region. The shooting method is applied to the inner region while standard finite difference scheme (FD) is applied for the outer region. Necessary error estimates are derived for the method. Computational efficiency and accuracy are verified through numerical examples. The method is easy to implement and suitable for parallel computing.     

A numerical solution of a two-dimensional IHCP

In this paper, we will first study the existence and uniqueness of the solution of a two-dimensional inverse heat conduction problem (IHCP) which is severely ill-posed, i.e., the solution does not depend continuously on the data. We propose a stable numerical approach based on the finite-difference method and the least-squares scheme to solve this problem in the presence of noisy data. We prove the convergence of the numerical solution, then to regularize the resultant ill-conditioned linear system of equations, we apply the Tikhonov regularization 0th, 1st and 2nd method to obtain the stable numerical approximation to the solution. The stability and accuracy of the scheme presented is evaluated by comparison with the Singular Value Decomposition (SVD) method.     

Numerical solution for the model of vibrating elastic membrane with moving boundary

In this work, we are interested in obtaining an approximated numerical solution for the model of vibrating elastic membranes with moving boundary. The model is an extension of Kirchhoff’s model, which takes into account the change of size during the vibration. We apply the finite element method with a finite difference method in time to obtain an approximated numerical solution. Some numerical experiments are presented to show the effect of moving boundary effects in vibrating elastic membranes.     

The Motion of a Viscous Fluid Past an Impulsively Started Semi-infinite Flat Plate

It is shown that by employing a suitable numerical method, theproblem of determining the motion of a viscous fluid past asemi-infinite flat plate which is started impulsively from restwith constant velocity parallel to itself can be solved in termsof similarity variables. The numerical solution is comparedwith previous numerical and theoretical work on the problem.The final decay to the steady-state solution described by theBlasius velocity profile is exponential in character and isfound to be substantially in agreement with theoretical predictions.