Method of fundamental solutions with optimal regularization techniques for the Cauchy problem of the Laplace equation with singular points

The purpose of this study is to propose a high-accuracy and fast numerical method for the Cauchy problem of the Laplace equation. Our problem is directly discretized by the method of fundamental solutions (MFS). The Tikhonov regularization method stabilizes a numerical solution of the problem for given Cauchy data with high noises. The accuracy of the numerical solution depends on a regularization parameter of the Tikhonov regularization technique and some parameters of the MFS. The L-curve determines a suitable regularization parameter for obtaining an accurate solution. Numerical experiments show that such a suitable regularization parameter coincides with the optimal one. Moreover, a better choice of the parameters of the MFS is numerically observed. It is noteworthy that a problem whose solution has singular points can successfully be solved. It is concluded that the numerical method proposed in this paper is effective for a problem with an irregular domain, singular points, and the Cauchy data with high noises.     

Three Boundary Meshless Methods for Heat Conduction Analysis in Nonlinear FGMs with Kirchhoff and Laplace Transformation

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.     

Burton–Miller-type singular boundary method for acoustic radiation and scattering

This paper proposes the singular boundary method (SBM) in conjunction with Burton and Miller?s formulation for acoustic radiation and scattering. The SBM is a strong-form collocation boundary discretization technique using the singular fundamental solutions, which is mathematically simple, easy-to-program, meshless and introduces the concept of source intensity factors (SIFs) to eliminate the singularities of the fundamental solutions. Therefore, it avoids singular numerical integrals in the boundary element method (BEM) and circumvents the troublesome placement of the fictitious boundary in the method of fundamental solutions (MFS). In the present method, we derive the SIFs of exterior Helmholtz equation by means of the SIFs of exterior Laplace equation owing to the same order of singularities between the Laplace and Helmholtz fundamental solutions. In conjunction with the Burton–Miller formulation, the SBM enhances the quality of the solution, particularly in the vicinity of the corresponding interior eigenfrequencies. Numerical illustrations demonstrate efficiency and accuracy of the present scheme on some benchmark examples under 2D and 3D unbounded domains in comparison with the analytical solutions, the boundary element solutions and Dirichlet-to-Neumann finite element solutions.     

基于局部基本解法的静电场仿真分析

将局部基本解方法应用于静电场问题的模拟与分析。局部基本解方法是利用控制方程的基本解,基于局部理论和移动最小二乘原理提出的一种无网格算法。相比于有限元和有限差分等传统网格类方法,该方法仅需离散节点,避免了复杂的网格剖分难题。作为一种半解析数值技术,物理问题的基本解被作为插值基函数建立数值离散模型,从而保证了算法的较高精度。此外,与具有全局离散格式的无网格方法相比,局部基本解法更适用于高维复杂几何和大尺度模拟。二维和三维数值试验表明,该方法具有实施方便灵活,计算精度高和计算速度快等优势。为静电场仿真研究开辟新的途径,拓展了局部基本解方法的应用领域。     

Singular meshless method using double layer potentials for exterior acoustics

Time-harmonic exterior acoustic problems are solved by using a singular meshless method in this paper. It is well known that the source points cannot be located on the real boundary, when the method of fundamental solutions (MFS) is used due to the singularity of the adopted kernel functions. Hence, if the source points are right on the boundary the diagonal terms of the influence matrices cannot be derived. Herein we present an approach to obtain the diagonal terms of the influence matrices of the MFS for the numerical treatment of exterior acoustics. By using the regularization technique to regularize the singularity and hypersingularity of the proposed kernel functions, the source points can be located on the real boundary and therefore the diagonal terms of influence matrices are determined. We also maintain the prominent features of the MFS, that it is free from mesh, singularity, and numerical integration. The normal derivative of the fundamental solution of the Helmholtz equation is composed of a two-point function, which is one of the radial basis functions. The solution of the problem is expressed in terms of a double-layer potential representation on the physical boundary based on the potential theory. The solutions of three selected examples are used to compare with the results of the exact solution, conventional MFS, boundary element method, and Dirichlet-to-Neumann finite element method. Good numerical performance is demonstrated by close agreement with other solutions.     

