“…The method enables simulations of the low-Reynolds-number hydrodynamics of simple fluids, of passive liquid crystal flows [65,39,1,2,14], and finally, of active polar viscous gels. Here, we present numerical results of the time-evolution of the polarization and velocity fields in a passive liquid crystal [1,2], in an active polar film [33,32], and in active polar viscous gels with topological defects in the polarization field [33,36].…”