• Proceedings of the Korean Society of Coastal and Ocean Engineers Conference
• /
• 1998.09a
• /
• pp.24-30
• /
• 1998

고도의 산업발달은 인간환경과 밀접한 하구와 해안의 오염을 심화시켜 최근 사회적 문제로 대두되고 있으며, 이에 오염물질의 이동을 예측하고 제어하려는 노력이 꾸준이 진행되어 왔다. 이 같은 노력은 수치모형을 이용한 수질관리 연구에도 많은 진전을 가지고 왔으며, 오일러적 모형, 라그랑쥐적 모형, 오일러-라그랑쥐 적 모형이 대표적이다. 오일러-라그랑쥐적 모형은 이류에 대해 라그랑쥐적 방법, 확산에 대해 오일러적 방법으로 해석함으로서 각 방법들의 장점을 취하여 수치적 진동, 수치적 확산이 적고 효율성이 뛰어나 최근 많이 연구되고 있다. (중략)

• Proceedings of the Korea Water Resources Association Conference
• /
• 2010.05a
• /
• pp.253-257
• /
• 2010

최근 들어 이상 강우 현상으로 인한 기록적인 집중 강우와 더불어 토지 이용 변화로 인해 탁수의 발생 빈도가 증가하고 있다. 강우 유출로 유입한 탁수는 저수지 내 장기간 체류하면서 하류 하천의 수질 및 수생태계 뿐만 아니라 저수지 내부의 영양단계에도 많은 영향을 준다. 특히, 성층화된 저수지의 경우 높은 영양염류 농도를 포함한 탁수는 밀도류 거동 특성을 보이면서 수평 및 수직 혼합 과정을 거치면서 국부적인 부영양화 현상과 조류의 수화 현상의 원인이 되고 있다. 따라서 대형 저수지의 수질관리에 있어 하천 유입 탁수의 밀도류 해석, 저수지의 수온 성층 구조 변화, 부유입자의 동력학적 해석이 중요한 요소로 부각되어 왔다. 본 연구에서는 소양호를 연구 대상 지역으로 선정하여 2005년과 2007년 수문 사상을 바탕으로 2차원 횡방향 수치 모형을 구축하였다. 수치모형을 통해 수온 성층 구조의 재현성을 확인하였으며, 다양한 탁수 거동 모형을 구축하여 적용성을 평가하였다. 유입수의 SS(Suspended Solid)를 단일 입경으로 가정한 TM-1 모형, SS의 입경분포에 따라 3개의 그룹(SSi)으로 구분한 TM-2 모형, 3개 그룹을 포함하면서 저수지내 탁수 장기화로 인한 탁수 저감 효과를 1차 반응상수로 매개 변수화(유기물 함량($a_0$) ${\times}$ 분해속도(${\lambda}_a$))하여 수정된 지배방정식을 적용한 TM-3 모형을 사용하였다. 각각의 탁수 거동 모형은 2005년과 2007년 수문 조건에서 수온 성층 구조를 잘 재현하였다. TM-1 모형과 TM-2 모형을 비교해보면, 탁수 중심축의 최고 탁도에 대한 예측 성능은 TM-2 모형이 우수한 결과를 나타냈었다. 하지만, 장기 탁수 모의 시 저수지 수중 잔류 SS가 지속적으로 높게 나타나 중층 탁도를 과대평가하는 경향을 보였다. TM-3 모형이 TM-2 모형에 비해 수심별 탁도 분포에 대한 중심축 탁도가 저평가되는 경향을 보였지만 저수지 내 잔류 탁도에 대한 영향 부분에서는 개선된 결과를 나타내었다. 본 연구 결과는 저수지 탁수 밀도류 해석 및 운영 시스템에 활용 될 수 있으며, 선택 취수 설비 등의 수리 구조물의 영향 평가에 활용할 수 있다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.26 no.5C
• /
• pp.351-358
• /
• 2006

The purpose of this study is to recommend the simulation method and procedure of behaviors of CFRD(Concrete Faced Rockfill Dam) concrete face slab with impoundment by centrifuge tests, to examine the effects of the flexural rigidity of the concrete face slab on the face slab deformation from the centrifuge tests, and to evaluate the effects of the stiffness of face supporting zone on the displacement and moment of face slab by numerical analysis which is verified by the centrifuge tests. In this study, the centrifuge tests on the two model dams with the concrete face slab of different flexural rigidity were carried out. Also, the centrifuge tests were simulated by numerical analysis of which input material properties were obtained by the triaxial tests on the model materials. The validity of numerical analysis was evaluated by comparison between the results of centrifuge tests and numerical simulation. The deformation pattern of the concrete face slab was examined with the various stiffness of the face supporting zone by numerical analysis. From the results of centrifuge tests, the effects of face slab thickness on the deformation of face slab were negligible. From the results of centrifuge tests and numerical analysis, it was found that the amplitude of the maximum displacement of face slab and the position where the maximum displacement was mobilized with impoundment were affected by the stiffness of face supporting zone rather than the flexural rigidity of concrete face slab.

