• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.2
• /
• pp.139-149
• /
• 1994

In this paper, a finite element technique is applied to the prediction of the wave resonance phenomena in harbors. The mild-slope equation is used with a partial reflection boundary condition introduced to model the energy dissipating effects on the solid boundary. For an efficient modeling of the radiation condition at infinity, a new infinite element is developed. The shape function of the infinite element is derived from the asymptotic behavior of the first kind of the Hankel's function in the analytical boundary series solutions. For the computational efficiency, the system matrices of the element are constructed by performing the relevant integrations in the infinite direction analytically. Comparisons with the results from experiments and other solution methods show that the present model gives fairly good results. Numerical experiments are also carried out to determine the proper distance to the infinite elements from the mouth of the halter, which directly affect the accuracy and efficiency of the solution.

• Journal of Korea Water Resources Association
• /
• v.47 no.11
• /
• pp.1061-1066
• /
• 2014

A stable second-order finite volume method was proposed to predict sediment transport under rapidly varied flow conditions such as transcritical flow. For the use under unsteady flow conditions, a sediment transport model was coupled with shallow water equations. HLLC approximate Riemann solver based on a monotone upstream-centered schemes for conservation laws (MUSCL) reconstruction was used for the computation of the flux terms. From the comparisons of dam break flow experiments on erodible beds in one- and two-dimensional channels, good agreements were obtained when proper parameters were provided. Lastly, dam surface erosion problem by overtopped water was simulated. Overall, the numerical solutions showed reasonable results, which demonstrated that the proposed numerical scheme could provide stable and physical results in the cases of subcritical and supercritical flow conditions.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2017.05a
• /
• pp.324-324
• /
• 2017

기후변화로 인하여 강우특성이 변화하여 기존 치수시설물의 설계능력을 초과하는 홍수사상이 빈발하고 있다. 특히, 도심지내의 하천변 저지대는 대표적인 침수취약지역으로 강우사상이 설계빈도를 상회하는 경우 침수피해를 입을 수 있는 지역이다. 기후변화로 인하여 강우특성이 변화하는 상황에서 기존 빗물펌프장의 치수능력 증대 또는 신규시설의 설치가 필요하지만, 부지확보의 어려움으로 신규사업의 진행이 어려움을 겪고 있는 상황이다. 이와 같은 문제점에 대비하여 기존 빗물펌프장 설계단계에서 고려하지 않고 있는 요소의 추가적인 검토가 필요하다. 빗물펌프장의 유수지 용량 설계단계에서는 유입되는 우수의 수문곡선과 펌프의 최대토출량을 비교하여 잔류하게되는 우량수문곡선을 적분하여 설계용량을 산정한다. 그러나, 빗물펌프장의 실제 운영단계에서는 유수지 내의 수위와 하천의 수위를 고려하여 단계적으로 펌프를 가동하기 때문에 설계단계의 가정과 차이가 발생하게 된다. 또한, 빗물펌프장은 넓은 부지에 유입구와 흡수정을 설치함에 있어서 배치에 따른 유수지내의 흐름특성을 고려하지 않고 있다. 흡수정 설계단계에서는 흡수정으로 유입되는 흐름 특성이 토출성능에도 영향을 미칠 수 있으므로 이에 대한 분석이 필요한 상황이다. 본 연구에서는 빗물펌프장의 유수지 및 흡수정 형상과 펌프운영 기준을 반영하여 2차원 수치해석을 수행하고 유수지내의 흐름특성을 재현하였다. 수치모의결과, 유수지 형상에 따라서 흐름특성이 상이하고 흡수정 인근의 흐름특성이 유수지와 흡수정 배치에 영향을 받음을 확인할 수 있었다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.2A
• /
• pp.155-162
• /
• 2009

