L'opérateur de Bhatnagar-Gross-Krook (en abrégé BGK) est un opérateur linéaire qui se substitue à l'opérateur de collision de l'équation de Boltzmann
la fonction de distribution statistique de la vitesse à l'instant au point , sa valeur d'équilibre donnée par la statistique de Maxwell et un temps caractéristique. Cette approximation a été introduite en 1954 par Prabhu Lal Bhatnagar, Eugene Gross et Max Krook[1]. Elle permet de substantielles simplifications de la résolution de l'équation de Boltzmann et est très utilisée dans la méthode de Boltzmann sur réseau.
Il s'agit d'un terme de relaxation vers l'équilibre beaucoup plus simple que l'opérateur exact mais qui respecte les propriétés fondamentales de celui-ci pour une interaction moléculaire :
Par contre il a l'inconvénient de conduire à un nombre de Prandtl égal à l'unité ainsi qu'on peut le voir en effectuant un développement de type Chapman-Enskog. Des modifications du modèle permettent de pallier cet inconvénient. Parmi celles-ci l'opérateur ES-BGK (Ellipsoidal Statistical BGK) proposé par Lowell H. Holway Jr.[2] où l'on remplace la solution d'équilibre par une distribution maxwellienne anisotrope permet d'obtenir un nombre de Prandtl égal à 2/3. Cette méthode a été étendue pour une plus grande généralité et permet d'obtenir des résultats proches de la solution exacte de l'équation de Boltzmann[3].