• Journal of the Korean Society for Marine Environment & Energy
• /
• v.13 no.3
• /
• pp.174-180
• /
• 2010

The inhomogeneous Helmholtz equation is introduced for variable water depth and potential function and separation of variables are introduced for the derivation. Only harmonic wave motions are considered. The governing equation composed of the potential function for irrotational flow is directly applied to the still water level, and the inhomogeneous Helmholtz equation for variable water depth is obtained. By introducing the wave amplitude and wave phase gradient the governing equation with complex potential function is transformed into two equations of real variables. The transformed equations are the first and second-order ordinary differential equations, respectively, and can be solved in a forward marching manner when proper boundary values are supplied, i.e. the wave amplitude, the wave amplitude gradient, and the wave phase gradient at a side boundary. Simple spatially-centered finite difference numerical schemes are adopted to solve the present set of equations. The equation set is applied to two test cases, Booij’ inclined plane slope profile, and Bragg’ wavy bed profile. The present equations set is satisfactorily verified against other theories including the full linear equation, Massel's modified mild-slope equation, and Berkhoff's mild-slope equation etc.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.34 no.2
• /
• pp.479-492
• /
• 2014

In this study, was proposed a numerical scheme imposing exact solutions as the internal boundary conditions for the shallow-water flows over a discontinuous transverse structure such as a step. The HLLL approximate Riemann solver with the MUSCL was used for the test of the proposed scheme. Very good agreement was obtained between simulations and exact solutions for various problems of the shallow-water flows over a step. In addition, results by the numerical model showed good agreement with those of dam-break experiments over a step and stepped chute one. Developed model can simulate the shallow-water flows over discontinuous bottom such as a drop structure without additional rating curve or topography smoothing. Given the proper evaluations for the flow resistance by a step and the energy loss by the nappe flow in the future, could be simulated flooding and drying of the shallow-water flows over discontinuous topography such as a weir or the river road with retaining wall.

• Journal of Korea Water Resources Association
• /
• v.41 no.8
• /
• pp.773-784
• /
• 2008

The present study was carried out to examine flow properties in linear drainage channels such as road surface drainage facilities. The finite difference formulation for the varied flow analysis was solved for flow profiles in the channels. Starting the first step at the control section, the Newton-Raphson method was applied for producing numerical solutions of the equation. We considered two types of linear drainage channels, a channel with one outlet at downstream end and a channel with two outlets at both ends. Moreover, the flow analysis for various channel slopes was performed. However, we considered channels with the two outlets of slopes satisfying the condition that the both ends are the control section. The maximum of those slopes was decided from the relation between the channel slope and the location of control section. The flow of a channel with one outlet was calculated upward and downward from the control section existing in channel or upward from the control section at downstream end. The flow of a channel with two outlets at both ends were calculated for upstream and downstream channel segments divided by the water dividend, respectively and the flow analysis was completed when the water depth at the water dividend calculated from upstream end was equal to that calculated from downstream end. If the slope was larger than the critical slope, the channel with two outlets was likely to behave like the channel with one outlet. The maximum water depth was investigated and compared with that calculated additionally from the uniform flow analysis. The uniform flow analysis was likely to lead a excessive design of a drainage channel with mild slope.

• Tunnel and Underground Space
• /
• v.31 no.4
• /
• pp.270-288
• /
• 2021

We proposed a numerical method to simulate the hydro-mechanical behavior of rock fracture using a grain-based distinct element model (GBDEM) in the paper. As a part of DECOVALEX-2023 Task G, we verified the method via benchmarks with analytical solutions. DECOVALEX-2023 Task G aims to develop a numerical method to estimate the coupled thermo-hydro-mechanical processes within the crystalline rock fracture network. We represented the rock sample as a group of tetrahedral grains and calculated the interaction of the grains and their interfaces using 3DEC. The micro-parameters of the grains and interfaces were determined by a new methodology based on an equivalent continuum approach. In benchmark modeling, a single fracture embedded in the rock was examined for the effects of fracture inclination and roughness, the boundary stress condition and the applied pressure. The simulation results showed that the developed numerical model reasonably reproduced the fracture slip induced by boundary stress condition, the fracture opening induced by fluid injection, the stress distribution variation with fracture inclination, and the fracture roughness effect. In addition, the fracture displacements associated with the opening and slip showed good agreement with the analytical solutions. We expect the numerical model to be enhanced by continuing collaboration and interaction with other research teams of DECOVALEX-2023 Task G and validated in further study experiments.

