分数阶微分方程边值问题解的存在性

应用Gteen函数将分数阶微分方程边值问题可转化为等价的积分方程．近来此方法被应用于讨论非线性分数阶微分方程边值问题解的存在性．讨论非线性分数阶微分方程边值问题,应用Green函数,将其转化为等价的积分方程,并设非线性项满足Caratheodory条件,利用非紧性测度的性质和M6nch’s不动点定理证明解的存在性．     

Some relationships between Green's functions and invariant imbedding

The method of invariant imbedding has been used to resolve the solution of linear two-point boundary-value problems into contributions associated with the homogeneous equation with homogeneous boundary conditions, with inhomogeneous boundary conditions, and with an inhomogeneous source term in the equation. The relationship between the Green's function and the invariant imbedding equations is described, and it is shown that the Green's function can be determined from an initial-value problem. Several numerical examples are given which illustrate the efficacy of the initial-value algorithm.This work was supported by the US Atomic Energy Commission.     

重调和椭圆边值问题的正则积分方程

我们熟知,利用位势理论或由Green公式及基本解出发区域内调和及重调和边值问题可归化为边界上的积分方程。近年来冯康又提出一种更自然而直接的归化,即从Green公式及Green函数出发将微分方程边值问题化为边界上的含有广义函数意义下发散积分有限部分的奇异积分方程,这种归化在各种边界归化中占有特殊地位,被称为正则边界归化,本文将这一理论应用于重调和椭圆边值问题,研究了其正则归化的性质,并通过利用Green函数、Fourier分析及复变函数论方法等不同途径求出了在上半平面、单位圆内部、单位圆外部三种区域的Poisson积分公式及正则积分方程,其离散化可用于实际计算。 本文是在导师冯康教授指导下完成的,作者谨在此对他表示衷心的感谢。     

一维热传导方程移动边界识别的计算方法

构造了一种正则化的积分方程方法来由Cauchy数据确定一维热传导方程的移动边界．在将区域延拓至规则区域后,通过Fourier方法将问题转化为一个第一类Volterra积分方程．然后分别用Lavrentiev正则化方法以及Tikhonov正则化方法将不稳定的第一类Volterra积分方程转化为适定的第二类积分方程,并分别将积分方程转化为常微分方程组,并用Runge—Kutta方法数值求解,以及直接离散来求解．最后通过自由边界上的条件得到数值的移动边界．通过一些数值试验表明此方法是有效可行的,并且给出的方法无需迭代,数值计算较简单．     

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

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

具两组高阶基本解系列的MRM边界积分方程

以双参数地基上板弯曲问题为模型,利用两组高阶基本解进行交替多重替换,得到边界积分方程,并证明该方程与边值问题常规的边界积分方程本质是一致的,具更便于计算。     

双空间指示函数方法在三维柱面对称波导中电磁波的反散射问题的推广

给出了双空间指示函数方法在三维柱面对称波导中电磁波的反散射问题的推广.基于这个观察:当Green函数的点源在障碍物内部时,那么远域数据的赋权积分可以很好地近似估计Green函数,但是当Green函数的点源在障碍物外部时,那么远域数据的赋权积分就不能很好地近似估计Green函数.建立一个积分方程:它的右边是点源在所重构区域的Green函数,那么我们可以知道这个积分方程的解的范数在未知障碍物的内部有极大值,而这些取得极大值的点所围成的区域恰好就是所重构的障碍物区域.方法最显著的优势在于它不依赖于未知障碍物的边界条件.     

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

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

Convergence of the Neumann Series in BEM for the Neumann Problem of the Stokes System

A weak solution of the Neumann problem for the Stokes system in Sobolev space is studied in a bounded Lipschitz domain with connected boundary. A solution is looked for in the form of a hydrodynamical single layer potential. It leads to an integral equation on the boundary of the domain. Necessary and sufficient conditions for the solvability of the problem are given. Moreover, it is shown that we can obtain a solution of this integral equation using the successive approximation method. Then the consequences for the direct boundary integral equation method are treated. A solution of the Neumann problem for the Stokes system is the sum of the hydrodynamical single layer potential corresponding to the boundary condition and the hydrodynamical double layer potential corresponding to the trace of the velocity part of the solution. Using boundary behavior of potentials we get an integral equation on the boundary of the domain where the trace of the velocity part of the solution is unknown. It is shown that we can obtain a solution of this integral equation using the successive approximation method.     

