Revisit on partial solutions in the Adomian decomposition method: Solving heat and wave equations

The Adomian decomposition method is considered in application to heat and wave equations. The so-called partial solution technique is used. It is shown that the fundamental equation of the method is well defined only for certain types of boundary conditions. In cases involving inhomogeneous boundary conditions, improper results may be obtained by former method. This paper presents a further insight into partial solutions in the decomposition method, and the resolution of such cases.     

板弯曲问题的具两组高阶基本解序列的MRM方法

讨论了双参数地基上薄板弯曲问题.利用两组高阶基本解序列,即调和及重调和基本解序列,采用多重替换方法(MRM方法),得到了板弯曲问题的MRM边界积分方程.证明了该方程与边值问题的常规边界积分方程是一致的.因此由常规边界积分方程的误差估计即可得到板弯曲问题MRM方法的收敛性分析.此外该方法还可推广到具多组高阶基本解序列的情形.     

A fundamental solution method for three-dimensional viscous flow problems with obstacles in a periodic array

In this paper, we propose a fundamental solution method for three-dimensional viscous flow problems with obstacles in a periodic array. Our problem is mathematically a boundary value problem of the Stokes equation with periodic boundary conditions, to which it is difficult to give a good approximation by the ordinary fundamental solution method. Our method gives an approximate solution by a linear combination of the periodic fundamental solutions. In addition, we can compute the drag forces on the obstacles by using the data obtained in our method. Numerical examples for the problems of flows past spheres show the effectiveness of our method.     

Pasternak地基上简支板振动问题的准格林函数方法

提出一种新的数值方法——准格林函数方法．以Pasternak地基上简支多边形薄板的振动问题为例,详细阐明了准格林函数方法的思想．即利用问题的基本解和边界方程构造一个准格林函数,这个函数满足了问题的齐次边界条件,采用格林公式将Pasternak地基上薄板自由振动问题的振型控制微分方程化为两个耦合的第二类Fredholm积分方程．边界方程有多种选择,在选定一种边界方程的基础上,可以通过建立一个新的边界方程来表示问题的边界,以克服积分核的奇异性,最后由积分方程的离散化方程组有非平凡解的条件,求得固有频率．数值方法表明,该方法具有较高的精度．     

Boundary element method and wave equation

In this paper the boundary integral expression for a one-dimensional wave equation with homogeneous boundary conditions is developed. This is done using the time dependent fundamental solution of the corresponding hyperbolic partial differential equation. The integral expression developed is a generalized function with the same form as the well-known D'Alembert formula. The derivatives of the solution and some useful invariants on the characteristics of the partial differential equation are also calculated.The boundary element method is applied to find the numerical solution. The results show excellent agreement with analytical solutions.A multi-step procedure for large time steps which can be used in the boundary element method is also described.In addition, the way in which boundary conditions are introduced during the time dependent process is explained in detail. In the Appendix the main properties of Dirac's delta function and the Heaviside unit step function are described.     

A method based on meshless approach for the numerical solution of the two‐space dimensional hyperbolic telegraph equation

In the current article, we investigate the RBF solution of second‐order two‐space dimensional linear hyperbolic telegraph equation. For this purpose, we use a combination of boundary knot method (BKM) and analog equation method (AEM). The BKM is a meshfree, boundary‐only and integration‐free technique. The BKM is an alternative to the method of fundamental solution to avoid the fictitious boundary and to deal with low accuracy, singular integration and mesh generation. Also, on the basis of the AEM, the governing operator is substituted by an equivalent nonhomogeneous linear one with known fundamental solution under the same boundary conditions. Finally, several numerical results and discussions are demonstrated to show the accuracy and efficiency of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Treatment of singularities in the method of fundamental solutions for two-dimensional Helmholtz-type equations

We investigate a meshless method for the accurate and non-oscillatory solution of problems associated with two-dimensional Helmholtz-type equations in the presence of boundary singularities. The governing equation and boundary conditions are approximated by the method of fundamental solutions (MFS). It is well known that the existence of boundary singularities affects adversely the accuracy and convergence of standard numerical methods. The solutions to such problems and/or their corresponding derivatives may have unbounded values in the vicinity of the singularity. This difficulty is overcome by subtracting from the original MFS solution the corresponding singular functions, without an appreciable increase in the computational effort and at the same time keeping the same MFS approximation. Four examples for both the Helmholtz and the modified Helmholtz equations are carefully investigated and the numerical results presented show an excellent performance of the approach developed.     

