In this work the fractional Hunter-Saxton equation applied in the study of diffusion of nematic liquid crystals was done involving partial operators with two fractional orders, \(\alpha\) and \(\beta\), via Atangana-Riemann and Atangana-Caputo with bi-order and via Riemann-Liouville, Caputo-Fabrizio-Riemann and Atangana-Baleanu-Riemann for the space domain. The mathematical equation underpinning this physical phenomenon was solved numerically using an iterative scheme where the numerical approximations for second order were developed. The new approach with two fractional orders is able to consider media with two different layers, scales and properties. The generalization of this equation exhibit different cases of anomalous behavior and the numerical solutions obtained describes the propagation of waves in a nematic liquid cristal.