• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.4
• /
• pp.58-65
• /
• 1992

In this paper, wave resistance for a submerged body is calculated by a higher order panel method. The Neumann-Kelvin problem is solved by the source or normal dipole distribution method. The body surface is represented by a bicubic B-spline and the singularity strengths are approximated by a bilinear form. The results calculated by the higher order panel method are compared with those by the lowest order panel method developed by Hess & Smith. The convergence rate of the higher order panel method is much better than the lowest order panel method. But the wave resistance calculated by the higher order panel method still shows discrepancy with an analytic solution at low Froude number like that by the lowest order panel method.

• Journal of Ocean Engineering and Technology
• /
• v.14 no.3
• /
• pp.41-45
• /
• 2000

In solving the unsteady potential flow problem of the ship in waves with the panel method, in general one can consider the basic flow as the free stream or double body solution. For the double body solution, the body boundary condition has the 2nd derivatives of the velocity potential. Low order panel methods are known to suffer from the significant error in the 2nd derivatives computed at the body surface. This paper analyzes the numerical error in the 2nd derivatives for a 2-D cylinder and a 3-D sphere problem, and an extrapolation method to obtain the correct derivatives on the body surface is suggested.

• Journal of the Society of Naval Architects of Korea
• /
• v.57 no.2
• /
• pp.114-123
• /
• 2020

A high-order potential-based panel method based on Green's theorem, with piecewise-linear dipole strength on triangular panels, is formulated for the analysis of potential flow around a three-dimensional wing. Previous low-order panel methods adopt square panels with piecewise-constant dipole strength, which results in inherent errors. Square panels can not represent a high curvature lifting body, such as propellers, since the four vertices of the square panel do not locate at the same flat plane. Moreover the piecewise-constant dipole strength induces inevitable errors due to the steps in dipole strength between adjacent panels. In this paper a high-order panel method is formulated to improve accuracy by adopting a piecewise linear dipole strength on triangular panels. Firstly, the square panels are replaced by triangular panels in order to increase the geometric accuracy in representing the shape of the object with large curvature. Next, the step difference of the dipole strength between adjacent panels is removed by adopting piecewise-linear dipole strength on the triangular panels. The calculated results by the present method is compared with analytical ones for simple non-lifting geometries, such as ellipsoid. The results for an elliptic wing with zero thickness at finite angle of attack are compared with Jordan's results. The comparison shows reasonable agrements for the both lifting and non-lifting bodies.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.1
• /
• pp.8-15
• /
• 2002

A two-dimensional higher order panel method using B-splines has been developed to overcome the disadvantages of the low order panel method and to obtain more accurate solution. The sources and the normal dipoles are distributed on both the body and the free surface. Instead of applying the upwind finite difference schemes to satisfy the linearized free surface and the radiation condition, the derivatives of the basis functions of the B-splines are directly applied to the linearized free surface condition. Numerical damping in the Dawson's method are avoided in the Present computations. In order to validate the present method, numerical computations are carried out for a submerged cylinder and a two-dimensional hydrofoil steadily moving beneath a free surface. The numerical results show that fast convergence and better accuracies have been achieved by the present method.

• Journal of Ship and Ocean Technology
• /
• v.3 no.4
• /
• pp.1-14
• /
• 1999

The improved Green integral equation using the Kelvin-type Green function in known free of irregular frequencies where the integral over the inner free surface integral is removed from the integral equation, resulting in an overdetermined integral equation. The solution of the overdetermined Green integral equation is shown identical with the solution of the improved Green integral equation Using the B-spline higher order panel method, the overdetermined equation is discretized in two different ways; one of the resulting linear system is square and the other is redundant. Numerical experiments show that the solutions of both are identical. Using the present methods, the exact values and higher derivatives of the potential at any place over the wetted surface of the body can be found with much fewer panels than low order panel method.

• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.325-334
• /
• 2018

The Nonlinear dynamic response of a sandwich plate subjected to the low velocity impact is theoretically and experimentally investigated. The Hertz law between the impactor and the plate is taken into account. Using the Extended High Order Sandwich Panel Theory (EHSAPT) and the Ritz energy method, the governing equations are derived. The skins follow the Third order shear deformation theory (TSDT) that has hitherto not reported in conventional EHSAPT. Besides, the three dimensional elasticity is used for the core. The nonlinear Von Karman relations for strains of skins and the core are adopted. Time domain solution of such equations is extracted by means of the well-known fourth-order Runge-Kutta method. The effects of core-to-skin thickness ratio, initial velocity of the impactor, the impactor mass and position of the impactor are studied in detail. It is found that these parameters play significant role in the impact force and dynamic response of the sandwich plate. Finally, some low velocity impact tests have been carried out by Drop Hammer Testing Machine. The results are compared with experimental data acquired by impact testing on sandwich plates as well as the results of finite element simulation.

• Journal of applied mathematics & informatics
• /
• v.17 no.1_2_3
• /
• pp.73-91
• /
• 2005

The application of Green's function in calculation of flow characteristics around submerged and floating bodies due to a regular wave is presented. It is assumed that the fluid is homogeneous, inviscid and incompressible, the flow is irrotational and all body motions are small. Two methods based on the boundary integral equation method (BIEM) are applied to solve associated problems. The first is a low order panel method with triangular flat patches and uniform distribution of velocity potential on each panel. The second method is a high order panel method in which the kernels of the integral equations are modified to make it nonsingular and amenable to solution by the Gaussian quadrature formula. The calculations are performed on a submerged sphere and some floating spheroids of different aspect ratios. The excellent level of agreement with the analytical solutions shows that the second method is more accurate and reliable.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.2 s.140
• /
• pp.113-123
• /
• 2005

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low order panel methods.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.3
• /
• pp.31-40
• /
• 2001

Low-order panel methods are being used to design marine propellers. Since the potential value over each panel for these methods is assumed to be a constant, the accuracy of prediction is known to be limited. Therefore, a higher order boundary element method(HOBEM) has been studied to enhance the accuracy of prediction. In this paper, a HOBEM representing the body boundary surfaces and physical quantities by a 9-node Lagrangian shape function is employed to analyse the flow around marine propellers in steady potential flow. First, the numerical results for a circular wing with thickness variations are compared with Jordan's linear solution. Then, the computational results of two propellers(DTRC 4119 & DTRC 4842 propeller) are compared with the experimental and numerical results published. The pressure distribution on the surface of the propeller is also compared with experimental data.

• Journal of Ship and Ocean Technology
• /
• v.4 no.2
• /
• pp.13-23
• /
• 2000

A low-order potential based boundary element method is applied to the prediction of the flow around the cavitating propeller in steady or in unsteady inflow. For given cavitation number, the cavity shape is determined in an iterative manner until the kinematic and the dynamic boundary conditions are both satisfied on the approximate cavity boundary. In order to improve the solution behavior near the tip region, a hyperboloidal panel geometry and a modified split panel method are applied. The method is then extended to include the analysis of time-varying cavitating flows around the propeller blades via a time-step algorithm in time domain. In the method, the steady state oscillatory solution is obtained by incremental stepping in the itme domain. Finally, the present method is validated through comparison with other numerical results and experimental data.