第二类混合变分不等式的MRM-边界元分析

本文以弹性力学中的摩擦问题为背景,采用多重互易方法(MRM方法),边界元方法,将摩擦问题中的第二类混合变分不等式化解为MRM-边界混合变分不等式,给出了MRM-边界混合变分不等式解的存在唯—性,通过引入变换将原MRM-边界混合变分不等式化解为标准的凸极值问题,采用正则化方法处理后,给出了MRM-边界混合变分不等式的迭代分解方法。文末给出了数值算例。     

一个扩散问题的自然边界元法与有限元法组合

本文讨论由Helmholtz方程描述的扩散问题的自然边界元法与有限元法的组合．取一个圆作为公共边界，用Fourier展开建立边界积分方程，将无界区域上的问题化为有界区域上的非局部边值问题．在变分方程中公共边界上的未知量只包含函数本身而不包含其法向导数，从而减少了未知数的数目，并且边界元剐度矩阵只有极少量不同的元素，有利于数值计算．这种组台方法优越于建立在直接边界元法基础上的组合方法．文中证明了变分解的唯一性，数值解的收敛性和误差估计．最后讨论了数值技术并给出一个算倒．     

椭圆外区域上Helmholtz问题的耦合法

以椭圆外区域上Helmholtz方程为例,研究一种带有椭圆人工边界的自然边界元与有限元耦合法,给出了耦合变分问题的适定性及误差分析并给出数值例子.理论分析及数值结果表明,用方法求解椭圆外问题是十分有效的.为求解具有长条型内边界外Helmholtz问题提供了一种很好的数值方法.     

偏微分方程Δ2u-sΔu+k2u=0边值问题的MRM边界变分方程

导出边值问题Δ2u-sΔu+k2u=o;x∈Ω∪Ω'(R2;u|г=uo;аu/аn|г=go的定解问题,MRM边界变分方程,全平面解的表达式.从中可以看出,MRM边界变分方程中只包含弱奇异积分核,并且自动消除了原第一、二MRM边界积分方程中出现的强奇异积分核.问题解的表达式后并不加任何多项式,因而也不需要引入Lagrange乘子求解该项,这给边界元数值求解过程带来极大的方便.数值分析结果表明该方法具有明显优势.     

边界元方法的数学分析

本文介绍边界元方法与其它数值计算方法的关系.用拟微分算子给边界元法中所遇到的各种类型的边界积分方程以统一的数学描述.由拟微分算子的强椭圆性,得出边界积分方程的可解性.本文还介绍边界元空间的建立和几种求解边界积分方程的离散化方法.对于边界积分方程的变分形式,给出边界元近似解的收敛性和渐近误差估计.     

有势场逆问题的边界元法

本文给出了位势方程逆问题的一种最小二乘边界元解法。控制方程为Ｌａｐｌａｃｅ方程，但一部分边界上未给出任何边值，而只在某些内点上给出了势函值。这一问题在数学上属不适定问题，但在一定条件下存在唯一解。本文同时给出了一种估计解的可靠性的方法。数值试验表明，这类逆问题采用边界元法是非常有效的。     

第二类四阶变分不等式解的存在唯一性

本文考虑弹性力学平板理论中的单侧稳定问题,讨论了重调和方程边值问题及第二类四阶变分不等式,并证明它们的等价性,从而可以把高阶偏微分方程转化为相应的变分不等式加以解决,同时也为求解重调和方程提供了更多的方法和理论依据。并且最后给出了该类变分不等式解的存在唯一性的证明。     

速多极边界元法解的存在唯一性

本文对求解3维弹性摩擦接触问题的快速多极边界元法(FM- BEM)在数学理论上作了深入探讨.首先,利用向量和子空间理论找出快速优化广义极小残余算法(GMRES(m) )求解边界元方程组所满足的代数条件,使对工程用FM- BEM解的研究转化为对代数问题的讨论,然后,分三步证明了FM- BEM解的存在唯一性,为FM- BEM求解弹性摩擦接触工程问题提供强有力的数学支撑.     

一种改进边界元法在弹性扭转问题中的应用

本文用一种改进边界元法分析与计算了椭圆截面等直杆的扭转问题.并与边界元法的解进行比较,其结果极为符合.然而,改进边界元法较边界元法所需要的数据量少得多,计算时间也将大大减少了.因此,本文方法对求解Poisson方程问题是一种经济而行之有效的数值计算方法.     

快速多极边界元法解的存在唯一性

