Skip to main content

Truncation Errors, Exact And Heuristic Stability Analysis Of Two-Relaxation-Times Lattice Boltzmann Schemes For Anisotropic Advection-Diffusion Equation

  • Irina Ginzburg (a1)

This paper establishes relations between the stability and the high-order truncated corrections for modeling of the mass conservation equation with the two-relaxation-times (TRT) collision operator. First we propose a simple method to derive the truncation errors from the exact, central-difference type, recurrence equations of the TRT scheme. They also supply its equivalent three-time-level discretization form. Two different relationships of the two relaxation rates nullify the third (advection) and fourth (pure diffusion) truncation errors, for any linear equilibrium and any velocity set. However, the two relaxation times alone cannot remove the leading-order advection-d if fusion error, because of the intrinsic fourth-order numerical diffusion. The truncation analysis is carefully verified for the evolution of concentration waves with the anisotropic diffusion tensors. The anisotropic equilibrium functions are presented in a simple but general form, suitable for the minimal velocity sets and the d2Q9, d3Q13, d3Q15 and d3Q19 velocity sets. All anisotropic schemes are complemented by their exact necessary von Neumann stability conditions and equivalent finite-difference stencils. The sufficient stability conditions are proposed for the most stable (OTRT) family, which enables modeling at any Peclet numbers with the same velocity amplitude. The heuristic stability analysis of the fourth-order truncated corrections extends the optimal stability to larger relationships of the two relaxation rates, in agreement with the exact (one-dimensional) and numerical (multi-dimensional) stability analysis. A special attention is put on the choice of the equilibrium weights. By combining accuracy and stability predictions, several strategies for selecting the relaxation and free-tunable equilibrium parameters are suggested and applied to the evolution of the Gaussian hill.