两平行圆板间径向层流进口段效应分析

本文首先将B.B.方法[1]推广至二平行圆板间的径向扩散流动,由边界层运动方程式同时推导出动量积分方程式和能量积分方程式,而后再用Picard逐次逼近法[2]解能量积分方程式,求得进口段通道长随边界层厚度而改变的二级近似显函数表达式.从而为进口段效应诸系数的直接解析分析提供了可能.特别是当圆板外径小于进口段长度时,更加突出地表现了本方法的优越性.由于采用了能量积分方程式,则压力损失系数的各项才得以从理论上独立地推导出来.本文所提供的压力损失系数计算值,在进口修正雷诺数Re<100时,和文献[3]比较与实验值更为接近.因此在该范围内本文的结果既可靠又简便.     

Null space distributions—A new approach to finite convolution equations with a Hankel kernel

In some earlier publications it has been shown that the solutions of the boundary integral equations for some mixed boundary value problems for the Helmholtz equation permit integral representations in terms of solutions of associated complicated singular algebraic ordinary differential equations. The solutions of these differential equations, however, are required to be known on some infinite interval on the real line, which is unsatisfactory from a practical point of view. In this paper, for the example of one specific boundary integral equation, the relevant solutions of the associated differential equation are expressed by integrals which contain only one unknown generalized function, the support of this generalized function is no longer unbounded but a compact subset of the real line. This generalized function is a distributional solution of the homogeneous boundary integral equation. By this null space distribution the boundary integral equation can be solved for arbitrary right-hand sides, this solution method can be considered of being analogous to the method of variation of parameters in the theory of ordinary differential equations. The nature of the singularities of the null space distribution is worked out and it is shown that the null space distribution itself can be expressed by solutions of the associated ordinary differential equation.     

On the asymptotic convergence of a collocation method for some boundary integral equations of the first kind

In this paper we present a certain collocation method for the numerical solution of a class of boundary integral equations of the first kind with logarithmic kernel as principle part. The transformation of the boundary value problem into boundary singular integral equation of the first kind via single-layer potential is discussed. A discretization and error representation for the numerical solution of boundary integral equations has been given. Quadrature formulae have been proposed and the error arising due to the quadrature formulae used has been estimated. The convergence of the solution with respect to the proposed numerical algorithm is shown and finally some numerical results have been presented.     

Ursell's approach to obtaining an a priori estimate for the solution of the Neumann problem for the Helmholtz equation

The geometric optics asymptotic expression for Green's function for the Helmholtz equation in the exterior of a bounded convex region in R³ is rigorously justified (it is assumed that the Neumann boundary condition is satisfied). An integral equation constructed by a three-dimensional analogue of Ursell's method from Green's function for a paraboloid of revolution is the basis of all considerations. Analysis of this integral equation also makes it possible to prove the exponential decay of Green's function for the Helmholtz equation in the zone of deep shadow.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 42, pp. 85–154, 1974.     

The exterior dirichlet problem for the Laplace equation

Using a standard application of Green's theorem, the exterior Dirichlet problem for the Laplace equation in three dimensions is reduced to a pair of integral equations. One integral equation is of the second kind and the other is of the first kind. It is known that the integral equation of the second kind is not uniquely solvable, however, it has been demonstrated that the pair of integral equations has a unique solution. The present approach is based on the observation that the known function appearing in the integral equation of the second kind lies in a certain Banach space E which is a proper subspace of the Banach space of continuous complex-valued functions equipped with the maximum norm. Furthermore, it can be shown that the related integral operator when restricted to E has spectral radius less than unity. Consequently, a particular solution to the integral equation of the second kind can be obtained by the method of successive approximations and the unique solution to the problem is then obtained by using the integral equation of the first kind. Comparisons are made between the present algorithm and other known constructive methods. Finally, an example is supplied to illustrate the method of this paper.     

The Behaviour of Elastic Fields and Boundary Integral Mellin Techniques Near Conical Points

