• Proceedings of the Korea Water Resources Association Conference
• /
• 2004.05b
• /
• pp.552-556
• /
• 2004

프랙탈 이송확산방정식은 정수 차수의 미분연산자로 구성된 고전적인 이송확산방정식과 비교하여 프랙탈 차수의 미분연산자로 구성된 보다 상위개념의 방정식으로써 정의된다. 지금까지의 프랙탈 이송확산방정식은 추계학적인 기법을 동원하여 푸리에-라플라스 공간에서 주로 해석되었으나, 본 연구에서는 실제 공간에서 유한차분개념을 도입하여 보다 직접적으희 하천에서의 오염물 이송확산에 관한 지배방정식을 유도하였다. 이러한 개념의 유도방법은 프랙탈 차수 및 관련 확산계수의 물리적인 추정에 관한 실마리를 제공할 수 있다. 고전적인 이송확산방정식과는 달리 프랙탈 이송확산방정식은 실제 하천에서 관측되는 오염물의 시간-농도 분포곡선의 왜곡현상과 분포곡선의 전후방부 농도를 보다 실제에 가깝게 모의할 수 있을 것으로 기대되어진다.

• Journal of Korea Water Resources Association
• /
• v.37 no.11
• /
• pp.889-896
• /
• 2004

The fractal advection-diffusion equation (ADE) is a generalization of the classical AdE in which the second-order derivative is replaced with a fractal order derivative. While the fractal ADE have been analyzed with a stochastic process In the Fourier and Laplace space so far, in this study a fractal ADE for describing solute transport in rivers is derived with a finite difference scheme in the real space. This derivation with a finite difference scheme gives the hint how the fractal derivative order and fractal diffusion coefficient can be estimated physically In contrast to the classical ADE, the fractal ADE is expected to be able to provide solutions that resemble the highly skewed and heavy-tailed time-concentration distribution curves of contaminant plumes observed in rivers.

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

기후변화 및 해안 구조물의 증가 등 여러 원인이 연안침식 및 해안선 변화와 같은 연안의 지형변화를 가속하고 있다. 빠르게 변화하는 연안의 지형변화예측 및 대응책 강구를 위해서는 연안의 유사이송 현상에 대한 신속한 예측이 필요하다. 본 연구에서는 GPU 엔진 기반 파랑해석모형인 Celeris Advent를 활용하여 실시간으로 연안의 유사이송 모의가 가능한 수치모형을 개발하였다. Celeris Advent는 GPU의 병렬코어를 활용해 실시간 연산과 GUI를 통한 사용자와의 실시간 상호작용이 가능한 모형이다. 지배방정식은 확장형 Boussinesq 방정식에 유사이송방정식을 양방향 결합하여 구성하였고, 지배방정식에는 하이브리드 유한체적-유한차분 수치기법을 적용하여 이송항은 유한체적법(Kurganov & Petrova, 2007), 소스항은 유한차분법을 통해 이산화하여 해석한다. 유사이송방정식은 수심적분형 이송확산방정식에 침식 및 퇴적 플럭스를 반영하는 소스항을 결합하여, 이송항 및 확산항을 통해 유사의 이송/확산을 고려함과 동시에 소스항을 통해 하상과의 상호작용을 고려하였다.

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

본 연구에서는 GPU 가속화 기반의 Boussinesq 모형인 Celeris Advecnt에 수심 적분된 2차원 이송-확산방정식을 추가하여 인터렉티브 시스템 기반의 추적자 이동 모형을 개발하였다. Celeris Advent는 최초로 개발된 인터렉티브 시스템을 갖춘 Boussinesq 모형으로, 시뮬레이션 중에 사용자가 모형의 파라미터뿐 아니라 모델 도메인 내 수위 및 수심을 바꿀 수 있다. 이를 통해 사용자는 모의가 진행되는 도중에 모델의 안정성 및 효율성을 위해 시간 간격을 조정할 수 있을 뿐 아니라 방파제 설치 등과 같은 지형 변화를 고려하기 위해 도메인 내 격자별 수심을 조정할 수 있다. 본 연구에서는 연안에서의 추적자 이동 모의를 위해 Boussinesq 방정식과 더불어 이송-확산방정식을 풀이하는 추적자 이동 모형을 개발하였다. 추적자의 확산항의 경우 분자 자체의 확산과 더불어 쇄파에 따른 난류 확산을 고려하였다. 난류 확산계수는 슈미트 수를 1로 두어 와동점성계수와 동일하게 두었으며, 와동점성은 단순화된 형태의 쇄파모형을 고려하여 계산하였다. 쇄파모형의 고려로 인해 이송-확산방정식과 더불어 운동량 방정식에서도 쇄파에 따른 운동량 소산이 고려되었다. 마지막으로, 추적자 농도에 대한 인터렉티브 시스템을 추가하여, 모델 구동 중에도 사용자가 수심적분된 추적자 농도를 조정할 수 있도록 하였다. 기수행된 2개의 수리실험 조건과 관측값을 이용하여 벤치마크 테스트를 수행하였으며, 관측값과 대체로 일치하는 것을 확인하였다.