New method for solving fractional partial integro-differential equations by combination of Laplace transform and resolvent kernel method

In this article, we obtained the approximate solution for a new class of Time-Fractional Partial Integro-Differential Equation (TFPIDE) of the Caputo-Volterra type in which the integral is not limited to the convolution type. This new class of TFPIDE is distinct from the common problem with the convolution integral kernel. The general expression of the analytical solution for this special type of TFPIDE was derived using a combination of Laplace transform and the resolvent kernel method. In the process, Laplace transform will transform the equation into a second kind Volterra integral equation in terms of the transformed function. Two main problems in deriving the approximate analytical solutions were identified as Case I and Case II problems. To obtain the approximate solutions for Case I and Case II problems, numerical methods were designed based on approximation of the resolvent kernel with truncated Neumann series as well as approximation of the Laplace transform based on truncated Taylor series. Several numerical examples are presented to indicate the plausibility, mechanism and performance of the proposed methods.     

Mathematical and numerical studies on meshless methods for exterior unbounded domain problems

The method of fundamental solutions (MFS) is an efficient meshless method for solving boundary value problems in an exterior unbounded domain. The numerical solution obtained by the MFS is accurate, while the corresponding matrix equation is ill-conditioned. A modified MFS (MMFS) with the proper basis functions is proposed by the introduction of the modified Trefftz method (MTM). The concrete expressions of the corresponding condition numbers are given in mathematical forms and the solvability by these methods is mathematically proven. Thereby, the optimal parameter minimizing the condition number is also mathematically given. Numerical experiments show that the condition numbers of the matrices corresponding to the MTM and the MMFS are reduced and that the numerical solution by the MMFS is more accurate than the one by the conventional method.     

The method of fundamental solutions for 2D and 3D Stokes problems

A numerical scheme based on the method of fundamental solutions (MFS) is proposed for the solution of 2D and 3D Stokes equations. The fundamental solutions of the Stokes equations, Stokeslets, are adopted as the sources to obtain flow field solutions. The present method is validated through other numerical schemes for lid-driven flows in a square cavity and a cubic cavity. Test results obtained for a rectangular cavity with wave-shaped bottom indicate that the MFS is computationally efficient than the finite element method (FEM) in dealing with irregular shaped domain. The paper also discusses the effects of number of source points and their locations on the numerical accuracy.     

An adaptive greedy technique for inverse boundary determination problem

In this paper, the method of fundamental solutions (MFS) is employed for determining an unknown portion of the boundary from the Cauchy data specified on parts of the boundary. We propose a new numerical method with adaptive placement of source points in the MFS to solve the inverse boundary determination problem. Since the MFS source points placement here is not trivial due to the unknown boundary, we employ an adaptive technique to choose a sub-optimal arrangement of source points on various fictitious boundaries. Afterwards, the standard Tikhonov regularization method is used to solve ill-conditional matrix equation, while the regularization parameter is chosen by the L-curve criterion. The numerical studies of both open and closed fictitious boundaries are considered. It is shown that the proposed method is effective and stable even for data with relatively high noise levels.     

A multilayer method of fundamental solutions for Stokes flow problems

