Physics-constrained polynomial chaos expansion for scientific machine learning and uncertainty quantification