A Radial Basis Function Method with Improved Accuracy for Fourth Order Boundary Value Problems

Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.

An improved local radial point interpolation method for transient heat conduction analysis

The smoothing thin plate spline （STPS） interpolation using the penalty function method according to the optimization theory is presented to deal with transient heat conduction problems. The smooth conditions of the shape functions and derivatives can be satisfied so that the distortions hardly occur. Local weak forms are developed using the weighted residual method locally from the partial differential equations of the transient heat conduction. Here the Heaviside step function is used as the test function in each sub-domain to avoid the need for a domain integral. Essential boundary conditions can be implemented like the finite element method （FEM） as the shape functions possess the Kronecker delta property. The traditional two-point difference method is selected for the time discretization scheme. Three selected numerical examples are presented in this paper to demonstrate the availability and accuracy of the present approach comparing with the traditional thin plate spline （TPS） radial basis functions.

Factors affecting accuracy of radial point interpolation meshfree method for 3-D solid mechanics

Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.

Performance of Compact Radial Basis Functions in the Direct Interpolation Boundary Element Method for Solving Potential Problems

This study evaluates the effectiveness of a new technique that transforms doma in integrals into boundary integrals that is applicable to the boundary element method.Si mulations were conducted in which two-dimensional surfaces were approximated by inter polation using radial basis functions with full and compact supports.Examples involving Poisson’s equation are presented using the boundary element method and the proposed te chnique with compact radial basis functions.The advantages and the disadvantages are e xamined through simulations.The effects of internal poles,the boundary mesh refinemen t and the value for the support of the radial basis functions on performance are assessed.

SOME PROBLEMS WITH THE METHOD OF FUNDAMENTAL SOLUTION USING RADIAL BASIS FUNCTIONS

The present work describes the application of the method of fundamental solutions （MFS） along with the analog equation method （AEM） and radial basis function （RBF） approximation for solving the 2D isotropic and anisotropic Helmholtz problems with different wave numbers. The AEM is used to convert the original governing equation into the classical Poisson＇s equation, and the MFS and RBF approximations are used to derive the homogeneous and particular solutions, respectively. Finally, the satisfaction of the solution consisting of the homogeneous and particular parts to the related governing equation and boundary conditions can produce a system of linear equations, which can be solved with the singular value decomposition （SVD） technique. In the computation, such crucial factors related to the MFS-RBF as the location of the virtual boundary, the differential and integrating strategies, and the variation of shape parameters in multi-quadric （MQ） are fully analyzed to provide useful reference.

Radial point collocation method (RPCM) for solving convection-diffusion problems

In this paper, Radial point collocation method （RPCM）, a kind of meshfree method, is applied to solve convectiondiffusion problem. The main feature of this approach is to use the interpolation schemes in local supported domains based on radial basis functions. As a result, this method is local and hence the system matrix is banded which is very attractive for practical engineering problems. In the numerical examination, RPCM is applied to solve non-linear convection-diffusion 2D Burgers equations. The results obtained by RPCM demonstrate the accuracy and efficiency of the proposed method for solving transient fluid dynamic problems. A fictitious point scheme is adopted to improve the solution accuracy while Neumann boundary conditions exist. The meshfree feature of the nresent method is verv attractive in solving comnutational fluid nroblems.

Gaussian Radial Basis Function interpolation in vertical deformation analysis

In many deformation analyses,the partial derivatives at the interpolated scattered data points are required.In this paper,the Gaussian Radial Basis Functions(GRBF)is proposed for the interpolation and differentiation of the scattered data in the vertical deformation analysis.For the optimal selection of the shape parameter,which is crucial in the GRBF interpolation,two methods are used:the Power Gaussian Radial Basis Functions(PGRBF)and Leave One Out Cross Validation(LOOCV)(LGRBF).We compared the PGRBF and LGRBF to the traditional interpolation methods such as the Finite Element Method(FEM),polynomials,Moving Least Squares(MLS),and the usual GRBF in both the simulated and actual Interferometric Synthetic Aperture Radar(InSAR)data.The estimated results showed that the surface interpolation accuracy was greatly improved by LGRBF and PGRBF methods in comparison withFEM,polynomial,and MLS methods.Finally,LGRBF and PGRBF interpolation methods are used to compute invariant vertical deformation parameters,i.e.,changes in Gaussian and mean Curvatures in the Groningen area in the North of Netherlands.

A numerical approach based on the meshless collocation method in elastodynamics

In this paper, a collocation technique with the modified equilibrium on line method （ELM） for imposition of Neumann （natural） boundary conditions is presented for solving the two-dimensional problems of linear elastic body vibrations. In the modified ELM, equilibrium over the lines on the natural boundary is satisfied as Neumann boundary condition equations. In other words, the natural boundary conditions are satisfied naturally by using the weak formulation. The performance of the modified version of the ELM is studied for collocation methods based on two different ways to construct meshless shape functions： moving least squares approximation and radial basis point interpolation. Numerical examples of two-dimensional free and forced vibration analyses show that by using the modified ELM, more stable and accurate results would be obtained in comparison with the direct collocation method.

