Numerical validation of first order theory in bimodal porous formations