直接边界元法中边界积分的解析处理

确立了平面位势和弹性力学问题的边界元直接法中边界积分的解析计算框架系统,从而避免了传统的主似求积分,数值算例表胆它具有较高的精度和效率,特别是在边界量和边界附近区域内点物理量的计算可获得较高的精度。     

平面弹性力学问题边界元直接法中边界积分的解析处理

确立平面位势和弹性力学问题的边界元直接法中边界积分的解析计算框架系统,从而避免了传统的高斯近似求积分,数值算例表明它具有较高的精度和效率,特别是在边界量和边界附近区域内点物理量的计算可获得较高的精度。     

边界元近奇异积分计算的迭代Sinh-Sigmoidal组合式变换法

精确有效地消除积分的近奇异性是三维边界元法在工程应用中的首要问题.当源点与三角形积分单元间的距离无限趋近于零时,会出现近奇异积分问题,积分单元的形状和投影点的位置都是影响近奇异积分计算精度的重要因素.现有的非线性变换法大多只关注径向上积分的近奇异性,而忽略了角度方向和积分单元形状的影响,在投影点接近三角形积分单元边界的情况下,无法获得令人满意的计算精度,并且对子三角形积分单元的形状非常敏感.因此提出了一种改进的基于自适应分块技术和不同坐标变换的迭代sinh sigmoidal组合式变换法,分别消除径向和角度方向积分的近奇异性,在确保计算精度的同时,大大减小了计算规模.数值算例验证了该方法的有效性.     

薄体结构温度场的高阶边界元分析

分析了二维问题边界元法3节点二次单元的几何特征,区分和定义了源点相对高阶单元的Ⅰ型和Ⅱ型接近度.针对二维位势问题高阶边界元中奇异积分核,构造出具有相同Ⅱ型几乎奇异性的近似核函数,在几乎奇异积分单元上分离出积分核中主导的奇异函数部分.原积分核扣除其近似核函数后消除几乎奇异性,成为正则积分核函数,并采用常规Gauss数值方法计算该正则积分;对奇异核函数的积分推导出解析公式,从而建立了一种新的边界元法高阶单元几乎奇异积分半解析算法.应用该算法计算了二维薄体结构温度场算例,计算结果表明高阶单元半解析算法能充分发挥边界元法优势,显著提高计算精度.     

具有代数精度的降维展开式与边界型求积公式

在著作[1]中曾研究了高维积分的边界型求积公式的构造法.本文主要研究具有代数精度的边界型求积公式的构造问题,针对较为一般类型的积分区域,我们给出了具有指定代数精度的边界型求积公式的一股构造原则,其中应用了具有较高代数精度的降维展开式,並对降维展开式的余项给出了估计.     

机械求积法及其外推算法求解Helmholtz非线性边界积分方程

当Helmholtz微分方程转化为非线性边界积分方程后,可以利用机械求积法求得近似解,此方法具有较高的收敛精度阶O(h3)和较低的计算复杂度.构造机械求积法时,一个非线性方程系统通过离散非线性积分方程得到.此外,每个矩阵元素的值都不需要计算任何奇异积分.根据渐近紧理论和Stepleman定理,整个系统的稳定性和收敛性得到了证明.利用h3-Richardson外推算法,收敛精度阶可以提高到O(h5).为了求解非线性方程组,利用Ostrowski不动点定理研究了Newton的解的收敛性.几个算例从数值上说明了本算法的有效性.     

位势平面问题的新的规则化边界积分方程

广泛实践集中在直接变量边界积分方程的规则化研究,其本质是利用简单解消除边界积分的奇异性．然而,至今关于平面位势问题的第一类边界积分方程的规则化研究尚未涉足．致力于间接变量边界积分方程的规则化方法研究,基于一种新的思想和观点,确立平面位势问题的间接变量规则边界积分方程,它不包含CPV强奇异积分和HFP超奇异积分．数值算例表明现在的方法可取得很好的精度和效率,特别是边界量的计算．     

角度-距离复合变换法消除边界积分方程近奇异性

针对薄型结构边界单元分析过程中出现的近奇异积分问题,研究了采用一种角度变换和距离变换相结合的方法,节省了计算量,提高了计算精度.研究发现,当积分单元上与配置点距离最近的点落在积分单元的边沿或者顶点附近时,经过基础变换后的积分在两个方向都表现出奇异性,因此,对两个方向同时使用一维非线性变换能够切实消除近奇异性.数值算例验证了复合变换对近奇异积分计算精度的提升效果.     

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

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

简支梯形底扁球壳自由振动问题的准Green函数方法

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

求解对流扩散方程的一致四阶紧致格式

目前,许多高精度差分格式,由于未成功地构造与其精度匹配的稳定的边界格式,不得不采用低精度的边界格式.本文针对对流扩散方程证明了存在一致四阶紧致格式,它的边界点的计算格式和内点的计算格式的截断误差主项保持一致,给出了具体内点和边界格式;并分析了此半离散格式的渐近稳定性.数值结果表明该格式是四阶精度;在对流占优情况下,本文边界格式的数值结果比四阶精度的显式差分格式的的数值结果的数值振荡小,取得了不错的效果,理论结果得到了数值验证;驱动方腔数值结果显示,本文对N-S方程的离散格式具有很好的可靠性,适合对复杂流体流动的数值模拟和研究.     

