Factorization method for linear and quasilinear singularly perturbed boundary value problems for ordinary differential equations

For linear singularly perturbed boundary value problems, we come up with a method that reduces solving a differential problem to a discrete (difference) problem. Difference equations, which are an exact analog of differential equations, are constructed by the factorization method. Coefficients of difference equations are calculated by solving Cauchy problems for first-order differential equations. In this case nonlinear Ricatti equations with a small parameter are solved by asymptotic methods, and solving linear equations reduces to computing quadratures. A solution for quasilinear singularly perturbed equations is obtained by means of an implicit relaxation method. A solution to a linearized problem is calculated by analogy with a linear problem at each iterative step. The method is tested against solutions to the known Lagerstrom-Cole problem.     

Numerical solution of nonlinear two-dimensional integral equations using rationalized Haar functions

Two-dimensional rationalized Haar (RH) functions are applied to the numerical solution of nonlinear second kind two-dimensional integral equations. Using bivariate collocation method and Newton–Cotes nodes, the numerical solution of these equations is reduced to solving a nonlinear system of algebraic equations. Also, some numerical examples are presented to demonstrate the efficiency and accuracy of the proposed method.     

Numerical solution for the weakly singular Fredholm integro-differential equations using Legendre multiwavelets

An effective method based upon Legendre multiwavelets is proposed for the solution of Fredholm weakly singular integro-differential equations. The properties of Legendre multiwavelets are first given and their operational matrices of integral are constructed. These wavelets are utilized to reduce the solution of the given integro-differential equation to the solution of a sparse linear system of algebraic equations. In order to save memory requirement and computational time, a threshold procedure is applied to obtain the solution to this system of algebraic equations. Through numerical examples, performance of the present method is investigated concerning the convergence and the sparseness of the resulted matrix equation.     

Study on convergence of hybrid functions method for solution of nonlinear integral equations

In this article, we present the uniform convergence analysis and accuracy estimation of hybrid functions (HFs) method for finding the solution of nonlinear Volterra and Fredholm integral equations. The properties of HFs which consist of block-pulse functions (BPFs) and Legendre polynomials are used to reduce the solution of nonlinear integral equations to the solution of algebraic equations. The superiority and accuracy of the HFs method to BPF and Legendre polynomial methods are illustrated through some numerical examples.     

A unified approach to teaching quadratic and cubic equations

A simple method for teaching the algebraic solution of cubic equations is presented. The approach is via ‘completion of the cube’. It is found that this method is readily accepted by students already familiar with completion of the square as a method for quadratic equations. Many people are surprised at how simple the algebraic solution of cubic equations turns out to be. It is not normally necessary to resort to numerical techniques for these equations.     

Triangular functions (TF) method for the solution of nonlinear Volterra–Fredholm integral equations

A numerical method based on an m-set of general, orthogonal triangular functions (TF) is proposed to approximate the solution of nonlinear Volterra–Fredholm integral equations. The orthogonal triangular functions are utilized as a basis in collocation method to reduce the solution of nonlinear Volterra–Fredholm integral equations to the solution of algebraic equations. Also a theorem is proved for convergence analysis. Some numerical examples illustrate the proposed method.     

对称正则长波方程的行波解分岔

通过平面动力系统的方法讨论了对称正则长波方程的分岔问题.得到了该方程的分岔条件,在一些参数的具体值的情况下给出相图并通过微分方程的数值模拟方法模拟出了该方程的周期行波解、孤立行波解及无界行波解.     

基于拟Shannon小波浅水长波近似方程组的数值解

本文研究了浅水长波近似方程组初边值问题的数值解.利用小波多尺度分析和区间拟Shannon小波,对浅水长波近似方程组空间导数实施空间离散,用时间步长自适应精细积分法对其变换所的非线性常微分方程组进行求解,得到了浅水长波近似方程组的数值解,并将此方法计算的结果与其解析解进行比较和验证.     

Boundary element solution of the two-dimensional groundwater flow equation with stochastic free surface boundary condition

An innovative approach to the approximate solution of stochastic partial differential equations in groundwater flow is presented. The method uses a formulation of the Ito's lemma in Hilbert spaces to derive partial differential equations satisfying the moments of the solution process. Since the moments equations are deterministic, they could be solved by any analytical or numerical method existing in the literature. This permits the analysis and solution of stochastic partial differential equations occurring in two-dimensional or three-dimensional domains of any geometrical shape. The method is tested for the first time in the present paper through a practical application in a sandy phreatic aquifer at the Chalk River Nuclear Laboratories, Ontario, Canada. The equation solved is the two-dimensional LaPlace equation with a dynamic, randomly perturbed, free surface boundary condition. The moments equations are derived and solved by using the boundary integral equation method. A comparison is made with a previous analytical solution obtained by applying the randomly forced one-dimensional Boussinesq equation, and some observations on modeling procedures are given.     

Differential quadrature solution of heat- and mass-transfer equations

The heat- and mass-transfer equations have an important role in various thermal and diffusion processes. These equations are nonlinear, due to the solution dependent diffusion coefficient and the source term. In this study, one- and two-dimensional nonlinear heat- and mass-transfer equations are solved numerically. To this end, the differential quadrature method is used to discretize the problem spatially and the resulting nonlinear system of ordinary differential equations in time are solved using the Runge–Kutta method. The solution is improved in time iteratively by solving considerably small sized linear system of resulting equations. To demonstrate its usefulness and accuracy, the proposed method is applied to four test problems, involving different nonlinearities.     

Modified Camassa–Holm and Degasperis–Procesi Equations Solved by Adomian’s Decomposition Method and Comparison with HPM and Exact Solutions

