二维多域弹性问题虚边界无网格伽辽金法分析

针对复杂的不同材料属性的多域组合问题（比如复合材料交界面上接触应力的计算）,虚边界无网格伽辽金法被进一步研究,提出了二维多域弹性问题虚边界无网格伽辽金法。简要介绍了多域组合思想、子域虚边界元法,详细推导了二维多域弹性问题分析的虚边界无网格伽辽金法,得到具体的离散格式,便于编程,推广研究。方程的加权系数为位移、面力、连续边界上的位移与面力关系式偏导,数值意义明确,公式具体。最后通过计算数值实例为复合材料交界面上接触应力的计算,给出了复合圆盘接触面上的法向、径向应力,分多种方案调整每个子域的虚边界半径值,所得结果与解析解、其他数值方法进行比较。结论是二维多域弹性问题虚边界无网格伽辽金法的方法计算可行、精确性与稳定性好。     

移动最小二乘法导数近似讨论

在移动最小二乘法（moving least squares method, MLS）构造无网格形函数的数值方法中,通常采用无单元伽辽金法（element-free Galerkin method, EFG）的建议,将系数向量a参与导数运算。为探讨这种导数近似算法在更一般无网格法中的适用性和合理性,针对系数向量a是否应参与运算的问题进行讨论和数值检验。结果表明：单纯从近似意义上讲,这种将系数向量代入导数运算的算法并不具有优势;从数值方法的应用意义上讲,这种导数近似算法对数值求解,特别是强式无网格法,会带来一系列潜在不稳定的问题。建议在MLS导数近似中,系数向量a不应当参与导数运算,并提出采用一种由核基函数代替普通基函数的核近似法。     

无网格伽辽金法节点分布求解精度研究

摘　要 基于移动最小二乘法的无网格伽辽金法（EFG）的计算精度的影响因素很多。为了分析无网格伽辽金法对节点布置的影响,文中通过分析经典的悬臂梁问题,定义了一个能量范数作为误差指标,用Matlab程序开发了无网格法计算程序。研究了离散节点的布置方法及节点个数对无网格法伽辽金法求解精度的影响。考察了均匀布点法及随机布点法情况下的无网格法的计算精度,最后提出了两种布点方法结合起来应用的混合布点方法,得出了一些有益的结论.。     

基于径向基函数与B样条的散乱数据拟合方法

针对散乱数据的曲面拟合问题,提出一种径向基函数与B样条插值结合使用的曲面拟合方法.通过分片径向基函数插值,三维散乱点,再从分片插值曲面上获取预先设定好的有序网格点的值,最后利用张量积B样条插值有序网格点,从而得到拟合曲面.该方法较好地解决散乱数据插值和拟合的计算不稳定性问题,最后给出算法实例.     

基于Matlab的径向基函数数值计算软件包

为推动基于无网格方法的计算软件的发展,介绍基于Matlab自主开发的径向基函数（Radial Basis Function,RBF）数值计算软件包,阐述软件的理论基础、设计思路以及该软件包的功能和特点,并结合边界节点法（Boundary Knot Method,BKM）的数值实例给出软件的使用过程．该软件包可以根据不同的数学物理模型选择合适的数值算法来求解多种实际物理问题,也可对不同数值算法得到的结果进行比较．最后,总结应用Matlab进行数值计算软件开发的优缺点．     

插值边界的四边网格离散极小曲面建模方法

如何实现极小曲面的快速三维建模,是几何设计与计算领域中的难点和热点问题.给定一条封闭的边界离散折线,本文研究如何构造以其为边界的四边网格离散极小曲面.首先从曲面的内蕴微分几何度量出发,给出了离散四边网格极小曲面的数学定义;然后利用保长度边界投影、四边网格生成、径向基函数插值映射和非线性优化技术,提出了由给定边界离散折线快速构造离散四边网格极小曲面的一般技术框架.最后通过若干建模实例验证了本文方法的有效性.该方法可实现四边网格极小曲面的高质量建模,在建筑几何领域具有一定的应用价值.     

