[hiwahiwa]

I have written a Python code recently on compressible LBM using both the DDF model and Hybrid but they both do not give me the correct results for the shock tube. The results of DDF are unstable and give me nan while the hybrid method gives no change in the density contour. I am confused where is the problem! Whether the nondimensionalization is not ok or the Lattice and physical units are not compatible. I thought the problem is with the shock waves but the point is that I can't the correct results for the Couette flow too. I will appreciate if anyone can help me in this regard. Details are available upon request