non Darcy flow simulation

Dear all,

Does the current code can simulate the non-Darcy flow in parallel fracture when the velocity get high? What I do now is simply increasing the Reynolds number which cannot reflect the non-Darcy flow character, perhaps it is wrong. any suggestions?

Thanks a lot