In this Equation, the density in the buoyancy force term conforms to the Boussinesq approximation, which can be considered to be a function of the temperature.
I have this information but i dont know I should use piecewise linear or another method for Boussinesq approximation.