مقايسه روش‌هاي شبكه بولتزمن و تفاضل محدود در حل هم‌زمان معادلات آب زيرزميني و انتقال آلودگي

Comparison of Lattice Boltzmann and Finite Difference Methods in the Solution of Coupled Groundwater and Contaminant Transport Equations

روش شبكه بولتزمن یك روش عددی قدرتمند برای شبیه‌سازی جریان سیال است. برای بیان روابط حاكم بر جریان آب زیرزمینی در محیط متخلخل دو رویكرد اصلی وجود دارد. رویكرد اول استفاده از معادله ناویر استوكس و رویكرد دوم استفاده از معادله انتشار برای جریان می‌باشد. در این پژوهش برای جریان آب زیرزمینی از رویكرد دوم استفاده شده است. همچنین معادله حاكم بر انتقال آلودگی در محیط متخلخل، معادله انتقال-پخش می‌باشد. حل هم‌زمان معادلات جریان و انتقال آلودگی به وسیله روش شبكه بولتزمن در این پژوهش برای اولین بار صورت پذیرفت. همچنین از روش تفاضل محدود صریح و كرنك نیكلسون برای مقایسه با روش شبكه بولتزمن استفاده شد. نتایج این پژوهش نشان داد كه روش شبكه بولتزمن قادر به حل هم‌زمان معادله آب زیرزمینی و انتقال آلودگی با دقت بالایی است، به طوریكه دقت آن با روش كرنك نیكلسون برابر می‌باشد، اما سرعت اجرای آن بسیار بیشتر از روش كرنك نیكلسون است. در ضمن در مقایسه با روش تفاضل محدود صریح، روش شبكه بولتزمن در حل معادلات فوق دارای محدوده پایداری و سازگاری بزرگتری می‌باشد. در این رابطه عدد پكلت شبكه كوچكتر از 7 برای روش شبكه بولتزمن D1Q2 پیشنهاد می‌شود.

Lattice Boltzmann method is a powerful numerical method for the simulation of fluid flow. There are two approaches for the description of groundwater flow in porous media. In the first approach, Navier-Stokes equations are used, while diffusion equation is applied to describe groundwater flow in the second approach. In this research, groundwater flow is modeled using the second approach. Furthermore, the governing equation of contaminant transport is called advection-dispersion equation. Here for the first time, the coupled solution of groundwater flow and contaminant transport using Lattice Boltzmann method is performed. The results indicated that Lattice Boltzmann method is capable of solving groundwater and contaminant transport equations simultaneously with high precision. Moreover, it was found that although the accuracy of Lattice Boltzmann and Crank-Nicolson method are the same, the speed of Lattice Boltzmann method is much higher than finite difference methods. In addition, Lattice Boltzmann method has higher range of stability and consistency in comparison with explicit finite difference method. Regarding this issue, grid Peclet number smaller than 7 is recommended for D1Q2 scheme of LBM.