• Proceedings of the Korea Water Resources Association Conference
• /
• 2006.05a
• /
• pp.1910-1913
• /
• 2006

To solve the interaction of incident monochromatic waves with a bottom-fixed vertical circular cylinder, a numerical analysis by boundary element method is developed using three-dimensional linear potential theory. A numerical analysis by boundary element method is based on Green's theorem and introduce to an integral equation for the fluid velocity potential around the vertical circular cylinder. These numerical results are compared with those of ManCamy and Fuchs(1954) and Williams and Mansour(2002), and it has shown good relationship with their results. This numerical analysis developed by boundary element method will be applied for various offshore structures to be constructed in coastal zones in the future.

• Journal of Navigation and Port Research
• /
• v.32 no.7
• /
• pp.543-552
• /
• 2008

The diffraction of waves by three bottom fixed vertical circular cylinders is investigated by using the boundary element method. This method has been successfully applied to the isolated vertical circular cylinder and now is used to study the interaction between waves and multiple vertical cylinders. In this paper, a numerical analysis by the boundary element method is developed by the linear potential theory. The numerical analysis by the boundary element method is based on Green's second theorem and introduced to an integral equation for the fluid velocity potential around the vertical circular cylinders. To verify this method, the results obtained in present study are compared with the results computed by the multiple scattering method. The results of the comparisons show strong agreement. Also in this paper, several numerical examples are given to illustrate the effects of various parameters on the wave exciting force such are the separation distance, the wave number and the incident wave angle. This numerical computation method might be used broadly for the design of various offshore structures to be constructed in the future.

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

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

Aerodynamic analysis of NACA wings moving with a constant speed over guideways are performed using an indirect boundary element method (potential-based panel method). An integral equation is obtained by applying Green's theorem on all surfaces of the fluid domain. The surfaces over the wing and the guideways are discretized as rectangular panel elements. Constant strength singularities are distributed over the panel elements. The viscous shear layer behind the wing is represented by constant strength dipoles. The unknown strengths of potentials are determined by inverting the aerodynamic influence coefficient matrices constructed by using the no penetration conditions on the surfaces and the Kutta condition at the trailing edge of the wing. The aerodynamic characteristics for the wings flying over nonplanar ground surfaces are investigated for several ground heights.

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

• Proceedings of the KSRS Conference
• /
• 2004.10a
• /
• pp.586-589
• /
• 2004

A numerical algorithm using the Method of Moments (MoM) is introduced to compute precisely the scattering matrices of very thin deciduous leaves in this paper. At first, a dyadic Green's function was formulated and an integral equation for a volumetric current distribution in a lossy dielectric body. Then, the MoM was applied to the scattering problem with a specific technique to handle the numerical poles. The accuracy of the numerical technique was verified by examining the technique with various ways, and used to examine the validity regions of the classical analytical models.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2001.05a
• /
• pp.202-207
• /
• 2001

The added resistance of a ship advancing in waves can be split into the resistance due to the radiation wave and the resistance due to the diffraction wave. In this study, the former has been calculated by a method based on Maruo's formula. The latter must be calculated by other methods. Ship motion is calculated by the usual strip method. The amplitude of two dimensional far-field waves is calculated using the improved Green integral equation. The present numerical method can be used for the estimation of the added resistance due to the radiation wave since the present numerical result is much smaller than other existing numerical results considered to be overestimated.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.05a
• /
• pp.185-191
• /
• 2003

A numerical study on the hydrodynamic properties of a floating flexible breakwater consisting of triple vertical porous membrane structures attached to a floating rigid pontoon restrained by moorings is carried out in the context of two-dimensional linear wave-flexible body interaction theory. The tensions in the triple membranes are achieved by hanging a clump weight from its lower ends. The clump weight is also restrained properly by moorings. The dynamic behavior of the breakwater was described through an appropriate Green function, and the fluid multi-domains are incorporated into the boundary integral equation. Numerical results are presented which illustrate the effects of the various wave and structural parameters on the efficiency of the breakwater as a barrier to wave action. It is found that the wave reflection and transmission properties of the structures depends strongly on the membrane length taking major fraction of water column, the magnitude of tensions on membrane achieving by the clump weight, proper mooring types and stiffness, the permeability on the membrane dissipating wave energy.

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

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.65-72
• /
• 2001

This paper demonstrates an explicit-implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for an existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory. Darcy's law is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full Newton-Raphson method is used for solving the nonlinear governing equations.