• Journal of Korea Water Resources Association
• /
• v.43 no.10
• /
• pp.921-932
• /
• 2010

In this study, the adjustment processes of the disturbed channels by sand or gravel mining were investigated by a two dimensional numerical model in the generalized coordinate system. As a numerical scheme, the CIP (cubic interpolated pseudoparticle method) method was used to calculate the advection term in the flow field and central difference method was used to the diffusion term in it. The pit of the channel was partially filled with sediment at the toe of the pit upstream. As time increased, the headcut erosion upstream in the pit was decreased due to the sediment inflow. The almost inflow sediment upstream was trapped into the pit and the sediment deposit wedge migrated downstream in the pit with the steep submerged angle of repose. The numerical model was reproduced well the evolution processes of the channel. The mining pit migrated with speed as the channel was steep, and the numerical results were in overall agreement with the experimental results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.2B
• /
• pp.165-173
• /
• 2011

In this study, FESWMS FST2DH model was used to analyze the change of flow characteristics after making the riparian forest. The additional flow resistance is calculated based on the drag-force concept acting on each tree and the lateral momentum transfer between planted and non-planted zone could be satisfactorily reproduced by parabolic turbulence model in this depth-averaged 2-D numerical model. For model validation, the simulated velocities were compared with the measured data, showing good agreement in both tree density cases of experiments. The previous method using a proper Manning's n coefficient gives reasonable solutions only to evaluate the conveyance, but the calculated approach velocity at each tree was different from realistic value. The proposed procedure could be widely used to evaluate hydraulic effects of riparian trees in practical engineering.

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

최근 우리나라는 집중호우로 인한 토석류의 발생이 현저하게 증가하고 있으며, 2002년 태풍 루사와 2003년 매미 그리고 2006년 7월 집중호우 등으로 인한 피해의 규모는 여러 조사와 문헌에서도 확인할 수 있다. 2011년 7월, 서울 우면산 일대에서는 집중호우로 인한 토석류가 발생하였으며, 16명의 인명손실을 포함한 큰 피해가 있었다. 우면산 토석류는 십여 개 지역에서 동시다발적으로 발생하였으며, 토석류 발생 유역 특성과 토석류의 유동특성을 분석하기 위한 현장 조사가 진행되었다. 조사한 자료에 따르면 우면산 일대의 일일 강우량은 서초 기상측정소를 기준으로 최대 24시간 누적 324mm 그리고 시간당 최대 68.5mm/hr를 기록하였다. 상업용 소프트웨어인 FLO-2D는 유사농도의 함수로서 점성(viscous)응력, 항복(yield)응력, 난류 및 분산(dispersive) 응력항을 포함하는 2차 유변학(quadratic rheology) 모델을 기본으로 사용하여 이류(mudflow)와 토석류(debris flow)를 모의할 수 있다. FLO-2D는 흐름의 운동량 및 에너지 보존을 고려하여 격자와 시간에 관계없이 유동심도, 속도, 압력을 예측할 수 있으며, 격자 기반의 모델로서 GIS 및 기타 응용 프로그램들과 연동이 쉽다는 장점이 있다. 그러나 하상침식에 의해 유발된 토석류의 체적 증가는 고려 할 수 없으므로 토석류의 전파 및 퇴적영역에서의 토석류 모의에만 사용할 수 있는 단점이 있다. 이 연구의 목적은 FLO-2D 소프트웨어를 이용하여 우면산에서의 토석류 현상을 재현하는 것이다. 우면산 일대에서 발생한 토석류 중 서초구 방배동의 래미안 아파트 부근에서 발생한 토석류에 대하여 수치지도(DEM)와 현장조사를 통해 얻은 지형자료, 해당 지역의 강우량 및 지질 특성 자료 등을 토대로 FLO-2D 모델을 적용하여 토석류의 흐름특성을 검토한다. 토석류 유동 및 퇴적에 대한 가용한 현장관측 자료와의 비교 분석을 통하여 토석류 특성 값을 산정하고, 모델의 적용성을 검증한다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.6 no.2
• /
• pp.93-101
• /
• 1986

