Effect of the nonuniform grid on the two-dimensional transport equation was investigated in terms of theoretical analysis and finite difference method (FDM). The nonuniform grid having a typical structure of the numerical weather forecast model was incorporated in the vertical direction, while the uniform grid was used in the zonal direction. The staggered and non-staggered grid were placed in the vertical and zonal direction, respectively. Time stepping was performed with the third-order Runge Kutta scheme. An error analysis of the spatial discretization on the nonuniform grid was carried out, which indicated that the combined effect of the nonuniform grid and advection velocity produced either numerical diffusion or numerical adverse-diffusion. An analytic function is used for the quantitative evaluation of the errors associated with the discretized transport equation. Numerical experiments with the non-uniformity of vertical grid were found to support the analysis.
범용유한요소프로그램인 Abaqus의 확장유한요소법(XFEM)의 사용성 검증을 위하여 2차원 모델에 적용하여 수치해석을 수행하였다. 기존의 연구에 많이 사용되었던 응집요소(cohesive element) 모델은 균열 경로를 예측하고 요소를 삽입하여야 하는 단점 때문에 실제 균열을 모사하는데 한계가 있다. 이러한 이유로 응력의 방향성 및 특이성을 바탕으로 균열의 경로를 예측하는 확장유한요소법(XFEM)이 균열 해석에 있어서 더 발전된 방법으로 이용되어 왔다. 이번 연구에서는 균열의 경로가 자명한 2차원 모델에 사용하여 응집요소해석과 XFEM에 응집요소의 물성을 적용한 해석을 비교하고 XFEM 적용의 타당성을 확인하였다. 수치해석으로 균열 발생 직전의 응력분포 및 응력 특이성을 확인하고 실제 균열 발생경로와의 비교를 한다. 본 연구를 바탕으로 몇 가지의 한계를 극복하면 실제 복잡한 모델의 실제 균열진전해석을 수행하여 균열을 모사할 수 있을 것으로 기대된다.
본 논문에서는 2개 이상의 기하형상이 순응 또는 비순응 경계면에서 접합된 멀티패치 문제에 대한 아이소-지오메트릭 해석에 대해서 연구하였다. 패치 경계면에서 응력의 연속성을 표현하는 방법으로 Nitsche 방법론과 마스터-슬레이브 방법에 기반한 방법론에 대해서 지배방정식을 유도하고 아이소-지오메트릭 이산화를 수행하였다. 멀티패치 문제에 대해서 두 방법 론의 차이점을 간단하게 비교하였으며, 후처리 과정에서 사용되는 NURBS 곡면 기반의 응력 복원법에 대해서 기술하였다. 수치예제에서 비순응 경계면을 가지는 멀티패치 빔 문제를 통해 Nitshce 방법론을 검증하였으며, 응력집중을 가지는 문제에서 소개된 두 방법론이 유사한 결과를 보이는 것을 확인하였다. 소개된 NURBS 곡면 기반의 응력 복원법을 후처리에서 도입할 경우 멀티패치 문제의 경계면에서 개선된 연속적인 응력을 보임을 알 수 있다.
월성원자력환경센터 중저준위방사성폐기물 처분시설은 서로 다른 방식의 처분시설이 혼재하고, 월성 원자력발전소와도 인 접해있다. 이와 같은 높은 복잡성으로 인해 처분시설 안전성 평가 시 보다 면밀한 현상이해가 필요하다. 기존 1단계 사일로 처분시설의 성능평가모델들에 포함된 불필요한 보수성을 줄이고 복합처분시설에 대한 보다 실제적인 성능을 파악하기 위해 서는 다차원 수리/핵종이동 모델이 필요하다. 이와 함께 향후 복합처분시스템의 특성에 기인한 다양한 불확실성을 관리하 고 파라미터의 중요도를 분석하기 위해 많은 계산이 필요할 것으로 예상하며, 이를 위해 보다 효율적인 성능평가 모델이 요 구된다. 본 논문에서는 두 요건을 충족시키기 위해 수리성능 모델과 핵종이동 모델을 연계한 2단계 천층처분시설의 근계영 역 2차원 통합성능평가 모델을 개발하였다. 수리 및 핵종이동은 PORFLOW와 GoldSim 전산 코드를 이용해 평가하였으며, GoldSim 핵종이동 모델은 PORFLOW 핵종이동 모델과의 벤치마크를 통해 검증하였다. GoldSim 모델은 계산효율이 뛰어났 으며 기존의 모델에 비해 핵종이동거동을 이해하는데 용이하였다.
The finite element model of various laminated composite structures with delamination is developed based on the higher-order shear deformation theory. In the finite element formulation for the delamination, the seven degrees of freedom per each node are used with transformations in order to fit the displacement continuity conditions at the delamination region. The numerical results obtained for various composite plates and shells with delaminations are in good agreement with those of other preceding investigations. The new results for laminated structures in this study mainly show the effect of the interactions between the geometries and other various parameters, for example, delamination size, the number of layer and location of delamination in the layer direction. Key observation points are discussed and a brief design guideline is given.
Effect of ball filling ratio on the behavior of balls motion and their collision characteristic in a tumbler-ball milling of metal powder are investigated by a computer simulation. The discrete element method and the extended Kelvin model composed of nonlinear spring and nonlinear dashpot were employed in the simulation. It can be possible that analysis of the individual balls motion in a three-dimensional actual mill by the two-dimensional model simulation, since the simulated trajectories of ball paths are in relatively good agreement with the actual ones. It knows that the balls motion in the tumbler-ball mill is strongly influenced by the surface conditions of the balls and mill container wall. The energy consumption of the individual balls during impact and the impact frequency of the individual balls increased with an increase in the ball filling ratio and showed maximum values at about 50-60% ball filling ratio, and then decreased.
본 연구에서는 4대강 살리기 사업 후 퇴적현상이 지배적으로 발생하는 남한강과 섬강 합류부 구간을 대상으로 2차원 수치모형인 CCHE2D 모형을 이용하여 하천의 흐름 및 하상변동에 대한 해석을 수행하였다. 대상지점 합류부는 남한강 본류의 만곡부에 지류 섬강이 유입되는 특성을 갖는다. CCHE2D 모형은 비평형 유사이송을 해석하며 소류사와 부유사 조정거리가 중요한 입력변수로 대상지점에서는 소류사 조정거리가 하상변동에 가장 큰 영향을 주는 것으로 나타났다. 수치모의 결과 유량비(Qr) 변화가 남한강과 섬강 합류부 지점에서 흐름 및 하상변동에 영향을 미쳤으며, Qr≤ 2.5인 경우에는 합류전 본류의 유속이 증가하여 흐름박리구역을 감소시켰으며 이로 인해 합류부 내측의 퇴적이 감소하였다. Qr>2.5이면 합류부 구간에 퇴적이 증가하여 사주가 형성될 가능성이 높은 것으로 나타났다. 수치모의를 통해 2013년에 발생한 유량비 변화에 의해 합류부에 고정사 주가 형성된 것을 알 수 있었다.
The stability of the river in the restored river is an important issue in maintenance and management. Bed elevation change simulation can be an effective way to predict the direction of river restoration by predicting long and short term bed elevation change of river. A 2D numerical model (CCHE2D) was implemented to simulate the long-term bed elevation change. The study area is located in the Cheongmi-Cheon Notap-ri and 1.2 km long. The flow scenario was constructed using the flow data that was measured at the water level observatory located upstream at the Janghowon Bridge. The bed elevation change pattern according to restoration of abandoned channel was analyzed and the stability of river was evaluated.
본 연구는 식재된 개수로에서 흐름특성을 모의할 수 있는 수심 적분된 2차원 수치모형을 이용하여 원형 식생역 주변의 흐름을 수치모의하였다. 식생영향을 고려하기 위해 식생항력 항을 지배방정식에 추가하였고 다양한 식생체적비율(SVF) 조건에 따른 수치모의를 수행하였다. 흐름이 원형 식생역을 통과하고 하류에 저유속 구간인 후류영역(wake region)을 형성하며 식생체적비율이 0.08 이상이면 재순환 영역이 발생하였다. 재순환 발생위치는 식생체적비율이 감소하면 식생역 에서 더욱 하류로 이동하였다. 후류영역을 지나 원형 식생역 양 측면에서 유발된 전단층들의 상호작용에 의해 von Kármán 와열이 발생하였다. 원형 식생역 하류에서 발생하는 와류는 식생체적비율이 0.08 이상이 되면 나타나기 시작하였고 발생위치는 난류운동에너지가 최대값을 보이는 위치와 일치하였다. 최대 난류운동에너지는 식생체적비율이 감소하면 줄어드는 것으로 나타났고 최대값의 발생위치는 점점 하류로 이동하였다.
본 연구에서는 SMS(Surface Modeling System)내의 Tuflow 모형을 이용하여 2차원 홍수범람모의를 수행하고 1차원 홍수범람모의 결과와의 단순비교를 실시하였다. 1차원 홍수범람모의는 HEC-RAS 모형을 통하여 얻어진 수위를 WMS 모형을 통해 홍수범람을 실시한 결과를 이용하였으며, 이를 2차원 홍수범람모의와 비교하였다. 홍수위를 중심으로 단순홍수위 비교를 수행하였으며, 비교결과 최대수심을 중심으로 일반적으로 2차원 홍수범람모의가 작게 산출되는 경향을 보였으나 일부구간에서는 1차원 홍수범람수위가 작게 나타나는 등 예측가능한 일반적 경향성은 보이지 않는 것으로 나타났다. 향후 홍수범람상황의 재현성 등에 대한 조사가 이루어졌을 경우 보다 신뢰할 만한 자료를 얻을 수 있을 것으로 기대되며, 추가적인 연구가 지속되어야 할 것이다.
기후변화로 인한 기상이변으로집중호우가 빈번히 발생함에 따라서 대규모 사면 재해의 발생 빈도 및 피해규모가 급증하고 있는 추세이다. 이러한 사면 재해 중, 토석류는 이동속도가 빠르고, 이동거리가 길어 피해 범위가 넓은 특징을 지니고 있어 인구와 재산이 밀집된 도심지에서 토석류가 발생시 대규모의 인명 및 재산피해를 예상하게된다. 이러한 토석류 재해로부터 시민의 생명을 보호하기 위해서는 토석류로 인한 하류부 피해 예상범위를 사전에 파악하는 일이 무엇보다 중요하다. 본 연구는 2차원 토석류 모의 수치모형인 Flo-2D와 RAMMS를 이용하여 토석류 흐름에 의한 대모산 하류부 피해 예상지역을 예측하였다. 수치모형의 검보정을 위하여 2011년에 발생한 우면산 산사태 자료를 이용하여 점성, 항복응력, 유사체적 농도 값을 보정한 후, 대모산에 100년 빈도 강우가 발생했을 때에 하류부 피해 영향 범위를 예측한 후 두 모형의 차이를 비교 분석하였다. 본 연구를 통하여 작성된 토석류 재해지도는 추후 토석류 예경보 및 주민대피계획 수립 등에 활용될 수 있을 것이다.
본 연구에서는 낙동강 본류의 강정고령보에 설치된 인공하도식 어도의 효율성을 2차원 물리서식처 모형인 River2D를 이용하여 평가해 보았다. 강정고령보 상하류단의 실측 수위를 이용하여 모형을 보정 및 검정하였다. 모의조건은 2012년 풍수량, 평수량, 저수량을 기준으로 어도 입구 및 어도 내부의 흐름장을 분석하였다. 특히 본 연구의 목표종인 피라미가 산란기인 봄철에 주로 이동을 하므로, 저수기를 기준으로 어도가 기능을 잘 유지하는지 살펴보았다. 분석결과, 풍수 및 평수기에는 본류 수위가 높아 어도로의 피라미 유입 효율성이 떨어져 보였다. 반면 저수기에는 본류의 낮은 수위로 인해 어도 입구에 상대적으로 적절한 유속이 발생하여 유입효율과 이동효율이 높아질 것으로 평가되었다. 어도 내 흐름은 어류가 소상하기에 무리가 없는 유속장이 발생하는 것으로 파악되었다.
본 연구에서는 2차원 수치모형을 이용하여 폭풍해일고에 따른 해안도시지역 내에서의 범람모의 및 분석을 수행하였으며, 적용된 대상지역은 2003년 태풍 "매미" 내습 시 가장 피해가 컸던 지역인 경남 창원시 마산 합포구 신도시 지역 일부이다. 적용된 수치모형의 검증차원에서 태풍 "매미" 발생기간 동안 관측된 폭풍해일고 자료와 대상지역 내의 주요 지점에 대한 침수흔적 자료를 이용하여 비교 분석을 수행하였다. 또한 향후 발생 가능한 수퍼 태풍에 의한 영향을 분
본 연구에서는 댐 붕괴파와 같이 연속 및 불연속 흐름해석에 적용되고 있는 HLLC 기법을 불규칙한 하상지형에서의 흐름해석에 적용할 때 생성항과 흐름률항의 사이의 수치적 불균형으로 인한 수치진동을 감소시키기 위해 well-balanced HLLC 기법과 천수방정식에 근거한 비구조적 유한체적모형을 개발하였으며, 이를 댐 붕괴파 문제에 적용하였다. 적용된 well-balanced HLLC 기법은 단순히 흐름률항을 계산할 때 하상지형경사를 직접 포함시키는 것으
본 연구에서는 자연하천과 같은 복잡한 지형의 처리에 효율적인 삼각형 및 사각형 혼합격자의 적용이 가능한 2차원 유한체적모형을 개발하였다. 이를 위해 계산격자의 인접격자를 찾는 알고리즘을 제안하고, 제안 기법을 개발모형에 적용하여 계산격자 및 인접격자의 경계면에서의 흐름률을 HLLC 근사 Riemann 해법을 이용하여 계산하였다. 또한 흐름률과 생성항사이의 균형에 중요한 영향을 주는 혼합격자의 하상경사 처리를 위해 삼각형 및 사각형 격자에 대해 각각 다른
본 연구의 목적은 유한요소법에 의한 이차원 오염물 이송에 대한 등수역학적 수치모형을 개발하는 것으로서 개발된 모형을 이용하여 지류와 폐수처리장으로부터의 오염물이 주하도로 유입되는 경우에 대한 수질해석을 수행하였다. 본 연구에서는, 합리적인 매개변수산출을 위한 보정과 모형의 신뢰성과 정확도를 위한 검증을 수행하였다. 또한, 실제 유역인 한강하류부에 본 모형을 적용한 결과 하도를 따라서 실제적인 오염물 분포양상을 재현하였으며 계산값과 실측값의 비교에서 잘
침식률 측정장치인 SEDFLYME을 사용하여 직접 측정한 침식률 자료를 이용하는 퇴적물 수송에 관한 2차원 모형이 개발되었다. 개발된 모형은 현장 적용성을 높이기 위하여 경계밀착좌표계를 수평방향에 대해 사용하며, 퇴적물 이동 모의시 침식률 산정에서 불확실성을 줄이기 위해 침식률 측정자료를 직접 사용한다. 개발된 모형은 부유사와 소유사 이동을 모두 고려한다. 모형의 정확성을 검토하기 위하여 1차원 수로에서 장갑화 현상을 모의하고 실험자료와 비교하였다. 비
홍수터 흐름의 모의를 위한 준2차원 계산모형을 수립하였다. 모형의 계산망으 2차원 홍수터 구획체계를 하도와 결합한 것으로서 일반적으로 폐합형 망으로 구성된다. 홍수터 흐름에 대해서는 각 구획에서의 수량보존에 관한 연속방정식 및 인접구획간 하도형 또는 월류형 수위-유량 관계식을, 하천 본류에 대해서는 1차원 부정류에 대한 St.Venant 방정식을 각각 지배방정식으로 하여 흐름을 모의하는 준2차원 계산모형으로서, 폐합형 하천수계에 대한 부정류 해석 수치기
충격파의 높이나 속도는 홍수제어조작이나 수로벽과 빠른 유속을 가지는 하천에서 교량의 설계에 중요한 자료가 된다. 따라서 광범위한 조건에서 흐름의 불연속면을 모의할 수 있는 수치모형이 요구된다. 본 연구에서는 천수방정식을 지배방정식으로 한 Godunov 형 유한체적법 모형을 개발하였다. Riemann 해법으로 Roe(1981)의 해법이 사용된다. 이 모형은 본 연구에서 비구조적격자(unstructured grids)를 사용하기 위해 개발된 수정 MUSCL
본 논문에서는 수심적분 2차원 모형을 하나 소개하고, 이 모형이 가질 수 있는 오차의 한계를 민감도분석을 통해 제시하였다. 민감도 분석 결과, Manning 조도계수, 혼합계수, 침강속도, 한계전단력 등이 부유사농도에 대해 모형에 미치는 영향은 매우 크게 나타났으며, 조위와 확산계수의 영향은 거의 없는 것으로 나타났다.