本对求解3维弹性摩擦接触问题的快速多极边界元法(FM—BEM)在数学理论上作了深入探讨．首先，利用向量和子空间理论找出快速优化广义极小残余算法(GMRES(m))求解边界元方程组所满足的代数条件．使对工程用FM—BEM解的研究转化为对代数问题的讨论，然后．分三步证明了FM-BEM解的存在唯一性，为FM-BEM求解弹性摩擦接触工程问题提供强有力的数学支撑．     

奇异杂交边界点方法求解泊松方程

As a boundary-type meshless method,the singular hybrid boundary node method(SHBNM)is based on the modified variational principle and the moving least square(MLS)approximation,so it has the advantages of both boundary element method(BEM)and meshless method.In this paper,the dual reciprocity method(DRM)is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution.The general solution is achieved by means of SHBNM,and the particular solution is approximated by using the radial basis function(RBF).Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain.The postprocess is very simple.Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.     

Decoupled Mixed Element Methods for Fourth Order Elliptic Optimal Control Problems with Control Constraints

In this paper, we study the finite element methods for distributed optimal control problems governed by the biharmonic operator. Motivated from reducing the regularity of solution space, we use the decoupled mixed element method which was used to approximate the solution of biharmonic equation to solve the fourth order optimal control problems. Two finite element schemes, i.e., Lagrange conforming element combined with full control discretization and the nonconforming Crouzeix-Raviart element combined with variational control discretization, are used to discretize the decoupled optimal control system. The corresponding a priori error estimates are derived under appropriate norms which are then verified by extensive numerical experiments.     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

Dual reciprocity hybrid radial boundary node method for the analysis of Kirchhoff plates

A meshless method of dual reciprocity hybrid radial boundary node method (DHRBNM) for the analysis of arbitrary Kirchhoff plates is presented, which combines the advantageous properties of meshless method, radial point interpolation method (RPIM) and BEM. The solution in present method comprises two parts, i.e., the complementary solution and the particular solution. The complementary solution is solved by hybrid radial boundary node method (HRBNM), in which a three-field interpolation scheme is employed, and the boundary variables are approximated by RPIM, which is applied instead of moving least square (MLS) and obtains the Kronecker’s delta property where the traditional HBNM does not satisfy. The internal variables are interpolated by two groups of symmetric fundamental solutions. Based on those, a hybrid displacement variational principle for Kirchhoff plates is developed, and a meshless method of HRBNM for solving biharmonic problems is obtained, by which the complementary solution can be solved.     

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

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

The solution of the non-homogeneous Helmholtz equation by means of the boundary element method

One of the most important advantages of the Boundary Element Method (BEM) is that no internal discretization of the domain is required. This advantage, however, is generally lost when source terms are present in the governing differential equation. It is shown here that for the non-homogeneous Helmholtz equation with a harmonic source term, it is possible to transform the volume integral into a surface integral thus retaining this feature. The transformation is achieved using the Green formula. The technique is applied to solve numerically a test problem with known simple analytical solution.     

On the Numerical Solution of Some Eikonal Equations: An Elliptic Solver Approach

The steady Eikonal equation is a prototypical first-order fully nonlinear equation. A numerical method based on elliptic solvers is presented here to solve two different kinds of steady Eikonal equations and compute solutions, which are maximal and minimal in the variational sense. The approach in this paper relies on a variational argument involving penalty, a biharmonic regularization, and an operator-splitting-based time-discretization scheme for the solution of an associated initial-value problem. This approach allows the decoupling of the nonlinearities and differential operators. Numerical experiments are performed to validate this approach and investigate its convergence properties from a numerical viewpoint.     

Multigrid technique for biharmonic interpolation with application to dual and multiple reciprocity method

A biharmonic-type interpolation method is presented to solve 2D and 3D scattered data interpolation problems. Unlike the methods based on radial basis functions, which produce a large linear system of equations with fully populated and often non-selfadjoint and ill-conditioned matrix, the presented method converts the interpolation problem to the solution of the biharmonic equation supplied with some non-usual boundary conditions at the interpolation points. To solve the biharmonic equation, fast multigrid techniques can be applied which are based on a non-uniform, non-equidistant but Cartesian grid generated by the quadtree/octtree algorithm. The biharmonic interpolation technique is applied to the multiple and dual reciprocity method of the BEM to convert domain integrals to the boundary. This makes it possible to significantly reduce the computational cost of the evaluation of the appearing domain integrals as well as the memory requirement of the procedure. The resulting method can be considered as a special grid-free technique, since it requires no domain discretisation. This revised version was published online in June 2006 with corrections to the Cover Date.     

Boundedness of the Hessian of a Biharmonic Function in a Convex Domain

We consider the Dirichlet problem for the biharmonic equation on an arbitrary convex domain and prove that the second derivatives of the variational solution are bounded in all dimensions.