• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.3
• /
• pp.71-79
• /
• 1992

In order to obtain the wave-making resistance of a ship, so-called the Neumann-Kelvin problem is solved numerically. For computing the Havelock source, which is the Green's function of the problem, we adopted the methods given by Newman(1987) for the term representing the local disturbance, and Baar and Price(1988) for the wave disturbance, respectively. In the numerical code we developed, the source strength is assumed as bilinear on each panel and continuous throughout the hull surface. The wave-making resistance is calculated using the algorithm of de Sendagorta and erases(1988), which makes use of the wave amplitude far downstream. The Wigley hull was chosen for the sample calculation, and our results showed a good agreement with other existing experimental and numerical results.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.12 no.1
• /
• pp.92-98
• /
• 2001

When analyzing the scatting problem of multi-layered planar structures using closed-form Green's function, one of the main difficulties is that the numerical integrations for the evaluation of diagonal matrix elements converge slowly and are not so stable. Accordingly, even when the integration fur the singularity of type $e^{-jkr}/{\gamma}$, corresponding to the source dipole itself, is performed using such a method, this difficulty persists in the integration corresponding to the finite number of complex images. In order to resolve this difficulty, a new technique based upon the Gaussian quadrature in polar coordinates for the evaluation of the two-dimensional generalized exponential integral is presented. Stability of the algorithm and convergence is discussed. Performance is demonstrated for the example of a microstrip patch antenna.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.5
• /
• pp.267-275
• /
• 2022

Using a nonconforming mesh in enrichment methods results in several numerical issues induced by discontinuities and singularities found within the solution spaces, including the computational overhead during integration. In this study, we present a novel enrichment technique based on the selective expansion technique of moment fitting (Düster and Allix, 2020). In particular, two modifications are proposed to address the inefficiency during the integration process. First, a feedforward artificial neural network is introduced to correlate the implicit functions and integration moments. Through numerical examples, it is shown that the efficiency of the method is greatly improved when compared with existing expansion techniques, whereas the solution accuracy is maintained. Additionally, the finite element and domain representation grids are separated, which in turn improves the solution accuracy even for coarse mesh conditions.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.