administrators

坚果云竟然挂了？？？？

单独从理论来看，这是通过固相应力（粒子应力）实现的，跟双流体模型类似，看Andrew and O'Rourke大牛1996年的文章，里面所有stress $\tau$的地方，固相应力是一个矢量，这个量用于修正接近极限的粒子速度

https://www.jianguoyun.com/p/DRX1P0wQ9s3ZBhiKwKoD

href是动网格用的，pimple各种参数比如rDeltaTSmoothingCoeff是LTS时间步用的

很有可能，cyclicAMI有可能把粒子算丢了

ddtCorr这个是个老梗了，一直没时间深究，不过今天发现在某些算例下，区别很大

     fvc::flux(HbyA)
  +  fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, Uf)

其中第二行加与不加的效果如下，左侧是不加，右侧是加

https://codepen.io/pkra/pen/EPeKjo

@xiaolin

那如果是预定变化轨迹的网格高度的动网格变化是不是应该参照movingConeTopoFvMesh呢

是的，这样简单的多

另外我发现方程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

感谢提醒！ 已更新

\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}

多谢老铁对我资料的认可