icoFoam中压力泊松方程



  • 在icoFoam中的压力泊松方程中
    f297e7ca-ac47-4411-b8ea-5df6e176af11-image.png

    可知(26)是通过(27)积分离散过来的

    下面进行了推导:
    3dc700d5-7329-4bbf-9514-2b1050a359a0-image.png

    我没弄明白的是
    051222da-7efa-479d-bb16-8c16f42dd79d-image.png



  • \begin{equation}
    \frac{1}{\Delta V}\sum p_f\bfS_f=\frac{1}{\Delta V}\int p \rd\bfS=\frac{1}{\Delta V}\int\nabla p\rd V
    \end{equation}



  • 另外我发现方程12-26都需要除掉网格体积$\frac{1}{\Delta V}$,https://www.cfd-china.com/topic/2639/请教-icofoam解析-中的问题 当时做这个改动的时候弄混乱了

    原来的方程13这样

    \begin{equation}\notag
    \left( {\frac{{{V_\mathrm{P}}}}{{\Delta t}} + \sum { {\frac{{F_f^t}}{2}} + \sum { {\nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} } } } \right)\mathbf{U}_\mathrm{P}^* = - \sum { {\left( {\frac{{F_f^t}}{2} - \nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} \right)\mathbf{U}_\mathrm{N}^*} } + \frac{{{V_\rP}}}{{\Delta t}}\mathbf{U}_\mathrm{P}^t-\sum \frac{{p_\mathrm{P}^t + p_\mathrm{N}^t}}{2}\bfS_f.
    \end{equation}

    要改为

    \begin{equation}\notag
    \left( {\frac{{{1}}}{{\Delta t}} + \frac{1}{V_\mathrm{P}}\sum { {\frac{{F_f^t}}{2}} + \frac{1}{V_\mathrm{P}}\sum { {\nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} } } } \right)\mathbf{U}_\mathrm{P}^* = - \sum \frac{1}{V_\mathrm{P}}{ {\left( {\frac{{F_f^t}}{2} - \nu \frac{{\left| \bfS_f \right|}}{{\left| \mathbf{d} \right|}}} \right)\mathbf{U}_\mathrm{N}^*} } + \frac{{{1}}}{{\Delta t}}\mathbf{U}_\mathrm{P}^t-\frac{1}{V_\mathrm{P}}\sum \frac{{p_\mathrm{P}^t + p_\mathrm{N}^t}}{2}\bfS_f.
    \end{equation}

    导致后面的都要做一些改动。应该是我写的一个bug

    感谢提醒!:xiexie: :xiexie: :xiexie: 已更新



  • @东岳 谢谢老师,还是我学艺不精~:xinxin:


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2