关于晃荡惯性力的植入

[李东岳]

这玩意挺有意思。不过你的方程里面epsilon、omega、Ur之类不知道是啥，有相关sci可以参考么

[shuo]

@李东岳 李老师您好！epsilon是液舱晃荡的角加速度，omega是液舱晃荡的角速度，ur是相对于液舱的速度，也就是以液舱为参考系的速度。

原文：
1-s2.0-S0360544223003298-main.pdf

[李东岳]

在循环前植入：

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

[李东岳]

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

这个方程里面非粗体的$\omega ,\epsilon$以及${\theta }_m,T$你知道是什么么

[shuo]

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

在循环前植入：

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

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

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

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

[shuo]

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

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

这个方程里面非粗体的$\omega ,\epsilon$以及${\theta }_m,T$你知道是什么么

老师我的理解是，$\omega ,\epsilon$是标量表示的角速度和角加速度，文章里是把惯性力以分量的形式添加的，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);

这时候有晃荡现象了

是因为重力作用抵消了惯性力的作用吗？
现在的模型没有表面张力系数sigema给的是0，使用的是层流laminar，和真实的横摇运动还有一定差距，请问李老师还有什么改进的建议吗 ？

[李东岳]

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

[李东岳]

怎么样 后来有进展么

[shuo]

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

[李东岳]

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