Мне кажется, что если если нам "понятно" распределение Больцмана, то объяснение конвекции можно свести к вопросу о том, почему возникает макроскопическое движение газа из области высокого давления в область более низкого давления.
Поясню, что я имею в виду. Предположим, что в объеме холодного газа с температурой

имеется область горячего газа с температурой

( окруженная со всех сторон холодным газом), и на некоторой высоте давления холодного и горячего газа равны:

(газ считаем идеальным). В холодной и горячей областях устанавливается распределение Больцмана с соответствующей температурой. Градиент давления на уровне, где давления равны (для холодного воздуха)

. Расстояние

отсчитываем от уровня, где давления равны. Назовем его нулевым для краткости. Аналогично для горячего воздуха получим

, где

- плотность газа,

- концентрация молекул.. Приращение давления при отходе от нулевого уровня равно:

. Очевидно

.
Таким образом, на границе областей выше нулевого уровня больше давление горячего газа, а ниже нулевого уровня - больше давление холодного газа. Поэтому объем горячей области выше нулевого уровня будет увеличиваться, а ниже нулевого уровня уменьшаться, это, собственно, и означает движение горячей области вверх.
Похожим образом можно объяснить неустойчивость горячего слоя газа, над которым находится слой холодного газа. Для этого нужно рассмотреть, как меняется небольшая выпуклость вверх (или вниз), образовавшаяся на границе раздела слоев.