Functional Relationship to Describe Drains with Entrance Resistance

Numerical flow models usually represent drains as a system dependent boundary condition. If soil is saturated, drains act as the Dirichlet boundary condition with pressure head set equal to zero, and if soil is unsaturated, drains act as the Neumann boundary condition with flow set equal to zero. The underlying assumption is that drains exhibit ideal behavior. In reality, however, this is generally not so, and the flow encounters additional resistances due to pipe slotting and clogging of the envelope material around the drains. To account for the resulting resistance, a Hooghoudt-type boundary condition was developed that prescribes drain flow in relation to the groundwater level at a reference position. The measured drain discharge in an old drainage system was compared with calculated discharge assuming an ideal drain. It was found that the ideal drain assumption led to large errors in simulated discharge. With a correctly formulated and calibrated Hooghoudt boundary condition, however, more accurate drain discharges were obtained.