• 제목/요약/키워드: 1차원 이송-확산 방정식

검색결과 17건 처리시간 0.029초

하이브리드 FVM/FDM 기반의 2차원 흐름 및 스칼라 이송 모형 개발 (Development of 2DH hydrodynamic and scalar transport model based on hybrid finite volume/finite difference method)

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

  • PDF

수치기법을 이용한 확산 - 이송의 2차원 분석 (Two-Dimensional Analysis of Convection-Dispersion Using Numerical Schme)

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

  • PDF

Eulerian-Lagrangian 방법을 이용한 1차원 종확산방정식의 수치모형 (Numerical Modeling of One-Dimensional Longitudinal Dispersion Equation using Eulerian-Lagrangian Method)

  • 서일원;김대근
    • 물과 미래
    • /
    • 제27권2호
    • /
    • pp.155-166
    • /
    • 1994
  • Eulerian-Lagrangian 방법을 이용하여 1차원 종확산방정식의 수치모형을 비교·분석하였다. 본 연구에서서 비교·분석한 모형은 지배방정식을 연산자 분리방법에 의해서 이송만을 지배하는 이송방정식과 확산만을 지배하는 확산방정식으로 분리한다. 이송방정식은 특성곡선을 따라서 유체입자를 추적하는 특성곡선법을 사용하여 해를 구하고, 그 결과를 고정된 Eulerian 격자상에 보간하였고, 확산방정식은 상기 고정격자상에서 Crank-Nicholson 유한차분법을 사용하여 해를 구하였다. 이송방정식의 풀이에서 다양한 보간방법이 적용되었는데, 일반적으로 Hermite 보간다항식을 사용한 경우가 Lagrange 보간다항식을 사용한 경우보다 수치확산 및 수치진동 등의 오차를 최소화할 수 있어서 더욱 우수한 것으로 밝혀졌다.

  • PDF

Boussinesq 방정식 기반 인터렉티브 추적자 이동 모형 개발 (Development of interactive tracer transport model coupled with Boussinesq equations)

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

  • PDF

유한체적법을 이용한 하천 오염물의 이송해석 (Analysis of Pollutant tranport in the natural river using FVM)

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

  • PDF

점착성 부유사 이동에 관한 수치모형 (A Numerical Model for Cohesive Suspended Load Movement)

  • 안수한;이상화
    • 물과 미래
    • /
    • 제23권1호
    • /
    • pp.119-127
    • /
    • 1990
  • 점착성 부유사 농도분포는 해수유동과 물질 확산에 의해서 결정되며 지배방정식으로는 2차원 수심적분된 Reynolds운동방정식, 연속방정식과 Fick의 확산법칙에 근거를 둔 대류-확산방정식이 사용되었다. 해수유동과 점성퇴적물 확산인 두개의 모형은 유한차분법을 이용하였고 유동모형은 양해법, 확산모형은 다증법을 사용하여 부유사 이동의 현상을 파악하였다. 해수유동방정식의 적용시 이송항의 포함여부에 대해서 조사하였으며 물질확산 방정식에 대해서는 한계전단응력값의 변화가 부유사농도에 영향을 주는가에 대해서 비교하였다.

  • PDF

연산자 분리기법에 의한 1차원 이송-확산방정식의 수치모형 (Numerical Model of One-Dimensional Advection-Diffusion Equation Applying Split-Operator Method)

  • 이정규;강창구;이종인
    • 한국수자원학회논문집
    • /
    • 제30권2호
    • /
    • pp.143-154
    • /
    • 1997
  • 이송항에는 5차 보간다항식을 사용하는 Holly-Pressmann 기법을, 확산항에는 Hobson 등이 제안한 양해법을 사용하는 연산자 분리기법을 사용하여 1차원 이송-확산방정식의 수치모형을 제안하였다. 제안된 모형을 검정하기 위하여 일정한 유속과 종확산계수를 갖는 순간적으로 부하된 오염원의 경우와 상류단에 연속적인 오염원을 갖는 경우에 대하여 본 모형의 해를 해석해와 기존의 모형으로부터 구한 해를 비교검토하였다. Courant 수와 Peclet 수를 가진 경우에 대한 수치해석을 통하여, 본 모형이 Courant 수가 1보다 큰 경우에 대해서도 안정된 해를 제공함을 알 수 있었으며, 해석해가 존재하는 경우에 본 모형을 적용하여 얻은 수치해와 비교한 바 전반적으로 잘 일치하였다. 본 모형의 확산항에 사용된 양해법에서는 일반적인 양해법의 단점인 계산시간간격의 제약이 상당히 완화되어 상대적으로 큰 계산시간간격에 대해서도 양호한 결과를 보였다.

  • PDF

