Variational formulation of active nematics: theory and simulation

TL;DR

This paper introduces a variational framework for active nematics based on Onsager's formalism, enabling systematic derivation of models and numerical schemes for complex biological systems like cytoskeletal gels and tissues.

Contribution

It develops a thermodynamically consistent variational approach for active nematics, including new models and structured-preserving numerical methods.

Findings

Successfully framed standard models within the variational formalism

Developed a new compressible model for density-dependent active nematics

Demonstrated the approach with biological examples like wound healing and cell colonies

Abstract

The structure and dynamics of important biological quasi-two-dimensional systems, ranging from cytoskeletal gels to tissues, are controlled by nematic order, flow, defects and activity. Continuum hydrodynamic descriptions combined with numerical simulations have been used to understand such complex systems. The development of thermodynamically consistent theories and numerical methods to model active nemato-hydrodynamics is eased by mathematical formalisms enabling systematic derivations and structured-preserving algorithms. Alternative to classical nonequilibrium thermodynamics and bracket formalisms, here we develop a theoretical and computational framework for active nematics based on Onsager's variational formalism to irreversible thermodynamics, according to which the dynamics result from the minimization of a Rayleighian functional capturing the competition between free-energy release, dissipation and activity. We show that two standard incompressible models of active nemato-hydrodynamics can be framed in the variational formalism, and develop a new compressible model for density-dependent active nemato-hydrodynamics relevant to model actomyosin gels. We show that the variational principle enables a direct and transparent derivation not only of the governing equations, but also of the finite element numerical scheme. We exercise this model in two representative examples of active nemato-hydrodynamics relevant to the actin cytoskeleton during wound healing and to the dynamics of confined colonies of elongated cells.