We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.

Solving potential-flow problems using the Boundary Element Method (BEM) is a strong tradition in marine hydrodynamics. An early example of the application of BEM is by Bai & Yeung [1]. The bottleneck of the conventional BEM in terms of CPU time and computer memory arises as the number of unknowns increases. Wu & Eatock Taylor [2] suggested that the Finite Element Method (FEM) field solver is much faster than the BEM based on their comparisons in a wave making problem. In this paper, we aim to find a highly efficient method to solve fully-nonlinear wave-body interaction problems based on potential-flow theory. We compare the efficiency and the accuracy of five different methods for the potential flows in two dimensions (2D), two of which are BEM-based while the other three are field solvers. The comparisons indicate that it is beneficial to use either an accelerated matrix-free BEM, e.g. Fast Multipole Method accelerated BEM (FMM-BEM), or any field solvers whose resulting matrix are sparse. Another highlight of this paper is that an efficient numerical potential-flow method named the harmonic polynomial cell (HPC) method is developed. The flow in each cell is described by a set of harmonic polynomials. The presented procedure has approximately 4th order accuracy, while its resulting matrix is sparse similarly as the other field solvers, e.g. Finite Element Method (FEM), Finite Difference Method (FDM) and Finite Volume Method (FVM). The method is verified by a linear wave making problem for which the steady-state analytical solution is available, and the forced oscillation of a semi-submerged circular cylinder for which the frequency-domain added mass and damping coefficients are compared. The fully-nonlinear wave making problem and nonlinear propagating waves over a submerged bar are also studied for validation purposes. Only 2D cases are studied in this paper.

In numerical simulations with moving bodies, and often with complex

Abstract:Further developments and applications of the 2D harmonic polynomial cell (HPC) method proposed by Shao and Faltinsen [22] are presented. First, a local potential flow solution coupled with the HPC method and based on the domain decomposition strategy is proposed to cope with singular potential flow characteristics at sharp corners fully submerged in a fluid. The results are verified by comparing them with the analytical added mass of a double-wedge in infinite fluid. The effect of the singular potential flow is not dominant for added mass and damping, but the error is non-negligible when calculating mean wave loads using direct pressure integration. Then, the double-layer nodes technique is used to simulate a thin free shear layer shed from lifting bodies, across which the velocity potential is discontinuous. The results are verified by comparing them with analytical results for steady and unsteady lifting problems of a flat plate in infinite fluid. The latter includes the Wagner problem and the Theodorsen functions. Satisfactory agreement with other numerical results is documented for steady linear flow past a foil and beneath the free surface.

No abstract