Leevan Ling¹, Tomoya Takeuchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 45-54, 2008, DOI:10.3970/cmes.2008.029.045

Abstract The method of fundamental solutions is coupled with the boundary control technique to solve the Cauchy problems of the Laplace Equations. The main idea of the proposed method is to solve a sequence of direct problems instead of solving the inverse problem directly. In particular, we use a boundary control technique to obtain an approximation of the missing Dirichlet boundary data; the Tikhonov regularization technique and the L-curve method are employed to achieve such goal stably. Once the boundary data on the whole boundary are known, the numerical solution to the Cauchy problem can be obtained by solving a direct… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.2, pp. 77-94, 2008, DOI:10.3970/cmes.2008.028.077

Abstract The method of fundamental solutions (MFS) is a truly meshless numerical method widely used in the elliptic type boundary value problems, of which the approximate solution is expressed as a linear combination of fundamental solutions and the unknown coefficients are determined from the boundary conditions by solving a linear equations system. However, the accuracy of MFS is severely limited by its ill-conditioning of the resulting linear equations system. This paper is motivated by the works of Chen, Wu, Lee and Chen (2007) and Liu (2007a). The first paper proved an equivalent relation of the Trefftz method and MFS for circular… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.1, pp. 53-66, 2007, DOI:10.3970/cmes.2007.021.053

Abstract A newly modified Trefftz method is developed to solve the exterior and interior Dirichlet problems for two-dimensional Laplace equation, which takes the characteristic length of problem domain into account. After introducing a circular artificial boundary which is uniquely determined by the physical problem domain, we can derive a Dirichlet to Dirichlet mapping equation, which is an exact boundary condition. By truncating the Fourier series expansion one can match the physical boundary condition as accurate as one desired. Then, we use the collocation method and the Galerkin method to derive linear equations system to determine the Fourier coefficients. Here, the factor… More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.2, pp. 111-122, 2007, DOI:10.3970/cmes.2007.020.111

Abstract A highly accurate new solver is developed to deal with interior and exterior mixed-boundary value problems for two-dimensional Laplace equation, including the singular ones. To promote the present study, we introduce a circular artificial boundary which is uniquely determined by the physical problem domain, and derive a Dirichlet to Robin mapping on that artificial circle, which is an exact boundary condition described by the first kind Fredholm integral equation. As a consequence, we obtain a modified Trefftz method equipped with a characteristic length factor, ensuring that the new solver is stable because the condition number can be greatly reduced. Then,… More >

D.L. Young¹, K.H. Chen², J.T. Chen³, J.H. Kao⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.3, pp. 197-222, 2007, DOI:10.3970/cmes.2007.019.197

Abstract A boundary-type method for solving the Laplace problems using the modified method of fundamental solutions (MMFS) is proposed. The present method (MMFS) implements the singular fundamental solutions to evaluate the solutions, and it can locate the source points on the real boundary as contrasted to the conventional MFS, where a fictitious boundary is needed to avoid the singularity of diagonal term of influence matrices. The diagonal term of influence matrices for arbitrary domain can be novelly determined by relating the MFS with the indirect BEM and are also solved for circular domain analytically by using separable kernels and circulants. The… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.2, pp. 145-162, 2007, DOI:10.3970/cmes.2007.019.145

Abstract A new method is developed to solve the Dirichlet problems for the two-dimensional Laplace equation in the doubly-connected domains, namely the meshless regularized integral equations method (MRIEM), which consists of three portions: Fourier series expansion, the Fredholm integral equations, and linear equations to determine the unknown boundary conditions onartificial circles. The boundary integral equations on artificial circles are singular-free and the kernels are degenerate. When boundary-type methods are inefficient to treat the problems with complicated domains, the new method can be applicable for such problems. The new method by using the Fourier series and the Fourier coefficients can be adopted… More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.19, No.1, pp. 99-110, 2007, DOI:10.3970/cmes.2007.019.099

Abstract A new meshless regularized integral equation method (MRIEM) is developed to solve the interior and exterior Dirichlet problems for the two-dimensional Laplace equation, which consists of three parts: Fourier series expansion, the second kind Fredholm integral equation and an analytically regularized solution of the unknown boundary condition on an artificial circle. We find that the new method is powerful even for the problem with complex boundary shape and with random noise disturbing the boundary data. More >

B. Jin^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 253-262, 2004, DOI:10.3970/cmes.2004.006.253

Abstract In this paper, we propose a new numerical scheme for the solution of the Laplace and biharmonic equations subjected to noisy boundary data. The equations are discretized by the method of fundamental solutions. Since the resulting matrix equation is highly ill-conditioned, a regularized solution is obtained using the truncated singular value decomposition, with the regularization parameter given by the L-curve method. Numerical experiments show that the method is stable with respect to the noise in the data, highly accurate and computationally very efficient. More >

Chein-Shan Liu^1,2

CMC-Computers, Materials & Continua, Vol.5, No.1, pp. 31-42, 2007, DOI:10.3970/cmc.2007.005.031

Abstract By using a meshless regularized integral equation method (MRIEM), the solution of elastic torsion problem of a uniform bar with arbitrary cross-section is presented by the first kind Fredholm integral equation on an artificial circle, which just encloses the bar's cross-section. The termwise separable property of kernel function allows us to obtain the semi-analytical solutions of conjugate warping function and shear stresses. A criterion is used to select the regularized parameter according to the minimum principle of Laplace equation. Numerical examples show the effectiveness of the new method in providing very accurate numerical solutions as compared with the exact ones. More >

Weichung Yeih¹, Chein-Shan Liu², Chung-Lun Kuo³, Satya N. Atluri⁴

CMC-Computers, Materials & Continua, Vol.17, No.3, pp. 275-302, 2010, DOI:10.3970/cmc.2010.017.275

Abstract In this paper, a multiple-source-point boundary-collocation Trefftz method, with characteristic lengths being introduced in the basis functions, is proposed to solve the direct, as well as inverse Cauchy problems of the Laplace equation for a multiply connected domain. When a multiply connected domain with genus p (p>1) is considered, the conventional Trefftz method (T-Trefftz method) will fail since it allows only one source point, but the representation of solution using only one source point is impossible. We propose to relax this constraint by allowing many source points in the formulation. To set up a complete set of basis functions, we… More >