For the computation of the local singular behaviour of an homogeneous anisotropic clastic field near the three-dimensional vertex subjected to displacement boundary conditions, one can use a boundary integral equation of the first kind whose unkown is the boundary stress. Mellin transformation yields a one - dimensional integral equation on the intersection curve 7 of the cone with the unit sphere. The Mellin transformed operator defines the singular exponents and Jordan chains, which provide via inverse Mellin transformation a local expansion of the solution near the vertex. Based on Kondratiev's technique which yields a holomorphic operator pencil of elliptic boundary value problems on the cross - sectional interior and exterior intersection of the unit sphere with the conical interior and exterior original cones, respectively, and using results by Maz'ya and Kozlov, it can be shown how the Jordan chains of the one-dimensional boundary integral equation are related to the corresponding Jordan chains of the operator pencil and their jumps across γ. This allows a new and detailed analysis of the asymptotic behaviour of the boundary integral equation solutions near the vertex of the cone. In particular, the integral equation method developed by Schmitz, Volk and Wendland for the special case of the elastic Dirichlet problem in isotropic homogeneous materials could be completed and generalized to the anisotropic case.     

An extrapolation method for a class of boundary integral equations

Boundary value problems of the third kind are converted into boundary integral equations of the second kind with periodic logarithmic kernels by using Green's formulas. For solving the induced boundary integral equations, a Nyström scheme and its extrapolation method are derived for periodic Fredholm integral equations of the second kind with logarithmic singularity. Asymptotic expansions for the approximate solutions obtained by the Nyström scheme are developed to analyze the extrapolation method. Some computational aspects of the methods are considered, and two numerical examples are given to illustrate the acceleration of convergence.

     

积分算子方程的连续型配置方法

1．引言由于对积分算子方程来说，配置法比Galerkin法具计算量小的优点（少算一重积分），故配置法更受人们重视．但已有的文献几乎都是将配置空间取作非连续的分片多项式样条空间，以得到某种超收敛结果（如［1，2］）．这种方法存在下列不足：（a）光滑核Volterra积分方程与光滑核Fredholm积分方程具完全不同的收敛性质［1］，且需用不同的方法获得其加速收敛结果（比较［31与［4］），尽管Volterra积分方程在理论上被看作是Fredholm积分方程的特殊情形；（b）光滑核Volterra积分方程的配置解不具任何超收敛性，其迭代配置解也只在结点…     

Wavelet Collocation Methods for a First Kind Boundary Integral Equation in Acoustic Scattering

In this paper we consider a wavelet algorithm for the piecewise constant collocation method applied to the boundary element solution of a first kind integral equation arising in acoustic scattering. The conventional stiffness matrix is transformed into the corresponding matrix with respect to wavelet bases, and it is approximated by a compressed matrix. Finally, the stiffness matrix is multiplied by diagonal preconditioners such that the resulting matrix of the system of linear equations is well conditioned and sparse. Using this matrix, the boundary integral equation can be solved effectively.     

A finite element method for the solution of a potential theory integral equation

This paper discusses a finite element approximation for an integral equation of the second kind deduced from a potential theory boundary value problem in two variables. The equation is shown to admit a unique solution, to be variational and coercive in the Hilbert space of functions σ ε H^1/2(Γ), f_rd γ = 0. The Galerkin method with finite elements as trial functions is shown to lead to an optimal rate of convergence.     

Generalized Dirichlet‐to‐Neumann Map in Time‐Dependent Domains

We study the heat, linear Schrödinger (LS), and linear KdV equations in the domain l(t) < x < ∞ , 0 < t < T , with prescribed initial and boundary conditions and with l(t) a given differentiable function. For the first two equations, we show that the unknown Neumann or Dirichlet boundary value can be computed as the solution of a linear Volterra integral equation with an explicit weakly singular kernel. This integral equation can be derived from the formal Fourier integral representation of the solution. For the linear KdV equation we show that the two unknown boundary values can be computed as the solution of a system of linear Volterra integral equations with explicit weakly singular kernels. The derivation in this case makes crucial use of analyticity and certain invariance properties in the complex spectral plane. The above Volterra equations are shown to admit a unique solution.