多元散乱数据的渐近正定径向基函数插值

本文以一元B样条和径向基函数为基础,构造了一种n元渐近正定径向基函数.因这种渐近正定径向基函数兼有一元B样条和已有径向基函数的优点,将其应用于n维欧氏空间的多元散乱数据插值,得到了一种新的高维散乱数据插值函数.数值例子表明,这种插值格式对n维欧氏空间多元散乱数据具有良好的逼近效果.     

求解半线性椭圆问题的牛顿-瀑布型两层网格法

选取一对合适的步长使用中心差分格式离散半线性椭圆问题形成粗网格和细网格,使用三次样条插值算子将粗网格上高精度近似解插值到细网格为其提供初始值,结合牛顿法提出了牛顿-瀑布型两层网格法．数值实验表明该算法具有稳健性强、计算效率高的优点．     

径向基/超限插值结构网格运动方法中的变形一致性问题

针对径向基函数插值和超限插值的结构网格运动方法(RBF_TFI)中变形一致性问题引起的空间网格棱线扭曲和网格单元交叉现象进行深入分析,并依据分析结果提出一种改进的RBF_TFI结构网格运动方`法.首先基于HIRENASD机翼模型计算网格绕y轴旋转30°的算例,从变形网格质量、棱线交叉点位移误差与空间棱线扭曲等方面分析了不同基点选择方案中变形一致性问题产生的原因;然后将与物面相连的空间棱线定义为边界棱线,并采用三次样条插值的方式更新边界棱线上网格点的坐标.使用HIRENASD算例对文中方法的效果进行了验证,结果表明,该方法能够消除变形一致性问题引起的棱线扭曲现象,有效地避免变形网格出现交叉和产生负体积单元.     

交互式棱柱网格生成方法

针对流动数值模拟中边界层区域棱柱网格生成与修改困难的问题,提出一种基于物面边界约束的交互式棱柱网格生成方法.该方法交互生成物面边界点的空间推进面网格;将边界点空间推进面网格的每层网格点看成边界网格点运动得到,计算对应的边界网格点位移;利用径向基函数插值将边界网格点位移光滑传递给内部网格点,获得内部网格点的空间推进网格;建立所有网格点之间的连接关系,构造棱柱网格单元.基于V型槽、飞机机头等外形对文中方法进行测试,并对该方法存在的一些问题进行了讨论;利用该方法生成一个大展弦比翼身组合体模型的棱柱/四面体混合计算网格,并在PHENGLEI计算平台上进行了数值模拟,结果表明,该方法能够生成高质量的粘性棱柱网格.     

Development of a meshless Galerkin boundary node method for viscous fluid flows

In this paper, a meshless Galerkin boundary node method is developed for boundary-only analysis of the interior and exterior incompressible viscous fluid flows, governed by the Stokes equations, in biharmonic stream function formulation. This method combines scattered points and boundary integral equations. Some of the novel features of this meshless scheme are boundary conditions can be enforced directly and easily despite the meshless shape functions lack the delta function property, and system matrices are symmetric and positive definite. The error analysis and convergence study of both velocity and pressure are presented in Sobolev spaces. The performance of this approach is illustrated and assessed through some numerical examples.     

A meshless Galerkin method for Stokes problems using boundary integral equations

A meshless Galerkin scheme for the simulation of two-dimensional incompressible viscous fluid flows in primitive variables is described in this paper. This method combines a boundary integral formulation for the Stokes equation with the moving least-squares (MLS) approximations for construction of trial and test functions for Galerkin approximations. Unlike the domain-type method, this scheme requires only a nodal structure on the bounding surface of a body for approximation of boundary unknowns, thus it is especially suitable for the exterior problems. Compared to other meshless methods such as the boundary node method and the element free Galerkin method, in which the MLS is also introduced, boundary conditions do not present any difficulty in using this meshless method. The convergence and error estimates of this approach are presented. Numerical examples are also given to show the efficiency of the method.     

An improved pseudospectral meshless radial point interpolation (PSMRPI) method for 3D wave equation with variable coefficients

In this paper, a pseudospectral meshless radial point interpolation (PSMRPI) technique is applied to the three-dimensional wave equation with variable coefficients subject to given appropriate initial and Dirichlet boundary conditions. The present method is a kind of combination of meshless methods and spectral collocation techniques. The point interpolation method along with the radial basis functions is used to construct the shape functions as the basis functions in the frame of the spectral collocation methods. These basis functions will have Kronecker delta function property, as well as unitary possession. In the proposed method, operational matrices of higher order derivatives are constructed and then applied. The merit of this innovative method is that, it does not require any kind of integration locally or globally over sub-domains, as it is essential in meshless methods based on Galerkin weak forms, such as element-free Galerkin and meshless local Petrov–Galerkin methods. Therefore, computational cost of PSMRPI method is low. Further, it is proved that the procedure is stable with respect to the time variable over some conditions on the 3D wave model, and the convergence of the technique is revealed. These latest claims are also shown in the numerical examples, which demonstrate that PSMRPI provides excellent rate of convergence.

     

基于无网格自然邻接点Petrov-Galerkin法求解带源参数瞬态热传导问题

基于无网格自然邻接点Petrov-Galerkin法,本文建立了一种求解带源参数瞬态热传导问题的新方法.为了克服移动最小二乘近似难以准确施加本质边界条件的缺点,采用了自然邻接点插值构造试函数.在局部多边形子域上采用局部Petrov-Galerkin方法建立瞬态热传导问题的积分弱形式.这些多边形子域可由Delaunay三角形创建.时间域则通过传统的两点差分法进行离散.最后通过算例验证了该数值算法的有效性和正确性.     

Application of NEM in seepage analysis with a free surface

As one kind of meshless methods, the natural element method (NEM) constructs shape functions based on the Voronoi diagrams, and it has advantages of both the conventional meshless method and the finite element method. Since the nodes are independent of the integral mesh, it is more suitable for the analysis of seepage with a free surface than the finite element method. In addition, its shape functions satisfy the Kronecker δ conditions, therefore, its boundary conditions can be dealt with much easily than those of such meshless method as element-free Galerkin method (EFGM). In this paper, the NEM was used in the seepage analysis of dams. The initial free surface was assumed first in the calculations, and the location of the free surface was adjusted according to the calculation results. The examples showed that the natural element method lead to satisfactory results.     

Pricing American options under jump-diffusion models using local weak form meshless techniques

Recently, several numerical methods have been proposed for pricing options under jump-diffusion models but very few studies have been conducted using meshless methods [R. Chan and S. Hubbert, A numerical study of radial basis function based methods for options pricing under the one dimension jump-diffusion model, Tech. Rep., 2010; A. Saib, D. Tangman, and M. Bhuruth, A new radial basis functions method for pricing American options under Merton's jump-diffusion model, Int. J. Comput. Math. 89 (2012), pp. 1164–1185]. Indeed, only a strong form of meshless methods have been employed in these lectures. We propose the local weak form meshless methods for option pricing under Merton and Kou jump-diffusion models. Predominantly in this work we will focus on meshless local Petrov–Galerkin, local boundary integral equation methods based on moving least square approximation and local radial point interpolation based on Wendland's compactly supported radial basis functions. The key feature of this paper is applying a Richardson extrapolation technique on American option which is a free boundary problem to obtain a fixed boundary problem. Also the implicit–explicit time stepping scheme is employed for the time derivative which allows us to obtain a spars and banded linear system of equations. Numerical experiments are presented showing that the presented approaches are extremely accurate and fast.     

Galerkin collocation for an improved boundary element method for a plane mixed boundary value problem

We present the numerical implementation of a boundary element approximation for the plane mixed boundary value problem of the Laplacian. The performed Galerkin procedure is based on the direct boundary integral method. Its accuracy is improved by using appropriate singularity functions as additional test and trial functions besides quadratic splines. We analize the consistency error of the used numerical integrations and present asymptotic error estimates for the fully discretized numerical scheme which are of the same optimal orders as the Galerkin errors.