• Journal of the Korean Geosynthetics Society
• /
• v.19 no.2
• /
• pp.35-46
• /
• 2020

In this study, the dynamic response characteristics of the embankment model were analyzed using shaking table experiments. Laminar shear box was used to minimize the boundary effect of the model. The ratio of the vertical length to horizontal length of the slopes were 1:1, 1:1.5, and 1:2. The sensor array which is consist of 12 accelerometers was used to measure acceleration time-histories at each location of the slope model. The dynamic response characteristics of the models were analyzed for sine wave, sinesweep wave, and artificial earthquake wave in this study. The experimental results show that the dynamic response of the embankment model is increased with the slope angle. Furthermore, the experimental setup used in this study was verified with the comparative analysis between experimental results and 1-D analytical simulation on the flat ground model.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2022.05a
• /
• pp.244-244
• /
• 2022

4차 산업혁명에 따라 유역 및 하천관리 사업, 각종 풍수해 예방사업 분야에 다양한 스마트기술이 도입되고 있으나 건설현장 침수 피해 사고는 지속적으로 발생하고 있다. 굴착공사 현장에서 발생할 수 있는 침수피해를 사전에 예측하기 위해서는 공정별로 변화하는 현장상황을 반영하여 다양한 강우 시나리오를 기반으로 침수 예측 모델링이 선행되어야한다. 따라서 본 연구에서는 2차원 동수역학 모형인 HDM-2D 모형을 기반으로 굴착공사 현장 침수피해 예측 모델을 개발하여 굴착공사현장 침수 예·경보 시스템에 탑재하고자 한다. 침수피해 예측 모델은 천수방정식을 Petrov-Galerkin stabilized scheme 으로 이산화하여 해석하는 수평 2차원 동수역학 흐름해석모델로서 수로 및 지표면 등 다양한 지형 상황에서의 물 흐름을 상세하게 해석할 수 있다. 지형자료 생성 이후 경계조건 부여를 통해 수행되며 침수발생지역의 유속, 수심, 수위를 취득할 수 있다. 배수지나 굴착공사 현장에 2차원 흐름해석을 적용하는 경우 지형의 경사나 배치가 공간에 따라 변화하므로 불연속적인 흐름을 유발하여 모의결과의 계산 오차를 검토해야 한다. 2차원 침수피해 예측 모형의 정확성을 확인하기 위해 지면 돌출부가 있는 흐름 문제와 테스트베드 대상지에 침수해석 모형을 적용하였다. 돌출부 흐름 문제의 경우 돌출부를 지나며 발생하는 유속과 수심 모의 결과를 상용모형과 비교검증 하였으며 테스트베드 대상지역에 침수피해예측모형을 적용했을때 지형 경사에 따른 흐름의 변화와 침수양상을 확인할 수 있었다.

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

하천은 여러 하천의 연결망으로 구성되어 있으며, 연결부는 절점으로 이루어져 있고, 이 절점은 하천이 합류하는 곳으로서 하천의 시스템을 구성하는 가장 중요한 요소이다. 최근에 친자연 하천 복원 등의 많은 사업을 하고 있음에도 불구하고 이에 대한 검토가 거의 이루어지지 않고 있으며, 이를 검토하고 이해하는 것이 하천을 복원하고 관리하는데 매우 중요하다. 본 연구는 본류와 지류의 하상경사 및 하상고가 다르고, 유사 유입조건이 다른 비평형 유입조건을 갖고 있으며, 하폭이 넓고 수심이 상대적으로 얕은 합류부에서 흐름 및 하상변동에 관하여 수치모의를 수행하여 그 특성을 파악하는데 그 목적이 있다. 적용 대상지역은 금강과 미호천의 합류부로서, 흐름 특성 및 하상변동을 모의하기 위하여 RMA-2와 SED2D를 이용하였다. 유량의 증가에 따라 흐름에 대한 유속이 증가하며, 전단층의 길이가 증가하고, 합류부의 우안의 정체수역에서 발생하는 와(vortex)의 크기도 증가하였다. 부유사의 농도는 실제 상황에 맞도록 비평형 유입조건을 사용하였으며, 미호천에서 유입되는 유사의 농도에 영향을 받아, 합류직후에도 유사의 농도가 높으며, 전단층을 경계로 하여 농도 띠가 하류로 유지되었으며, 합류 후에는 우안에서는 하상이 상승되었다.

