关于晃荡惯性力的植入

李东岳

在循环前植入：

我忽然注意到你写的这个。你要在循环中更新一下这个量。你更新没有。

李东岳

$$
\begin{array}{l}
\mathbf{F}_y(t)=-\rho\left(g \sin \theta(t)+y \omega(t)^2+z \varepsilon(t)+2 \omega(t) u_y(t)\right) \boldsymbol{j} \\
\mathbf{F}_z(t)=-\rho\left(g \cos \theta(t)+z \omega(t)^2-y \varepsilon(t)-2 \omega(t) u_z(t)\right) \boldsymbol{k}
\end{array}
$$

这个方程里面非粗体的$\omega,\varepsilon$以及$\theta_m,T$你知道是什么么

shuo

@李东岳在关于晃荡惯性力的植入中说：

在循环前植入：

我忽然注意到你写的这个。你要在循环中更新一下这个量。你更新没有。

之前没有更新，我在UEqn.H文件的开头又加了一行：

bf=(angulara^rotateR) + (angularv^(angularv^rotateR))+(2*angularv^U);

但是算出来后的结果还是像上面的图那样

shuo

@李东岳在关于晃荡惯性力的植入中说：

$$
\begin{array}{l}
\mathbf{F}_y(t)=-\rho\left(g \sin \theta(t)+y \omega(t)^2+z \varepsilon(t)+2 \omega(t) u_y(t)\right) \boldsymbol{j} \\
\mathbf{F}_z(t)=-\rho\left(g \cos \theta(t)+z \omega(t)^2-y \varepsilon(t)-2 \omega(t) u_z(t)\right) \boldsymbol{k}
\end{array}
$$

这个方程里面非粗体的$\omega,\varepsilon$以及$\theta_m,T$你知道是什么么

老师我的理解是，$\omega,\varepsilon$是标量表示的角速度和角加速度，文章里是把惯性力以分量的形式添加的，openfoam里的动量方程是矢量方程，我就直接把角速度和角加速度写成矢量，把惯性力直接加到方程里了，没有将其分解到y轴和z轴上。

李东岳

你那个bf带密度没有，把rho乘上看看

shuo

@李东岳李老师我重新检查了一下，之前bf是在求解UEqn时加入的，现在改成把rho*bf放在了方程UEqn里面。

bf=(angulara^rotateR) + (angularv^(angularv^rotateR))+(2*angularv^U);
fvVectorMatrix UEqn
(
    fvm::ddt(rho, U)
  + fvm::div(rhoPhi, U)
  - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
  + turbulence->divDevRhoReff(rho, U)
  -rho*bf     //加在这里了                                                       
);

直接算还是不行，然后尝试在算例中关掉了重力

dimensions   [0 1 -2 0 0 0 0];
value             (0 0 0);

这时候有晃荡现象了
output_嗨格式压缩副本.gif
是因为重力作用抵消了惯性力的作用吗？
现在的模型没有表面张力系数sigema给的是0，使用的是层流laminar，和真实的横摇运动还有一定差距，请问李老师还有什么改进的建议吗？

李东岳

应该是重力作用大于惯性力的作用，
你得植入效果已经出来了，你核实一下这个乘的rho应该是水的固定值密度，还是水跟空气混合的密度

李东岳

怎么样后来有进展么

shuo

@李东岳不好意思李老师，中间去做了一些别的工作，没有怎么学习of，最近又回头看了一下这个代码，用这个方法和of的动网格oscillatingRotatingMotion计算了同一个工况，结果还是显示在相同的振幅和频率下，惯性力不能像动网格一样使液舱出现明显的晃荡现象。感觉在液舱尺度是实验缩尺尺度下，惯性力也不是很大

李东岳

我觉得你这个挺有意思的。初学可以练手用。但是个人觉得算法应该没有问题，出现不一致的情况可能还是哪里出现了小bug。

CFD中文网

关于晃荡惯性力的植入