Диссипативная функция многофазной гетерогенной среды
       2. Вывод выражения для изменения энтропии системы

     Подставим в формулу (3.2) соотношения Гиббса для каждой из фаз (3.1), а в них - уравнения изменения внутренней энергии (2.8), (2.9) и концентраций компонентов (2.10) и, приведя подобные члены, получим:
  (3.3)
где j12 - поток массы между сплошной фазой и r-фазой (для процесса кристаллизации );  q12 - поток тепла между сплошной фазой и r-фазой, не связанный с фазовыми переходами; J1k - диффузионный поток k-го компонента в сплошной фазе; q1 - поток тепла в сплошной фазе за счёт процесса теплопроводности; f12 - сила взаимодействия между сплошной фазой и включениями, возникающая вследствие действия сил трения при контакте фаз; - тензор вязких напряжений; - тензор скоростей деформаций несущей фазы; i2 - энтальпия r-фазы.

     Проинтегрируем соотношение (3.3) по объёму V, занимаемому смесью, используя формулу Гаусса Остроградского:
 
где А - вектор; Аn - проекция вектора А на нормаль к поверхности FS..
     При этом получим выражение для изменения энтропии всего объёма V рассматриваемой системы:
  (3.4)