On the Application of BEM to Boundary-fitted Conformal Grid Generation and Associated Singularity Problems

博士 === 國立臺灣大學 === 土木工程學研究所 === 93 === In this thesis, the process of a conformal grid generation by the boundary element method and the complex mapping technique is presented in details. To elaborate the grid generation method, the associated singularity problems of mathematical si...

Full description

Bibliographic Details
Main Authors: Jennhann Wang, 王鄭翰
Other Authors: T. K. Tsay
Format: Others
Language:en_US
Published: 2004
Online Access:http://ndltd.ncl.edu.tw/handle/4t7ka9
Description
Summary:博士 === 國立臺灣大學 === 土木工程學研究所 === 93 === In this thesis, the process of a conformal grid generation by the boundary element method and the complex mapping technique is presented in details. To elaborate the grid generation method, the associated singularity problems of mathematical singularity and geometrical singularity, which are related to the BEM and the complex mapping, are studied comprehensively. It is found that the improvement in eliminating the numerical boundary layer leads to the ability of present grid generation method producing grids near the boundary in $10^{-9}$ scale when the boundary element method is applied. The problem of the grids'' overlapping at a sharp corner in a conformal grid system is dealt with a complex mapping technique, which can manage such a geometrical singularity adequately. However, there are still some limitations on the application of the complex mapping. In this dissertation, discussions are made to illustrate the practical way in generating a conformal grid system that was once thought difficult to construct. In the end of this thesis, the boundary-fitted curvilinear grid system is applied to calculate a case of computational fluid dynamics in connection with storm surge. In the present study, the fundamental researches on Boundary Element Method (BEM) are performed at first. Through the studies on the mathematical singularity and the geometrical singularity, the resolution of the conformal grids near the boundary and the overlapping of grids at a sharp corner can be improved. In Chapter ef{chap2}, detailed comparisons of the contour method and the direct method to evaluate the accuracy of singular boundary integrals are made. In contrast to conventional numerical integration methods which suffer from the numerical boundary layer, the so-called boundary layer can be shown to vanish when the contour and the direct methods are applied. The singularity problems, including mathematical singularities and geometrical singularities, are discussed respectively in Chapter ef{chap2} and Chapter ef{chap3}. Furthermore, in Chapter ef{chap3}, the geometrical singularities are overcome by applying a complex mapping technique. Several benchmarks are proposed to offer detailed discussions. In Chapter ef{chap4}, based on the BEM and the complex mapping method, a powerful numerical grid generation method is presented. By using the correct evaluation technique in the integrals of the BEM, the new grid generation method can avoid the effect of mathematical singularity mentioned in Chapter ef{chap2}. This enables us to produce conformal grids near the boundary. Moreover, by using a complex mapping technique in Chapter ef{chap3}, the methods serve to avert the overlapping of the grids at a sharp corner, which occurred in the past studies. The new grid generation can produce a boundary-fitted orthogonal grids. In addition, it maintains conformal properties, either angles or length ratios, so as to map a domain of Laplace equation to another domain of still the Laplace equation. This makes the governing equation most concise after the conformal transformation. Nevertheless, some defects of the grid generation arising from the mathematical singularity and the geometrical singularity are found. The grid generation method appears to be destined to have multi-value and the branch-cut lines'' problems, thus limiting its application in natural geometry. In this study, angle constraints are presented to serve as a criterion in checking the overlapping of grids. It helps in the trial and error operation during the step-by-step complex mappings. In addition, various numerical examples are employed to examine the validity of this grid generation system. Finally, in Chapter ef{chap5}, further applications of the orthogonal grid system are made on the storm surge around Taiwan; the trend of the results is the same as the observed data. Possible ways for improvement are discussed.