Algorithm composition scheme for problems in composite domains based on the difference potential method

An algorithm composition scheme for the numerical solution of boundary value problems in composite domains is proposed and illustrated using an example. The scheme requires neither difference approximations of the boundary conditions nor matching conditions on the boundary between the subdomains. The scheme is suited for multiprocessor computers.     

A local iteration scheme for nonlinear two-dimensional steady-state heat conduction: a BEM approach

An efficient algorithm is proposed to solve the steady-state nonlinear heat conduction equation using the boundary element method (BEM). Nonlinearity of the heat conduction equation arises from nonlinear boundary conditions and temperature dependence of thermal conductivity. Using Kirchhoff's transformation, the case of temperature dependence of thermal conductivity can be transformed to the nonlinear boundary conditions case. Applying the BEM technique, the resulting matrix equation becomes nonlinear. The nonlinearity, however, only involves the boundary nodes that have nonlinearboundary conditions. The proposed local iterative scheme reduces the entire BEM matrix equation to a smaller matrix equation whose rank is the same as the number of boundary nodes with nonlinear boundary conditions. The Newton-Raphson iteration scheme is used to solve the reduced nonlinear matrix equation. The local iterative scheme is first applied to two one-dimensional problems (analytical solutions are possible) with different nonlinear boundary conditions. It is then applied to a two-region problem. Finally, the local iterative scheme is applied to two cavity problems in which radiation plays a role in the heat transfer.     

A fourth-order method of the convection-diffusion equations with Neumann boundary conditions

In this paper, we have developed a fourth-order compact finite difference scheme for solving the convection-diffusion equation with Neumann boundary conditions. Firstly, we apply the compact finite difference scheme of fourth-order to discrete spatial derivatives at the interior points. Then, we present a new compact finite difference scheme for the boundary points, which is also fourth-order accurate. Finally, we use a Padé approximation method for the resulting linear system of ordinary differential equations. The presented scheme has fifth-order accuracy in the time direction and fourth-order accuracy in the space direction. It is shown through analysis that the scheme is unconditionally stable. Numerical results show that the compact finite difference scheme gives an efficient method for solving the convection-diffusion equations with Neumann boundary conditions.     

Viscosity-splitting scheme for the Navier-Stokes equations

A viscosity-splitting scheme for the initial boundary value problems of the Navier-Stokes equations is considered. In the scheme, the Stokes equation is solved in conjunction with a nonhomogeneous boundary condition which connects the tangent flow with a no-slip flow. Convergence is proved.     

A non-iterative immersed boundary-lattice Boltzmann method with boundary condition enforced for fluid–solid flows

A non-iterative immersed boundary lattice Boltzmann method (IB-LBM) is proposed in this work for the simulation of fluid–solid flows. In the scheme, the interface is implemented by the correction of the neighboring distribution functions, similar to that of the LBM. Such treatment of the boundary is contrary to the traditional methods, where the interface is usually modeled as a generator of external force. Therefore, an advantage of the present method is to remove the efforts to evaluate the IB force and then incorporate it into the governing equation. Furthermore, an adjustment parameter is introduced to the immersed boundary scheme, which ensures the interpolated distribution functions derive the desired velocity at the boundary. Compared with the solution of a large boundary matrix and the multiple force correction that generally used in the previous studies, the present method is simpler and efficient without any iterative procedures. Those above-mentioned features make the present scheme based on the correction of the distribution function, with the enforcement of no-slip boundary condition. Simulation of flow past a fixed cylinder shows that there is no penetration of streamlines to the cylinder surface, indicating a well enforcement of the no-slip boundary condition. This scheme is further validated in the flows of a cylinder oscillating in a quiescent fluid, circular and elliptical particles settling in a channel. The results have good agreement with those data available in the literature.     

The new alternating direction implicit difference methods for solving three-dimensional parabolic equations

A new alternating direction implicit (ADI) scheme for solving three-dimensional parabolic equations with nonhomogeneous boundary conditions is presented. The scheme is also extended to high-order compact difference scheme. Both of them have the advantages of unconditional stability and being convenient to compute the boundary values of the intermediates. Besides this, the compact scheme has high-order accuracy and uses less computational time. Numerical examples are presented and the results are very satisfactory.     

Fourth Order Symmetric Finite Difference Schemes for the Acoustic Wave Equation

We present an explicit, symmetric finite difference scheme for the acoustic wave equation on a rectangle with Neumann and/or Dirichlet boundary conditions. The scheme is fourth order accurate both in time and space. It is obtained by mass lumping of a finite element scheme. The accuracy and the difference approximations at the boundary are analyzed in terms of local and global errors. AMS subject classification (2000) 65M10