STURM-LIOUVILLE PROBLEMS WITH EIGENDEPENDENT BOUNDARY AND TRANSMISSIONS CONDITIONS

1IntroductionIt is well-known that the Sturmian Theory is an important aid in solving many problemsin mathematical physics.Therefore this theory is one of the most actual and extensivelydeveloped field in spectral analysis of boundary-value problems of St…     

Singular boundary method for 3D time-harmonic electromagnetic scattering problems

A frequency domain singular boundary method is presented for solving 3D time-harmonic electromagnetic scattering problem from perfect electric conductors. To avoid solving the coupled partial differential equations with fundamental solutions involving hypersingular terms, we decompose the governing equation into a system of independent Helmholtz equations with mutually coupled boundary conditions. Then the singular boundary method employs the fundamental solutions of the Helmholtz equations to approximate the scattered electric field variables. To desingularize the source singularity in the fundamental solutions, the origin intensity factors are introduced. In the novel formulation, only the origin intensity factors for fundamental solutions of 3D Helmholtz equations and its derivatives need to be considered which have been derived in the paper. Several numerical examples involving various perfectly conducting obstacles are carried out to demonstrate the validity and accuracy of the present method.     

The method of fundamental solutions for the calculation of the eigenvalues of the Helmholtz equation

We investigate the application of the method of fundamental solutions (MFS) for the calculation of the eigenvalues of the Helmholtz equation in the plane subject to homogeneous Dirichlet boundary conditions. We present results for circular and rectangular geometries.     

Analysis of vibrations in large flexible hybrid systems

The mathematical model of a real flexible elastic system with distributed and discrete parameters is considered. It is a partial differential equation with non-classical boundary conditions. Complexity of the boundary conditions makes it impossible to find exact analytical solutions. To address the problem, we use the asymptotical method of small parameters together with the numerical method of normal fundamental systems of solutions. These methods allow us to investigate vibrations, and a technique for determination of complex eigenvalues of the considered boundary value problem is developed. The conditions, at which vibration processes of different characteristics take place, are defined. The dependence of the vibration frequencies on the physical parameters of the hybrid system is studied. We show that introduction of different feedbacks into the system allows one to control the frequency spectrum, in which excitation of vibrations is possible.     

求解对流扩散方程的Haar小波方法

本文用Haar小波求解对流扩散方程，将满足初始和边界条件的常系数偏微分方程简化为较简单的代数方程组进行求解．实例说明了这种方法具有收敛速度快和计算容易的特点，同时又避免了用Daubechies小波求解微分方程需要计算相关系数的麻烦．本文所使用的方法可以求解一般的微（积）分方程．     

Numerical simulation of reaction‐diffusion telegraph systems in unbounded domains

We present a boundary element method for computing numerical solutions of the reaction‐diffusion telegraph equation in unbounded domains. This technique does not need artificial boundary conditions at the computational domain and uses a new algorithm to compute the Fourier transform, the convolution theorem, and the fact that the exact solution of the telegraph equation can be written as an integral transform in terms of the fundamental solution. We use the logistic growth model to find how the population of an organism evolves according to its growth rate. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 326–335, 2015     

Solution of Poisson's equation with global,local and nonlocal boundary conditions

Boundary value problems (BVPs) for partial differential equations are common in mathematical physics. The differential equation is often considered in simple and symmetric regions, such as a circle, cube, cylinder, etc., with global and separable boundary conditions. In this paper and other works of the authors, a general method is used for the investigation of BVPs which is more powerful than existing methods, so that BVPs investigated by the method can be considered in anti-symmetric and arbitrary regions surrounded by smooth curves and surfaces. Moreover boundary conditions can be local, non-local and global. The BVP is expanded in a convex and bounded region D in a plane. First, by generalized solution of the adjoint of the Poisson equation, the necessary boundary conditions are obtained. The BVP is then reduced to the second kind of Fredholm integral equation with regularized singularities.     

万因斯坦-钱伟长法的应用——简支与固支混合边界矩形板基频上下限

本文将边界条件放松法[1][2]应用于简支与固支混合边界矩形板,求得此类板的基频下限值.文中还设计了一种能满足位移边界条件的多项式作为振型试函数,从而用里兹法得到相应的上限值.具体算例得到了满意的结果.最后本文指出,通常用力法迭加法得到的此类板的所谓精确解,若考虑实际计算的级数截断误差,本质上是一种下限解.     

A modification of the invariant imbedding method for a singular boundary value problem

