“…, J E , J E presents the number of elements). The coordinates x m (m = 1, 2, 3), pressure p, normal velocity v n of one point on an element can be interpolated by those values of nodal points of this element with a set of shape functions in a local coordinate system ξ 1 oξ 2 : (8) where N l (ξ 1 , ξ 2 ) (l = 1 ∼ L) is the shape function, L is the number of nodal points on the element, −1 ≤ ξ 1 , ξ 2 ≤ 1 are the local coordinates. For node I, the boundary integral equation is discretized as…”