Numerical flow solutions on a backward-facing step using the lattice Boltzmann equation method

Numerical solutions of 2-D laminar flow over a backward-facing step using the lattice Boltzmann equation method (LBEM) are presented in this article. Unlike conventional numerical schemes based on macroscopic continuum equation (mass conservation and Navier-Stokes) discretisation, the LBEM is based on microscopic models and mesoscopic kinetic equations. The simulations were validated for a wide range of Reynolds numbers (100 £ Re £ 1,000), comparing them to previous studies. Several flow features, such as primary and secondary vortex location at the bottom and top of the wall, respectively, were investigated regarding Reynolds number. Two typical classes of boundary condition were implemented in the LBEM model: the Drichlet condition at the inlet flow (parabolic speed profile) and the Newman condition at the outlet flow (zero gradient speed). The results showed that the LBEM gave accurate results over a wide range of Reynolds number; these were compared with other numerical methods and experimental data. Se presenta una solución numérica del flujo sobre un escalón en dos dimensiones utilizando el método de la ecuación reticular de Boltzmann (LBEM). A diferencia de los métodos numéricos tradicionales basados en la discretización de las ecuaciones macroscópicas del continuo (conservación de la masa y Navier-Stokes), los LBEM se fundamentan en modelos microscópicos y mesoscópicos de las ecuaciones cinéticas. Se muestran los resultados obtenidos para este flujo en el estado estacionario y para un amplio rango de números de Reynolds (100 £ Re £ 1000), y se han comparado con estudios previos. Se ha investigado la aparición y localización de los principales vórtices en el flujo, tanto en la pared inferior como en la superior, y su comportamiento en función del número de Re. Se han implementado al modelo LBEM dos tipos comunes de condiciones de frontera: condición de Drichlet a la entrada (perfil de velocidad parabólico) y condición de Newman a la salida (derivada nula de la velocidad). Los resultados obtenidos muestran gran exactitud del método utilizado para un amplio rango de números de Reynolds, al ser comparados con resultados experimentales y numéricos de otros autores.