圆柱绕流问题

symon90

palabos软件算的，给出几个前几个速度云图。这是属于反射吗？请问怎么解决？

[ 本帖最后由 symon90 于 2014-5-5 15:17 编辑 ]

luo@odu.edu

This is typical for any method of artificial compressibility, and the LBE belongs to this category (it does not solve the Poisson equation). The acoustic waves in the system of is hard to kill, and that is what you are seeing.

The MRT-LBE can adjust the bulk viscosity to dissipate the acoustic waves. I suspect palabos use the LBGK algorithm, thus it cannot dissipate acoustics effectively.

onesupeng

恩，好像LBM启动计算时常常会遇见这种情况，演化一段时间大部分就会被耗散掉。

@罗老师，以后我用LBM可能会经常请教你，但是到时候用email沟通，你不介意吧？

symon90

谢谢罗老师的回复~您所提到用MRT-LBE，能够避免上面提到的波，也就是多松弛因子，那么，palabos里的LES模型，也是修正松弛因子，请问，使用palabos的大涡模拟能够和MRT-LBE一样避免上述情况嘛？刚接触LBM，不太懂，问题不专业还希望罗老师别见外~

symon90

您好，我越往后算，流场越乱。这个波在流场里来回“跑”~势必会对计算结果造成严重的影响，还请多多指教啊。

symon90

罗老师，您好，我今天用MRT试了下，还是和BGK出现的情况一样，描述下具体原因：palabos里给的圆柱绕流例子是泊肃叶流中的例子，我是只是将其中的速度改成了一个固定的值，就出现了上述情况~还请您多多指教

flowinflow

Is this high order result? Is there a non-reflecting b.c.?

symon90

您好，不是高阶格式，固面边界处理采用反弹边界处理。边界条件：速度边界条件，左侧速度入口，右侧出口。上下 freeslip或是默认都出现上述情况。

luo@odu.edu

I did not know you were doing LES. Truly non-reflective BCs at the outlet are difficult to achieve -- I don't know a perfect one.   But, something has to be done at outlet, i.e., using free-stream BCs.

symon90

谢谢罗老师的回答，我一直定义的出口为:boundary:utflow  .我再尝试下别的定义。

luo@odu.edu

Both inlet and outlet are important.

symon90

好的，谢谢罗老师的建议。我发现，这个波是启动计算时从圆柱传出来的，出不去，和圆柱形成的卡门涡街不同，计算过程中，形成的涡街是可以出去的。

[ 本帖最后由 symon90 于 2014-5-7 21:29 编辑 ]

symon90

就像上图 的 naca0012 绕流（高RE，200W）脱落的涡可以从出口边界跑出去的，但是由于启动时在翼型周围生成的“波”在流场中来回跑，出不去，等到计算一定时间后，整个流场就崩溃了，如图下。

luo@odu.edu

The acoustic waves are generated at inlet, it appears. They interact with the cylinder then go wild. There are a few ways to reduce (but not completely eliminate) this kind of effect.

1) correct BCs (inlet, outlet, at fluid-solid boundaries);
2) adjusting dissipation to dump out the acoustics in the initial stage;
3) use odd-number grid size.

MRT-LBE will help greatly.

通流

看起来像个不错的非定常计算。波出不去，只能是远方边界条件的问题。