• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.1
• /
• pp.26-32
• /
• 2012

The 3-D unsteady viscous flow around an impulsively started rotating marine propeller is simulated using VIC(Vortex-In-Cell) method which is adequate to analyze the strong vortical flow around complicatedly-shaped body. The computational procedure is governed by the vorticity transport equation in Lagrangian form. In order to solve the equation, a regular grid which is independent to the shape of a body is introduced and each term of the equation is evaluated numerically on the grid by applying immersed boundary concept. In this paper, the overall algorithm including the formulation of governing equations and boundary conditions is described and some computational results are presented with discussing their physical validity.

• The Journal of Korea Institute of Information, Electronics, and Communication Technology
• /
• v.3 no.2
• /
• pp.31-38
• /
• 2010

It is difficult to design accurately the bandgap of metamaterial depending on metamaterial pattern and array configuration. In this paper, we propose a design method for the wanted bandgap frequency using any metamaterial pattern in 2 dimensional array. Metamaterial structure is consisted of periodic array. Therefore the calculation area in FDTD(finite difference time domain) method can be reduced by applying the periodic boundary condition to 2-D metamaterial array. The method for design and calculation the L and C values by using 2-D is also considered. So it can be designed more accurately and rapidly. For example, we designed metamaterial square pattern array in 5 GHz, and compared with the 1-D metamaterial pattern using analysis method in microstrip line. As a result, we found that the accuracy of this proposed method can be incresed to 14.7%.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.5 no.4
• /
• pp.546-558
• /
• 2013

The research performed in this paper was carried out to investigate the numerical analysis of the sheet cavitation on marine propeller. The method is boundary element method (BEM). Using the Green's theorem, the velocity potential is expressed as an integral equation on the surface of the propeller by hyperboloid-shaped elements. Employing the boundary conditions, the potential is determined via solving the resulting system of equations. For the case study, a DTMB4119 propeller is analyzed with and without cavitating conditions. The pressure distribution and hydrodynamic performance curves of the propellers as well as cavity thickness obtained by numerical method are calculated and compared by the experimental results. Specifically in this article cavitation changes are investigate in both the radial and chord direction. Thus, cross flow variation has been studied in the formation and growth of sheet cavitation. According to the data obtained it can be seen that there is a better agreement and less error between the numerical results gained from the present method and Fluent results than Hong Sun method. This confirms the accurate estimation of the detachment point and the cavity change in radial direction.

• Journal of Ocean Engineering and Technology
• /
• v.20 no.6 s.73
• /
• pp.67-74
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time-domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with large amplitude under the free~surface are solved in time-domain. Through the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.

• Journal of applied mathematics & informatics
• /
• v.26 no.5_6
• /
• pp.1057-1069
• /
• 2008

We consider singularly perturbed Boundary Value Problems (BVPs) for third and fourth order Ordinary Differential Equations(ODEs) of convection-diffusion type with discontinuous source term and a small positive parameter multiplying the highest derivative. Because of the type of Boundary Conditions(BCs) imposed on these equations these problems can be transformed into weakly coupled systems. In this system, the first equation does not have the small parameter but the second contains it. In this paper a computational method named as 'An asymptotic finite element method' for solving these systems is presented. In this method we first find an zero order asymptotic approximation to the solution and then the system is decoupled by replacing the first component of the solution by this approximation in the second equation. Then the second equation is independently solved by a fitted mesh Finite Element Method (FEM). Numerical experiments support our theoritical results.

• Journal of the Optical Society of Korea
• /
• v.14 no.4
• /
• pp.444-450
• /
• 2010

The traditional approach to inverting aerosol extinction makes use of the assumption of a constant LIDAR ratio in the entire Mie-LIDAR signal profile using the Fernald method. For the large uncertainty in the cloud optical depth caused by the assumed constant LIDAR ratio, an not negligible error of the retrieved aerosol extinction below the cloud will be caused in the backward integration of the Fernald method. A new algorithm to determine aerosol extinction below a cirrus cloud from Mie-LIDAR signals, based on a new cloud boundary detection method and a Mie-LIDAR signal modification method, combined with the backward integration of the Fernald method is developed. The result shows that the cloud boundary detection method is reliable, and the aerosol extinction below the cirrus cloud found by inverting from the modified signal is more efficacious than the one from the measured signal including the cloud-layer. The error due to modification is less than 10% taken in our present example.

• Journal of the Korean Society for Precision Engineering
• /
• v.25 no.4
• /
• pp.127-133
• /
• 2008

As the application of the MEMS parts increases, the structural safety of MEMS appears importantly. A lot of MEMS parts are made by a multi-grain silicon wafer, which is an orthotropic material. Moreover directions of the materials on each grain are distributed randomly. The stress analysis for the multi-grain is important factor in order to apply the MEMS parts to industrial applications. The finite element method (FEM) is commonly used by a stress analysis method but the boundary element method (BEM) is known as the result of the BEM is more accurate than that of the FEM since the fundamental solution are used. In this study, we derived the boundary integration equation for the orthotropic material by applying fundamental solutions with complex variables. The multi-region analysis procedure for the BEM and the multi-grain generation procedure by a random process technique are developed in order to apply the analysis of the multi-grain orthotropic material. The discontinuous element is used in order to remove the comer problem in the BEM. The results of the present method are compared with those of the finite element method in order to verify the present procedure.

• Journal of Power System Engineering
• /
• v.14 no.6
• /
• pp.75-82
• /
• 2010

In this study, the topics for free-surface wave simulation, nonlinear hydrodynamic force, and the critical resonance frequency of so-called ${\tau}=U{\omega}/g$=1/4 are discussed. A high-order spectral/boundary element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with forward speed under the free-surface are solved in time domain.

• The Journal of the Acoustical Society of Korea
• /
• v.14 no.2E
• /
• pp.37-42
• /
• 1995

The polar parabolic equation method (Polar PE) which introduces a series of "cascaded" boundary fitting coordinates into the parabolic equation method has been verified as a good numerical method for atmospheric sound propagation over a curved surface and hills. Polar PE is applied here to underwater sound propagation over a sea mountain assuming locally reacting boundary sea bottom and pressure release water surface for the boundary conditions. Calculations are presented for underwater propagation over a 450 m high sea mountain. Feasibility of Polar PE application for underwater sound propagation over a smooth mountain is discussed.

• Journal of Ocean Engineering and Technology
• /
• v.37 no.2
• /
• pp.49-57
• /
• 2023

Since Chappelear developed a Fourier approximation method, considerable research efforts have been made. On the other hand, Fourier approximations are unsuitable for deep water waves. The purpose of this study is to provide a Fourier approximation suitable even for deep water waves and a numerical method to determine the Fourier coefficients and the wave properties. In addition, the convergence of the solution was tested in terms of its order. This paper presents a velocity potential satisfying the Laplace equation and the bottom boundary condition (BBC) with a truncated Fourier series. Two wave profiles were derived by applying the potential to the kinematic free surface boundary condition (KFSBC) and the dynamic free surface boundary condition (DFSBC). A set of nonlinear equations was represented to determine the Fourier coefficients, which were derived so that the two profiles are identical at specified phases. The set of equations was solved using Newton's method. This study proved that there is a limit to the series order, i.e., the maximum series order is N=12, and that there is a height limitation of this method which is slightly lower than the Michell theory. The reason why the other Fourier approximations are not suitable for deep water waves is discussed.