The behavior of a plane buoyant jet discharged vertically upward into a stagnant uniform environment is analyzed by continuity, momentum transport equation by numerical scheme. The governing equations are solved by finite difference method employing stream function and vorticity transport and Prandtl's turbulent model. Results for centerline velocities and temperatures, temperature distribution and flow pattern in receiving environment due to buoyant jet in the range of discharge densimetric Froude number of 4 to 32 show good agreement with published data. Spreading rate and dispersion ratio, which are required in integral type analysis of whole range of buoyant jet and have not been obtained yet, are derived in terms of discharge densimetric Froude number and vertical distance from source.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.46 no.1
• /
• pp.41-51
• /
• 2018

Downwash from helicopter rotor blades and external winds from various maneuvering make an unguided rocket change its trajectory and range. For the prediction of the trajectory and range, it is essential to consider the downwash effect. In this study, an algorithm was developed to calculate 6-Degree-Of-Freedom(6 DOF) forces and moments exerting on the rocket, and total flight trajectory of a 2.75-inch unguided rocket in a helicopter downwash flow field. Using Actuator Disk Model(ADM) analysis result, the algorithm could analyze the entire trajectory in various initial launch condition such as launch angle, launch velocity, and external wind. The algorithm that considered the interference between a fuselage and external winds could predict the trajectory change more precisely than inflow model analysis. Using the developed algorithm, the attitude and trajectory change mechanism by the downwash effect were investigated analyzing the effective angle of attack change and characteristics of pitching stability of the unguided rocket. Also, the trajectory and range changes were analyzed by considering the downwash effect with external winds. As a result, it was concluded that the key factors of the rocket range change were downwash area and magnitude which effect on the rocket, and the secondary factors were the dynamic pressure of the rocket and the interference between a fuselage and external winds. In tailwind case which was much influential on the range characteristics than other wind cases, the range of the rocket rose as increasing the tailwind velocity. However, there was a limit that the range of the rocket did not increase more than the specific tailwind velocity.

• Water for future
• /
• v.29 no.1
• /
• pp.163-179
• /
• 1996

Flood flow characteristics of the Han River (from Goan to Indo Bridge) are analyzed using 1-dimensional St. Venant equations. NETWORK, a finite difference model, is used to calculate each term (local acceleration term, convective acceleration term, pressure force term, gravity force term, and friction force term) of the momentum equation and to analyze the flow characteristics. By the result of the study, as the general characteristics of flow in a channel that acceleration terms are very small and the other three terms are much greater is presented, three terms(pressure force term, gravity force term, friction force term) are to be main terms which decide the characteristics of flow. Specially in this region the acceleration term is noted so large that it cannot be ignored according to the shape of hydrograph and the region.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2023.05a
• /
• pp.130-130
• /
• 2023

제트류는 다양한 크기와 운동량의 에디가 복잡하게 혼합되어 이루어져 있으며, 이를 정확하게 모델링하고 이해하기 위해서는 제트류의 다양한 특성들을 잘 반영하여 연구를 수행해야 한다. 다양한 연구 수행 방법 중 수치해석 방법은 상대적으로 공간 및 시간적 비용이 적게 들어서 널리 사용되고 있다. 이러한 수치해석 방법에는 DNS(Direct Numerical Simulation), LES(Large Eddy Simulation), RANS(Reynolds Averaged Navier Stokes) 등이 있으며, 그중 LES는 난류 모델링을 사용하는 RANS 방법에 비해 더욱 정확한 흐름 모델링을 제공하는 장점이 있다. 이러한 LES는 대규모 에디는 직접 해석하면서, 일정 크기 이하의 에디는 모델링을 사용해 해석하는 것이 특징이다. 하지만, LES를 사용하기 위해서는 적절한 그리드 크기를 결정하는 것이 중요하며, 이는 모델의 정확성과 연산 비용에 큰 영향을 미친다. 하지만, 여전히 적절한 그리드 크기를 결정하는 것은 어려운 문제이다. 이러한 LES 모델링을 사용할 때 적절한 그리드 크기를 결정하기 위해서는 정확한 시간 평균 속도 변동을 연구하는 것이 앞서 선행되어야 한다. 따라서, 본 연구에서는 기계학습 기반 접근 방식을 사용하여 난류 제트 내 시간 평균 속도 변동을 예측하는 연구를 진행하였다. 즉, 난류 제트 역학을 이해하는 데 중요한 파라미터인 시간 평균 유속을 이용하여 시간 평균 속도 변동을 예측하는 데 초점을 맞추었다. 모델의 성능은 평균 제곱 오차와 R-제곱 등 다양한 지표를 사용하여 평가되었다.

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

