Reply To: channel3d in olb1.6 and olb1.5

[Fuxin He]

Dear Fedor

Are you still setting physRefL = 1.0, so charPhysU = 0.0415? I thank you for providing me with the proper resolution for charPhysU = 1. Also I set these as you said:
const int rhoMethod = 2; // method for density reconstruction
// 0: Zou-He
// 1: extrapolation
// 2: constant
const int fneqMethod = 3; // method for fneq reconstruction
// 0: regularized NEBB (Latt)
// 1: extrapolation NEQ (Guo Zhaoli)
// 2: regularized second order finite Differnce
// 3: equilibrium scheme
const int wallProfile = 1; // wallfunction profile
// 0: Musker profile
// 1: power law profil

The processed pressure still shows the presence of negative values, and based on the relationship between pressure and density, I’m guessing it’s because of the calculated negative density, which is not the same as if I followed the parameters:
const T latticeWallDistance = 0.5; // lattice distance to boundary
const int rhoMethod = 0; // method for density reconstruction
// 0: Zou-He
// 1: extrapolation
// 2: constant
const int fneqMethod = 1; // method for fneq reconstruction
// 0: regularized NEBB (Latt)
// 1: extrapolation NEQ (Guo Zhaoli)
// 2: regularized second order finite Differnce
// 3: equilibrium scheme
const int wallProfile = 0; // wallfunction profile
// 0: Musker profile
// 1: power law profile

The two wall model parameters I named 231 and 010 and the processed pressures are as follows, I think the pressure field is unphysical.
https://postimg.cc/gallery/c3PSgbY

sincerely,
Fuxin He