Physically-constrained emulation of partial differential equations