Efficient Numerical Schemes for a Two-Phase Hydrodynamical Model of Active Liquid Crystals and Solids
This work provides efficient numerical methods for simulating active nematic fluids interacting with solids, which is important for understanding soft confined active matter systems.
The paper proposes linear, fully decoupled numerical schemes for a two-phase hydrodynamic model of active liquid crystals and solids, proving thermodynamic consistency and energy stability. Numerical simulations investigate how activity and parameters affect active nematic fluid-solid interactions.
We propose several linear, fully decoupled numerical schemes with first- and second-order temporal accuracy for a novel Q-tensor-based two-phase hydrodynamic model describing the coupling of active nematic liquid crystal solutions with isotropic solid substrates. The model is derived from the generalized Onsager principle and includes nontrivial terms that contribute zero to the total free-energy dissipation. We prove that the proposed decoupled linear schemes are thermodynamically consistent at the discrete level. In the passive limit, the SGE-BDF1 and SGE-PDG schemes are unconditionally energy stable, while the SGE-BDF2 scheme is energy stable with respect to a modified energy under a standard boundedness assumption and a sufficiently large stabilization parameter. We perform extensive numerical simulations to investigate how activity and other model parameters affect active nematic fluid-solid interactions. Finally, we analyze the physical mechanisms underlying the observed behaviors, providing deeper insight into the dynamics of soft confined active nematic fluids.