Numerical Assessment of Nanofluid Natural Convection Using Local RBF Method Coupled with an Artificial Compressibility Model

In this paper,natural heat convection inside square and equilateral triangular cavities was studied using a meshless method based on collocation local radial basis function(RBF).The nanofluids used were Cu-water or Al_(2)O_(3)-water mixture with nanoparticle volume fractions range of 0≤φ≤0.2.A system of continuity,momentum,and energy partial differential equations was used in modeling the flow and temperature behavior of the fluids.Partial derivatives in the governing equations were approximated using the RBF method.The artificial compressibility model was implemented to overcome the pressure velocity coupling problem that occurs in such equations.Themain goal of this work was to present a simple and efficient method to deal with complex geometries for a variety of problem conditions.To assess the accuracy of the proposed method,several test cases of natural convection in square and triangular cavities were selected.For Rayleigh numbers ranging from 103 to 105,a validation test of natural convection of Cu-water in a square cavity was used.The numerical investigation was then extended to Rayleigh number 106,as well as Al_(2)O_(3)-water nanofluid with a volume fraction range of 0≤φ≤0.2.In a second investigation,the same nanofluids were used in a triangular cavitywith varying volume fractions to test the proposed meshless approach on non-rectangular geometries.The numerical results appear to be in agreement with those from earlier investigations.Furthermore,the suggested meshless method was found to be stable and accurate,demonstrating that it may be a viable alternative for solving natural heat transfer equations of nanofluids in enclosures with irregular geometries.

Meshless Method of Lines for Numerical Solution of Kawahara Type Equations

In this work, an algorithm based on method of lines coupled with radial basis functions namely meshless method of lines (MMOL) is presented for the numerical solution of Kawahara, modified Kawahara and KdV Kawahara equations. The motion of a single solitary wave, interaction of two and three solitons and the phenomena of wave generation is discussed. The results are compared with the exact solution and with the results in the relevant literature to show the efficiency of the method.

An Efficient Meshless Method for Hyperbolic Telegraph Equations in (1 + 1) Dimensions

Numerical solutions of the second-order one-dimensional hyperbolic telegraph equations are presented using the radial basis functions.The purpose of this paper is to propose a simple novel direct meshless scheme for solving hyperbolic telegraph equations.This is fulfilled by considering time variable as normal space variable.Under this scheme,there is no need to remove time-dependent variable during the whole solution process.Since the numerical solution accuracy depends on the condition of coefficient matrix derived from the radial basis function method.We propose a simple shifted domain method,which can avoid the full-coefficient interpolation matrix easily.Numerical experiments performed with the proposed numerical scheme for several second-order hyperbolic telegraph equations are presented with some discussions.

Simulation of Oil-Water Flow in a Shale Reservoir Using a Radial Basis Function

Due to the difficulties associated with preprocessing activities and poor grid convergence when simulating shale reservoirs in the context of traditional grid methods,in this study an innovative two-phase oil-water seepage model is elaborated.The modes is based on the radial basis meshless approach and is used to determine the pressure and water saturation in a sample reservoir.Two-dimensional examples demonstrate that,when compared to the finite difference method,the radial basis function method produces less errors and is more accurate in predicting daily oil production.The radial basis function and finite difference methods provide errors of 5.78 percent and 7.5 percent,respectively,when estimating the daily oil production data for a sample well.A sensitivity analysis of the key parameters that affect the radial basis function’s computation outcomes is also presented.

A Local Meshless Method for Two Classes of Parabolic Inverse Problems

A local meshless method is applied to find the numerical solutions of two classes of inverse problems in parabolic equations. The problem is reconstructing the source term using a solution specified at some internal points;one class is that the source term is time dependent, and the other class is that the source term is time and space dependent. Some numerical experiments are presented and discussed.

Three-dimensional magnetotelluric forward modeling using structural-point EFGM and FEM coupling method

In this paper,the authors propose a method of three-dimensional(3D)magnetotelluric(MT)forward modeling algorithm based on the meshfree and finite element coupling method.The model is discretized by regular nodes in the central area,and the radial point interpolation method(RPIM)based on the global weakness is utilized to construct the meshfree shape function.The Governing equations in each background gird are solved by Gaussian integration.In the extended area where the points are sparsely distributed,to avoid the instability of the meshfree method,finite element method(FEM)with regular grids is used to solve the governing equation.Finally,the meshfree and finite element governing equations are coupled by the continuity of the field at the interfaces,and the direct solution technique is used to realize the 3D MT forward modeling.Numerical experiments of several typical electrical models are used to verify the effectiveness of the method.

极限下限分析的区域光滑径向点插值法