• Journal of Korean Port Research
• /
• v.11 no.1
• /
• pp.121-128
• /
• 1997

To analyze the ocean wave more efficiently, more fine grids are used with relatively less computer memory. Each element of free surface is discretized into more fine grids because the ocean waves are much influenced by the mesh used in the finite difference scheme. According to the flow analysis, remarkable improvements could be seen in the free surface generation. The multi grid is applied to confirm the validity of scheme. The Baldwin Lomax turbulence model is used for the analysis of S103 Inuid ship. Finally some discussion on experiments was made for the physical phenomena of the viscous

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.43-43
• /
• 2011

본 연구에서는 한강제방 취약구간인 성산대교와 잠실철교 부근 지역에 대한 홍수범람을 수치모의 하였다. 사용된 수치모형은 스위스의 Beffa에 의해 개발된 FLUMEN(FLUvial Modeling ENgine)으로서 스위스, 독일, 오스트리아 등에서 홍수범람해석에 사용된 바 있는 모형이다. 제방 취약지역은 한강 하천정비기본계획(2002)과 대학과 연계한 하천관리에 관한 연구용역(2단계 4차년)에 제시된 HEC-RAS 부등류 해석에 의해 계산된 홍수위와 기존 제방의 높이를 비교하여 산정하였다. 범람모의를 위해 HEC-RAS 부정류 해석을 통해 경계조건을 산정하고, FLUMEN을 이용하여 한강제방의 취약지역에 대한 범람모의와 파제 시나리오를 작성하여 파제에 따른 범람모의를 실시하고 적용성을 검토하였다. 제방 취약지역에서 월류로 인한 범람이 발생하였고, 이로 인해 일부 주거지역이 침수되었다. 가상 파제시나리오를 통한 수치모의 결과에서 여의도에서 $2.179km^2$의 넓은 지역에 침수현상이 계산되었으며, 최대 침수심은 5.054m로 성산대교 남단의 가상 파제 시나리오에서 계산되었다. 수치모의 결과 FLUMEN은 한강유역의 범람모의에 적합하다고 판단되며, 본 연구는 홍수 방어대책을 수립함에 있어 홍수 취약지역의 선정과 수공구조물의 설치 방향을 결정하는데 기초자료로 활용이 기대 된다. 정확도를 높이려면 보다 정밀한 제내지, 제외지의 측량자료의 적용, 내수 침수모형과의 연계, 조밀한 불규칙 삼각망의 작성이 필요하며, 국내실정에 맞는 정확한 적용기준마련과 유역의 수리학적 특성이 고려된 홍수해석 모형의 개발이 필요하고 보다 정확한 조도계수의 산정을 위한 지속적인 연구가 필요하다.

• Water for future
• /
• v.29 no.2
• /
• pp.221-233
• /
• 1996

To investigate the hydrodynamic characteristics of a two-dimensional buoyant surface jet, the most important factors of the numerical analysis are the evaluation of the free surface and the turbulence transportation under the stratification. In present study, a numeriacal simulation model used with the semi-implicit method for pressure-linked equations (SIMPLE), the non-hydrostatic approximation and the algebraic stress model (ASM) is applied to investigate the vertical structure of internal flow hydrodynamically. The ASM enables to take account of anisotropy of turbulence, the damping effects of the density interface, and the free surface on the turbulence structure accurately. The ASM tested produces better agreement than the $\kappa-\varepsilon$ model with measurements by Nakatsuji (1984) on the flow development and turbulence structure. Applicability of the ASM to a two-dimensional buoyant surface jet is examined through comparison with experimental data.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2007.05a
• /
• pp.1268-1271
• /
• 2007

