By setting exergy destruction rate as an optimization objective and fluid power consumption as a constraint condition, a momentum equation with an additional volume force is constructed through functional variation to numerically simulate
convective heat transfer in coupling with energy equation