Many physical and scientific phenomena are modeled by nonlinear partial differential equations (NPDEs); it is difficult to handle nonlinear part of these equations. Recently some analytical methods are applied to solve such equations. In this work, modified Camassa–Holm and Degasperis–Procesi equation is studied. Adomian’s decomposition method (ADM) is applied to obtain solution of this equation. The results are compared to those of homotopy perturbation method (HPM) and exact solution. The study highlights the significant features of the employed method and its ability to handle nonlinear partial differential equations.     

Application to differential transformation method for solving systems of differential equations

In this paper, we present an analytical solution for different systems of differential equations by using the differential transformation method. The convergence of this method has been discussed with some examples which are presented to show the ability of the method for linear and non-linear systems of differential equations. We begin by showing how the differential transformation method applies to a non-linear system of differential equations and give two examples to illustrate the sufficiency of the method for linear and non-linear stiff systems of differential equations. The results obtained are in good agreement with the exact solution and Runge–Kutta method. These results show that the technique introduced here is accurate and easy to apply.     

Stein implicit Runge-Kutta methods with high stage order for large-scale ordinary differential equations

The Runge-Kutta method is one of the most popular implicit methods for the solution of stiff ordinary differential equations. For large problems, the main drawback of such methods is the cost required at each integration step for computing the solution of a nonlinear system of equations. In this paper, we propose to reduce the cost of the computation by transforming the linear systems arising in the application of Newton's method to Stein matrix equations. We propose an iterative projection method onto block Krylov subspaces for solving numerically such Stein matrix equations. Numerical examples are given to illustrate the performance of our proposed method.     

Numerical solution of some partial differential equations by means of a deterministic method of approximate functional integration

A numerical method of solution of some partial differential equations is presented. The method is based on representation of Green functions of the equations in the form of functional integrals and subsequent approximate calculation of the integrals with the help of a deterministic approach. In this case the solution of the equations is reduced to evaluation of usual (Riemann) integrals of relatively low multiplicity. A procedure allowing one to increase accuracy of the solutions is suggested. The features of the method are investigated on examples of numerical solution of the Schrödinger equation and related diffusion equation.     

Conditions for partitioning a system of differential algebraic equations into weakly coupled subsystems

Systems of differential algebraic equations are examined. A method is proposed for transforming the rectangular matrix of algebraic equations to block diagonal form. This method ensures the prescribed accuracy of the solution with respect to the original system of equations.     

The Elementary Solution of Triple Integral Equations and The Solution of Triple Series Equations Involving Associated Legendre Polynomials and their Application

A method is developed for the formal solution of an important class of triple integral equations involving Bessel functions. The solution of the triple integral equations is reduced to two simultaneous Fredholm integral equations and the results obtained are simpler than those of other authors and also superior for the purposes of solution by iteration. In the same manner the formal solution of triple series equations involving associated Legendre polynomials is presented. The solution of the problem is reduced to that of solving a Fredholm integral equation of the first kind. Finally to illustrate the application of the results an electrostatic problem is discussed.     

广义Benjamin-Bona-Mahony方程Fourier谱方法的大时间误差估计

§ 1　ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＢｅｎｊａｍｉｎ Ｂｏｎａ Ｍａｈｏｎｙｅｑｕａｔｉｏｎｕｔ+ｕｘ+ｕｕｘ -ｕｘｘ-ｕｘｘｔ =0 ( 1 .1 )ｉｎｃｏｒｐｏｒａｔｅｓｎｏｎｌｉｎｅａｒｄｉｓｐｅｒｓｉｖｅａｎｄｄｉｓｓｉｐａｔｉｖｅｅｆｆｅｃｔｓ ,ａｎｄｈａｓｂｅｅｎｐｒｏｐｏｓｅｄａｓａｍｏｄｅｌｆｏｒｂｏｔｈｔｈｅｂｏｒｅｐｒｏｐａｇａｔｉｏｎａｎｄｔｈｅｗａｔｅｒｗａｖｅｓ[1,2 ] .Ｔｈｅｅｘｉｓｔｅｎｃｅａｎｄｕｎｉｑｕｅｎｅｓｓｏｆｓｏｌｕｔｉｏｎｓｆｏｒｔｈｉｓｅｑｕａｔｉｏｎｈａｖｅｂｅｅ…     

The primal mixed finite element method and the LBB condition

Spurious or kinematic modes have posed a major obstacle to the implementation of the mixed finite element method. This research shows that spurious modes resulting from the approximation spaces not satisfying the LBB condition do not prevent a well posed problem. When the LBB condition is not satisfied, the resulting matrix equations are singular. A direct solution method is presented for the efficient solution of the possibly singular equations. Orthogonal flux basis functions are introduced to simplify the problem. Then the solution procedure is based on nested domain decomposition. This solution procedure is shown to be competitive with direct solution methods for the displacement finite element method. Examples are included to demonstrate various aspects of the LBB condition and the solution procedure.     

Banach空间半线性发展方程的周期解

把上、下解方法引入到有序Ｂａｎａｃｈ空间中的半线性发展方程周期解问题，利用正算子半群特征与单调迭代程序，获得了最大周期解与最小周期解的存在性．所得的结果概括和推广了常微分方程与偏微分方程中的有关结论     

Existence of solutions and Gevrey class regularity for Leray-alpha equations

In this paper we consider some equations similar to Navier-Stokes equations, the three-dimensional Leray-alpha equations with space periodic boundary conditions. We establish the regularity of the equations by using the classical Faedo-Galerkin method. Our argument shows that there exist an unique weak solution and an unique strong solution for all the time for the Leray-alpha equations, furthermore, the strong solutions are analytic in time with values in the Gevrey class of functions (for the space variable). The relations between the Leray-alpha equations and the Navier-Stokes equations are also considered.