Abstract

We propose a two-population lattice Boltzmann model on standard lattices for the simulation of compressible flows. The model is fully on-lattice and uses the single relaxation time Bhatnagar–Gross–Krook kinetic equations along with appropriate correction terms to recover the Navier–Stokes–Fourier equations. The accuracy and performance of the model are analyzed through simulations of compressible benchmark cases including Sod shock tube, sound generation in shock–vortex interaction, and compressible decaying turbulence in a box with eddy shocklets. It is demonstrated that the present model provides an accurate representation of compressible flows, even in the presence of turbulence and shock waves.

Highlights

  • The development of accurate and efficient numerical methods for the simulation of compressible fluid flows remains a highly active research field in computational fluid dynamics (CFD), and is of great importance to many natural phenomena and engineering applications

  • We propose a two-population lattice Boltzmann model on standard lattices for the simulation of compressible flows

  • It is impossible to use a grid size fine enough to resolve the physical shock structure defined by the molecular viscosity, most numerical schemes rely on some numerical dissipation to stabilize the simulation and capture the shock over a few grid points.[1,2]

Read more

Summary

INTRODUCTION

The development of accurate and efficient numerical methods for the simulation of compressible fluid flows remains a highly active research field in computational fluid dynamics (CFD), and is of great importance to many natural phenomena and engineering applications. Scitation.org/journal/phf extent, they increase significantly the computational cost and suffer from a limited temperature range,[16] as well Another approach, which has received considerable attention in recent years, maintains the simplicity and efficiency of the standard lattices and employs correction terms in order to remove the aforementioned spurious terms in the stress tensor.[17,18] Due to intrinsic nonuniqueness of the correction term, different implementations exist in the literature, all recover the same equations in the hydrodynamic limit.[19,20,21] See Hosseini et al.[22] for a detailed review of different implementations.

Kinetic equations
Discrete velocities and factorization
Extended f-equilibrium
Hydrodynamic limit
À1 x1 2
Equations of gas dynamics
Summary of the lattice Boltzmann model
General implementation issues
Sod’s shock tube
Shock–vortex interaction
Decaying of compressible isotropic turbulence
Low turbulent Mach number
EðjÞdj pffiffiffiffiffi 2p : j0
Effect of deviation discretization on the accuracy
High turbulent Mach number
CONCLUSION
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call