• Journal of the Korean GEO-environmental Society
• /
• v.19 no.12
• /
• pp.25-32
• /
• 2018

This paper investigates the effect of rock discontinuities on a shear stress wave that is induced by earthquake or blasting and provides the result of numerical parametric studies. The numerical tests of different conditions of rock and discontinuity have been carried out after confirming that the numerical approach is valid throughout a verification analysis from which the test results were compared with a theoretical solution. In-situ stress condition was considered as a rock condition and internal friction angle and cohesive value, which are the shear strength parameters, were considered as discontinuities condition. The joint inclination angle was also taken into account as a parameter. With the various conditions of different parameters, the test results showed that a shear stress wave propagating through a mass is highly influenced by the shear strength of discontinuities and the condition of joint inclination angle as well as in-situ stress. The study results indicate that when earthquake or blasting-induced dynamic loading propagates through a jointed rock mass or a stratified soil ground the effect of in-situ stress and discontinuities including a stratum boundary should be taken into account when evaluating the dynamic effect on nearby facilities and structures.

• Journal of Biosystems Engineering
• /
• v.10 no.1
• /
• pp.24-38
• /
• 1985

경사진 밀폐 공간에서 마주 보는 두 벽면의 온도 차로 인하여 발생되는 자연 대류 현상은 여러 공학 분야에서 볼 수 있는 중요한 열전달 현상으로서, 최근 들어 평판형 태양열 집열기를 설계하려는 사람들에게 많은 관심의 대상이 되고 있다. 평판형 태양열 집열기의 경우 덮개판으로 부터의 대류 열손실을 감소시킴으로서 집열 효율을 높일 수 있을 뿐만 아니라 사용목적에 따라 소형 집열기를 제작할 수 있어 경제적으로 유리하게 될 것이다. 밀폐된 공간에서 최초에 정지 상태에 있는 얇은 유체층을 하부에서 가열시켜 주면 열팽창 현상이 일어나고, 이것에 의한 부력이 점도나 열전도도 등의 안정화 요인을 극복할 수 있을 정도로 커지면 System이 불안정하게 되어 자연 대류 현상이 수반되며 이 때문에 열전달율이 급격히 증가하게 된다. 이러한 현상의 지배 방정식은 연립 비선형 편미분 방정식으로 특수한 경계 조건외에는 일반적으로 해석적 해를 구하기가 어렵기 때문에 실험적 연구가 많이 실시되어 왔고 이들 결과의 대부분은 전반적인 열전달 특성치만을 구하는데 집중되어 왔다. 본 연구에서는 수치 해석법의 하나인 유한 차분법을 도입하여 이차원으로 가정한 경사진 평판형 밀폐 공간에서의 자연 대류 현상의 지배 방정식을 유한 차분화시켜, $$2.74{\times}10^3\leq_-Gr\leq_-2.0{\times}10^6$$, Pr=0.73, $$15^{\circ}\leq_-a\leq_-150^{\circ}$$, 종횡비는 1, 2, 3, 5, 9에 대하여 정상 상태에서의 해를 구하면서 수치적으로 실험하였다. 본 연구에서 얻어진 결론을 요약하면 다음과 같다. (1) 해석적으로 구하기 어려운 경사진 밀폐 공간에서 자연대류현상의 지배 방정식을 유한 차분법으로 해결할 수 있으며, 대류항과 확산항을 따로 고려한 유한차분법이 효과적임을 확인하였다. (2) 저온과 고온 벽면에서의 온도를 각각 균일하게 놓고 단변으로 이루어진 벽면은 완전히 절연되어 있는 경우에 대하여 수치해를 구한결과, 이전의 해석적 및 실험적 결과와 일치하였으며, 시간의 경과에 따른 온도 및 유선의 변화를 현상학적으로 관찰할 수 있었다. (3) 평균 열전달 계수에 미치는 경사각의 효과를 살펴본 결과 종횡비가 1 인 경우 경사각이 $45^{\circ}$에서, 종횡비가 2, 3, 5, 9인 경우 경사각이 $60^{\circ}$에서 각각 평균 열전달 계수 최대치가 나타났다. (4) Ra수(Rayleigh number) 가 증가될수록, 경사각에 상관없이 평균 열전달 계수도 증가되었다. 한편 Ra수 및 경사각의 변화에 따라 종횡비가 증가될수록 평균 열전달 계수는 경사각이 $90^{\circ}$인 경우를 제외하고는 감소됨을 볼 수 있었다. 경사각이 $90^{\circ}$인 경우, 평균 열전달 계수는 종횡비가 2 인 곳에서 최대치를 얻을 수 있었으며, 종횡비가 계속 증가될수록 평균 열전달 계수는 점차 감소되어짐을 볼 수 있었다.

