Hard Constraint Guided Flow Matching for Gradient-Free Generation of PDE Solutions