scholarly journals Theory, Analysis, and Applications of the Entropic Lattice Boltzmann Model for Compressible Flows

Entropy ◽  
2020 ◽  
Vol 22 (3) ◽  
pp. 370 ◽  
Author(s):  
Nicolò Frapolli ◽  
Shyam Chikatamarla ◽  
Ilya Karlin

The entropic lattice Boltzmann method for the simulation of compressible flows is studied in detail and new opportunities for extending operating range are explored. We address limitations on the maximum Mach number and temperature range allowed for a given lattice. Solutions to both these problems are presented by modifying the original lattices without increasing the number of discrete velocities and without altering the numerical algorithm. In order to increase the Mach number, we employ shifted lattices while the magnitude of lattice speeds is increased in order to extend the temperature range. Accuracy and efficiency of the shifted lattices are demonstrated with simulations of the supersonic flow field around a diamond-shaped and NACA0012 airfoil, the subsonic, transonic, and supersonic flow field around the Busemann biplane, and the interaction of vortices with a planar shock wave. For the lattices with extended temperature range, the model is validated with the simulation of the Richtmyer–Meshkov instability. We also discuss some key ideas of how to reduce the number of discrete speeds in three-dimensional simulations by pruning of the higher-order lattices, and introduce a new construction of the corresponding guided equilibrium by entropy minimization.

2015 ◽  
Vol 26 (01) ◽  
pp. 1550006 ◽  
Author(s):  
Mohsen Ghadyani ◽  
Vahid Esfahanian ◽  
Mohammad Taeibi-Rahni

Attempts to simulate compressible flows with moderate Mach number to relatively high ones using Lattice Boltzmann Method (LBM) have been made by numerous researchers in the recent decade. The stability of the LBM is a challenging problem in the simulation of compressible flows with different types of embedded discontinuities. The present study proposes an approach for simulation of inviscid flows by a compressible LB model in order to enhance the robustness using a combination of Essentially NonOscillatory (ENO) scheme and Shock-Detecting Sensor (SDS) procedure. A sensor is introduced with adjustable parameters which is active near the discontinuities and affects less on smooth regions. The validity of the improved model to capture shocks and to resolve contact discontinuity and rarefaction waves in the well-known benchmarks such as, Riemann problem, and shock reflection is investigated. In addition, the problem of supersonic flow in a channel with ramp is simulated using a skewed rectangular grid generated by an algebraic grid generation method. The numerical results are compared with analytical ones and those obtained by solving the original model. The numerical results show that the presented scheme is capable of generating more robust solutions in the simulation of compressible flows and is almost free of oscillations for high Mach numbers. Good agreements are obtained for all problems.


Author(s):  
Venkata Ravishankar Kasibhotla ◽  
Danesh Tafti

The paper is concerned with the prediction and analysis of dynamic stall of flow past a pitching NACA0012 airfoil at 1 million Reynolds number based on the chord length of the airfoil and at reduced frequency of 0.25 in a three dimensional flow field. The turbulence in the flow field is resolved using large eddy simulations with the dynamic Smagorinsky model at the sub grid scale. The development of dynamic stall vortex, shedding and reattachment as predicted by the present study are discussed in detail. This study has shown that the downstroke phase of the pitching motion is strongly three dimensional and is highly complex, whereas the flow is practically two dimensional during the upstroke. The lift coefficient agrees well with the measurements during the upstroke. However, there are differences during the downstroke. The computed lift coefficient undergoes a sharp drop during the start of the downstroke as the convected leading edge vortex moves away from the airfoil surface. This is followed by a recovery of the lift coefficient with the formation of a secondary trailing edge vortex. While these dynamics are clearly reflected in the predicted lift coefficient, the experimental evolution of lift during the downstroke maintains a fairly smooth and monotonic decrease in the lift coefficient with no lift recovery. The simulations also show that the reattachment process of the stalled airfoil is completed before the start of the upstroke in the subsequent cycle due to the high reduced frequency of the pitching cycle.


Sign in / Sign up

Export Citation Format

Share Document