트렌치 수로에서 동수압을 고려한 자유수면 흐름을 해석하고 난류 모형의 적용성을 제시하였다. 본 연구에 사용된 지배방정식으로 비정상 상태의 비압축성 유체에 대한 연속방정식과 비점성 Navier-Stokes 방정식을 사용하였다. 난류완결문제를 해결하기 위해서 $\kappa-\varepsilon$방정식을 사용하여서 난류 와점성을 구할 수 있다. 자유수면과 동수압을 고려하기 위하여 3단계로 나누어서 해석하였다. 제 1단계에서는 운동량방정식을 연직방향에 대해 음해적으로 차분하였다. 제2단계에서는 유속과 동수압 보정항으로 이루어진 식을 연속방정식에 대입하여 타원형 방정식인 동수압-포와송 방정식을 해석하여 얻어지는 유속은 질량보존법칙을 만족하게 된다. 마지막으로 자유수면과 최종유속을 보정 및 계산하였다. 본 연구에서 제시한 수치모형을 검증하기 위해서, 트렌치 수로에서 난류의 흐름에 변화를 고려하기 위해 수치모의를 하였다. 전반적으로 수치모의에 의한 결과와 실험 자료가 일치하는 경향을 보였다.

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

보나 댐에 설치된 배수관이나 지하에 매설된 상하수도관과 같이 물을 운반하기 위한 관로 뿐만 아니라 유류를 운송하는 관, LPG와 같은 기체를 운송하는 관처럼 유체를 운송할 때 다양한 관로를 사용한다. 그 중 사용범위나 제작에 대한 용이성 등에 의해서 원형관이 주로 사용된다. 따라서 위와 같은 원형관 내의 관수로 흐름분포에 대한 연구는 아주 중요하며 필요하다. 원형관 내 흐름분포는 관에 연결된 펌프 혹은 수조나 저수지의 수위에 의한 관내 압력에 의해 지배되어지며, 관 내부 표면의 거칠기나 관의 꺾인 정도 등 다양한 요인에 영향을 받는다. 본 연구에서는 일반적인 실험결과를 도출하기 위하여 직선의 원형관을 대상으로 실험과 수치모의를 동시에 수행한다. 실험은 높이 3m, 길이 4.5m, 폭 1.5m 수조의 바닥에서부터 0.45m 위에 위치한 1.8m 길이의 아크릴 재질의 원형관이 설치된 수조에서 진행되었으며, 수치모의 또한 동일한 규모에서 수행되었다. 수조의 수위를 변화시켜 여러가지 레이놀즈 수에 따른 관 내 흐름의 변화에 대하여 분석하는 것이 목적이며, 실험결과와 수치모의 계산결과 간의 비교검증을 통해 분석한다. 이 연구에서는 난류의 영향을 RANS(Reynolds-averaged Navier-Stokes)와 LES(large eddy simulation)을 혼합한 형태인 Hybrid RANS/LES 모형 중 하나인 DES(detacged-eddy simulation)기법을 이용하여 해석함으로써, 빠른 유속으로 인한 관 내의 강한 와(eddy)에 대한 효과를 예측하는데 한계가 있는 RANS 모형과 벽면근처에서의 흐름 해석을 위해 굉장히 높은 격자해상도가 요구되어 계산적 비용 측면에서 효율이 떨어지는 LES 모형의 한계를 극복하고자 한다. DES 모형은 벽에서 떨어진 와에 대하여 LES로 직접해석하고, 벽 근처에서의 흐름에 대해선 난류모형을 통해 모델링함으로써 벽 근처 계산격자와 계산량을 줄이면서 LES와 비슷한 결과를 얻을 수 있다. 실험결과와 수치모의 계산결과 사이의 비교검증을 통하여 관 내 흐름에 대한 수치모의의 적용성을 평가하고 실험에서 측정하기 어려운 난류강도와 압력변동의 상세한 특성을 수치모의를 통해 분석함으로서 관 내 흐름특성에 대하여 자세히 분석하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2020.06a
• /
• pp.149-149
• /
• 2020

