I'm Trying to model the natural convection of a PCM. Since i have to take account of the density inversion of this PCM, i cannot use Boussinesq Model in ANSYS fluent with a constant thermal expansion coefficient, because then i cannot use the function for the density. Besides this, using the temperature dependent properties for all terms of the momentum equation is too complicated for my case. I can't simply use a temperature dependent density for the whole momentum equation.
Is there a way or a UDF that can help me to use a specific function for density which only affects the buoyancy term In ANSYS FLUENT?