FvMatrix::H

Материал из ru.wiki.laduga.ru
Перейти к: навигация, поиск


Реализация оператора H operator, недиагональной функции, используемой при линеаризации конвективного члена. Уравнение импульса может быть выражено следующим образом:

$ A_p \mathbf{U}_p = \mathbf{H}(\mathbf{U}) - \Delta p $

Где:

  • p - индекс объёма ячейки;
  • A - коэффициенты дискретизации;
  • U - скорость;
  • p - давление;
  • H - задаётся следующим образом:

По индексу грани:

$ \mathbf{H}_{o,i} = A_{n,i} \boldsymbol \psi_{n,i} + \mathbf{S}_i $

$ \mathbf{H}_{n,i} = A_{o,i} \boldsymbol \psi_{o,i} + \mathbf{S}_i $

Где:

  • A - матричные коэффициенты matrix coefficients;
  • o и n - индексы "владельца" и "соседа" owner and neighbour indices;
  • i - индекс грани;
  • $ \boldsymbol\psi $ - переменная поля fvMatrix (обычно скорость, U);
  • S - дискретизация слагаемого источников.

По индексу объёма ячейки:

$ \mathbf{H}_p = \sum\limits_r A_r \boldsymbol \psi_r + \mathbf{S}_p $

Где:

  • r связанные ячейки[1].

Вычисление суммы связанных ячеек выполняется lduMatrix::H.

Ссылка на оригинал - http://openfoamwiki.net/index.php/OpenFOAM_guide/FvMatrix;;H

Примечание

  1. То, что мы называем связанные ячейки условно имеем в виду соседние ячейки. Но OpenFOAM имеет другое значение для соседей (neighbours), поэтому мы используем термин связанные ячейки.