使用Breakup Model下無法顯示破碎的拉格朗日粒子(child droplets)
-
@r07522121 你的意思是液滴分布和你预期的不同吗？应该分布更加散而广？
-
@r07522121 你可以增加初始液滴的入射速率parcelsPerSecond，拉格朗日方法本身就是靠离散的液滴入注来模拟连续的液体入注，只有入注的速率足够快才能模拟连续液体入注，很多文献都会研究不同的入注速率，然后说达到某个速率就行了，就像网格无关性一样
-
@r07522121
2.1.3.1 Statistical Convergence Study
In order to ensure a sufficient number of computational parcels are injected into
the simulation to statistically represent the spray, a statistical convergence study is
performed.
来自文章MODELING THE INFLUENCE OF
NOZZLE-GENERATED TURBULENCE ON DIESEL
SPRAYS
我的理解是必须要有足够多的粒子才能代表连续的的液滴，调节 parcelsPerSecond是最简单的方法 -
@yfclark 在 使用Breakup Model下無法顯示破碎的拉格朗日粒子(child droplets) 中说：
MODELING THE INFLUENCE OF
NOZZLE-GENERATED TURBULENCE ON DIESEL
SPRAYS謝謝!這篇好像蠻實用。不過還是疑惑如果要是改變ER的話要動哪一個值會比較恰當?
-
@r07522121 你好 你使用的是哪个破碎模型？KHRT还是TAB还是什么，不是所有的模型都会添加childparcel的
-
@r07522121
这你看一下code KHRT 是可以添加child的，其他的需要具体看破碎模型。
另外你的injection 应该 parcelBasisType 是什么？ -
@r07522121 在 使用Breakup Model下無法顯示破碎的拉格朗日粒子(child droplets) 中说：
另外在 lagrangian/sprayCloud/nParticle中能看到每一粒的值為1~74000，是否代表我在paraview看到的只是Parcels而不是全部的particle?
是的 paraview只能看到parcels 不是particles。
可以根据你的parcelBasisType
还有massflow rate 以及nparticle 算出流量什么的我现在用的是fixed，就是在injection model 中需要设定每个parcel有固定数量的particle，然后控制parcelpersecond 来控制最终射入的parcel number 以及根据设定的distribution model 可以手动计算出 massflow rate。 这样的好处是你可以设定固定的 nparticle 就是计算麻烦点。
-
@星星星星晴 謝謝回覆，
剛看了一下KHRT的確是有addParcel這一項，可是沒有在後處理看到有child parcel的出現目前是使用mass, 然後用flowRateAndDischarge，調高parcelPerSecond之後，nParticle有明顯減少到MAX~500左右，結果還不錯。可是還是覺得沒有解決根本的問題...
fixed->nParticle 5,我有嘗試過，把parcelPerSecond調很高之後再把flowrateprofile調到很低。可是massflowRate會固定被提高到爆掉(發散)
補充一下我是用RosinRammler從最大5e-4,最小1e-6,d 5e-4, spread parameter 3.
-
@r07522121 嗯 所以不知道你的code 到底运行到那边没有 可以试试在你的code中输出一些到terminal put banana，是不是没有到addparcel的那一步。
你有没有解 y equation
-
@r07522121 所以要手动算一下先，我一直用flowprofile 1，这部分到底有没有影响我也不太了解 哈哈哈哈
-
@r07522121 另外建议你先搞一下sprayparcel.c 这块，add parcel 就是一个pointer 把child的data传过去 然后 td.cloud.addparcel（）差不多这样，添加child 是在sprayparcel中处理的
-
@星星星星晴 好的，我試試Info一下到底哪裡沒跑到，因為看完幾個檔案及sprayparcel.C裡面是有植入addparcel的代碼。貌似calcBreakup是由bool update()決定而update()是跟bool addParcels而決定，最後addParcels是由KHRT.C 裡面的
if (nChildDrops >= nParticle)而決定，也許是criterion沒到的原因，可是我調到1還是沒有幫助哭哭，我再研究研究看看。謝謝你!關於flowRateProfile，我也不確定那個值是不是有用的，之前在CFDONLINE看到有人說設為1 是initialize的作用。
-
我想起来之前做过的一个Crossflow 的case，我隐约记得是可以显示child的。
不过我没达到你的超音速的flow。
而且你计算y方程了没有，y和yDot。
我建议你可以单独只喷射一个parcel，然后看破碎的情况~
祝你好运~
-
@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; }
16/27