• 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)를 산정하여 수치모형의 성능을 정량적으로 평가하였으며, 비소산함수를 적용함에 따라 스칼라의 수치확산이 감소하게 되었음을 확인하였다. 또한, 세 차례의 벤치마크 테스트 결과는 공통적으로 수치모형에 의해 계산된 결과값이 비소산함수를 고려함에 따라 해석해와 잘 일치함을 확인하였다.

• Water for future
• /
• v.27 no.2
• /
• pp.155-166
• /
• 1994

Various Eulerian-Lagrangian numerical models for the one-dimensional longitudinal dispersion equation are studied comparatively. In the model studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing adveciton and the other dispersion. The advection equation has been solved using the method of characteristics following fluid particles along the characteristic line and the results are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpolation polynomials are superior to Lagrange interpolation polynomials in reducing dissipation and dispersion errors in the simulation.

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

하천으로 유입된 오염물질의 거동을 정확하게 예측하는 것은 하천 시스템의 수질 유지관리에 매우 중요하다. 본 연구에서는 1차 감쇠율(decay rate)을 가진 비보존성 오염물의 비정상 이송해석 방정식의 해를 위해 유한체적기법이 개발되어졌다. 하천 흐름 해석을 위해 자연형 단면에서의 마름-젖음 해석이 가능한 기법이 도입되었다. 이 기법은 2차-정도의 정확성와 Courant 수가 1 까지 안정함을 보장한다. 도입된 기법은 Godnov 형의 유한체적기법을 이용하여 St. Venant 방정식들을 해석하였고 질량 및 운동량 플럭스는 Roe 형의 Riemann Solver 를 사용하여 연산하였다. 오염물의 이송 해석은 추가적인 이송-확산 방정식을 도입을 통해 기존의 St. Venant 방정식과 함께 풀려질 수 있다. 추가된 방정식과 St. Venant 식은 3×3 eigenstructure를 구성하였고 이는 2차원 흐름해석 기법과 유사하게 해석될 수 있었다. 본 연구 모형의 검증을 위해 오염물의 계속적 주입을 가정한 가상 및 실제 하천에 적용되었다. 연구된 기법은 모든 적용에서 합리적 정확도를 가지고 오염물질의 연속적인 특성을 잘 모의하고 있었다.

• Water for future
• /
• v.29 no.4
• /
• pp.187-198
• /
• 1996

A numerical model for unsteady dispersion of horizontal line source in turbulent shear flow is developed. A fractional step finite difference method is used which splits the unsteady two-dimensional advective diffusion equation into the longitudinal advection and the vertical diffusion equations, and solves them alternately for half time intervals by the Holly-Preissmann scheme and the Crank-Nicholson scheme, respectively. The developed numerical model is verified using a semi-analytic solution for steady dispersion in turbulent shear flow. Dispersion of an instantaneous plane source in turbulent shear flow is analyzed using the model. The degree of mixing at the same dimensionless time is almost the same regardless of the friction factor, and the travel distance required to reach a certain degree of mixing is inversely proportional to the square root of the friction factor.

• Proceedings of the Korea Water Resources Association Conference
• /
• 1987.07a
• /
• pp.201-214
• /
• 1987

하천에서 종, 횡방향이 고려된 확산-이송 현상을 수치기법인 유한요소법을 이용하여 2차원으로 해석하였다. 유한요소법으로는 Galerkin의 가중잔차 방법을 수십에 대해 적분을 취한 연속, 운동량 및 확산-이송방정식에 적용하였고, 선형보간함수와 선형삼각형요소가 이용되었다. 모형의 타당성을 입증하기 위해 단순화된 1차원 수로에서 수차례 검정한 결과 정확해와 거의 일치하는 만족할만한 결과가 도출되었다. 개발된 모형의 실험이 2차원수로에서 행하여져 지류의 유입에 따른 확산-이송현상이 모의되었으며, 실험적용은 개발사업후의 한강본류 9km 구간에 적용되어 탄천과 중량천의 지천 영향을 받는 오염 농도가 2차원적으로 도시되었다.

• Water for future
• /
• v.28 no.5
• /
• pp.151-162
• /
• 1995

A two-dimensional stream tube dispersion model is developed to simulate accurately transverse dispersion processes of pollutants in natural channels. Two distinct features of the stream tube dispersion model derived herein are that it employs the transverse cumulative discharge as an independent variable replacing the transverse distance and that it is developed in a natural coordinate system which follows the general direction of the channel flow. In the model studied, Eulerian-Lagrangian method is used to solve the stream tube dispersion equation. The stream tube dispersion equation is decoupled into two components by the operator-splitting approach; one is governing advection and the other is governing dispersion. The advection equation has been solved using the method of characteristics and the results are interpolated onto Eulerian grid on which the dispersion equation is solved by centered difference method. In solving the advection equation, cubic spline interpolating polynomials is used. In the present study, the results of the application of this model to a natural channel are compared with a steady-state flow measurements. Simulation results are in good accordance with measured data.