The main goal of this paper is to solve a class of Darboux problems by converting them into the two-dimensional nonlinear Volterra integral equation of the second kind. The scheme approximates the solution of these integral equations using the discrete Galerkin method together with local radial basis functions, which use a small set of data instead of all points in the solution domain. We also employ the Gauss–Legendre integration rule on the influence domains of shape functions to compute the local integrals appearing in the method. Since the scheme is constructed on a set of scattered points and does not require any background meshes, it is meshless. The error bound and the convergence rate of the presented method are provided. Some illustrative examples are included to show the validity and efficiency of the new technique. Furthermore, the results obtained demonstrate that this method uses much less computer memory than the method established using global radial basis functions.

This article addresses the resolution of the inverse problem for the parameter identification in orthotropic materials with a number of measurements merely on the boundaries. The inverse problem is formulated as an optimization problem of a residual functional which evaluates the differences between the experimental and predicted displacements. The singular boundary method, an integration-free, mathematically simple and boundary-only meshless method, is employed to numerically determine the predicted displacements. The residual functional is minimized by the Levenberg-Marquardt method. Three numerical examples are carried out to illustrate the robustness, efficiency, and accuracy of the proposed scheme. In addition, different levels of noise are added into the boundary conditions to verify the stability of the present methodology.

A novel meshless scheme is proposed for inverse source identification problems of Helmholtz-type equations. It is formulated by the non-singular general solutions of the Helmholtz-type equations augmented with radial basis functions. Under this meshless scheme, we can determine smooth source terms from partially accessible boundary measurements with accurate results. Numerical examples are presented to verify validity and accuracy of the present scheme. It is demonstrated that the present scheme is simple, accurate, stable and computationally efficient for inverse smooth source identification problems.

In this paper, we study a class of contact smoothed particle hydrodynamics (SPH) by introducing Riemann solvers and using high-order limiters. In particular, a promising concept of WENO interpolation as limiter is presented in the reconstruction process. The physical values relating interactional particles used as the initial values of the Riemann problem can be reconstructed by the Taylor series expansion. The contact solvers of the Riemann problem at contact points are incorporated in SPH approximations. In order to keep the fluid density at the wall rows to be consistent with that of the inner fluid wall boundaries, several lines of dummy particles are placed outside of the solid walls, which are assigned according to the initial configuration. At last, the method is applied to compressible flows with sharp discontinuities such as the collision of two strong shocks and the interaction of two blast waves and so on. The numerical results indicate that the method is capable of handling sharp discontinuity and efficiently reducing unphysical oscillations.

Kernel-based methods are popular in computer graphics, machine learning, and statistics, among other fields; because they do not require meshing of the domain under consideration, higher dimensions and complicated domains can be managed with reasonable effort. Traditionally, the high order of accuracy associated with these methods has been tempered by ill-conditioning, which arises when highly smooth kernels are used to conduct the approximation. Recent advances in representing Gaussians using eigenfunctions have proven successful at avoiding this destabilization in scattered data approximation problems. This paper will extend these techniques to the solution of boundary value problems using collocation. The method of particular solutions will also be considered for elliptic problems, using Gaussian eigenfunctions to stably produce an approximate particular solution.

This paper presents three boundary meshless methods for solving problems of steady-state and transient heat conduction in nonlinear functionally graded materials (FGMs). The three methods are, respectively, the method of fundamental solution (MFS), the boundary knot method (BKM), and the collocation Trefftz method (CTM) in conjunction with Kirchhoff transformation and various variable transformations. In the analysis, Laplace transform technique is employed to handle the time variable in transient heat conduction problem and the Stehfest numerical Laplace inversion is applied to retrieve the corresponding time-dependent solutions. The proposed MFS, BKM and CTM are mathematically simple, easy-to-programming, meshless, highly accurate and integration-free. Three numerical examples of steady state and transient heat conduction in nonlinear FGMs are considered, and the results are compared with those from meshless local boundary integral equation method (LBIEM) and analytical solutions to demonstrate the efficiency of the present schemes.

An application of the meshless Local Radial Basis Function Collocation Method (LRBFCM) in solution of incompressible turbulent combined forced and natural convection is for the first time explored in the present paper. The turbulent flow equations are described by the low-Re number к – ε model with Launder and Sharma and Abe et al. closure coefficients. The involved temperature, velocity, pressure, turbulent kinetic energy and dissipation fields are represented on overlapping 5-noded sub-domains through the collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second order partial derivatives of the fields are calculated from the respective derivatives of the RBF’s. The involved equations are solved through the explicit time stepping. The pressure-velocity coupling is based on Chorin’s fractional step method. The adaptive upwinding technique, proposed by Lin and Atluri , is used because of the convection dominated situation. The solution procedure is represented for a 2D upward channel flow with differentially heated walls. The results have been assessed by achieving a reasonable agreement with the direct numerical simulation of Kasagi and Nishimura for Reynolds number 4494, based on the channel width, and Grash of number 9.6×105. The advantages of the represented mesh-free approach are its simplicity, accuracy, similar coding in 2D and 3D, and straightforward applicability in non-uniform node arrangements.