Two-step discretization method for 2D/3D Allen-Cahn equation base on RBF-FD scheme

A hybrid version of the strong form meshless Radial Basis Function-Finite Difference (RBF-FD) method is introduced for solving thermo-mechanics. The thermal model is spatially discretised with RBF-FD, where trial functions are polyharmonic splines augmented with polynomials. For time discretisation, the explicit Euler method is employed. An extension of RBF-FD, the hybrid RBF-FD, is introduced for solving mechanical problems. The model is one-way coupled, where temperature affects displacements. The thermo-elastoplastic material response is considered where the stress field is generally non-smooth. The hybrid RBF-FD, where the finite difference method is used to discretise the divergence operator from the balance equation, is shown to be successful when dealing with such problems. The mechanical model is introduced in a plane strain and in a generalised plane strain (GPS) assumption. For the first time, this work presents a strong form RBF-FD for GPS problems subjected to integral form constraints. The proposed method is assessed regarding h-convergence and accuracy on the benchmark with heating an elastoplastic square. It is proven to be successful at solving one-way coupled thermo-elastoplastic problems. The proposed novel meshless approach is efficient, accurate, and robust. Its use in an industrial situation is provided in Part 2 of this paper.

Numerical study of Fisher’s equation by the RBF-FD method

RBF-FD method for the high dimensional time fractional convection-diffusion equation

A local radial basis function-finite difference (RBF-FD) method for solving 1D and 2D coupled Schrödinger-Boussinesq (SBq) equations

This study introduces a deep neural network approach that utilizes radial basis functions (RBFs) to solve forward and inverse problems in the process of diffusion. The input layer incorporates multiquadric (MQ) RBFs, symbolizing the radial distance between the boundary points on the spacetime boundary and the source points positioned outside the spacetime boundary. The output layer is the initial and boundary data given by analytical solutions of the diffusion equation. Utilizing the concept of the spacetime coordinates, the approximations for forward and backward diffusion problems involve assigning initial data on the bottom or top spacetime boundaries, respectively. As the need for discretization of the governing equation is eliminated, our straightforward approach uses only the provided boundary data and MQ RBFs. To validate the proposed method, various diffusion scenarios, including forward, backward, and inverse problems with noise, are examined. Results indicate that the method can achieve high-precision numerical solutions for solving diffusion problems. Notably, only 1/4 of the initial and boundary conditions are known, yet the method still yields precise results.

On the static, vibration, and transient responses of micro-plates made of materials with different microstructures

In this investigation, we concentrate on solving the regularized long-wave (RLW) and extended Fisher–Kolmogorov (EFK) equations in one-, two-, and three-dimensional cases by a local meshless method called radial basis function (RBF)–finite-difference (FD) method. This method at each stencil approximates differential operators such as finite-difference method. In each stencil, it is necessary to solve a small-sized linear system with conditionally positive definite coefficient matrix. This method is relatively efficient and has low computational cost. How to choose the shape parameter is a fundamental subject in this method, since it has a palpable effect on coefficient matrix. We will employ the optimal shape parameter which results from algorithm of Sarra (Appl Math Comput 218:9853–9865, 2012). Then, we compare the approximate solutions acquired by RBF–FD method with theoretical solution and also with results obtained from other methods. The numerical results show that the RBF–FD method is suitable and robust for solving the RLW and EFK equations.

Local radial basis function-based differential quadrature (RBF-DQ) method was recently proposed by us. The method is a natural mesh-free approach. Like the conventional differential quadrature (DQ) method, it discretizes any derivative at a knot by a weighted linear sum of functional values at its neighbouring knots, which may be distributed randomly. However, different from the conventional DQ method, the weighting coefficients in present method are determined by taking the radial basis functions (RBFs) instead of high order polynomials as the test functions. The method works in a similar fashion as conventional finite difference schemes but with truly mesh-free property. In this paper, we mainly concentrate on the multiquadric (MQ) radial basis functions since they have exponential convergence. The effects of shape parameter c on the accuracy of numerical solution of linear and nonlinear partial differential equations are studied, and how the value of optimal c varies with the number of local support knots is also numerically demonstrated. The proposed method is validated by its application to solve incompressible Navier-Stokes equations. Excellent numerical results are obtained on an irregular knot distribution.

Semi-discretized numerical solution for time fractional convection–diffusion equation by RBF-FD

Stability analysis for a maximum principle preserving explicit scheme of the Allen–Cahn equation

