  In the first part of this study, we compared the performances of two
categories of no-slip boundary treatments, i.e., the interpolated bounce-back
schemes and the immersed boundary methods in a series of laminar flow
simulations within the lattice Boltzmann method. In this second part, these
boundary treatments are further compared in the simulations of turbulent flows
with complex geometry to provide a next-level assessment of these schemes. Two
non-trivial turbulent flow problems, a fully developed turbulent pipe flow at a
low Reynolds number, and a decaying homogeneous isotropic turbulent flow laden
with a large number of resolved spherical particles are considered. The major
problem of the immersed boundary method revealed by the present study is its
incapability in computing the local velocity gradients inside the diffused
interface, which can result in significantly underestimated dissipation rate
and viscous diffusion locally near the particle surfaces. Otherwise, both
categories of the no-slip boundary treatments are able to provide accurate
results for most of turbulent statistics in both the carrier and dispersed
phases, provided that sufficient grid resolutions are used.