The method of fundamental solutions (MFS) is a meshless method for the solution of boundary value problems and has recently been proposed as a simple and efficient method for the solution of Stokes flow problems. The MFS approximates the solution by an expansion of fundamental solutions whose singularities are located outside the flow domain. Typically, the source points (i.e. the singularities of the fundamental solutions) are confined to a smooth source layer embracing the flow domain. This monolayer implementation of the MFS (monolayer MFS) depends strongly on the location of the user-defined source points: On the one hand, increasing the distance of the source points from the boundary tends to increase the convergence rate. On the other hand, this may limit the achievable accuracy. This often results in an unfavorable compromise between the convergence rate and the achievable accuracy of the MFS. The idea behind the present work is that a multilayer implementation of the MFS (multilayer MFS) can improve the robustness of the MFS by efficiently resolving different scales of the solution by source layers at different distances from the boundary. We propose a block greedy-QR algorithm (BGQRa) which exploits this property in a multilevel fashion. The proposed multilayer MFS is much more robust than the monolayer MFS and can compute Stokes flows on general two- and three-dimensional domains. It converges rapidly and yields high levels of accuracy by combining the properties of distant and close source points. The block algorithm alleviates the overhead of multiple source layers and allows the multilayer MFS to outperform the monolayer MFS.     

Dirac方程的数值解法及其在原子总能量计算中的应用

中心势近似下径向Dirac方程的求解是相对论性原子(离子)结构计算的基础.本文通过相对论性方程中径向波函数大分量与非相对论方程径向波函数的类比,提出了径向Dirac方程的一种数值解法.为了验证数值解法的精度和可靠性,首先将数值结果与类氢势作用下的解析解进行比较.然后,将这种算法扩展到基于解析势的相对论性原子结构计算中,并将计算出的总能量与实验结果和其他方法得到的结果进行对比.     

Least-squares collocation meshless approach for radiative heat transfer in absorbing and scattering media

A least-squares collocation meshless method is employed for solving the radiative heat transfer in absorbing, emitting and scattering media. The least-squares collocation meshless method for radiative transfer is based on the discrete ordinates equation. A moving least-squares approximation is applied to construct the trial functions. Except for the collocation points which are used to construct the trial functions, a number of auxiliary points are also adopted to form the total residuals of the problem. The least-squares technique is used to obtain the solution of the problem by minimizing the summation of residuals of all collocation and auxiliary points. Three numerical examples are studied to illustrate the performance of this new solution method. The numerical results are compared with the other benchmark approximate solutions. By comparison, the results show that the least-squares collocation meshless method is efficient, accurate and stable, and can be used for solving the radiative heat transfer in absorbing, emitting and scattering media.     

Solution procedure of residue harmonic balance method and its applications

This paper presents a simple and rigorous solution procedure of residue harmonic balance for predicting the accurate approximation of certain autonomous ordinary differential systems. In this solution procedure, no small parameter is assumed. The harmonic residue of balance equation is separated in two parts at each step. The first part has the same number of Fourier terms as the present order of approximation and the remaining part is used in the subsequent improvement. The corrections are governed by linear ordinary differential equation so that they can be solved easily by means of harmonic balance method again. Three kinds of different differential equations involving general, fractional and delay ordinary differential systems are given as numerical examples respectively. Highly accurate limited cycle frequency and amplitude are captured. The results match well with the exact solutions or numerical solutions for a wide range of control parameters. Comparison with those available shows that the residue harmonic balance solution procedure is very effective for these autonomous differential systems. Moreover, the present method works not only in predicting the amplitude but also the frequency of bifurcated period solution for delay ordinary differential equation.     

A review on the solution of Grad–Shafranov equation in the cylindrical coordinates based on the Chebyshev collocation technique

Equilibrium reconstruction consists of identifying, from experimental measurements, a distribution of the plasma current density that satisfies the pressure balance constraint. Numerous methods exist to solve the Grad–Shafranov equation, describing the equilibrium of plasma confined by an axisymmetric magnetic field. In this paper, we have proposed a new numerical solution to the Grad–Shafranov equation (an axisymmetric, magnetic field transformed in cylindrical coordinates solved with the Chebyshev collocation method) when the source term (current density function) on the right-hand side is linear. The Chebyshev collocation method is a method for computing highly accurate numerical solutions of differential equations. We describe a circular cross-section of the tokamak and present numerical result of magnetic surfaces on the IR-T1 tokamak and then compare the results with an analytical solution.     

Collocation Methods for Hyperbolic Partial Differential Equations with Singular Sources

