A three-field based optimization formulation for flow simulations in networks of fractures on non-conforming meshes