This paper presents and compares two different strategies in the numerical simulation of passive microfluidic mixers based on chaotic advection. In addition to flow velocity field calculations, concentration distributions of molecules and trajectories of microscale particles were determined and compared to evaluate the performance of the applied modeling approaches in the proposed geometries. A staggered herringbone type micromixer (SHM) was selected and studied in order to demonstrate finite element modeling issues. The selected microstructures were fabricated by a soft lithography technique, utilizing multilayer SU-8 epoxy-based photoresist as a molding replica for polydimethylsiloxane (PDMS) casting. The mixing processes in the microfluidic systems were characterized by applying molecular and particle (cell) solutions and adequate microscopic visualization techniques. We proved that modeling of the molecular concentration field is more costly, in regards to computational time, than the particle trajectory based method. However, both approaches showed adequate qualitative agreement with the experimental results.