数学物理中的总体和局部场GL方法

为了求解物理化学生物材料和金融中的微分方程,提出了一种总体(Global)和局部(Local)场方法.微分方程的求解区域可以是有限域,无限域,或具曲面边界的部分无限域.其无限域包括有限有界不均匀介质区域.其不均匀介质区域被分划为若干子区域之和.在这含非均匀介质的无限区域,将微分方程的解显式地表示为在若干非均匀介质子区域上和局部子曲面的积分的递归和.把正反算的非线性关系递归地显式化.在无限均匀区域,微分方程的解析解被称为初始总体场.微分方程解的总体场相继地被各个非均匀介质子区域的局部散射场所修正.这种修正过程是一个子域接着另个子域逐步相继地进行的.一旦所有非均匀介质子区域被散射扫描和有限步更新过程全部完成后,微分方程的解就获得了.称其为总体和局部场的方法,简称为GL方法.GL方法完全地不同于有限元及有限差方法,GL方法直接地逐子域地组装逆矩阵而获得解.GL方法无需求解大型矩阵方程,它克服了有限元大型矩阵解的困难.用有限元及有限差方法求解无限域上的微分方程时,人为边界及其上的吸收边界条件是必需的和困难的,人为边界上的吸收边界条件的不精确的反射会降低解的精确度和毁坏反算过程.GL方法又克服了有限元和有限差方法的人为边界的困难.GL方法既不需要任何人为边界又不需要任何吸收边界条件就可以子域接子域逐步精确地求解无限域上的微分方程.有限元和有限差方法都仅仅是数值的方法,GL方法将解析解和数值方法相容地结合起来.提出和证明了三角的格林函数积分方程公式.证明了当子域的直经趋于零时,波动方程的GL方法的数值解收敛于精确解.GL方法解波动方程的误差估计也获得了.求解椭圆型,抛物线型,双曲线型方程的GL模拟计算结果显示出我们的GL方法具有准确,快速,稳定的许多优点.GL方法可以是有网,无网和半网算法.GL方法可广泛应用在三维电磁场,三维弹塑性力学场,地震波场,声波场,流场,量子场等方面.上述三维电磁场等应用领域的GL方法的软件已经由作者研制和发展了。     

弹性力学平面问题的等价的边界积分方程

严格地推导出与弹性力学平面问题的微分方程边值问题等价的间接和直接未知量边界积分方程,用实例指出某些习用的边界积分方程有时不必要或不充分,并进行了数值比较.     

边界元方法的抽象误差估计及其应用

§1.引言 边界元方法是近二十年来发展的一种求解偏微分方程的数值方法,其基本思想是:先利用Green公式或位势将区域上的偏微分方程转化成边界上的积分方程,此时偏微分方程的解由边界积分方程的解表出;然后数值求解边界积分方程,进而求得偏微分方程的近     

三维非局部弹性场中裂纹问题的分析方法

通过求解得到了三维非局部弹性力学对称情形的单位集中不连续位移基本解·基于该基本解和三维局部（经典）弹性力学的不连续位移边界积分方程———边界元方法·提出了三维非局部弹性力学中的平片裂纹Ⅰ型问题的通用解法,并给出了算例     

关于薄板的无网格局部边界积分方程方法中的友解

无网格局部边界积分方程方法是最近发展起来的一种新的数值方法,这种方法综合了伽辽金有限元、边界元和无单元伽辽金法的优点,是一种具有广阔应用前景的、真正的无网格方法．把无网格局部边界积分方程方法应用于求解薄板问题,给出了薄板无网格局部边界积分方程方法所需要的友解及其全部公式．     

固体弹性三维问题统一解

依据三维弹性力学问题的Kelvin解,用三维虚边界元法来建立积分方程,从而使三维实体和各类板、壳等问题的求解思想得到统一.对各类三维问题采用统一的思想建立数学模型,更有利于程序模块化,增强了程序的通用性.另外,建立积分方程时直接引用Kelvin解,而未引入任何其它假设,使该方法的解更偏于实际,且使应用范围拓宽.再者,与边界元直接法相比,该方法的优点在于无需处理奇异积分,且系数阵是对称的.最后给出部分算例,以证明方法的有效性和计算精度.     

非齐次弹性力学问题双互易边界元方法研究北大核心CSCD

基于弹性力学边界元方法理论,将边界元法与双互易法结合,采用指数型基函数对非齐次项进行插值得到双互易边界积分方程.将边界积分方程离散为代数方程组,利用已知边界条件和方程特解求解方程组,得出域内位移和边界面力.指数型基函数的形状参数是由插值点最近距离的最小值决定,采用这种形状参数变化方案,分析径向基函数(RBF)插值精度以及插值稳定性.再次将指数型基函数应用到双互易边界元法中,分析双互易边界元方法下计算精度及稳定性,验证了指数型插值函数作为双互易边界元方法的径向基函数解决弹性力学域内体力项问题的有效性.     

Black-Scholes方程的Jacobi有理谱方法

1 引言无界区域问题的有理谱方法已经得到广泛地应用．它有很多优点,特别是我们不需要添加任何人工边界以及作任何变量变换就可以直接逼近微分方程．此外,Jacobi 有理谱方法可以用来数值求解变系数的微分方程,如金融数学中的基本方程-Black-     

抛物型初边值问题的边界积分-微分方程及其边界元方法

本文提出和研究了抛物型方程Neumann初边值问题的一个新的边界归化方法。它将原始初边值问题归化成一类新的边界积分-微分方程。由此导出一种新的既保持原始问题的自伴性,又具有可积弱奇性积分核的边界变分方程和边界元方法,给出了近似解在各种范数意义下的先验误差估计。     

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

GL method for solving equations in math-physics and engineering

