适用范围：用于不可压，层流，牛顿流体的瞬态求解器
控制方程：
$$
\nabla\cdot\mathbf u=0\tag{1}
$$
$$
\rho\frac{\partial \mathbf u}{\partial t}+\rho\nabla\cdot{\partial (\mathbf u\mathbf u)}=-\nabla p+\mu\nabla\cdot(\nabla\mathbf u)\tag{2}
$$
在不可压缩流体的计算中，通常将压力定义为$p/\rho$
之后便是采用之前提到的PISO算法计算，此处不再赘述
$$
{ c } \frac { ( \rho \mathbf { v } ) ^ { n + 1 } - ( \rho \mathbf { v } ) ^ { n } } { \Delta t } + C \left( \mathbf { v } ^ { n + 1 } \right) = L \left( \mathbf { v } ^ { n + 1 } \right) - G \left( p ^ { n + 1 } \right) \tag{1}
$$
$$
\frac { 3 ( \rho \mathbf { v } ) ^ { n + 1 } - 4 ( \rho \mathbf { v } ) ^ { n } + ( \rho \mathbf { v } ) ^ { n - 1 } } { 2 \Delta t } + C \left( \mathbf { v } ^ { n + 1 } \right) = L \left( \mathbf { v } ^ { n + 1 } \right) - G \left( p ^ { n + 1 } \right)\tag{2}
$$
第一个方程：隐式Euler方案；第二个方程：时间导数近似
备注：本文主体内容为[https://krasjet.github.io/quaternion/quaternion.pdf] 的归纳，如有侵权，请联系本人删除。
本文介绍了标量运输求解器scalarTransPortFoam
$$
\frac {\partial(\rho Y)}{\partial t}= \omega_i(Y_i,T)
$$
$$
h=u_0+\frac p \rho+\int_0^t\frac q \rho d\tau=u+\frac p \rho
$$
$$
p = \sum p _ { i } = \sum n _ { i } \frac { R T } { V } = \sum \frac { m _ { i } } { M _ { i } } \frac { R T } { V } = \sum \frac { y _ { i } } { M _ { i } } \underbrace { \frac { m _ { t o t } } { V } } R T
$$
转载自：C++ 中的 mutable 关键字
此篇介绍 C++ 中的mutable关键字。