使用Breakup Model下無法顯示破碎的拉格朗日粒子(child droplets)
-
Foamers您們好，
最近在做超音速下液體噴注在Ma>1 的crossflow 的penetration及SMD分佈validation("Numerical simulation of the gas-liquid interaction of a liquid jet in supersonic crossflow", Acta Astronautica 134 (2017))，目前卡在無論使用任何的breakup model都沒有辦法在paraview(VTK)顯示child droplets(如圖)
且打開lagrangian/sprayCloud/d 的粒子數量也沒有因為breakup而提升，唯能一確定的是drops有的確進行breakup和mass是有符合我給的初始條件(mass flow rate)，請問各位大大我的驗證計算裡面是哪裡設定不對或者有問題？謝謝！
Best regards,
Rick -
另外在 lagrangian/sprayCloud/nParticle中能看到每一粒的值為1~74000，是否代表我在paraview看到的只是Parcels而不是全部的particle?
-
@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
9/27