• Journal of Ocean Engineering and Technology
• /
• v.32 no.6
• /
• pp.447-457
• /
• 2018

In this study, a two-dimensional fully nonlinear transient wave numerical tank was developed using a desingularized indirect boundary integral equation method. The desingularized indirect boundary integral equation method is simpler and faster than the conventional boundary element method because special treatment is not required to compute the boundary integral. Numerical simulations were carried out in the time domain using the fourth order Runge-Kutta method. A mixed Eulerian-Lagrangian approach was adapted to reconstruct the free surface at each time step. A numerical damping zone was used to minimize the reflective wave in the downstream region. The interpolating method of a Gaussian radial basis function-type artificial neural network was used to calculate the gradient of the free surface elevation without element connectivity. The desingularized indirect boundary integral equation using an isolated point source and radial basis function has no need for information about the element connectivity and is a meshless method that is numerically more flexible. In order to validate the accuracy of the numerical wave tank based on the desingularized indirect boundary integral equation method and meshless technique, several numerical simulations were carried out. First, a comparison with numerical results according to the type of desingularized source was carried out and confirmed that continuous line sources can be replaced by simply isolated sources. In addition, a propagation simulation of a $2^{nd}$-order Stokes wave was carried out and compared with an analytical solution. Finally, simulations of propagating waves in shallow water and propagating waves over a submerged bar were also carried and compared with published data.

• Proceedings of the KIEE Conference
• /
• 1998.11a
• /
• pp.88-90
• /
• 1998

In this paper, the indirect BIEM(boundary integral equation method) is adopted to analyze 3-D eddy currents in cover plate of power transformer. In indirect BIEM, the equivalent magnetic surface charge density and the eqivalent magnanetic surface current density are the unknowns. Using triangular constant elements, the integral equations are discretized into boundary element equations of minimum order. Eddy currents are obtained in terms of euqivalent magnetic surface sources. And the locad overheating can be predicted using the eddy currents distribution in cover plate of power transformer.

• Journal of KSNVE
• /
• v.7 no.6
• /
• pp.975-984
• /
• 1997

A direct boundary element method (DBEM) is developed for thin bodies whose surfaces are rigid or compliant. The Helmholtz integral equation and its normal derivative integral equation are adoped simultaneously to calculate the pressure on both sides of the thin body, instead of the jump values across it, to account for the different surface conditions of each side. Unlike the usual assumption, the normal velocity is assumed to be discontinuous across the thin body. In this approach, only the neutral surface of the thin body has to be discretized. The method is validated by comparison with analytic and/or numerical results for acoustic scattering and radiation from several surface conditions of the thin body; the surfaces are rigid when stationary or vibrating, and part of the interior surface is lined with a sound-absoring material.

• Proceedings of the KIEE Conference
• /
• 1998.07a
• /
• pp.113-115
• /
• 1998

This paper presents the indirect boundary integral equation method(BIEM) to analyze 3D moving conductor problem. Instead of an artificial upwind algothm, the proposed method uses a fundamental Green's function which is a particular solution of diffusion equation. Therefore, this method yields a stable and accurate solution regardless of the Peclet number. The indirect BIEM is compared with 3D upwind FEM for a numerical model which has analytic solutions.

• 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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2167-2173
• /
• 1996

An approximate weight function technique using the indirect boundary integral equation has been presented for the analysis of stress intensity foactors(SIFs) of a thick pipe. One-term boundary integral was introduced to represent the crack surface displacement field for the displacement based weight function technique. An explicit closed-form SIF solution applicable to symmetric cracked pipes without any modification of the solution including both circumferential and radial cracks has been derived. The necessary information in the analysis is two or three reference SIFs. In most cases the SIF solution were in good agreement with those available in the literature.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.4
• /
• pp.981-989
• /
• 1995

A new weight function approach to determine SIF(stress intensity factor) using single-layer potential has been presented. The crack surface displacement field was represented by one boundary integral term whose kernel was modified from Kelvin's fundamental solution. The proposed method enables the calculation of SIF using only one SIF solution without any modification for the crack geometries symmetric in two-dimensional plane such as a center crack in a plate with or without an internal hole, double edge cracks, circumferential crack or radial cracks in a pipe. The application procedure to those crack problems is very simple and straightforward with only one SIF solution. The necessary information in the analysis is two reference SIFs. The analysis results using present closed-form solution were in good agreement with those of the literature.

• International Journal of Aeronautical and Space Sciences
• /
• v.3 no.2
• /
• pp.76-81
• /
• 2002

A numerical simulation of wake behavior behind three-dimensional wings in ground effect is done using an indirect boundary element method (Panel Method). An integral equation is obtained by applying Green's 2nd Identity on all surfaces of the flow domain. The AIC is constructed by imposing the no penetration condition on solid surfaces, and the Kutta at the wing's trailing edge. The ground effect is included using an image method. At each time step, a row of wake panels from wings' trailing edge are convected downstream following the force-free condition. The roll-up of wake vortices behind wings in close proximity is simulated.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1997.04a
• /
• pp.128-135
• /
• 1997

The indirect boundary integral equation is formulated to analyze the behavior of a cavity in a multilayered half-plane subjected to earthquake waves. This formulation uses the fundamental solutions that are numerically calculated by the generalized transmission and reflection coefficient method. The free surface of the cavity without external excitation influences the behavior of the half-plane. Consequently this analysis adds the consideration of scattering-field into the analysis and the total motion field of the cavity is decomposed into the free-field and scattering-field motions. The free-field motion is obtained from the modification of the transmission and reflection coefficient method. The scattering-field motion is calculated is calculated by the indirect boundary value problem which has the ficticious boundaries and sources. In this study, P wave, SV wave, SH wave, and Rayleigh wave are analyzed respectively.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.790-795
• /
• 2002

The object of this research is limited to the reduction of compression process noise only among the main sources of compressor noise such as motor noise, compression process noise, and valve port flow noise. Thus the research is focused on the wave motion rather than the particle motion of sound wave travels. A muffler is a commonly used device to reduce the compression process noise, generated by the pressure pulsations caused by the cyclic compression process. In this research, the acoustic characteristics of the muffler are analyzed by using the normal gradient integral equation proposed by Wu and Wan. Moreover, a commercial code SYSNOISE developed by indirect variational boundary integral equation is also used to validate the results. For the noise reduction, the topology optimization technique using a genetic algorithm is used. The number, size and position of the muffler holes are considered as design variables. Compared with original design, the optimized design has very improved acoustic characteristics. Both numerical and experimental analyses are used to evaluate new design.