• Journal of the Society of Naval Architects of Korea
• /
• v.56 no.1
• /
• pp.11-22
• /
• 2019

In this paper, a Rankine source method is applied and validated to analyze the hydrodynamic response of a three-dimensional floating structure in the frequency domain. The boundary value problems for radiation and diffraction problem are solved by using a desingularized indirect boundary integral equation method (DIBIEM). The DIBIEM is simpler and faster than conventional methods based on the numerical surface integration of Green's function because the singularities of Green's function are located outside of fluid regions. In case of floating structure with complex geometry, it is difficult to desingularize the singularities of Green's function consistently. Therefore a mixed approach is carried out in this study. The mixed approach is partially desingularized except singularities of the body. Wave drift loads are calculated by the middle-field formulation method that is mathematically simple and has fast convergence. In order to validate the accuracy of the developed program, various numerical simulations are carried out and these results are analyzed and compared with previously published calculations and experiments.

• The Bulletin of The Korean Astronomical Society
• /
• v.44 no.1
• /
• pp.46.1-46.1
• /
• 2019

We present an accurate and efficient method to calculate the gravitational potential of an isolated system in three-dimensional Cartesian and cylindrical coordinates subject to vacuum (open) boundary conditions. Our method consists of two parts: an interior solver and a boundary solver. The interior solver adopts an eigenfunction expansion method together with a tridiagonal matrix solver to solve the Poisson equation subject to the zero boundary condition. The boundary solver employs James's method to calculate the boundary potential due to the screening charges required to keep the zero boundary condition for the interior solver. A full computation of gravitational potential requires running the interior solver twice and the boundary solver once. We develop a method to compute the discrete Green's function in cylindrical coordinates, which is an integral part of the James algorithm to maintain second-order accuracy. We implement our method in the {\tt Athena++} magnetohydrodynamics code, and perform various tests to check that our solver is second-order accurate and exhibits good parallel performance.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.14 no.4
• /
• pp.303-309
• /
• 2004

Absorptive material arrangement method for effective interior noise control is proposed. Sound field with arbitrary boundary condition is formulated by Kirchhoff-Helmholtz integral equation. A simple example such as a rectangular cavity will present physical meaning between changing boundary condition and control of sound field. The effect of changing boundary condition is expressed in modal admittance. From this formulation. an admittance map is presented. The admittance map is the figure to represent position where absorptive material is attached for effective interior noise control. The admittance map can be assigned to each resonant frequency. There. however, may be common area of those maps. Then, frequency robust arrangement of absorptive material in noise control will be presented.

• Structural Engineering and Mechanics
• /
• v.56 no.4
• /
• pp.605-623
• /
• 2015

A formulation of the boundary element method (BEM) based on Kirchhoff's hypothesis to analyse stiffened plates composed by beams and slabs with different materials is proposed. The stiffened plate is modelled by a zoned plate, where different values of thickness, Poisson ration and Young's modulus can be defined for each sub-region. The proposed integral representations can be used to analyze the coupled stretching-bending problem, where the membrane effects are taken into account, or to analyze the bending and stretching problems separately. To solve the domain integrals of the integral representation of in-plane displacements, the beams and slabs domains are discretized into cells where the displacements have to be approximated. As the beams cells nodes are adopted coincident to the elements nodes, new independent values arise only in the slabs domain. Some numerical examples are presented and compared to a wellknown finite element code to show the accuracy of the proposed model.

• Journal of the Korean Mathematical Society
• /
• v.41 no.2
• /
• pp.345-368
• /
• 2004

A preconditioning algorithm is developed in this paper for the iterative solution of the linear system of equations resulting from the p-version boundary element approximation of the three dimensional integral equation with hypersingular operators. The preconditioner is derived by first making the nodal and side basis functions locally orthogonal to the element internal bases, and then by decoupling the nodal and side bases from the internal bases. Its implementation consists of solving a global problem on the wire-basket and a series of local problems defined on a single element. Moreover, the condition number of the preconditioned system is shown to be of order $O((1＋ln/p)^{7})$. This technique can be applied to discretization with triangular elements and with general basis functions.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1702-1707
• /
• 2003

Absorptive material arrangement method for effective interior noise control is proposed. Sound field with arbitrary boundary condition is formulated by Kirchhoff-Helmholtz integral equation. A simple example such as a rectangular cavity will present physical meaning between changing boundary condition and control of sound field. The effect of changing boundary condition is expressed in modal admittance. From this formulation, an admittance map is presented. The admittance map is the figure to represent position where absorptive material is attached. The admittance map can be assigned to each resonant frequency. There, however, may be common area of those maps. Then, frequency robust arrangement of absorptive material in noise control will be presented.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.8 no.3
• /
• pp.232-240
• /
• 1984

The plane elastostatic boundary value problem with the sharp V-notched singularity is formulated by a contour integral method for determining numerically the stress intensity factors. The integral formula is based on Somigliana type of reciprocal work in terms of displacement and traction vectors on the plate boundary. The characteristic singular solutions can be identified on the basis of traction free boundary conditions of two radial notch edges. Two numerical example examples are treated in detail; a symmetric mode-I type of notched plate with various interior angles and a mixed mode type of cantilever subjected to end shear.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.5
• /
• pp.635-644
• /
• 1986

A Compuressible turbulent boundary layer effect of the high temperature, accelerating gas flow through the De-Laval nozzle on combustion chamber pressure is numerically investigated. For this purpose, the coupled momentum integral equation and energy integral equation are solved by the Bartz method, and 1/7 power law for both the turbulent boundary layer velocity distribution and temperature distribution is assumed. As far as the boundary layer thicknesses are concerned, we can obtain reasonable solutions even if relatively simple approximations to the skin friction coefficient and stanton number have been used. The effects of nozzle wall cooling and/or mass flow rate on the boundary layer thicknesses and the combustion chamber pressure are studied. Specifically, negative displacement thickness is appeared as the ratio of the nozzle wall temperature to the stagnation temperature of the free stream decreases, and, consequently, it makes the combustion chamber pressure low.

• Journal of Ocean Engineering and Technology
• /
• v.36 no.4
• /
• pp.232-242
• /
• 2022

This study presents a method for level ice-structure interaction analysis to estimate the fatigue damage of arctic structures by applying plate theory to the behavior of level ice. The boundary element method (BEM), which incurs a lower computational cost than the finite element method (FEM), was introduced to solve the plate bending problem. The BEM formulation was performed by applying the BEM to plate theory. Finally, to check the validity of the proposed method, the BEM results and FEM results obtained using the ABAQUS commercial software were compared. The response results of the BEM analysis agreed well with those of the FEM analysis. Based on the results of the analysis, the BEM approach is considered to be very powerful in level ice-structure interaction analysis for estimating level ice-induced fatigue damage. Further work is being conducted to perform level ice fracture analysis based on the stress field calculated using the boundary element method.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.08a
• /
• pp.39.4-39.4
• /
• 2008