• 대한조선학회논문집
• /
• 제29권4호
• /
• pp.58-65
• /
• 1992

본 논문에서는 물체표면을 겹 3차 B-Spline 방법에 의하여 표현하고, 특이점세기를 겹선형으로 근사하는 고차판요소법을 이용하여 몰수체에 대한 조파저항을 계산하였다. 여기서 Neumann-kelvin 문제는 쏘오스분포법과 법선다이폴분포법에 의하여 해석되었다. 고차 판요소법에 의하여 계산된 결과는 Hess & Smith가 사용한 최저차 판요소법 결과와 비교하였으며, 고차 판요소법의 수렴도는 보통 판요소법보다 훨씬 좋은 것으로 나타났다. 그러나, 고차 판요소법에 의하여 계산된 조파저항도 최저차 판요소법에 의한 것과 마찬가지로 낮은 Froude 수에서는 해석해와의 차이를 보이고 있다.

• 한국해양공학회지
• /
• 제14권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.

• 대한조선학회논문집
• /
• 제57권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.

• 대한조선학회논문집
• /
• 제39권1호
• /
• pp.8-15
• /
• 2002

B-spline을 이용하여 물체의 형상과 포텐셜을 표현함으로써, 저차경계요소법의 단점들을 극복하고 수치계산의 정도를 높이기 위한 고차 패널법을 개발하였다. 물체표면과 자유표면에 법선 다이폴과 쏘스를 분포시켰으며, 자유표면 및 방사조건을 만족시키기 위해 상류차분식을 사용하는 대신 B-spline의 기저함수를 미분하여 선형화된 자유표면 경계조건에 직접 적용하였다. 이 방법을 적용함으로써 Dawson 방법에서 문제가 되었던 수치감쇠 문제를 피할 수 있었다. 수치계산 프로그램을 검증하기 위해 2차원 원주주위의 유동계산과 날개면 주위의 유동해석을 수행하였으며, B-spline 기저 고차패널법에 의한 수치계산 결과가 저차패널법에 비해 빠른 수렴성과 정확성을 보였고 계산에 필요한 패널 수가 현저히 줄어드는 대단히 만족스러운 결과를 얻었다.

• Journal of Ship and Ocean Technology
• /
• 제3권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
• /
• 제68권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
• /
• 제17권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.

• 대한조선학회논문집
• /
• 제42권2호
• /
• 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.

• 대한조선학회논문집
• /
• 제38권3호
• /
• pp.31-40
• /
• 2001

일반적으로 특이점들을 분포하여 물체주위 유동장을 해석할 때 저차패널법(Low-Order Panel Method)이 유용하게 사용되어져 왔다. 저차패널법과는 다르게 적분방정식의 계산정도와 해의 수렴성을 높일 수 있는 방법으로서, 물체경계의 형상을 2차 이상의 고차곡면요소로 표현하고 각 곡면요소 내에서 물리량의 변화를 동일한 차수를 갖는 고차경계요소법을 사용할 수 있다. 본 연구에서는 물체표면의 곡면요소와 물리량의 변화를 9절점 라그란지안(Lagrangian) 형상함수를 사용하여 정상 포텐셜유동 중에 작동하는 프로펠러 주위의 유동을 해석하였다. 개발된 프로그램을 프로펠러의 해석에 앞서 원형 날개에 대하여 Jordan의 선형 해석해와 비교하였으며, 두께의 영향을 관찰하였다. 프로펠러의 해석을 위해 DTRC 4119 프로펠러와 DTRC 4842 프로펠러에 적용하여 실험치와 계산치를 비교 검토하였다.

• Journal of Ship and Ocean Technology
• /
• 제4권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.