CFD中文网

@yhdthu 因为用了moving reference frame (MRF) 来调整气泡纵向的速度，让它保持在空间域的中心，所以没用interIsoFoam。不过听了你说的，我会用interIsoFoam直接试试，谢谢:happy:

@yhdthu fvSolution里对alpha是如下这样设置的：

   alpha.water
    {
     interfaceMethod "isoAdvector";
     isoFaceTol  1e-8;
     surfCellTol 1e-8;
     nAlphaBounds 3;
     snapAlphaTol 1e-12; 
     clip   false; 
     writeIsoFaces false;
     gradAlphaNormal false;
     nAlphaCorr      2;
     nAlphaSubCycles 1;
     cAlpha          1;
    }

对PISO的设置是

{
	pRefPoint (0 0 0);
	pRefValue 0;
    momentumPredictor true;
    nCorrectors     3;
    nOuterCorrectors	1;
    nNonOrthogonalCorrectors 0;
}

controlDict里

maxCo           0.05;
maxAlphaCo      0.05;
maxDeltaT		1e-5;

求解器其实没用interIsoFoam，是以interFoam为基础进行了修改，在src里加了isoAdvection；压力速度耦合一开始想用PIMPLE算，但有问题，就改成了PISO。您在改isoAdvection时遇到过哪些影响很大的设置呀？

@李东岳 模拟的是一个二维算例，80mm*440mm,一共有30660个网格，感觉这个并不大，算着算着气泡就消失了:sad: ...不知道如果我增加nAlphaSubCycles能不能对alpha的运算结果有点好的影响
哎，三维算例很大，运算很慢，alpha有类似的问题...

Emmm....正在尝试模拟在水中上升的气泡，处理边界使用了isoAdvection，压力速度耦合使用了PISO，求解器以InterFoam为基础。图片是alpha.water的截图，还想请问为什么使用isoAdvection时气泡会发散的这么严重?前辈们有没有什么解决的建议？
P.S. 在使用MULES是气泡的边界会稳定很多