The collocation Trefftz method (CTM) proposed in [36] is employed to annular shaped domains, and new error analysis is made to yield the optimal convergence rates. This popular method is then applied to the special case: the Dirichlet problems on circular domains with circular holes, and comparisons are made with the null field method (NFM) proposed , and new interior field method (IFM) proposed in [35], to find out that both errors and condition numbers are smaller. Recently, for circular domains with circular holes, the null fields method (NFM) is proposed by Chen and his groups. In NFM, the fundamental solutions (FS) with the source nodes Q outside of the solution domains are used in the Green formulas, and the FS are replaced by their series expansions. The Fourier expansions of the known or the unknown Dirichlet and Neumann boundary conditions on the circular boundaries are chosen, so that the explicit discrete equations can be easily obtained by means of orthogonality of Fourier functions. The NFM has been applied to elliptic equations and eigenvalue problems in circular domains with multiple holes, reported in many papers; here we cite those for Laplace¡¦s equation only (see [18, 19, 20]). For the boundary integral equation (BIE) of the first kind, the trigonometric functions are used in Arnold [4, 5], and error analysis is made for infinite smooth solutions, to derive the exponential convergence rates. In Cheng¡¦s Dissertation [21, 22], for BIE of the first kind, the source nodes are located outside of the solution domain, the linear combination of fundamental solutions are used, and error analysis is made only for circular domains. This fact implies that not only can the CTM be applied to arbitrary domains, but also a better numerical performance is provided. Since the algorithms of the CTM is simple and its programming is easy, the CTM is strongly recommended to replace the NFM for circular domains with circular holes in engineering problems.

In this thesis, the boundary errors are defined for the NFM to explore the convergence rates, and the condition numbers are derived for simple cases to explore numerical stability. The optimal convergence (or exponential) rates are discovered numerically. This thesis is also devoted to seek better choice of locations for the field nodes of the FS expansions. It is found that the location of field nodes Q does not affect much on convergence rates, but do have influence on stability. Let £_ denote the distance of Q to ∂S. The larger £_ is chosen, the worse the instability of the NFM occurs. As a result, £_ = 0 (i.e., Q ∈ ∂S) is the best for stability. However, when £_ > 0, the errors are slightly smaller. Therefore, small £_ is a favorable choice for both high accuracy and good stability. This new discovery enhances the proper application of the NFM. However, even for the Dirichlet problem of Laplace¡¦s equation, when the logarithmic capacity (transfinite diameter) C_£F = 1, the solutions may not exist, or not unique if existing, to cause a singularity of the discrete algebraic equations. The problem with C_£F = 1 in the BEM is called the degenerate scale problems. The original explicit algebraic equations do not satisfy the conservative law, and may fall into the degenerate scale problem discussed in Chen et al. [15, 14, 16], Christiansen [35] and Tomlinson [42]. An analysis is explored in this thesis for the degenerate scale problem of the NFM. In this thesis, the new conservative schemes are derived, where an equation between two unknown variables must satisfy, so that one of them is removed from the unknowns, to yield the conservative schemes. The conservative schemes always bypasses the degenerate scale problem; but it causes a severe instability. To restore the good stability, the overdetermined system and truncated singular value decomposition (TSVD) are proposed. Moreover, the overdetermined system is more advantageous due to simpler algorithms and the slightly better performance in error and stability. More importantly, such numerical techniques can also be used, to deal with the degenerate scale problems of the original NFM in [15, 14, 16]. For the boundary integral equation (BIE) of the first kind, the trigonometric functions are used in Arnold [3], and error analysis is made for infinite smooth solutions, to derive the exponential convergence rates. In Cheng¡¦s Ph. Dissertation [18], for BIE of the first kind the source nodes are located outside of the solution domain, the linear combination of fundamental solutions are used, error analysis is made only for circular domains. So far it seems to exist no error analysis for the new NFM of Chen, which is one of the goal of this thesis. First, the solution of the NFM is equivalent to that of the Galerkin method involving the trapezoidal rule, and the renovated analysis can be found from the finite element theory. In this thesis, the error boundary are derived for the Dirichlet, the Neumann problems and its mixed types. For certain regularity of the solutions, the optimal convergence rates are derived under certain circumstances. Numerical experiments are carried out, to support the error made.