Dear OpenLb developer,
First of all I should thank you for sharing this great opensource code. I have some questions regarding the LES SGS model that you developed in the source code. In the aorta_example, there is a function that compute y_plus at the boundaries, it seems that the value is less than one all the time, considering that the (a) wall boundary in LBM is like staircase and (b) LES wall BL resolving is very tricky (even in FVM which uses unstructured mesh) how is it possible to reach that value? I was wondering if it is possible to resolve boundary layer in LBM. I will be appreciated if you would please let me know whether there is any recent development in capturing BL in LBM, I am specifically interested in turbulet wallshearstress and the way you resolve BL in the boundaries.
Many thanks again
Ryan