• Transactions of the Korean Society of Mechanical Engineers
• /
• v.7 no.2
• /
• pp.175-185
• /
• 1983

When a circular cylinder near the plane fluid-interface of different viscosities is in parallel and normal motion, solutions of the Oseen equation are obtained. Classical image method with Faxen's integral form is used to satisfy the boundary conditions on the plane interface. Coefficients of drag and lift increase as a cylinder approaches to the interface. But drag-coefficients of parallel motions with viscosity-ratio less than unity are decreased slightly. They show monotonic increase with Reynolds number in case of parallel motion, but minimum values of drag coefficients in normal motion are appeared. On the other hand Stokes' solution are obtained by taking limits of low Reynolds number except the case of parallel motion with viscosity-ratio not equal to infinity.

• Proceedings of the KIEE Conference
• /
• 2001.07b
• /
• pp.687-689
• /
• 2001

In this study, BIEM(Boundary Integral Equation Method) using Green's function is applied in order to analyze 3D electric field precisely. The algorithm is developed with equivalent electric surface charge as a variable, which promises less unknown variables and higher accuracy of electric field analysis. The validity of the developed program is varified by applying it to a coaxial cylinder mode and 3-phase GIS model.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.13 no.5
• /
• pp.437-443
• /
• 1988

We calculate the distribution of the current on the strip by the incident waves on the infinite conducting strip line. The boundary equations represented as the spatial domain function become very complicated equations including convolution integral. Transformed it to the spectral domain, we have a very simple equation is composed by some algebraic multiplication of the current density function and Green's function. the acceleration of iteration procedure is achieved by Kastner's method. The result of iteration gives us the optimum value when it satisfies the iteration stop condition presented in this paper. We confirmed that the induced current density distribution on the stripline has been changed as variaties of the width.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.11 no.5
• /
• pp.323-331
• /
• 1986

By solving the current distribution, natual frequencies, and the amplitude of natural modes about the perfect conducting scatterer of revolution are presented. Singularity expansion method is applied to slove the time-domain response based on the simplified E-field integral equation by the application of the extended boundary conditin for the perfect conductin scatterer of revolution.

• Journal of KSNVE
• /
• v.7 no.2
• /
• pp.215-222
• /
• 1997

The flow and acoustic fields due to a vortex ring interaction with a rigid sphere are simulated numerically. The flow field is regarded as three-dimensional inviscid and incompressible. The vorticity is assumed to be concentrated inside the finite core of vortex filament. The vortex filament curve, described by parabolic blending curve function, is used to effectively solve the modified Biot-Savart equation. The interaction between a vortex ring and a rigid sphere using the parabolic blending curve is calculated. The trajectory of the vortex ring is obtained with several different initial positions between the ring and the sphere. The force variations acting on the sphere are calculated by using the boundary integral method. Finally, we can also obtain the acoustic signals at the far field observation positions from the force variations acting on the rigid surface. We can find that the dipole axis of the directivity patterns are rotated during the interacting phenomena.

• Interaction and multiscale mechanics
• /
• v.1 no.3
• /
• pp.329-337
• /
• 2008

This paper presents the derivation of the generalized governing equations in theory of elasticity, their weak forms and the some applications in the numerical analysis of structural mechanics. Unlike the differential equations in classical elasticity theory, the generalized equations of the equilibrium and compatibility equations presented here take the form of integral equations, and the generalized equilibrium equations contain the classical differential equations and the boundary conditions in a single equation. By using appropriate test functions, the weak forms of these generalized governing equations can be established. It can be shown that various variational principles in structural analysis are merely the special cases of these weak forms of generalized governing equations in elasticity. The present weak forms of elasticity equations extend greatly the choices of the trial functions for approximate solutions in the numerical analysis of various engineering problems. Therefore, the weak forms of generalized governing equations in elasticity provide a powerful modeling tool in the computational structural mechanics.

• Journal of Korea Water Resources Association
• /
• v.32 no.6
• /
• pp.649-655
• /
• 1999

The maximum run-up heights of periodic waves are numerically investigated in this study. Incident waves are sinusoidal and enoidal waves. The maximum run-up height of enoidal wave approaches that of sinusoidal wave as the wave length decreases, while it approaches that of solitary wave as the wave length increases. If wave height is fixed, the maximum run up heights of enoidal waves are always greater than those of sinusoidal waves but smaller than those of solitary waves.

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

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