Refactored thermal_buoyancy as thermal_expansion_coefficient