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

댐 제방 등의 붕괴로 인하여 발생하는 급격한 유량의 변화와 흐름영역의 변화로 인한 천이류 및 도수의 발생, 불규칙한 하천단면에서 갈수기 저수기의 흐름해석은 기존의 수치해법의 한계로 인하여 수리모형실험 및 경험식 또는 단면의 단순화 등에 의존하고 있는 실정이다. 본 연구에서는 자연하천에서 비선형 흐름율 계산에 불연속초기조건의 해석해인 Riemann 근사해법을 사용하여 수치적으로 안정되고 정확한 1차원 모형을 개발하고자 한다. 이를 위하여 유한체적법을 사용하였고, 수위와 유량의 계산을 위하여 요구되는 유한체적을 유출입하는 흐름율의 계산에 HLL Riemann 해법을 사용하였으며, MUSCL 기법으로 2차 정확도기법으로 확장하였다. Riemann 해법을 통하여 계산된 비선형의 흐름율과 보존 특성을 만족시켜줄 수 있는 하상 및 하폭변화로 인한 생성항을 처리하는 기법을 제안함으로서 새로운 1차원 수치해석모형을 개발하였다. 개발된 모형의 실제하천의 적용성을 확인하기 위하여 하상과 하폭이 변화하는 부정류 흐름에 적용하여 모형의 적용성 및 정확성을 검증하였다.

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

본 연구에서는 2차원 비선형 천수모형과 수심평균된 스칼라 이송모형을 해석하는 수치모형에 대해 기술하였다. 수치모형의 정확성을 보장함과 동시에 안정성을 높이기 위해 유한체적법, 플럭스 재구성 및 minmod 제한자를 사용하였다. 비선형 천수방정식의 이송항과 바닥 경사항은 계산된 수심의 양수 보존과 흐름의 정상 상태를 보장하기 위한 second order well-balanced positivity preserving central-upwind method를 이용하여 수치적으로 이산화되었다. 마찬가지로, 이송-확산 방정식 내 이송항은 동일한 2차 풍상차분법을 통해 수치적으로 풀이하였다. 격자점 경계면에서의 불연속으로 인한 수치진동을 방지하기 위해 이송항의 계산에 포함된 보존항의 차이로 인해 발생하는 스칼라의 수치확산을 최소화하기 위해 무차원의 비소산함수를 도입하였다. 또한, 확산항은 유한차분법을 이용하여 이산화하였다. 제안된 수치모형은 시간미분항의 계산을 위해 오일러 기법을 적용하여 계산된 수심 및 스칼라의 양수 보존여부와 함께 정지된 흐름의 정상 상태의 보존여부를 확인하였다. 제안된 수치모형의 해석 정확성을 평가하기 위해 1, 2차원 공간 내 다양한 흐름 조건에서의 해석해를 이용한 3개의 벤치마크 테스트를 수행하였다. 평균 제곱근 오차(Root Mean Squared Error, RMSE)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2011.05a
• /
• pp.175-175
• /
• 2011

하천을 횡단하는 구조물은 흐름을 지체시켜 상 하류에 급격한 수위차를 발생시키며, 이에 동반되는 국부유속의 증가는 제방 및 하상을 세굴시킨다. 제방 및 하상 세굴은 하천횡단구조물의 구조적 안정성에 영향을 주어 불안정한 상태에 이르게 할 수 있다. 이와 같은 현상을 방지하기 위하여 하상유지시설 및 제방 보호공 등의 대책을 하천설계기준에 제시하고 있다. 하천설계기준에는 하천횡단구조물 보호시설의 종류, 설치길이 등이 명시되어 있으나, 보호공 제원 선정에 대한 구체적인 설계공식을 제시하지 않아 보완이 필요한 상황이다. 보호공 제원 선정 설계기준은 설계유속을 1차원 단면평균 유속으로 활용하는 국외의 설계공식을 활용할 수 있다. 그러나, 하천횡단구조물 주변의 흐름은 매우 복잡하므로 국외 설계공식을 그대로 적용하기에 앞서서 이에 대한 검토가 필요하다. 본 연구에서는 하천횡단구조물 주변의 흐름을 2차원 수치모형으로 재현하고 그 결과를 1차원 수치모의 결과 및 수리실험 관측결과와 비교하였다. 실험수로에 각기 다른 제원의 보호공을 설치하고 관측한 임계상태의 수리상황을 1, 2차원 수치모형으로 재현하였다. 각각의 모형에서 도출된 수리상황 결과를 비교 분석하여 보호공 제원 결정 공식을 검토하였다. 또한, 실제 하천구간에 임의의 보와 수제를 설치하고, 계획홍수량에서의 흐름특성을 비교하였다. 1, 2차원 모형에 의하여 도출된 관계를 분석함으로써, 하천횡단구조물 설계시 적용할 수 있는 보호공 제원 결정 설계공식을 제안하였다.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.1
• /
• pp.97-103
• /
• 1993

This paper presents a quasi-analytical method using integral equations to obtain head distributions in unsaturated porous media with different hydrogeologic properties. One-dimensional soultion algorithms were developed for two cases of boundary conditions at the top: 1) constant head and 2) constant flux. Water table elevation at the bottom was assumed known for both cases. The methodology was applied to a fly ash disposal facility in an alluvium area. The results showed that the pressure head distributions had high nonlinearity with large gradients slightly above the interface of two media which made preliminary numerical solutions implausible. The developed method helped to structure finite element grids for improving convergence and accuracy.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.2 s.25
• /
• pp.65-71
• /
• 2007

