C.C. Tsai¹, Y.C. Lin², D.L. Young^2,3, S.N. Atluri⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.2, pp. 103-114, 2006, DOI:10.3970/cmes.2006.016.103

Abstract In the applications of the method of fundamental solutions, locations of sources are treated either as variables or a priori known constants. In which, the former results in a nonlinear optimization problem and the other has to face the problem of locating sources. Theoretically, farther sources results in worse conditioning and better accuracy. In this paper, a practical procedure is provided to locate the sources for various time-independent operators, including Laplacian, Helmholtz operator, modified Helmholtz operator, and biharmonic operator. Wherein, the procedure is developed through systematic numerical experiments for relations among the accuracy, condition number, and source positions in different… More >

K.H. Chen¹, J.T. Chen², J.H. Kao³

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.1, pp. 27-40, 2006, DOI:10.3970/cmes.2006.016.027

Abstract In this paper, we employ the regularized meshless method (RMM) to search for eigenfrequency of two-dimension acoustics with multiply-connected domain. The solution is represented by using the double layer potentials. The source points can be located on the physical boundary not alike method of fundamental solutions (MFS) after using the proposed technique to regularize the singularity and hypersingularity of the kernel functions. The troublesome singularity in the MFS methods is desingularized and the diagonal terms of influence matrices are determined by employing the subtracting and adding-back technique. Spurious eigenvalues are filtered out by using singular value decomposition (SVD) updating term… More >

W.J. Santos¹ , J.A.F. Santiago¹, J.C.F Telles¹

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.1, pp. 23-40, 2012, DOI:10.3970/cmes.2012.087.023

Abstract The aim of this paper is to present numerical simulations of Cathodic Protection (CP) Systems using a Genetic Algorithm (GA) and the Method of Fundamental Solutions (MFS). MFS is used to obtain the solution of the associated homogeneous equation with the non-homogeneous equation subject to nonlinear boundary conditions defined as polarization curves. The adopted GA minimizes a nonlinear error function, whose design variables are the coefficients of the linear superposition of fundamental solutions and the positions of the source points, located outside the problem domain. In this work, the anodes added to the CP system are considered as point sources… More >

C.H. Tsai¹, D.L. Young², J. Kolibal³

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.1, pp. 53-72, 2010, DOI:10.3970/cmes.2010.066.053

Abstract The time evolution method of fundamental solutions (MFS) is proposed to solve backward heat conduction problems (BHCPs). The time evolution MFS belongs to one of the mesh-free numerical methods and is essentially composed of a sequence of diffusion fundamental solutions which exactly satisfy the heat conduction equations. Through correct treatment of temporal evolution, the resulting system of the time evolution MFS is smaller, and effectively decreases the possibility of ill-conditioning induced by such strongly ill-posed problems. Both one-dimensional and two-dimensional BHCPs are examined in this study, and the numerical results demonstrate the accuracy and stability of the MFS, especially for… More >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.3, pp. 221-246, 2010, DOI:10.3970/cmes.2010.060.221

Abstract We investigate the implementation of the method of fundamental solutions (MFS), in an iterative manner, for the algorithm of Kozlov, Maz'ya and Fomin (1991) in the case of the Cauchy problem in two-dimensional isotropic linear elasticity. At every iteration, two mixed well-posed and direct problems are solved using the Tikhonov regularization method, while the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion is also presented. The iterative MFS algorithm is tested for Cauchy problems for isotropic linear elastic materials to confirm the numerical convergence, stability and accuracy of… More >

Chia-Cheng Tsai¹, Chein-Shan Liu², Wei-Chung Yeih³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.2, pp. 131-152, 2010, DOI:10.3970/cmes.2010.056.131

Abstract The fictitious time integration method (FTIM) previously developed by Liu and Atluri (2008a) is combined with the method of fundamental solutions and the Chebyshev polynomials to solve Poisson-type nonlinear PDEs. The method of fundamental solutions with Chebyshev polynomials (MFS-CP) is an exponentially-convergent meshless numerical method which is able to solving nonhomogeneous partial differential equations if the fundamental solution and the analytical particular solutions of the considered operator are known. In this study, the MFS-CP is extended to solve Poisson-type nonlinear PDEs by using the FTIM. In the solution procedure, the FTIM is introduced to convert a Poisson-type nonlinear PDE into… More >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.48, No.2, pp. 121-154, 2009, DOI:10.3970/cmes.2009.048.121

Abstract We investigate the numerical implementation of the alternating iterative algorithm originally proposed by ` 12 ` 12 `$12 `&12 `#12 `^12 `_12 `%12 `~12 *Kozlov91 in the case of the Cauchy problem for the two-dimensional Laplace equation using a meshless method. The two mixed, well-posed and direct problems corresponding to every iteration of the numerical procedure are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases… More >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.46, No.3, pp. 221-254, 2009, DOI:10.3970/cmes.2009.046.221

Abstract We investigate the stable numerical reconstruction of an unknown portion of the boundary of a two-dimensional domain occupied by a functionally graded material (FGM) from a given boundary condition on this part of the boundary and additional Cauchy data on the remaining known portion of the boundary. The aforementioned inverse geometric problem is approached using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. The optimal value of the regularization parameter is chosen according to Hansen's L-curve criterion. Various examples are considered in order to show that the proposed method is numerically stable with respect to… More >

Chia-Cheng Tsai¹

CMES-Computer Modeling in Engineering & Sciences, Vol.45, No.3, pp. 249-272, 2009, DOI:10.3970/cmes.2009.045.249

Abstract Analytical particular solutions of Chebyshev polynomials are obtained for problems of Reissner plates under arbitrary loadings, which are governed by three coupled second-ordered partial differential equation (PDEs). Our solutions can be written explicitly in terms of monomials. By using these formulas, we can obtain the approximate particular solution when the arbitrary loadings have been represented by a truncated series of Chebyshev polynomials. In the derivations of particular solutions, the three coupled second-ordered PDE are first transformed into a single six-ordered PDE through the Hörmander operator decomposition technique. Then the particular solutions of this six-ordered PDE can be found in the… More >

D. L. Young^1,3, K. H. Chen², T. Y. Liu³, L. H. Shen³, C. S. Wu³

CMES-Computer Modeling in Engineering & Sciences, Vol.40, No.3, pp. 225-270, 2009, DOI:10.3970/cmes.2009.040.225

Abstract In this article, a hypersingular meshless method (HMM) is extended to solve 3D potential problems for arbitrary domains after a 2D model was successfully developed (Young et al. 2005a). The solutions are represented by a distribution of the double layer potentials instead of the single layer potentials as generally used in the conventional method of fundamental solutions (MFS). By using the desingularization technique to regularize the singularity and hypersingularity of the double layer potentials, the source points can be located exactly on the real boundary to avoid the sensitivity of locating fictitious boundary for putting the singularity outside the computational… More >