مدل‌سازی جریان دانه ‌ای در سیلو به روش نیمه ضمنی ذرات متحرک با تراک ضعیف (wcmps) با مدل رئولوژیکی هرشل بالکی نمایی

سیلوها در بسیاری از فرآیندهای صنعتی و مهندسی به طور گسترده مورد استفاده قرار می گیرند، بنابراین تخلیه مواد دانه ای و غلات از سیلوها به عنوان یکی از مهم ترین مسائل مورد توجه قرار گرفته است. با توجه به هزینه بالای مطالعات آزمایشگاهی بر روی مواد و شرایط مختلف، روش های محاسباتی به عنوان یک رویکرد جایگزین با هزینه بسیار کمتر مورد استفاده قرار می گیرند. به دلیل توانایی روش‌های لاگرانژی برای مدل‌سازی تغییر شکل‌ها و ناپیوستگی‌های بزرگ، در این مطالعه یک مدل لاگرانژی بدون شبکه بر اساس یک فرمول mps با تراکم ضعیف برای شبیه‌سازی تخلیه سیلو از مواد دانه‌ای توسعه و ارزیابی شده است. در این مطالعه ماده دانه‌ای به‌عنوان یک سیال ویسکو پلاستیک غیرنیوتنی در نظر گرفته می‌شود و یک مدل رئولوژیکی هرشل بالکی (h b) تنظیم‌شده نمایی برای مدل‌سازی رفتار دانه‌های استفاده شده است. قابلیت روش عددی توسعه‌یافته برای تخلیه مواد دانه ای از سیلوها ارزیابی شده است و با نتایج آزمایشگاهی و روش dem مقایسه شده است. مقایسه نتایج پروفیل سطح و پروفیل سرعت روش عددی توسعه‌یافته برای تخلیه دانه ها از سیلو با اندازه‌گیری‌های تجربی موجود و روش dem، قابلیت‌های مدل پیشنهادی را برای پیش‌بینی دقیق پروفیل‌های سطح و سرعت در این مسئله نمونه را نشان می‌دهد. 

modeling of granular flow in a silo using the weaklycompressible moving particle semi-implicit method (wcmps) with the exponential herschel-balky rheological model

silos are widely used in many industrial and engineering processes, so discharge grain from the silos is considered as one of the most important issues. due to the high cost of laboratory studies on different materials and conditions, computational methods are used as a substitute approach for much less cost. because of the ability of nonlinear lagrangian methods to model large deformations and discontinuities, this study develops and evaluates a mesh free lagrangian model based on a weakly compressible mps formulation to simulate of the discharge of the granular silo. in lagrangian methods, unlike the eulerian method, instead of networking the solution field and breaking the equations on the nodes, the solution field is divided into a number of particles and the broken equations are solved on these particles. in fact, the governing equations are transformed into particle interaction equations using different operators. in the meantime, the particles that are closer to the particle under study will have a greater effect on that particle. in such a way that the effect of relatively distant particles can be ignored in comparison with closer particles and the interaction between particles can be limited to a specific domain called the radius of effect. the effect of each particle on the calculated particle is measured by a weight function. in the wc mps method, the system is considered as a system with weakly compressibility and calculates the pressure of each particle using the state equation. in this study, the tait’s state equation is used, which is used for high pressure water flow. the mps method uses particle density to track the free surface. because there are no particles outside the free surface, the density of the particles on the free surface decreases sharply. a particle is known as a free surface particle whose density is somewhat lower than the standard particle density. the value of this limit may be selected from 80% to 99% depending on the problem therefore, the pressure of this particle on the free surface will be set to zero in each time step and in the mps method, and there is no need to apply any additional condition for the free surface. for solid (impermeable) boundaries, such as walls or beds, this boundary condition is applied. in the vicinity of solid boundaries, the particle density decreases, which can lead to computational disturbances. therefore, a number of ghost particles are located outside the boundaries to prevent this density reduction.  in this method the granular material is considered as a non newtonian visco plastic fluid and an exponentially herschel balky (h b) rheological model in combination with pressure dependent yield criteria model is employed to model non cohesive grain behavior. the ability of the developed numerical method to discharge grain from silos has been evaluated and compared with the experimental results and the dem method. comparison of the results of the developed numerical method for the discharge of grains from the silo with the available experimental measurements and the dem method shows the capabilities of the proposed model to accurately predict the surface and velocity profiles in this sample problem.