Divergence-free neural operators for stress field modeling in polycrystalline materials