• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.2
• /
• pp.121-129
• /
• 2008

In this paper, we present a fast analysis of multilayer microstrip fractal structure by using the fast multipole method (FMM). In the analysis, accurate spatial green's functions from the real-axis integration method(RAIM) are employed to solve the mixed potential integral equation(MPIE) with FMM algorithm. MoM's iteration and memory requirement is $O(N^2)$ in case of calculation using the green function. the problem is the unknown number N can be extremely large for calculation of large scale objects and high accuracy. To improve these problem is fast algorithm FMM. FMM use the addition theorem of green function. So, it reduce the complexity of a matrix-vector multiplication and reduce the cost of calculation to the order of $O(N^{1.5})$, The efficiency is proved from comparing calculation results of the moment method and Fast algorithm.

• Economic and Environmental Geology
• /
• v.27 no.2
• /
• pp.191-199
• /
• 1994

We have developed an algorithm based on the method of integral equations to simulate the magnetotelluric (MT) responses of three-dimensional (3-D) bodies in a layered half-space. The inhomogeneities are divided into a number of cells and are replaced by an equivalent current distribution which is approximated by pulse basis functions. A matrix equation is constructed using the electric Green's tensor function appropriate to a layered earth, and is solved for the vector current in each cell. Subsequently, scattered fields are found by integrating electric and magnetic Green's tensor functions over the scattering current About a 3-D conductive body near the earth's surface, interpretation using 2-D transverse electric modeling schemes can imply highly erratic low resistivities at depth. This is why these routines do not account for the effect of boundary charges. However, centrally located profiles across elongate 3-D prisms may be modeled accurately with a 2-D transverse magnetic algorithm, which implicitly includes boundary charges in its formulation. Multifrequency calculations show that apparent resistivity and impedance phase are really two complementary parameters. Hence, they should be treated simultaneously in broadband MT interpretation.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.10 no.3
• /
• pp.173-180
• /
• 2007

Oscillating wave amplitude in a bottom-mounted OWC chamber designed for wave energy converter is investigated by applying characteristic wave conditions in Korean coastal water. The effects of shape parameters of OWC chamber in a view of wave energy absorbing capability are analyzed. Both experimental and numerical approaches are adopted and their results are compared to optimize the shape parameters which can result in a maximum power production under given wave distribution. The experiment was carried out in a wave flume under 2-D assumption of OWC chamber. The numerical scheme employed a hybrid Green integral equation which adopts the Rankine Green function inside chamber to take account of fluctuating air pressure, while it uses the Kelvin Green function in outer domain. Air duct diameter, chamber width, and submerged depths of front skirt and back wall of chamber changes the magnitude and peak frequency of wave absorption significantly.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.4
• /
• pp.73-81
• /
• 1994

The source and source/dipole distribution methods using 3-dimensional panel method have been widely used for ship motion analysis. When these methods are used, large errors in the predicted hydrodynamic coefficients are introduced around the irregular frequencies caused by the resonance of imaginary internal flow. Therefore, the irregular frequencies need to be removed for an accurate prediction of ship motion. This paper adopts 3-dimensional translating and oscillating Green function derived by Wu. The adaptive integration method, stretching transform and stationary phase method are used for the calculation of the calculation of Green function and the integral equation is derived by distributing the Green function n ship surface and inner free-surface. The condition of zero normal velocity, that is, wall condition on inner free-surface has been successfully used for the removal of irregular frequencies in oscillating problems. The calculations are carried out for series 60($C_B=0.7$) vessel and the results are compared with those of other theoretical analyses and experiment.

• Economic and Environmental Geology
• /
• v.27 no.5
• /
• pp.491-498
• /
• 1994

The integral equation method is used for magnetotelluric (MT) modeling of a finite inhomogeneity in a two-layered earth. An integral equation relates the incident plane-wave field and the scattering currents in the three-dimensional (3-D) inhomogeneity through the electric tensor Green's function appropriate to a layered earth. This paper describes an effect of overburden and basement on MT responses of 3-D body. The effect of overburden is to reduce the detectability of target, and the reduction of detectability is more apparent for conductive overburden than for resistive one. The effect of basement, on the other hand, may enhance the anomaly due to 3-D body in the upper layer. In case of the resistive basement current perturbations about the body tend to be confined to the more conductive upper layer.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.4
• /
• pp.16-21
• /
• 2003

The very large container ships have been built recently and those ships have very small structural rigidity compared with the other conventional ships. As a result, the destruction of ship hull is occurred by the springing including to warping phenomena due to encounter waves. In this study, the solutions of hydrodynamic coefficients are obtained by solving the three dimensional source distribution method and the forward speed Green function representing a translating and pulsating source potential for infinite water depth is used to calculating the integral equation. The vessel is longitudinally divided into various sections and the added mass, wave damping and wave exciting forces of each section is calculated by integrating the dynamic pressures over the mean wetted section surface. The equations for six degree freedom of motions is obtained for each section in the frequency domain and stiffness matrix is calculated by Euler beam theory. The computations are carried out for very large ship and effects of bending and torsional ridigity on the wave frequency and angle are investigated.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.12 no.6
• /
• pp.931-939
• /
• 2001

The slot coupling characteristics was analysed in a radial-line slot antenna for its design. The previously proposed waveguide model with a periodic boundary condition on its narrow walls and periodically arranged slots on its wide wall was used. The magnetic field integral equation and two dyadic Green\`s functions for respective regions was derived and the method of moments was used. To maximize the efficiency of numerical analysis and to extract singularities, two different kinds of basis functions, the entire domain basis function and the sub-domain one, are used. In addition, the Ewald sum technique for the rectangular waveguide and the Shanks transform for the half space were used to accelerate the computation of the slowly convergent potential Green\`s functions. Simulation results expressed the effects of the various design parameters on the slot coupling.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.11 no.7
• /
• pp.1204-1210
• /
• 2000

An analysis method for obtaining transient response of electromagnetic pulse(EMP) coupled via a slot aperture into a parallel-plate waveguide, is considered. The mixed-potential integral equation(MPIE), in which the vector and scalar potential Green's functions for the unknown equivalent magnetic surface current of the aperture in parallel-plate region are expressed in closed forms derived by use of the improved complex image method, is formulated. When the method of moments(MoM) is used to solve the MPIE, the matrix-fill time is significantly reduced because of closed-form Green's functions. In order to check the validity of the present method, the numerical results obtained by the present method are compared with those obtained by the previous method. Fairly good agreements between them are observed.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.10-17
• /
• 2002

A study on the dynamic ground effect on three-dimensional wings is done using an indirect boundary element method(unsteady panel method). An integral equation is obtained by applying Green's theorem on all surfaces of the fluid domain. Constant strength dipole and source panels arc distributed on a wing's surface. The wake sheet is represented by constant strength dipoles. At each time step, a row of wake panels is assumed to be convected from the trailing edge of the wing. The tip vortex behind wings in dynamic ground effect moves outward. The amplitudes of the aerodynamic coefficients for the wings in dynamic ground effect are augmented much more comparing to the case in static ground effect.