Publication
Active nematic fluids confined in narrow channels are known to generate spontaneous flows when the activity is sufficiently intense. Recently, it was demonstrated [R. Green, J. Toner, and V. Vitelli, Phys. Rev. Fluids 2, 104201 (2017)] that if the molecular anchoring at the channel walls is conflicting, i.e., perpendicular on one plate and parallel on the other, flows are initiated even in the zero activity limit. An analytical laminar velocity profile for this specific configuration was derived within a simplified nematohydrodynamic model in which the nematic order parameter is a fixed-magnitude unit vector n. The solution holds in a regime where the flow does not perturb the nematic order imposed by the walls. In this study, we explore systematically active flows in this confined geometry with a more general theoretical model that uses a second-rank tensor order parameter Q to express both the magnitude and orientation of the nematic phase. The Q-model allows for the presence of defects and biaxial, in addition to uniaxial, molecular arrangements. Our aim is to provide a unified picture, beyond the limiting regime explored previously, to serve as a guide for potential microfluidic applications that exploit the coupling between the orientational order of the molecules and the velocity field to finely control the flow and overcome the intrinsic difficulties of directing and pumping fluids at the microscale. We reveal how the nematic-flow coupling is not only dependent on geometrical constraints, but is also highly sensitive to material and flow parameters. We specifically stress the key role played by the activity and the flow aligning parameter, and we show that solutions mostly depend on two dimensionless parameters. We find that for large values of the activity parameter, the flow is suppressed for contractile particles while it is either sustained or suppressed for extensile particles depending on whether they tend to align or tumble when subject to shear. We explain these distinct behaviors by an argument based on the results of the stability analysis applied to two simpler configurations: active flows confined between parallel plates with either orthogonal or perpendicular alignment at both walls. We show that the analytical laminar solution derived for the n model in the low activity limit is found also in the Q model, both analytically and numerically. This result is valid for both contractile and extensile particles and for a flow-tumbling as well as aligning nematics. We remark that this velocity profile can be derived for generic boundary conditions. To stress the more general nature of the Q model, we conclude by providing a numerical example of a biaxial three-dimensional thresholdless active flow for which we show that biaxiality is especially relevant for a weakly first-order isotropic-nematic phase transition.