求问有什么思路，可以给流体添加一个持续的，自定义分布的外力场？

[hotairballon1997]

现在我得到了一个不均匀分布的外力场，想查看流体在这个外力场持续的作用下的流向和流速变化情况。

我有尝试过将外力场转化为压力场p，并保持速度场U不变，然后进行模拟。但是因为这样的p和U并不是稳态，压力场p很快就变化掉了，不能持续。所以这个办法应该不可行。

想问问各位前辈，在fvOption里面能不能添加自定义分布的，并且持续的外力场？或者有其他什么思路可以实现？

[bestucan]

这个操作即使实现了，也很容易发散，
P=FV
恒定的力作用下，速度不断增加，功率不断增加，能超光速。

能量转换总有功率，功率一定，力保持不住才是正常的，速度上去了，力就下来了。

非要搞，最相近的大概是磁场了，磁流体，在磁场中受力。磁场不变。给拉格朗日粒子设定磁极强度，每次迭代按不均匀分布的力刷新拉格朗日粒子的磁极强度。

[hotairballon1997]

@bestucan 感谢回复，这么一说确实是，就没法达到稳定状态。

不过如果按照您说的，用拉格朗日粒子，如果我设置粒子为passive scalar，用icoUncoupledKinematicParcelFoam对一个不稳定的场来求解粒子，这样粒子速度也会不断增加，最终导致发散吗？

[bestucan]

emmm……超出我的理解范围了，不懂了。

我觉的还会发散，但是会可控一些。比如，这个是线性发散，直接给固定力，感觉是指数增长发散。我猜的

[OItoCFD]

有点不懂，难道不是直接加给体积力项不就可以了吗？

[船只]

有一个思路是可行的：找一个求解器，在其动量方程的中填加一个力F项（当做常数，不进行求解）进行编译；然后在算例0文件中加一个力F的边界条件设置，这里可以用每个网格的值来设置。根据你的要求，可以将外力作用边界面处的网格的F值设置为你想要的值，其它不相关网格的F值设为（0,0,0）就行了。
这样你的力在该位置就是持续作用的，而且不会随着流场的变化而改变（求解器中没有进行求解）
这样的难点在于，你需要知道各位置处流场网格的编号及其规律，然后根据网格编号来进行设置网格的F值

[Samuel-Tu]

@船只 如果网格不同可以用MapField先映射，然后读取F场

[船只]

@samuel-tu 这个可以有。如果可以的话用结构网格是最方便的了，直接在该区域定义一个block，单独进行网格划分