A numerical study is given on the spectral methods and the high order WENO finite difference scheme for the solution of linear and nonlinear hyperbolic partial differential equations with stationary and non-stationary singular sources. The singular source term is represented by the $δ$-function. For the approximation of the $δ$-function, the direct projection method is used that was proposed in [6]. The $δ$-function is constructed in a consistent way to the derivative operator. Nonlinear sine-Gordon equation with a stationary singular source was solved with the Chebyshev collocation method. The $δ$-function with the spectral method is highly oscillatory but yields good results with small number of collocation points. The results are compared with those computed by the second order finite difference method. In modeling general hyperbolic equations with a non-stationary singular source, however, the solution of the linear scalar wave equation with the non-stationary singular source using the direct projection method yields non-physical oscillations for both the spectral method and the WENO scheme. The numerical artifacts arising when the non-stationary singular source term is considered on the discrete grids are explained.     

A comparison of differential and integral equations of radiative transfer

The equations of radiative transfer and of statistical equilibrium of a two-level atom are solved by means of differential and integral equations for a one-dimensional medium. The numerical solutions are compared to the analytic solution. It is found that the integral equation for piecewise quadratic source functions gives more accurate results than does the differential equation.     

An Improved Formulation of Singular Boundary Method

This study proposes a new formulation of singular boundary method (SBM) to solve the 2D potential problems, while retaining its original merits being free of integration and mesh, easy-to-program, accurate and mathematically simple without the requirement of a fictitious boundary as in the method of fundamental solutions (MFS). The key idea of the SBM is to introduce the concept of the origin intensity factor to isolate the singularity of fundamental solution so that the source points can be placed directly on the physical boundary. This paper presents a new approach to derive the analytical solution of the origin intensity factor based on the proposed subtracting and adding-back techniques. And the troublesome sample nodes in the ordinary SBM are avoided and the sample solution is also not necessary for the Neumann boundary condition. Three benchmark problems are tested to demonstrate the feasibility and accuracy of the new formulation through detailed comparisons with the boundary element method (BEM), MFS, regularized meshless method (RMM) and boundary distributed source (BDS) method.     

Numerical study of paraxial beam formation

A method for numerical study of different integral problems with sufficiently strongly localized solutions is considered. Approximations of the solution with different degrees of accuracy were constructed and, using the existing approximation, an approach to an increase of the solution accuracy is proposed. Numerical solutions are obtained for the problem of free oscillations of laser cavities the optical surface of which was approximated by polynomials of the second and fourth degree inclusive. The obtained numerical solutions were stable.     

A Jacobi collocation approximation for nonlinear coupled viscous Burgers’ equation

This article presents a numerical approximation of the initial-boundary nonlinear coupled viscous Burgers’ equation based on spectral methods. A Jacobi-Gauss-Lobatto collocation (J-GL-C) scheme in combination with the implicit Runge-Kutta-Nyström (IRKN) scheme are employed to obtain highly accurate approximations to the mentioned problem. This J-GL-C method, based on Jacobi polynomials and Gauss-Lobatto quadrature integration, reduces solving the nonlinear coupled viscous Burgers’ equation to a system of nonlinear ordinary differential equation which is far easier to solve. The given examples show, by selecting relatively few J-GL-C points, the accuracy of the approximations and the utility of the approach over other analytical or numerical methods. The illustrative examples demonstrate the accuracy, efficiency, and versatility of the proposed algorithm.     

Axisymmetric solutions to fractional diffusion-wave equation in a cylinder under Robin boundary condition

The axisymmetric time-fractional diffusion-wave equation with the Caputo derivative of the order 0 < α ≤ 2 is considered in a cylinder under the prescribed linear combination of the values of the sought function and the values of its normal derivative at the boundary. The fundamental solutions to the Cauchy, source, and boundary problems are investigated. The Laplace transform with respect to time and finite Hankel transform with respect to the radial coordinate are used. The solutions are obtained in terms of Mittag-Leffler functions. The numerical results are illustrated graphically.