• Proceedings of the Korean Association of Geographic Inforamtion Studies Conference
• /
• 2003.04a
• /
• pp.62-62
• /
• 2003

Landsat TM 영상을 이용, 명암차가 높은 산악 지역에 적용해왔던 알고리즘을 개선하여 비교적 명암차가 낮고 넓게 분포하는 충적층 지역의 선구조를 추출하는 알고리즘을 개발하였다. 수치지형모델(OEM)에 대하여 Local Enhancement 를 이용해서 평탄한 지역을 선정하여 이로부터 충적층을 추출하였다. Zevenbergen & Thorno's Method를 3×3 moving windowing을 통해서 최대 경사방향과 경사를 구해서 충적층을 지나는 선구조 요소를 추출하고 다시 Hough 변환을 이용해서 1차 선구조를 추출하였다 이를 이용하여 충적층의 직각방향의 지형단면의 경사를 유추해서 spline 보간법을 이용해 단면의 최저점을 구하고 이 구해진 점들을 다시 Hough 변환을 이용해서 최종 선구조를 추출하였다. 본 연구에서 사용한 알고리즘은 기존 알고리즘에서 사용된 소창문보다 훨씬 큰 충적층이 분포하는 지역의 지형 경사가 수렴하는 부분에 선구조가 뚜렷이 나타남을 볼 수 있다. 최대경사방향과 경사를 구해서 얻어진 1 차선구조와 최종 선구조에서 선구조 방향이 다소 차이를 보인다. 1 차 선구조의 수직방향 지형단면의 자료를 이용함에 있어, 지형 단면의 시작정과 끝지점을 임의적으로 결정하는 것이 아니라, 충적층을 가로질러 최고점까지 또는 다음 충적층이 나을 때까지의 자료를 이용해서 보간법을 사용하였고, 충적층의 넓이에 따라 보간할 자료량의 차이에 의한 오차가 발생할 수 있다. 넓은 충적층에서 선구조가 잘 추출되는 반면에 좁은 충적층이 분포하거나 계곡에 해당하는 지역에l서는 경사수렴부와 일치하지 않는 선구조가 추출되었다. 이는 향후 계속적으로 연구해서 보완되어야 할 것으로 사료된다.페클 잡영 제거 정도에 있어 다른 필터들과 큰 차이가 없지만 경계선보존지수는 다른 필터들에 비하여 가장 우수함을 확인할 수 있었다.rbon 탐식효율을 조사한 결과 B, D 및 E 분획에서 유의적인 효과를 나타내었다. 이상의 결과를 종합해볼 때, ${\beta}$-glucan은 고용량일 때 직접적으로 또는 $IFN-{\gamma}$ 존재시에는 저용량에서도 복강 큰 포식세로를 활성화시킬 뿐 아니라, 탐식효율도 높임으로써 면역기능을 증진 시키는 것으로 나타났고, 그 효과는 crude ${\beta}$-glucan의 추출조건에 따라 달라지는 것을 알 수 있었다.eveloped. Design concepts and control methods of a new crane will be introduced in this paper.and momentum balance was applied to the fluid field of bundle. while the movement of′ individual material was taken into account. The constitutive model relating the surface force and the deformation of bundle was introduced by considering a representative prodedure that stands for the bundle movement. Then a fundamental equations system could be simplified considering a steady state of the process. On the basi