In this paper, we propose a GL method for solving the ordinary and the partial differential equation in mathematical physics and chemics and engineering. These equations govern the acustic, heat, electromagnetic, elastic, plastic, flow, and quantum etc. macro and micro wave field in time domain and frequency domain. The space domain of the differential equation is infinite domain which includes a finite inhomogeneous domain. The inhomogeneous domain is divided into finite sub domains. We present the solution of the differential equation as an explicit recursive sum of the integrals in the inhomogeneous sub domains. Actualy, we propose an explicit representation of the inhomogeneous parameter nonlinear inversion. The analytical solution of the equation in the infinite homogeneous domain is called as an initial global field. The global field is updated by local scattering field successively subdomaln by subdomain. Once all subdomains are scattered and the updating process is finished in all the sub domains, the solution of the equation is obtained. We call our method as Global and Local field method, in short , GL method. It is different from FEM method, the GL method directly assemble inverse matrix and gets solution. There is no big matrix equation needs to solve in the GL method. There is no needed artificial boundary and no absorption boundary condition for infinite domain in the GL method. We proved several theorems on relationships between the field solution and Green＇s function that is the theoretical base of our GL method. The numerical discretization of the GL method is presented. We proved that the numerical solution of the GL method convergence to the exact solution when the size of the sub domain is going to zero. The error estimation of the GL method for solving wave equation is presented. The simulations show that the GL method is accurate, fast, and stable for solving elliptic, parabolic, and hyperbolic equations. The GL method has advantages and wide applications in the 3D electromagnetic （EM）     

On a variational problem with unknown boundaries and the determination of optimal shapes of elastic bodies : PMM vol. 39, n 6, 1975, pp. 1082–1092

The optimization problem is considered for a partial differential equation of elliptic type. The boundary of the domain in which the equation is given emerges as the control function and is to be determined from the condition of the extremum of the integral of the solution of the boundary value problem. Seeking the extremals is reduced to solving a va national problem without differential constraints. Necessary conditions for optimality are obtained, and shapes of elastic bars possessing the maximum stiffness under torsion are found with their aid.     

An improved boundary element formulation for the motion of spherical bubbles

We suggest a modified boundary element method for modeling the potential flow caused by the motion of many spherical bubbles. The problem for the fluid velocity potential is reduced to a Fredholm integral equation of the second kind. The kernel of that integral equation has no singularity; as a result, its numerical solution does not encounter any difficulties. The matrix representation of the integral equation is diagonally dominant and is well suited to handle multiple bubble system.     

The Robin problem for bending of elastic plates

The solution of the Robin problem in a finite domain for the system of equations modeling the bending of elastic plates with transverse shear deformation is approximated by means of a generalized Fourier series method closely connected to the structure of the boundary integral equation treatment of the problem. The theory is exemplified by numerical computation that shows a high degree of accuracy and efficiency.     

Universal boundary value problem for equations of mathematical physics

A new statement of a boundary value problem for partial differential equations is discussed. An arbitrary solution to a linear elliptic, hyperbolic, or parabolic second-order differential equation is considered in a given domain of Euclidean space without any constraints imposed on the boundary values of the solution or its derivatives. The following question is studied: What conditions should hold for the boundary values of a function and its normal derivative if this function is a solution to the linear differential equation under consideration? A linear integral equation is defined for the boundary values of a solution and its normal derivative; this equation is called a universal boundary value equation. A universal boundary value problem is a linear differential equation together with a universal boundary value equation. In this paper, the universal boundary value problem is studied for equations of mathematical physics such as the Laplace equation, wave equation, and heat equation. Applications of the analysis of the universal boundary value problem to problems of cosmology and quantum mechanics are pointed out.     

An iterative penalty method for a nonlinear problem of equilibrium of a Timoshenko-type plate with a crack

A variational problem of equilibrium of an elastic Timoshenko-type plate in a domain with a slit is considered. A nonpenetration boundary condition in the form of an inequality is specified on the edges of the slit. A penalized equation and an iterative linear equation in integral and differential forms are constructed. Some results on solution convergence and an error estimate are obtained.     

The link between integral equations and higher order ODEs

We study a class of partial integro-differential equations defined on a spatially extended domain that arise in the modeling of neuronal networks. In previous studies the Fourier transform was used to derive associated fourth order ordinary differential equations when the solution of the time independent integral equation is a homoclinic orbit. This gives rise to the question of whether solutions of the ODE whose Fourier transform is not well defined are also solutions of the time independent integral equation. We address this question and show that any solution of the ODE that satisfies a fairly relaxed growth condition is also a solution of the integral equation. Applications to two specific examples are given.     

An effective boundary element method for inhomogeneous partial differential equations

A method for removing the domain or volume integral arising in boundary integral formulations for linear inhomogeneous partial differential equations is presented. The technique removes the integral by considering a particular solution to the homogeneous partial differential equation which approximates the inhomogeneity in terms of radial basis functions. The remainder of the solution will then satisfy a homogeneous partial differential equation and hence lead to an integral equation with only boundary contributions. Some results for the inhomogeneous Poisson equation and for linear elastostatics with known body forces are presented.     

应用边界元法模拟纤维增强复合材料平面弹性问题

将含有随机分布多种夹杂相复合材料的二维弹性力学问题归结为复连通区域的边界积分方程,进而转化成矩阵方程进行求解和分析．根据同类夹杂相外在边界上的面力与位移之间关系矩阵完全相同的特点,使得最后的矩阵方程阶数得到大规模减少,这正是此处提出改进的边界元方法的主要思路．数值算例表明,对于此类问题,与常规的边界元分域解法相比更加有效．以该方法为基础,可以详细给出纤维增强复合材料二维条件下的宏观等效力学性质．