This paper deals with the strongest beams with the solid regular polygon cross-section, whose volumes are always held constant. The differential equation of the elastic deflection curve of such beam subjected to the concentrated and trapezoidal distributed loads are derived and solved numerically. The Runge-Kutta method and shooting method are used to integrate the differential equation and to determine the unknown initial boundary condition of the given beam. In the numerical examples, the simple beams are considered as the end constraint and also, the linear, parabolic and sinusoidal tapers are considered as the shape function of cross sectional depth. As the numerical results, the configurations, i.e. section ratios, of the strongest beams are determined by reading the section ratios from the numerical data related with the static behaviors, under which static maximum behaviors become to be minimum.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.3B
• /
• pp.261-270
• /
• 2008

Even though there is a great deal of progress in a numerical method of high caliber like SPH, it is very rarely deployed in a water resources community. Despite the great stride in computing environment, depth averaged approach like a nonlinear shallow equation is still efficient tool for flood routing in large watershed, but it can give some misleading information like the inundation height of flood. In this rationale, we numerically simulate the flow into the dry channel, dry channel with an obstacle triggered by the collapse of a two dimensional water column using SPH (Smoothed Particle Hydrodynamics) in order to boost the application of numerical method of high caliber like SPH in a water resources community. As a most severe test of the robustness of SPH, we also carry out the simulation of the flow through a clearance into the wet channel driven by the rapid removal of a water gate. As a hydrodynamic model, we used the Navier-Stokes equation, a numerical integration of which was carried out using SPH. To verify the validity of newly proposed numerical model, we compare the numerically simulated flow with the others in the literature mainly from VOF and MAC, and hydraulic experiments by Martin and Moyce (1952), Koshizuka et al. (1995) and Janosi et al. (2004). It was shown that agreements between the numerical results in this study and hydraulic experiments are remarkable.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.14 no.4
• /
• pp.435-444
• /
• 2001

Even todays, accurate and efficient algorithms for the large deformation analysis of elastoplastic frame structures lack due to the complexities of kinematics, material nonlinearities and numerical methods to cater for. The author suggests appropriate beam element based upon the incremental formulation from the 3D rod theory where Cauchy stress and engineering strain are variables to incorporate plasticity equations so that objectivity may be satisfied. A rectum mapping methods which can integrate and satisfy yield criteria efficiently is suggested and a continuation method which has global convergency and quadratic speed is developed as well. leading-unloading example problems are tested and the ideas are proved to be valuable.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.321-327
• /
• 2007

This study presents a new gridless finite difference method for solving elastic crack problems. The method constructs the Taylor expansion based on the MLS(Moving Least Squares) method and effectively calculates the approximation and its derivatives without differentiation process. Since no connectivity between nodes is required, the modeling of discontinuity embedded in the domain is very convenient and discontinuity effect due to crack is naturally implemented in the construction of difference equations. Direct discretization of the governing partial differential equations makes solution process faster than other numerical schemes using numerical integration. Numerical results for mode I and II crack problems demonstrates that the proposed method accurately and efficiently evaluates the stress intensity factors.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.2
• /
• pp.57-69
• /
• 2000

A higher order panel method based on representation for both the geometry and the velocity potential is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then the singular part is integrated analytically whereas the non-singular part 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 bodies with much fewer panels than existing low order panel methods.

• Journal of the Korean Institute of Intelligent Systems
• /
• v.5 no.2
• /
• pp.86-96
• /
• 1995

In this paper the fuzzy extension for the classical engineering mechanics problems is studied. The governing differential equation is derived for the buckling loads of the columns with uncertain mediums: the their own weight and the flexural rigidity. The columns with one typical end constraint(hinged1 clarnped/free) and the other finite rotational spring with fuzzy constant are considered in numerical examples. The vertex method is used to evaluate the fuzzy functions. The Runge-Kutta method and Determinant Search method are used to solve the differential equation and determine the buckling loads, respectively. The membership functions of the buckling load are calculated. The index of fuzziness to quantitatively describe the propagation of fuzziness is defined. According to the fuzziness of governing factors, the varlation of index of fuzziness for buckling load is investigated, and the sensitivity for the end constraints is analyzed.