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