Transient analysis of composite plates by a local radial basis functions-finite difference techique

Radial Basis Function Finite Difference method for solving the generalized time-fractional Burgers equation with three types of boundary conditions

ABSTRACTA framework is proposed to bridge coarse and fine meshes in a single simulation within the phase field method for fracture. Fine meshes are used in the vicinity of localized defects to accurately capture crack initiation, while coarse meshes are used away from initial defects and include only crack propagation paths. This reduces the prohibitive computational times associated with uniformly fine meshes over the entire domain, or with the use of complex adaptive meshes in the phase field method. The coupling between the two overlapping meshes is achieved using a variational formulation in which the energies of the models associated with the fine and coarse meshes are weighted in the superposition zone. Two situations are considered. The first includes the resolution of the phase field problem only in the fine mesh, while the coarse mesh is limited to the undamaged elastic problem. In the second situation, cracks can propagate in both the fine and coarse mesh. Variational formulations and associated finite element implementations are detailed. Numerical examples are presented, showing the potential of this approach to significantly reduce computational costs in the phase field method for cracking without affecting the accuracy.

A high-order accurate flow solution method for quadrilateral elements is developed and tested. The method employs the flux reconstruction (FR) approach proposed by Huynh. Ideal orders of accuracy are achieved for several problems on structured-based grids while some deterioration is observed on purely unstructured grids. Cares for curved wall boundary are also shown. The localized artificial diffusivity (LAD) method is incorporated to capture discontinuities that is crucial for compressible flow simulations. The FR+LAD works very well for 1D test cases but multi-dimensional shock problems requires further modifications to confine the artificial diffusivity to the shock location. Hybrid scheme of the FR and finite difference method taking advantage of 1D-feature alleviates the instability. I. Introduction Expanding demands for the computational fluid dynamics (CFD) requires more reliable numerical schemes. Specifically, turbulent flow simulations, acoustic simulation, combustion simulations and their combined problems. The key words in such simulations may be high accuracy and high resolution. In structured grid CFD, a compact finite difference scheme satisfies both requirements at higher level than existing other methods. If we put more emphasis on easier handling of complex geometries and flow-data dependent grid generation (grid adaptation), then unstructured grid method should also be explored and it should depart from 2nd-order scheme with limiter. Recent efforts in this area are development of the discontinuous Galerkin (DG) method, 1,2 the spectral finite volume (SV) method 3 and the staggered-grid multi-domain method spectral finite difference (SG) method 4,5 and its evolutional method called a spectral difference method. 6 The common and maybe the only way for those unstructured methods is to add degrees of freedom in each cell to achieve high order since the wide stencil of the original grid is not possible due to unstructured nature. The ENO-type method that uses original degrees of freedom in the grids and extends adaptive stencils to the flow solutions is also explored but the difficulty of extending stencils due to undirectional nature of unstructured grids and the unsuitability for parallel computing prohibits wide applications. A flux reconstruction (FR) scheme 7,8 studied in the paper is the new one among the former group. It uses differential form of conservation laws like the SD method but it is more general and includes equivalent counter part of the schemes in the group. In this paper, the performance of the FR is tested by a newly developed Navier-Stokes solver. Realistic flow applications of the FR in the Navier-Stokes equations are seldom seen since it is relatively new scheme. The Localized artificial diffusivity (LAD) 9–13 is incorporated for the shock capturing in the FR. The last item is motivated by the challenge and success of combining the SD and LAD in the recent study. 11,12 Wang 14 developed the lifting collocation penalty (LCP) method to extend the FR to triangular or tetrahedral grids and applied the LCP-FR combined method for hybrid unstructured grids. In the present study, only quadrilateral cells are treated to directly apply the LAD method that was developed on structured grids. The shock capturing by the FR+LAD high order scheme is tried.

The radial basis function (RBF) method, especially the multiquadric (MQ) function, was introduced in solving linear integral equations. The procedure of MQ method includes that the unknown function was firstly expressed in linear combination forms of RBFs, then the integral equation was transformed into collocation matrix of RBFs, and finally, solving the matrix equation and an approximation solution was obtained. Because of the superior interpolation performance of MQ, the method can acquire higher precision with fewer nodes and low computations which takes obvious advantages over thin plate splines (TPS) method. In implementation, two types of integration schemes as the Gauss quadrature formula and regional split technique were put forward. Numerical results showed that the MQ solution can achieve accuracy of1E-5. So, the MQ method is suitable and promising for integral equations.