Sensitivity Analysis of Infiltration using a Mass Conservative Numerical Solution of Richards Equation

Water flow into unsaturated soils is most often modeled by Richards equation consisting of the mass conservation law and Darcy's law. Three standard forms of Richards equation are presented as the head (${\Psi}$)-based form, the moisture content (${\theta}$) based form, and the mixed form. Numerical solutions of these partial differential equations with highly nonlinear terms can cause poor results along with significant mass balance errors. The numerical solution based on the mixed form of Richards equation is known that the mass is perfectly conserved without any additional computational efforts. The aim of this study is to develop fully implicit numerical scheme of Richards equation for one-dimensional vertical unsaturated flow in homogeneous soils using the finite difference approximation, and then to perform sensitivity analysis of infiltration to the variations in the unsaturated soil properties and to different soil types.

Richairds 방정식의 질량보존적 수치해석 해법에 의한 침투량의 민감도분석

질량보존의 법칙과 Darcy의 법칙으로 표현되는 Richards 방정식은 비포화대의 토양수분흐름을 모의하는데 널리 사용되어 왔다. Richards 방정식은 압력수두의 항으로 표현되는 방정식, 토양수분의 항으로 표현되는 방정식, 그리고 이 둘을 혼합한 형태의 방정식 등, 세가지 형태로 표현할 수 있다. 고차의 비선형 항들을 포함하는 이 편미분방정식들을 수치해석방법으로 풀 때, 질량 비보존을 수반하는 오류의 결과가 초래될 수 있다. 세가지 방정식들 중 혼합형 Richards 방정식이, 다른 추가적인 계산없이 질량을 온전히 보존하는 것으로 알려져 있다. 이 연구의 목적은 동질성 토양에서의 1차원적 연직방향 비포화수 흐름모의를 위해, Richards 방정식의 질량보존적 수치해석법을 완전음해 유한차분법으로 개발하고, 이를 통해 민감도 분석을 실시하여 토양특성인자들과 토양종류에 따른 침투율의 변화를 살펴보는 데 있다.