Dissipative partial differential equations for fluid flows of various types share a common mathematical structure derived from the generalized Onsager principle and variational principle.The combination of these two principles allows one to develop a general framework to design numerical algorithms for the dissipative partial differential equations.In this talk,I will discuss some recent development in energy stable numerical algorithms for various multiphase and complex fluid flow models.