add thermaldynamic_equations class

prepare to calculate mass flow of each station, but still a lot of troubles