In this research, we made a one and two-dimensional analysis of numerical data collected from the bend curvature of a bended river section. According to the result from the numerical analysis, the inflow & output angle caused a water level deviation which increased with an increase of the flood discharge. From the water level deviation of our two-dimensional numerical model, we obtained the maximum slope of 6,67% when the inflow and output angle was 105 degrees and the flood discharge was 500 CMS. As for the right side, the differences with the one-dimensional numerical model were reduced when the angle was more than $90^{\circ}$. As for the left side the differences were reduced when the angle was more than $105^{\circ}$. For a river with more than 90 degrees bend curvature, a hydraulic experiment would be more appropriate than a numerical analysis.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.5 no.4
• /
• pp.336-339
• /
• 2004

In this paper the two-dimensional scattering rates were calculated in pseudomorphic Al/sub x//Ga/sub 1-x//As/Ga/sub y/In/sub l -y//As/GaAs heterostructure systems. The electronic states of the square quantum well were determined by the numerical self-consistent solution of Poisson's and Schrodinger's equations. The numerically obtained wave functions and energy levels were used to obtain the major two-dimensional scattering rates in this structure. Polar optical- and acoustic-phonon scattering, piezoelectric, ionized impurity and alloy scattering were considered for the first two sub-bands. The results were compared to the three-dimensional scattering rates also calculated in the same region.

• 한국방재학회:학술대회논문집
• /
• 2010.02a
• /
• pp.63.1-63.1
• /
• 2010

수중구조물에 의한 파랑의 변형을 예측하기 위해 3차원 수치모형을 도입하여 수치모형 실험을 수행하였다. 본 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(Spatially Averaged Navier-Stokes) 방정식의 해를 구하는 LES(large-eddy-simulation) 기반의 수치모형이다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two Step projection 기법을 사용하여 SANS 방정식을 계산하였으며, Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 선형파를 일정 수심상에서 조파시켜 해석해와 비교한 후 수중구조물이 설치된 지형에 적용하여 파랑의 변형을 수치모의하여 수리모형 실험 결과와 비교 및 분석하였다.

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

1차원 수치해석망은 수위관측소와 우량관측소를 경계로 광역모형을 구성하고 HEC-RAS를 이용하여 부정류 해석을 실시하였다. 그 결과로 얻은 유량과 수위자료는 2차원 해석 모형인 RMA-2의 경계조건으로 적용하였다. 또한 수치지형도를 이용하여 2차원 격자망을 구성하고 RMA-2를 적용하여 수리특성을 분석하였다. 광역모형의 구성에 의한 HEC-RAS 부정류 해석 결과 제방설치 전과 후의 수리특성의 변화는 그리 크게 나타나지 않았다. 또한 RMA-2 해석 결과 계획홍수량 초과하는 홍수에서 대상유역에서 제방설치 전과 후의 변화는 통수단면적의 감소에 따른 빠른 흐름과 수위의 변화를 가져왔으나 그 크기는 현저하지 않음을 보여준다.

• Journal of Wetlands Research
• /
• v.15 no.1
• /
• pp.91-103
• /
• 2013

In this Research, in order to improve problems of not enough technical validation and structural and hydraulic stability evaluation when nature-friendly revetment block is applied to field, hydraulic stability evaluation according to hydraulic behavior change of porosity soil block for vegetation reinforcement that secures ecological function was reviewed. By selecting object section, numerical analysis and hydraulic model experiments were performed; for numerical analysis, by using 1-dimensional numerical analysis model HEC-RAS and 2-dimensional numerical analysis RMA-2, one-dimensional(1D) and two-dimensional(2D) numerical analysis were performed; by applying Froude's similarity law, reduced-scale hydraulic model experiments according to vegetation existence were performed. In hydraulic model experiment, for validity of experiment result, the result of velocity and tractive force of reduced-scale hydraulic model experiments was converted to prototype so that it can be compared and reviewed under the same condition of one-dimensional(1D) and two-dimensional(2D) numerical analysis result; as a result, it was confirmed that comparatively united result appeared, and by comparing prototype-converted tractive force result with revetment's allowable tractive force coming from an existing research, block's hydraulic stability was suggested.

• Journal of Korean Tunnelling and Underground Space Association
• /
• v.25 no.1
• /
• pp.1-11
• /
• 2023

In this study, the stability of the room-and-pillar underground method was investigated using numerical analysis method. In-situ geotechnical investigation was conducted, and a supporting pattern was selected based on the geotechnical investigation data. For the supporting pattern, Type-1, 2, 3 were selected for each ground condition. A 3D numerical analysis model was developed for effective simulation as the room-and-pillar underground method consist of a pillar and room. As a review of numerical analysis, it was confirmed that the crown settlement, convergence, shotcrete and rock bolt were all stable in all supporting patterns. As a result of the analysis by the construction stage, it was confirmed that excessive stress was generated in the room when the construction stage of forming pillar. So, precise construction is required during the actual construction stage of the pillar formation.