Optimal operation of large district heating networks through fast fluid-dynamic simulation