Hi, this is my first thread and nice to meet you all!
I want to use the D2Q25ZOT (zeroonethree) lattices introduced by Chikatamarla and Karlin on Physical Review E 79,046701(2009) for 2 dimensional LidDriven Cavity problem.
At my first attempt to solve this, I was using the Halfway bounce back BC just like when we use the D2Q9 lattices. Unfortunately, the results is kinda weird and unphysical even for moderate Reynolds number.
Now, I realize that the problem is on the BC, one have to use the BC such as in Physical Review E 90,043306(2014) rather than the bounceback which is an extension of the TammMottSmith (TMS) BC along with the Malaspinas’ “Regularized Boundary Condition”. From those papers, I conclude (with a little doubt) that the steps are:
 evaluate the derivatives of velocity and temperature on target nodes with finitedifference.
 evaluate the nonequilibrium and equilibrium distribution function at target nodes.
 compute fi’ (equation 93 on PRE 90 043306)
 replace the distribution function fi such as fi = fi + fi’(target)  fi’(local)
Hence, I will summarize my problem in two questions,

Did I missed something here? because this is a new thing for me as I am used to implement the common boundary conditions only(bounceback, halfway bounceback, periodic) and I am pretty confused when choosing which are target nodes, and which are local nodes particularly on D2Q5ZOT lattices facing the north, south, west, and east walls.

How to implement those steps exactly in my code? the reason is similar to question number one.
Feel free if you want to ask me about details and want to read those papers. You can write down your email here or contact me personally at pradipto.academics@gmail.com and I will give you everything you need.
Any kind of help would be appreciated.
Regards,
Pradipto