We consider a dynamically-consistent analytical model of a 3D topographic vortex. The model is governed by equations derived from the classical problem of the axisymmetric Taylor–Couette flow. Using linear expansions, these equations can be reduced to a differential sixth-order equation with variable coefficients. For this differential equation, we formulate a boundary value problem, which has a number of issues for numerical solving. To avoid these issues and find the eigenvalues and eigenfunctions of the boundary value problem, we suggest a modification of the invariant imbedding method (the Riccati equation method). In this paper, we show that such a modification is necessary since the boundary conditions possess singular matrices, which sufficiently complicate the derivation of the Riccati equation. We suggest algebraic manipulations, which permit the initial problem to be reduced to a problem with regular boundary conditions. Also, we propose a method for obtaining a numerical solution of the matrix Riccati equation by means of recurrence relations, which allow us to obtain a matrizer converging to the required eigenfunction. The suggested method is tested by calculating the corresponding eigenvalues and eigenfunctions, and then, by constructing fluid particle trajectories on the basis of the eigenfunctions.     

Boundary element method for spatially-periodic potential problems

In this paper, we propose a boundary element method for two-dimensional potential problems with one-dimensional spatial periodicity, which have been difficult to be solved by the ordinary boundary element method. In the presented method, we reduce the potential problems with Dirichlet and Neumann boundary conditions to integral equation problems with the periodic fundamental solution of the Laplace operator and, then, obtain approximate solutions by solving linear systems given by discretizing the integral equations. Numerical examples are also included. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Singular boundary characteristics of the Hamilton–Jacobi equation

Situations exist in boundary value problems for first order partial differential equations arising in physics (the Hamilton–Jacobi equation), optimal control theory (the Bellman equation) and the theory of differential games (the Isaacs equation) when the value of the required function is not given on a part of the boundary or not at all, or it is not the limit of the (generalized) solution of the problem. Nevertheless, such conditions are required for constructing the solution (by the method of characteristics, for example). It is shown that the required boundary values can be exposed as a specific continuation of the conditions that are known in the boundary submanifolds of the given part of the boundary. This extension of the conditions is accomplished using the characteristic curves starting in a known submanifold of the boundary and running along the boundary. The characteristics are a generalization of the classical characteristics associated with a partial differential equation. They are called singular characteristics, and the theory of these has been developed in a number of the author's papers. After obtaining these “natural” boundary conditions, the solution is constructed using the conventional method of integrating the equations of the classical characteristics. Conditions of the Dirichlet and Neumann type are considered. The technique is illustrated using a numerical example from the theory of differential games containing a number of parameters.     

任意厚度具有自由边叠层板的精确解析解

自由边问题一直是三维弹性力学中的难题,通常很难满足自由边上一个正应力和两个剪应力都等于0．基于三维弹性力学基本方程和状态空间方法,引入自由边界位移函数并考虑全部弹性常数,建立了正交异性具有自由边单层和叠层板的状态方程．对状态方程中的变量以级数形式展开,通过边界条件的满足精确求解任意厚度具有自由边叠层板的位移和应力,此解满足层间应力和位移的连续条件．算例计算表明,采用引入的位移函数形式,简化了计算过程并且采用较少的级数项可以获得收敛解．与有限元方法计算结果进行了对比,可以得到较高精度的数值结果．其解可以作为其它数值方法和半解析方法的参考解．     

Lax pairs and a new spectral method for linear and integrable nonlinear PDEs

A new transform method for solving initial-boundary value problems for linear and integrable nonlinear PDEs in two independent variables has been recently introduced in [1]. For linear PDEs this method involves: (a) formulating the given PDE as the compatibility condition of two linear equations which, by analogy with the nonlinear theory, we call a Lax pair; (b) formulating a classical mathematical problem, the so-called Riemann-Hilbert problem, by performing a simultaneous spectral analysis of both equations defining the Lax pair; (c) deriving certain global relations satisfied by the boundary values of the solution of the given PDE. Here this method is used to solve certain problems for the heat equation, the linearized Korteweg-deVries equation and the Laplace equation. Some of these problems illustrate that the new method can be effectively used for problems with complicated boundary conditions such as changing type as well as nonseparable boundary conditions. It is shown that for simple boundary conditions the global relations (c) can be analyzed using only algebraic manipulations, while for complicated boundary conditions, one needs to solve an additional Riemann-Hilbert problem. The relationship of this problem with the classical Wiener-Hopf technique is pointed out. The extension of the above results to integrable nonlinear equations is also discussed. In particular, the Korteweg-deVries equation in the quarter plane is linearized.