• Proceedings of the Korea Water Resources Association Conference
• /
• 2016.05a
• /
• pp.470-474
• /
• 2016

Madsen et al. (2002)이 제안한 일차원 고차 Boussinesq 방정식에 대하여 불연속갤러킨 유한요소법(Discontinuous Galerkin Finite Element Method)을 적용하였다. 연속적인 Boussinesq 방정식에서 각 요소경계에 불연속을 허용할 수 있도록 공간차분하고, 시간방향으로 4차 Runge-Kutta 시간적분법, 각 요소사이에는 Lax-Friedrichs 수치흐름률을 사용하였다. 계산영역의 양쪽에 불필요한 파랑의 반사를 억제하도록 흡수층을 설치하였으며, 영역 내부에서 조파할 수 있도록 하였다. Luth et al.(1994)의 수중잠제 실험에 적용하여 관측값과 잘 일치함을 확인하였다.

• Proceedings of the SAREK Conference
• /
• 2006.06a
• /
• pp.116-121
• /
• 2006

Numerical simulations for dendritic growth of crystals are conducted in this study by the level set method. The effect of order of difference is tested for reinitialization error in simple problems and authors founded in case of 1st order of difference that very fine grids have to be used to minimize the error and higher order of difference is desirable to minimize the reinitialization error The 2nd and 4th order Runge-Kutta scheme in time and 3rd and 5th order of WENO schemes with Godunov scheme are applied for space discretization. Numerical results are compared with the analytical theory, phase-field method and other researcher's level set method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.2
• /
• pp.87-95
• /
• 2018

This paper introduces a new explicit spectral element method for the simulation of SH-waves in semi-infinite domains. To simulate the wave motion in unbounded domains, it is necessary to reduce the infinite extent to a finite computational domain of interest. To prevent the wave reflection from the trunctated boundaries, perfectly matched layer(PML) wave-absorbing boundary is introduced. The forward problem for simulating SH-waves in PML-truncated domains can be formulated as second-order PDEs. The second-order semi-discrete form of the governing PDEs is constructed by using a mixed spectral elements with Legendre-gauss-Lobatto quadrature method, which results in a diagonalized mass matrix. Then the second-order semi-discrete form is transformed to a first-order, whose solutions are calculated by the fourth-order Runge-Kutta method. Numerical examples showed that solutions of SH-wave in the two-dimensional analysis domain resulted in stable and accurate, and reflections from truncated boundaries could be reduced by using PML boundaries. Elastic wave propagation analysis using explicit time integration method may be apt for solving larger domain problems such as three-dimensional elastic wave problem more efficiently.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.4
• /
• pp.413-420
• /
• 1994

Dynamics of the river plume is a very complicated non-linear problem with the free boundary changing in time and space. Mixing with the ambient water through the boundary makes the problem more complicated. In this paper we reduced 3-dimensional problem into 1-dimensional one by using the integral analysis method. Basic equations have been integrated over the lateral and vertical variations. For these integrations we adopted the well-established assumption that the flow-axis component of plume velocity and the density difference of the plume with the ambient water have Gaussian distributions in directions which are perpendicular to the flow-axis of the plume. We also used the result of our previous study on the lateral spreading velocity of the plume derived under the same assumption. And entrainment was included as a mixing process. The resultant 1-dimensional equations were solved by Runge-Kutta numerical method. Consequently, comparatively easy method of numerical analysis is presented for the 3-dimensional river plume. The method can also be used for the analysis of the thermal plume of cooling water of power plants.

• Journal of Korea Water Resources Association
• /
• v.41 no.3
• /
• pp.243-249
• /
• 2008

Continuous stirred tank reactor(CSTR) model which can be applied to control and management of the surface flow wetland is developed to simulate the water quality in this research. The model solution is obtained from the optimization model using the least-squares and 4th-order Runge-Kutta methods. The model is applied to simulate BOD and TSS in the wetland database of U.S. EPA, in which the hydraulic and water quality data are enough and the number of pond is just one for simple analysis of running results. The model is tested in two different cases, one constant volume case and another constant volume and flow rate case considering only reaction term, mass flux term and both reaction and mass flux terms respectively. It is found that the model simulates the real water quality very well with both reaction and mass flux terms rather than only reaction term and the settling velocity of TSS becomes $0.3{\sim}0.4\;m/d$. The model can be applied in wetlands treatment efficiently.

• Journal of the Society of Naval Architects of Korea
• /
• v.57 no.5
• /
• pp.271-277
• /
• 2020

Recoverability and vulnerability of navy ships under underwater explosion are critical verification factors in the acquisition phase of navy ships. This paper aims to establish numerical analysis techniques for the underwater explosion of navy ships. Doubly Asymptotic Approach (DAA) Equation of Motion (EOM) of primary shock wave and secondary bubble pulse proposed by Geers-Hunter was introduced. Assuming a non-compressive fluid, reference solution of the DAA EOM of Geers-Hunter using Runge-Kutta method was derived for the secondary bubble pulse phase with an assumed charge conditions. Convergence analyses to determine fluid element size were performed, suggesting that the minimum fluid element size for underwater explosion analysis was 0.1 m. The spherical and cylindrical fluid domains were found to be appropriate for the underwater explosion analyses from the fluid domain shape study. Because the element size of 0.1 m was too small to be applied to the actual navy ships, a very slender beam with the square solid section was selected for the study of fluid domain existence effect. The two underwater explosion models with/without fluid domain provided very similar results in terms of the displacement and stress processes.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.13 no.3
• /
• pp.259-259
• /
• 2003

Genetic algorithms using a Process of genetic evolution of an organism are appropriate for hard problems that have not been solved by any deterministic method. Up to now, the controller design method has been made with the frequency dependent specification but the design method with the time specification has gotten little progress. In this paper, we study the controller design to satisfy the performance of a plant using the generalized Manabe standard form. When dealing with a controller design in the case of two parameter configurations, there are some situations that neither a known pseudo inverse technique nor the inverse method can be applicable. In this case, we propose two methods of designing a controller by the gradient algorithm and the new pseudo inverse method so that the desired closed polynomials are either equalized to or approximated to the designed polynomial. Design methods of the proposed controller are implemented in Java.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.15 no.3
• /
• pp.1-9
• /
• 2011

The purpose of this study is to examine the dynamic characteristics of low, medium and high speed Maglev trains and guideways through dynamic interaction analysis. The coupled dynamic equations of motion for a vehicle of 10-dof and the associated guideway girders are developed by superposing vibration modes of the girder itself. The controller used in the UTM-01 Maglev vehicle is adopted to control the air gap between the bogie and guideway in this study. The effect of roughness, the guideway deflection-ratio and vehicle speed on the dynamic response of the maglev vehicle and guideway are then investigated using the 4th Runge-Kutta method. From the numerical simulation, it is found that the air gap increases with an increase of vehicle speed and the roughness condition. In particular, the dynamic magnification factor of the guideway girder is small at low and medium speeds, but the factor is noticeable at super-high speeds.