Corresponding author
Hide All
[1]Adhikari R. and Succi S., Duality in matrix lattice Boltzmann methods, Phys. Rev. E, 78:0667011-9, 2008.
[2]Ahrenholz B., Toölke J., and Krafczyk M., Lattice-Boltzmann simulations in reconstructed parametrized porous media. Int. J. Comput. Fluid Dyn., 20(6):369377, 2006.
[3]Ancona M. G., Fully-Lagrangian and Lattice-Boltzmann Methods for solving systems of conservation equations. J. Comput. Phys., 115:107120, 1994.
[4]Benzi R., Succi S., and Vergassola M.The lattice Boltzmann equation: theory and applications. Physics Reports, 222:145197, 1992.
[5]Chun B. and Ladd A. J. C., Interpolated boundary conditions for Lattice Boltzmann simulations of flows in narrow gaps. Phys. Rev. E, 75:0667051-11, 2007.
[6]Dong Y., Zhang J., and Yan G., A higher-order moment method of the lattice Boltzmann model for the conservation law equation, Appl. Math. Model., 34:481494, 2010.
[7]Dubois F., Equivalent partial differential equations of a lattice Boltzmann scheme. Comp. Math. Appl., 55:14411449,2008.
[8]Dubois F. and Lallemand P., Towards higher order lattice Boltzmann schemes. J. Stat. Mech.: Theory and Experiment, PO6:P06006, 2009. doi: 10.1088/1742-5468/2009/06/006006.
[9]Dubois F., Lallemand P., and Tekitek M., On a superconvergent lattice Boltzmann boundary scheme. Comp. Math. Appl., 59(7):21412149, 2010.
[10]Du Fort E. C. and Frankel S. P., Stability conditions in the numerical treatment of parabolic differential equations. M.T.A.C., 7:135, 1953.
[11]Ginzburg I. and d’Humières D., Multi-reflection boundary conditions for lattice Boltzmann models. Phys. Rev. E, 68:0666141-30,2003.
[12]Ginzburg I., Equilibrium-type and Link-type Lattice Boltzmann models for generic advection and anisotropic-dispersion equation. Adv. Wat. Res., 28:11711195, 2005.
[13]Ginzburg I., Generic boundary conditions for Lattice Boltzmann models and their application to advection and anisotropic-dispersion equations. Adv. Wat. Res., 28:11961216, 2005.
[14]Ginzburg I., Variably saturated flow described with the anisotropic Lattice Boltzmann methods. J. Comp. Fluids, 25:831848, 2006.
[15]Ginzburg I. and d’Humières D., Lattice Boltzmann and analytical modeling of flow processes in anisotropic and heterogeneous stratified aquifers. Adv. Wat. Res., 30:22022234, 2007.
[16]Ginzburg I., Lattice Boltzmann modeling with discontinuous collision components. Hydro-dynamic and advection-diffusion equations. J. Stat. Phys., 126:157203, 2007.
[17]Ginzburg I., Consistent Lattice Boltzmann schemes for the Brinkman model of porous flow and infinite Chapman-Enskog expansion. Phys. Rev. E, 77:0666704:1-12, 2008.
[18]Ginzburg I., Verhaeghe F., and d’Humières D., Two-relaxation-time Lattice Boltzmann scheme: about parametrization, velocity, pressure and mixed boundary conditions. Commun. Comput. Phys, 3:427478, 2008.
[19]Ginzburg I., Verhaeghe F., and d’Humières D., Study of simple hydrodynamic solutions with the two-relaxation-times Lattice Boltzmann scheme. Commun. Comput. Phys, 3:519581, 2008.
[20]Ginzburg I., d’Humières D., and Kuzmin A., Optimal stability of advection-diffusion Lattice Boltzmann models with two relaxation times for positive/negative equilibrium. J. Stat. Phys., 139(6):10901143, 2010. DOI 10.1007/s10955-010-9969-9.
[21]Hammou H., Ginzburg I., and Boulerhcha M., Two-relaxation-times Lattice Boltz-mann schemes for solute transport in unsaturated water flow, with a focus on stability. Adv. Wat. Res., 34(6):779793, 2011.
[22]Higuera F. J., Succi S., and Benzi R., Lattice gas dynamics with enhanced collisions. Euro-phys. Lett., 9:345349, 1989.
[23]Higuera F. J. and Jiménez J., Boltzmann approach to lattice gas simulations. Europhys. Lett., 9:663668, 1989.
[24]Hindmarsch A. C., Grescho P. M., and Griffiths D. F., The stability of explicit time-integration for certain finite difference approximation of the multi-dimensional advection-diffusion equation. Int. J. Num. Methods Fluids, 4:853897, 1984.
[25]Hirt C. W., Heuristic stability theory for finite-difference equations. J. Comput. Phys., 2:339355, 1968.
[26]Holdych D. J., Noble D. R., Georgiadis J. G., and Buckius R. O., Truncation error analysis of Lattice Boltzmann methods. J. Comput. Phys., 193(2):595619, 2004.
[27]d’Humières D., Generalized Lattice-Boltzmann equations. AIAA Rarefied Gas Dynamics: Theory and Simulations. Progress in Astronautics and Aeronautics, 59:450548, 1992.
[28]d’Humières D., Bouzidi M., and Lallemand P., Thirteen-velocity three-dimensional Lattice Boltzmann model. Phys. Rev. E, 63:0667021-7, 2001.
[29]d’Humières D., Ginzburg I., Krafczyk M., Lallemand P., and Luo L.-S., Multiple-relaxation-time Lattice Boltzmann models in three dimensions. Phil. Trans. R. Soc. Lond. A, 360:437451, 2002.
[30]d’Humières D. and Ginzburg I., Viscosity independent numerical errors for Lattice Boltz-mann models: from recurrence equations to “magic” collision numbers. Comp. Math. Appl., 58(5):823840, 2009.
[31]d’Humières D., Some stability results, in preparation, 2010.
[32]Junk M. and Yang Z., Analysis and invariant properties of a Lattice Boltzmann method. Adv. Appl. Math. Mech., 2/5:640669, 2010.
[33]Koivu V., Decain M., Geindreau C., Mattila K., Bloch J.-F., and Kataja M., Transport properties of heterogeneous materials. Combining computerised X-ray micro-tomography and direct numerical simulations. Int. J. Comput. Fluid Dyn., 23(10):713721, 2009.
[34]Kuzmin A., Ginzburg I., and Mohamad A. A., A role of the kinetic parameter on the stability of two-relaxation-times advection-diffusion Lattice Boltzmann scheme. Comp. Math. Appl., 61:34173442, 2011.
[35]Kwok Y.-K. and Tam K.-K., Stability analysis of three-level difference schemes for initial boundary problems for multi-dimensional convective-diffusion equations. Commun. Numer. Methods Engr., 9:595605, 1993.
[36]Lallemand P. and Luo L.-S., Theory of the lattice Boltzmann method: dispersion, dissipation, isotropy, Galilean invariance, and stability. Phys. Rev. E, 61:65466562, 2000.
[37]McCracken M. E. and Abraham J., Multiple-relaxation-time lattice-Boltzmann model for multiphase flow. Phys. Rev. E, 71:0367011-9, 2005.
[38]Li Y. and Huang P., A coupled lattice Boltzmann model for advection and anisotropic dispersion problem in shallow water. Adv. Wat. Res., 31(12):17191730, 2008.
[39]Miller J. J. H., On the location of zeros of certaines classes of polynomials with application to numerical analysis. J. Inst. Math. Appl., 8:397406, 1971.
[40]Pan C., Luo L.-S., and Miller C. T., An evaluation of lattice Boltzmann schemes for porous media simulation. J. Comp. Fluids, 35(4):898909, 2006.
[41]Qian Y., d’Humières D., and Lallemand P., Lattice BGK models for Navier-Stokes equation. Europhys. Lett., 17:479484, 1992.
[42]Rasin I., Saucci S. and Miller W., A multi-relaxation lattice kinetic method for passive scalar diffusion, J. Comput. Phys., 206:453462, 2005.
[43]Rheinländer M., Stability and multiscale analysis of an advective Lattice Boltzmann scheme. Progress Comp. Fluid Dynamics, 8(1-4):5668, 2008.
[44]Servan-Camas B. and Tsai F. T. C., Lattice Boltzmann method for two relaxation times for advection-diffusion equation: third order analysis and stability analysis. Adv. Wat. Res., 31:11131126, 2008.
[45]Servan-Camas B. and Tsai F. T. C., Saltwater intrusion modeling in heterogeneous confined aquifers using two-relaxation-time Lattice Boltzmann method. J. Comput. Phys., 228:236256, 2009.
[46]Servan-Camas B. and Tsai F. T. C., Two-relaxation-time lattice Boltzmann method for the anisotropic dispersive Henry problem. Water Resources Research, 46:W0251525, 2009.
[47]van der Sman R. G. M., MRT Lattice Boltzmann schemes for confined suspension flows. Comp. Phys. Commun., 181:15621569,2010.
[48]Suga S., Stability and accuracy of Lattice Boltzmann Schemes for anisotropic advection-diffusion equations. Int. J. Mod. Phys. C, 20(4):633650, 2009.
[49]Suga S., An accurate multi-level finite difference scheme for 1D diffusion equations derived from the Lattice Boltzmann method. J. Stat. Phys., 140:494503.
[50]Wolfram S., software package “Mathematica 6”.
[51]Yoshida H. and Nagaoka M., Multiple-relaxation-time lattice Boltzmann model for the convection and anisotropic diffusion equation. J. Comput. Phys., 229:77747795, 2010.
[52]Zhang X., Bengough A. G., Deeks L. K., Crawford J. W., and Young I. M.A lattice BGK model for advection and anisotropic dispersion equation. Adv. Wat. Res., 25:18, 2002.
Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Communications in Computational Physics
  • ISSN: 1815-2406
  • EISSN: 1991-7120
  • URL: /core/journals/communications-in-computational-physics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 38 *
Loading metrics...

Abstract views

Total abstract views: 189 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 16th January 2018. This data will be updated every 24 hours.