하천에서의 유사이동의 영향을 고려한 홍수위험도 분석 (Flood Risk Analysis Considering the Sediment Transport in a River)

  • 손인호;김병현;한건연
    • 한국수자원학회:학술대회논문집
    • /
    • 한국수자원학회 2015년도 학술발표회
    • /
    • pp.302-302
    • /
    • 2015
  • 본 연구에서는 PMP 조건하에서 농업용 목적의 시례저수지의 가상붕괴에 따른 홍수파 해석을 통해서 하류부에 대한 수리학적 분석을 실시하였다. 그리고, 사력댐으로 구성된 댐 제체의 붕괴로 인한 유사 이송을 모의하고 이로 인한 하류부에서의 영향을 분석하였다. 이를 위하여 경상남도 김해에 위치한 시례저수지 유역을 대상으로 수리 수문자료 및 하천현황을 상세히 조사하고 댐의 가상 붕괴로 인한 하류부 홍수파 해석 및 댐 본체의 유실에 따른 유사의 이송으로 인해 발생 가능한 하류부에서의 영향을 분석하였다. 시례저수지의 댐 형식은 흙댐으로써, 구성방식은 침하량이 적은 중심코아형이며 댐 높이는 34.3 m, 유역면적은 3.13 km이다. 본 연구에서는 하천에서의 흐름과 유사이송 계산이 가능한 1차원 수치 모형인 CCHE1D를 이용하였다. 적용 모형은 확산파 및 동역학적 방정식을 모두 고려한다. 또한, 비평형 이송 모형을 사용하여 유사의 발생과 퇴적을 산정하고 비균질 상태의 유사 이송을 계산 하는데 있어서 하천 단면형의 변화와 하상물질의 구성 상태, 제방침식, 수로 확폭의 과정과 함께 산정된다. CCHE1D는 하상물질의 공극률과 non-equilibrium adaptation length, mixing layer thickness 등의 여러 변수들을 현재까지 개발된 식들을 제공하며, 토사이동해석을 위한 공식은 SEDTRA 모듈, Wu et al(2000) 공식, 수정 Ackers and White(1973) 공식, 수정 Engelund and Hansen(1967)이 사용된다. 비균질 유사이송과 bed deformation, bed-material sorting은 완전연계과정(coupling)로, 유사와 흐름방정식의 계산은 반연계과정(semi coupling)을 적용하여 계산하였다.

  • PDF

토양 오염원의 이동에 관한 연구 (감쇠항이 있는 3차원 이송-확산 방정식의 수치모형 개발) (A Study on the Transport of Soil Contaminant (A Development of FDM Model for 3-D Advection-Diffusion Equation with Decay Term))

  • 김상준
    • 한국수자원학회논문집
    • /
    • 제45권2호
    • /
    • pp.179-189
    • /
    • 2012
  • 오염물질의 이동 현상을 모의하기 위하여, 감쇠항이 있는 3차원 이송-확산 방정식의 수치모형이 개발되었다. 개발된 모형은 유한차분 모형으로서 시간단계의 가중치 ${\alpha}$를 포함하는 음해법(implicit finite difference method)과, 반복법인 Gauss-Seidel SOR(successive over relaxation)이 사용되었다. 모형은 보다 단순화된 가정 하에서 존재하는 두 가지의 해석적인 해와 비교되었다. 그 결과 Peclet number가 5~20 이하에서는 수치 분산의 영향이 크지 않았고 작은 오차범위 내에서 해석적인 해와 동일하였다. 또한 가중치 ${\alpha}$의 변화에 대한 모형의 거동은 Crank-Nicolson 모형(${\alpha}$=0.5)이 fully-implicit 모형(${\alpha}$=1)보다 해석적인 해에 접근함을 보여주었다. 모형의 검증과 실효성 제고를 위하여, mass balance를 검토하였다. 즉, 이송, 확산 및 감쇠항 각각에 대한 질량 이동을 산출하였으며, 그 결과 질량 이동의 계산 오차는 약 3% 이내였다. 본 모형은 감쇠 과정이 수반되는 3차원 이송-확산의 농도분포와 질량이동을 산출할 수 있으며 다양한 경계조건을 설정함으로서 현장조건을 반영할 수 있다. 그러나본 모형은 고정격자를 기반으로하는 유한차분 모형이므로 Peclet number가 비교적 작게 나타날 수 있는 토양 및 지하수계의 오염물질 이동 등의 문제에서 유용하게 적용될 수 있을 것으로 사료된다.