使用Breakup Model下無法顯示破碎的拉格朗日粒子(child droplets)

[星星星星晴]

@r07522121 所以要手动算一下先，我一直用flowprofile 1，这部分到底有没有影响我也不太了解 哈哈哈哈

[星星星星晴]

@r07522121 另外建议你先搞一下sprayparcel.c 这块，add parcel 就是一个pointer 把child的data传过去 然后 td.cloud.addparcel（）差不多这样，添加child 是在sprayparcel中处理的

[r07522121]

@星星星星晴 好的，我試試Info一下到底哪裡沒跑到，因為看完幾個檔案及sprayparcel.C裡面是有植入addparcel的代碼。貌似calcBreakup是由bool update()決定而update()是跟bool addParcels而決定，最後addParcels是由KHRT.C 裡面的if (nChildDrops >= nParticle)而決定，也許是criterion沒到的原因，可是我調到1還是沒有幫助哭哭，我再研究研究看看。謝謝你!

關於flowRateProfile，我也不確定那個值是不是有用的，之前在CFDONLINE看到有人說設為1 是initialize的作用。

[星星星星晴]

@r07522121

我想起来之前做过的一个Crossflow 的case，我隐约记得是可以显示child的。
不过我没达到你的超音速的flow。
而且你计算y方程了没有，y和yDot。
我建议你可以单独只喷射一个parcel，然后看破碎的情况~
祝你好运~

[r07522121]

@星星星星晴 你好，發現是y跟yDot=1e-38，所以原來沒有計算到...可是奇怪我有把solvOscillationEqn調到yes。還是那兩項是只有TAB在算?
前輩的案例看起來的確有breakup，有用primary breakup model嗎? 謝謝分享。

[星星星星晴]

@r07522121
不是只有TAB算，y和ydot是计算droplet的形变的，我这也是好久之前的case了，我当时是根据TAB扩展了一个TAB模型，后来就没有再继续做了

Dahms R N, Oefelein J C. The significance of drop non-sphericity in sprays[J]. International Journal of Multiphase Flow, 2016, 86: 67-85.
你可以看看这个文章，我当时就是根据这个做的

  ReitzDiwakarCoeffs
    {
        solveOscillationEq yes;
        Cbag            6;
        Cb              0.785;
        Cstrip          0.5;
        Cs              10;
    }


    ReitzKHRTCoeffs
    {
        solveOscillationEq yes;
        B0              0.61;
        B1              40;
        Ctau            1;
        CRT             0.1;
        msLimit         0.2;
        WeberLimit      6;
    }

    TABCoeffs
    {
		solveOscillationEq true;
        y0              0;
        yDot0           0;
        Cmu             10;
        Comega          8;
        WeCrit          12;
		Cb				0.5;
		SMDCalculationMethod method2;
		NChild			16;
    }

[星星星星晴]

@r07522121 你最好看一下模型的.H，看看都有什么东西可以定义