하천에서 물 흐름이 보와 댐과 같은 수공구조물을 지날 때 일반적으로 흐름상태에 다양하고 급진적인 변화가 발생한다. 특히 흐름이 구조물을 지나면서 사류(supercritical flow)로 변하고 다시 상류(subcritical flow)로 복원되면서 일어나는 도수(hydraulic jump) 현상은 수위의 급변화, 흐름 에너지 소산, 변동성이 강한 압력 분포 등이 특징이다. 이러한 흐름 특성들은 보나 여수로와 같은 수공구조물 자체의 성능뿐만 아니라 이들 수공구조물의 하류에서 발생하는 국부세굴로 인해 구조물의 안정성에 부정적인 영향을 줄 수 있다. 따라서 수공구조물을 설계할 때는 이들 구조물을 통과하는 흐름의 비정상 난류 흐름 특성을 정확하게 해석하여 반영하여야 한다. 이 연구에서는 k-omega SST 난류 모형과 자유수면의 급격한 변동을 해석하기 위한 하이브리드-VOF(hybrid volume of fluid)기법을 이용하여 도수현상을 수치적으로 재현하고자 한다. 기존 CFD(computational fluid Dynamics) 모델링에서는 자유수면 변동의 영향을 고려하기 위해 VOF 기법을 많이 사용하였다. 하지면 전통적인 VOF 기법은 다상흐름(multiphase flow)을 오직 부피분율(volume fraction)의 함수로만 고려하며 모의함으로써 강한 수면변동뿐만 아니라 공기연행(air entrainment)를 동반하는 난류 흐름을 모의하는데는 한계가 있다. 이 연구에서 이용하는 Eulerian 기법인 하이브리드 VOF 기법은 물과 공기의 각 상에 대하여 흐름 특성들을 개별적으로 계산하기 때문에 공기연행을 포함한 급변류 흐름에서 전통적인 VOF 기법보다 적용성이 우수하다. 이와 같은 난류모형과 자유수면 포착기법을 이용하여 3차원 비정상 난류 흐름 수치모형을 구축하여 수공구조물 주변에서 발생하는 강한 공기연행과 난류 특성를 보이는 급변류를 수치적으로 재현한다. 이 연구는 계산된 수치해석 결과를 기존 수리실험 결과와 비교하여 수치모형의 적용성을 평가하고 도수 현상에서 발생하는 독특한 흐름 특성을 제시한다.

• 한국방재학회:학술대회논문집
• /
• 2011.02a
• /
• pp.59-59
• /
• 2011

파랑의 전파와 변형에 대한 연구에는 수심방향으로 적분한 2차원방정식인 완경사방정식과 Boussinesq 방정식을 기반으로 한 수치모형을 이용한 연구가 최근까지 가장 활발하게 진행되어 오고 있다. 그러나 실제 구조물의 설계에는 2차원 수치모형에서 고려할 수 없는 수심방향 유속에 기인한 정확도의 문제로 인해 구조물의 형상과 재원을 설계하기 위한 정교한 수치모형실험이 어려워 주로 수리모형실험에 의존해 왔다. 수리모형실험은 실제 현상을 가장 잘 재현해낼 수 있어 신뢰성이 매우 높지만 다양한 실험을 수행하기가 어렵고 많은 시간과 비용이 소요되는 단점이 있다. 이에 따라 최근 수심방향으로 완전한 운동방정식인 Navier-Stokes 방정식을 푸는 3차원 수치모형에 대한 연구가 활발히 진행되고 있다. 이론적으로 매우 우수한 모형이긴 하나 정확도 높은 결과를 얻기 위해서는 매우 조밀한 격자를 필요로 하기 때문에 아직까지 막대한 계산시간이 필요하다는 단점이 있으나 컴퓨터 기술이 급격한 속도로 발전하고 있어 Navier-Stokes 방정식 모형의 적용 가능성은 계속 높아지고 있다. 파랑변형을 다루는 수치모형실험을 수행할 때 외부조파를 사용할 경우 구조물이나 지형에 의해 반사되어 나온 파랑이 조파지점에 도달할 때 실험영역으로 재 반사되는 문제가 발생한다. 이를 해결하기 위해 내부조파기법의 개발에 대한 연구가 필수적이었으며, 자유수면변위를 변수로 사용하는 모형의 경우 그 연구가 매우 활발하게 진행되어 왔다. 한편, Navier-Stokes 방정식 모형의 경우 자유수면변위를 변수로 사용하는 2차원 모형에 비해 상대적으로 연구가 미흡하였다. 본 연구에서는 기존의 연직 2차원 Navier-Stokes 방정식 모형에 사용된 연속방정식에 질량 원천항을 추가하는 내부조파기법을 도입하여 3차원 수치모형에서 고립파를 내부조파하고, 급경사에서의 고립파의 처오름 및 처내림 현상을 수리모형 실험결과와 비교 및 분석하였다. 수치모형은 Navier-Stokes 방정식을 엇갈림 격자체계에서 계산하는 동수압 모형으로서, Two-step projection 기법을 사용하는 유한차분모형을 사용하였다. 본 수치모형은 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(spatially averaged Navier-Stokes) 방정식을 계산하는 LES(large-eddy-simulation) 기반의 수치모형으로, 난류 모델링을 위해 Smagorinsky LES 모형을 사용한다. 또한, 압력장의 계산을 위해 Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하였으며, 자유수면 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 수치모형실험이 전체적으로 수리모형실험에서 관측한 파랑의 처오름 및 처내림 현상을 잘 재현하고 있는 것으로 나타났으며, 정량적인 비교를 통해 수치모형의 성능을 검증하였다.