شبیه‌سازی سقوط دو ذره صلب دایروی در کانال عمودی: ترکیب روش مرز غوطه‌ور-شبکه بولتزمن و روش اجرا گسسته

در این مطالعه، روش مرزغوط هور – شبکه بولتزمن به عنوان حلگر ناحیه سیال با روش اجزا گسسته به عنوان شبیه ساز برخورد ذرات ترکیب می‌شود. حاصل این ترکیب ایجاد مدل عددی توانمند خواهد بود که قادر به تحلیل جریانهای ذرهای با دقت مرتبه دو می‌باشد. در روش اجزا گسسته جهت تحلیل نیروی برخورد، به دو ذره اجازه داده می‌شود که در یکدیگر نفوذ کنند که این امر باعث ایجاد یک مدل واقعیتری از برخورد نسبت به مدلهای قبلی می‌شود. روش ترکیبی ذکر شده به طور جداگانه با شبیه سازی برخورد خشک (عدم حضور سیال) دو ذره صلب و همچنین شبیه سازی سقوط تک ذره صلب دایروی در یک کانال عمودی صحت سنجی می‌شود. در نهایت سقوط دو ذره صلب دایروی در کانال عمودی مورد مطالعه قرار خواهد گرفت و تاثیر پارامترهای فیزیکی موثر بر رفتار ذرات در حین برخورد، نظیر میزان صلبیت، ضریب بازگشت و اصطکاک بر نتایج بررسی می‌شود. در فصل نتایج نشان داده شده است که ضریب اصطکاک باعث فزایش مدت زمان کیسینگ می‌شود که در نتیجه باعث ایجاد تغییرات قابل توجهی در مسیر حرکت ذرات خواهد شد. همین‌طور نشان داده شد، که در این فیزیک خاص، ضریب بازگشت تاثیر چندانی بر رفتار ذرات حین برخورد نخواهد داشت.

Simulation of Two Circular Particles Falling in Vertical Channel: Combination of Immersed Boundary Lattice Boltzmann Method and Discrete Element Method

In this study, Immersed BoundaryLattice Boltzmann Method (IBLBM) as a fluid solver is combined with Discrete Element Method (DEM) as a collision model. The consequences of this arrangement go to a numerical great model (IBLBDEM) which is capable to simulate particulate flows with secondorder accuracy. To apply nonslip boundary condition, Eulerian velocities are interpolated in Lagrangian nodes using diffuse delta function. In DEM, two particles can penetrate to each other which this approach generates more realistic model rather other collision rules. Generally, in this model, the most important parameter is overlap distance between two particles which is directly related to amount of particles rigidity. The mentioned hybrid method is validated by simulation of drycontact of two particles and sedimentation of single particle in vertical channel, individually. Finally, sedimentation of two circular particles in vertical channel is studied and effects of physical parameters such as rigidity, restitution coefficient and friction coefficient in particles behavior has been investigated. Finally, it is shown that increasing friction coefficient leads to increasing in kissing time that causes a change in particles path. For this particular model, It is also dedicated that restitution coefficient does not have significant effect in particles behavior.