极限分析的高效数值计算方法在结构设计和安全评定中具有非常重要的作用.为了更加有效地求解极限分析问题,将区域光滑径向点插值法与二阶锥规划相结合,提出了理想弹塑性结构极限下限分析的一种新方法.将问题域离散为简单的三角形背景单元,每个单元进一步划分成若干个光滑域.为了将复杂的域积分转化为简单的边界积分,并且避免计算形函数的导数,采用广义梯度光滑技术对每个光滑域进行应变光滑处理.由于径向点插值法构造的形函数满足Kronecker delta性质,本质边界条件可以直接施加.依据下限定理,在满足以等效积分弱形式表达的自平衡应力场平衡条件的基础上,用二阶锥规划成功构建了极限分析下限法的计算模型,从而可方便地通过基于原始-对偶内点法的数学规划求解器MOSEK直接求解该问题.数值算例结果表明,本文所提方法有效地克服了维数障碍问题,具有较高的计算精度,并且计算结果对网格畸变十分不敏感.

二维水下声辐射问题的改进径向基点插值法研究

径向基点插值法是一种典型的无网格数值计算方法,在分析声学问题时,相比于传统有限元法能更好地抑制频散误差,且在相同的节点分布下通常可以得到更精确的数值解。本文提出一种改进的节点选取方案用于构造插值形函数,即改进径向基点插值法。该方案采取一个简单而直接的格式,可确保在进行数值积分时同一背景积分单元中的被积函数是连续可微的,从而减小数值积分误差,得到比原始径向基点插值法更精确的数值解。同时,为了处理外声场问题,本文采用DtN映射技术将无限域截断为有界计算域,满足索默菲尔德辐射条件。数值试验表明,相比于传统有限元法和原始径向基点插值法,本文改进方法具有更高的计算精度和计算效率,在研究水下声辐射问题时具有良好的应用前景。

无人机DEM测点插值算法的矿区沉陷测绘技术

为实现对矿区沉陷深度的精准测量,完善测绘技术的具体实施方案,针对无人机DEM测点插值算法的矿区沉陷测绘技术展开研究;根据DEM测点获取结果,实施分维值测定处理,再以此为基础,确定插值拐点所处区域,实现对无人机DEM测点插值算法应用流程的完善;联合像控点布设原则,定义空中三角区域,并通过推导径向插值基函数的方式,建立矿区沉陷区域的地理模型;在空间坐标系转换条件的基础上,确定矢量数据叠加的处理结果,再结合像片倾角与旋偏角的实际取值,求解测量精度评价指标,完成无人机DEM测点插值算法的矿区沉陷测绘技术实施方法的设计;实验结果表明,上述测绘技术方法作用下,沉陷深度测量结果与矿区真实沉陷深度之间的差值不超过1 m,符合精准测量的应用需求,对于测绘实施方案的完善可以起到一定的促进性影响作用。

某关停铅蓄电池企业土壤重金属污染特征及生态风险评价

为探明江苏省某铅蓄电池企业生产6年关停后对该地块土壤是否造成污染,对土壤中Pb、Cu、Zn、Hg、As、Cd、Ni、Cr这8种重金属采样分析,基于主成分分析和径向基函数(RBF)插值法分析地块土壤重金属的污染分布特征。利用单因子指数法和Nemerow综合指数法对污染地块土壤环境质量分析,对于涉及人体健康生态风险评价采用Hakanson潜在生态危害指数法开展综合评价。结果表明:研究区域土壤已遭受明显的人为外源重金属污染,且局部重金属富集程度较高,纵向污染区域主要集中在100cm以上土壤中,区域整体污染程度达到轻度污染,其中危废库局部土壤受污染程度已为中度污染,建议对该单元土壤开展进一步风险评估和修复;研究区域土壤尚不存在潜在生态风险,但是需重点关注危废库构筑物单元的土壤环境质量,主要威胁来自土壤中的镉,其EiCd潜在生态风险指数等级已达到强的级别。建议该地块后续开发利用时不作为农用地使用或进行局部区域土壤修复治理。

FREE VIBRATION ANALYSIS OF ARBITRARY SHAPED PLATES BY BOUNDARY KNOT METHOD

The boundary knot method （BKM） is a truly meshless boundary-type radial basis function （RBF） collocation scheme, where the general solution is employed instead of the fundamental solution to avoid the fictitious outside boundary of the physical domain of interest. In this study, the BKM is first used to calculate the free vibration of free and simply-upported thin plates. Compared with the analytical solution and ANSYS （a commercial FEM code） results, the present BKM is highly accurate and fast convergent.

ELASTIC DYNAMIC ANALYSIS OF MODERATELY THICK PLATE USING MESHLESS LRPIM

A meshless local radial point interpolation method （LRPIM） for solving elastic dy-namic problems of moderately thick plates is presented in this paper. The discretized system equation of the plate is obtained using a locally weighted residual method. It uses a radial basis function （RBF） coupled with a polynomial basis function as a trial function,and uses the quartic spline function as a test function of the weighted residual method. The shape function has the properties of the Kronecker delta function,and no additional treatment is done to impose essen-tial boundary conditions. The Newmark method for solving the dynamic problem is adopted in computation. Effects of sizes of the quadrature sub-domain and influence domain on the dynamic properties are investigated. The numerical results show that the presented method can give quite accurate results for the elastic dynamic problem of the moderately thick plate.