• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.1376-1380
• /
• 2004

본 연구에서는 홍수범람지도 제작을 위해 사용되는 수치모형의 검증을 수행하고자 한다. 검증할 수치모형은 스위스의 Beffa에 의해 개발된 FLUMEN(FLUvial Modeling ENgine)으로서 수심 적분된 2차원 비선형 천수방정식에 불규칙 심각망을 이용한 유한체적법(finite volume method)이 적용된 수치모형이며 스위스, 독일, 오스트리아 등에서 홍수범람해석에 사용된 바 있는 모델이다. FLUMEN 모형의 검증을 위하여 범람해석시 가장 중요한 문제인 이동경계조건(moving boundary condition)을 포함하고 있는 원형섬에서의 고립파에 처오름높이를 계산하여 수리모형실험 결과와 비교한다. 수리모형실험은 미국 육군 공병단 산하의 해안공학연구소(CERC, Coastal Engineering Research Center)에서 수행되었으며(Liu 등, 1995) 수조의 중앙에 높이 0.625m, 해저지름 7.2m, 경사각 $14.04^{\circ}$의 원형섬이 위치한다. 본 연구 결과, 수치해석으로 계산된 섬에서의 실제 처오름높이와 입사파의 파고의 비(R/A)는 수리모형실험의 결과와 어느 정도 일치하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2012.05a
• /
• pp.711-711
• /
• 2012

이 연구에서는 수심적분된 이차원 비정상 수치모형인 KU-RLMS 모형의 잠김/드러남 기법을 개선하고, 기존의 WAD 기법과 새로 도입한 기법의 정확도 평가를 수행하였다. 기존의 WAD 기법은 수학적으로는 다소 불완전하지만 수치적으로는 손쉬운 방법으로써, 각 시간 단계에서 잠긴 격자 또는 드러난 격자를 시험하고, 각 격자의 경계에서 플럭스에 대한 개폐 조건을 적용하는 방법을 사용하였다. 개선된 기법은 WAD 경계에서의 질량보존을 고려하는 방법이다. 두 가지 잠김/드러남 처리 과정의 정확도 비교는 포물형 수조에 대한 해석해와 수치모형의 결과를 비교하는 방법을 사용하였다. 기존에 WAD 기법은 수치해와 해석해의 위상차가 발생하는 것을 확인할 수 있으며, 진폭은 조금씩 감소하는 현상이 나타났지만, 개선된 기법을 사용할 경우 위상차와 진폭감소 현상이 현저히 개선됨을 확인할 수 있었다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2015.05a
• /
• pp.324-328
• /
• 2015

본 연구에서는 수심적분 2차원 수치모형인 Nays2D를 이용하여 수공구조물(낙차공)이 사주의 거동과 하천의 평면변화에 미치는 영향을 분석하였다. 길이 20m, 저수로 폭이 0.48m, 좌우안 홍수터 폭이 0.26m인 복단면 하도에서 낙차공에 의한 사주의 변화를 파악하기 위한 실험결과와 수치모의 결과를 비교하여 그 적용성을 검증하였다. 최심 하상고의 변화와 교호사주의 발생을 모의하였으며, 실험결과와 잘 일치하였다. 하폭이 상대적으로 넓고 수심이 얕은 복렬사주가 발생하는 조건에서 사주의 거동을 파악하기 위한 수치실험을 수행하였으며, 상류에서 발생한 사주의 형상이 낙차공 하류의 하도변화에 영향을 미치는 것을 파악하였다. 낙차공 상류와 하류에서 교호사주가 발달하였으며, 낙차공 상류에서 발달한 사주의 형상을 반영하여 낙차공 하류에서 교호사주가 형성되었다. 사주의 파장은 낙차공 상류의 파장보다 낙차공하류에서 짧고, 사주의 파고는 낙차공 하류에서 파고보다 상대적으로 높다. 또한 사주의 이동속도는 낙차공 상류보다 하류에서 빠른 특성을 보여주고 있으며, 낙차공에 의하여 흐름이 강하게 교란되면서 하류에서 사주가 형성되는 영향을 주기 때문으로 판단된다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.4
• /
• pp.223-234
• /
• 1991

In this paper, the harbor tranquility problem is analysed by an improved boundary integral equation method. The effect of the diffracted wave transformation induced by the breakwater and structures located at a harbor mouth is considered. Partial reflection concept is also used to consider energy dissipation effects. The present model is verified by comparing the results of the model for rectangular and semi-circular harbors with the analytic solutions. they show a reasonable agreement. Also the wave height distribution of the HUPO harbor computed using the present model agree well to those from the previous hydraulic model tests. It also shows a good agreement with the results from the time-dependent mild slope equation model.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.465-474
• /
• 1999

A recently developed numerical method based on a volume integral formulation is applied to calculate the accurate stress intensity factors at the crack tips in unbounded isotropic solids in the presence of multiple anisotropic inclusions and cracks subject to external loads. In this paper, a detailed analysis of the stress intensity factors are carried out for an unbounded isotropic matrix containing an orthotropic cylindrical inclusion and a crack. The accuracy and effectiveness of the new method are examined through comparison with results obtained from analytical method and finite element method using ANSYS. It is demonstrated that this new method is very accurate and effective for solving plane elastostatic problems in unbounded solids containing anisotropic inclusions and cracks.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.7 no.4
• /
• pp.679-683
• /
• 2006

Numerical integration-based program which simulates motion of pile driven by vibratory pile driver was developed for predicting rate of penetration of pile. Rate of penetration of pile calculated from developed program was compared with those of field test. As pile penetration depth increases, the difference between predicted rate of penetration and measured rate of penetration decreases. It was concluded that the reason for large difference between the predicted value and the measured value at shallower depths was attributed to decrease of vertical compressive force caused from relatively larger flexural and torsional motion of sheet pile.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.4 s.70
• /
• pp.385-393
• /
• 2005

This paper established the dynamic model of a flexible Timoshenko beam capable of geometrical nonlinearities subject to large overall motions by using the finite element method. Equations of motion are derived by using Hamilton principle and are formulated in terms of finite elements in which the nonlinear constraint equations are adjoined to the system using Lagrange multipliers. The Newmark direct integration method and the Newton-Raphson iteration are employed here for the numerical study which is to demonstrate the efficiency of the proposed formulation.

• Water for future
• /
• v.23 no.4
• /
• pp.487-497
• /
• 1990

The qualitative and quantitative prediction for the dispersion of thermal discharge from nuclear / fossil power plant, steel works etc. has significant roles for the cooling system. Design and environmental management. In this study, the several important physical properties for the behavior of a thermal discharge with strong turbulent and buoyant effects are described. The comparative evaluation between MIT and PDS models is carried out, which have the different model structures. In general, MIT and PDS models are commonly used to calculate the thermal discharge behavior with considering the ambient current and the angle of jet in an unstratified water body. The simulated results by these models have great discrepancies due to the different assumptions in modling.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.4
• /
• pp.9-20
• /
• 1999

A higher order panel method based on B-spline representation for both the geometry and the velocity potential is developed for the solution of the flow around two-dimensional lifting bodies. The self-influence functions due to the normal dipole and the source are separated into the singular and nonsingular parts, and then the former is integrated analytically whereas the latter is integrated using Gaussian quadrature. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution around lifting foils with much fewer panels than existing low order panel methods.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.