  The LBM is combined with the Volume Penalization (VP LBM) approach to
simulate flows in the presence of obstacles. The single relaxation time LBM and
the multiple relaxation time LBM are used. For cases where the fluid motion is
enhanced by moving boundaries, the comparison of our results with analytical
results is satisfactory. The flow around a cylinder at a relatively high
Reynolds number is also calculated. For that case,the VP LBM method was
parallelized. A comparison of our results with those found in the literature,
demonstrates the reliability of the VP LBM to treat complex flows.
