Numerical and Analytical Modeling of Flow Partitioning in Partially Saturated Fracture Networks