The object of this study is to develop the model that solves the numerically difficult problems in hydraulic engineering and to demonstrate the applicability of this model by means of various test examples, such as, verification in the gradually varied unsteady condition, three steady flow problems with the change of bottom slope with exact solution, and frictional bed with analytical solution. The governing equation of this model is the integral form of the Saint-Venant equation satisfying the conservation laws, and finite volume method with the Riemann solver is used. The evaluation of the mass and momentum flux with the HLL Riemann approximate solver is executed. MUSCL-Hancock scheme is used to achieve the second order accuracy in space and time. This study introduce the new and simple technique to discretize the source terms of gravity and hydrostatic pressure force due to longitudinal width variation for the balance of quantity between nonlinear flux and source terms. The results show that the developed model's implementation is accurate, robust and highly stable in various flow conditions with source terms, and this model is reliable for one-dimensional applications in hydraulic engineering.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.34-34
• /
• 2021

본 연구에서는 격자기반 분포형 수문모델링을 통해 하천갈수량을 추정하고자 한다. 분포형 수문모형은 단방향흐름 알고리즘에 의한 토양 물수지식을 기반으로 개발되었으며 운동파(kinematic wave) 이론을 적용하여 지표 및 지표하 유출을 모의한다. 또한, 격자별로 수문학적 물수지요소인 차단량, 증발산량, 침투 및 침루량, 지하수충전량 등을 계산하며, 댐·보 방류량을 해당 지점 격자의 물수지에 적용할 수 있도록 개발하였다. 본 모형은 2개의 다목적댐과 3개의 다기능보가 위치한 금강유역(9,645.5 km²)에 적용하였으며, 유역 면적과 하천 유속을 고려하여 1 km × 1 km 격자를 구성하고 10분 간격으로 2013년부터 2020년까지 수문모의를 진행하였다. 모형의 입력자료로 유역 인근의 12개 기상관측소로부터 시단위 기상자료를 구축하였으며, 모형의 검보정은 일단위 관측유량(Q), 플럭스 타워 증발산량, 실측 토양수분 및 지하수위 자료를 구축하여 활용하였다. 댐 및 보 지점에 대해 Q와 1/Q로 검보정을 수행한 결과, 평균 결정계수(R²)는 댐 지점에서 0.53~0.65, 보 지점에서 0.46~0.69의 값을 나타냈으며, Nash-Shtcliffe efficiency(NSE)는 댐 지점에서 0.46~0.55, 보 지점에서 0.31~0.65의 값을 나타냈다. 공간 보정을 위해 증발산량, 토양수분, 지하수위에 대한 검보정을 수행할 예정이며, 유황곡선을 활용하여 하천차수, 토양속성 및 토지이용에 따른 하천갈수량을 분석할 예정이다.

• Korean Journal of Materials Research
• /
• v.6 no.1
• /
• pp.99-105
• /
• 1996

본 연구에서는 TIMA를 전구체로 하는 수직형 MOCVD 반응기를 대상으로 수학적 모델을 세우고 컴퓨터에 의한 수치모사를 수행하여 반응기 설계 변수 및 공정조건이 AI의 증착속도와 증착두께 분포에 미치는 영향을 알아보았다. 수학적 모델은 수직형 반응기를 축대칭으로 보아 2차원으로 수립하였으며 반응기내의 운동량전달, 열전달, 물질전달을 포함한다. 이 수학적 모델의 지배 방정식들에 대하여 Galerkin 유한요소법을 적용하여 수치적으로 반응기 내의 유체 흐름 구조, 온도분포와 반응물의 농도 분포를 구하였다. 수치모사 결과 AI의 증착속도는 반응기 압력이 0.47torr, 기판온도가 25$0^{\circ}C$, 유량이 7.5sccm일 경우, 190-230$\AA$/min로 나타났다.