CFD中文网

@Samuel-Tu
記得好像是用巨集 #calc

試試U,p參考airfoil2D算例的邊界條件才給?

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

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

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

@星星星星晴 謝謝回覆，
剛看了一下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.

@星星星星晴 你好，我使用跟文獻一樣的KH-RT模型。詢問一下用哪一個模型會添加子粒子?

@yfclark 在 使用Breakup Model下無法顯示破碎的拉格朗日粒子(child droplets) 中说：

MODELING THE INFLUENCE OF
NOZZLE-GENERATED TURBULENCE ON DIESEL
SPRAYS

謝謝!這篇好像蠻實用。不過還是疑惑如果要是改變ER的話要動哪一個值會比較恰當?

@yfclark 感謝指點，這點對我來講真的是新知識。

想確認一點，假如我之後要處理的問題是liquid燃燒，那入注速率parcelsPerSecond應該是dependent of Equivalence ratio所以不能隨意增加，這樣說法是沒錯嗎？

@yfclark 是的，這樣描述更貼切。所以我很好奇lagrangian/sprayCloud/nParticle值的意思是不是代表真實particle數量。

@东岳 目前的嘗試也是打算增加parcelsPerSecond, 但這樣會跟文獻的ER多許多。

不好意思是輸入法問題，我是對岸的研究生。

@yfclark 你好，從最小的粒子直徑上來說粒子是有破碎，但是後處理的SMD只有沿著主粒子路徑有值，並沒有像文獻中的這樣。

另外我沒有使用primary atomization model, 只有使用secondary breakup(Pilch-Erdman), 粒子初始直徑為0.5mm(assumed跟噴口直徑一樣)

@东岳 你好，我在paraFoam看所顯示的粒子數量一樣少，只是d沒辦法像foamToVTK用paraview的scale粒子大小。我懷疑是設定parcelsPerSecond太小，可是透過換算出來的質量流率卻是如此(大概每秒10萬個parcels)，現在很不知所措...

另外在 lagrangian/sprayCloud/nParticle中能看到每一粒的值為1~74000，是否代表我在paraview看到的只是Parcels而不是全部的particle?

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

@刘雄国 您好，

是的因為粒子源在某很小區域噴發出來，所以通常只會被一個或兩個processor domain所計算,拖慢整體計算效率。

請問一下你說的是reconstructMesh嗎? 可是我沒有使用Dynamic mesh。請問該如可implenment?

謝謝

大家好，先簡單說一下我的狀況：
最近在仿真超音速下sprayFoam的問題，固定Equivalence ratio情況下, 粒子在每個timeStep下的數量都非常巨量。但平行運算下即使我使用64核計算機計算情況下運算速度覺得不理想。後來發現如果parcel injection並非均勻的話(i.e. 某位置injection), 即使Mesh 被decompose成n個domain，大量什至全部的拉格朗日粒子仍然只會被一個processor所處理，Eularian-phase早就算完，處理器都在等某processorN計算粒子。

目前想到是手動decompose 整個computational domain(例如injection的附多切好幾個subdomain)，看到decomposeParDict有MANUAL的type，但需手動一個一個cell的label放在input file,而且對於cell number非常多的問題非常不實際。 想請問各位大大有解決方法或者想法嗎？

Regards,
Rick

顆粒和顆粒的碰撞次數應該要自己編，顆粒跟patch的碰撞次數倒是有