Abstract
In this paper, we described numerically several scenarios of formation of vortex flows (VF) in microsized hybrid-oriented liquid crystal (HOLC) channels with orientation defects using a nonlinear generalization of the classical Ericksen–Leslie theory that allows taking into account termomechanical contribution, both in the expression for the shear stress and in the entropy balance equation. An analysis of the numerical results showed that there are two or one vortices in the HOLC channel although two vortices directed towards each other are generated at the initial stage of the VT formation Thermomechanically Excited Vortical Flow.
Similar content being viewed by others
REFERENCES
H. A. Stone, A. D. Stroock, and A. Ajdari, Ann. Rev. Fluid Mech. 36, 381 (2004).
T. M. Squires and S. R. Quake, Rev. Mod. Phys. 77, 977 (2005).
R. B. Shoch, J. Y. Han, and P. Renaud, Rev. Mod. Phys. 80, 839 (2008).
P. Hanggi and F. Marchesoni, Rev. Mod. Phys. 81, 387 (2009).
W. Sparreboom, A. van den Berg, and J. C. T. Eijkel, New J. Phys. 12, 0115004 (2010).
P. Popov, L. W. Hanaker, M. Mirheydari, E. K. Mann, and A. Jakli, Sci. Rep. 7, 1603 (2017).
J. Beekman, I. Nys, O. Willekens, and K. Neyts, J. Appl. Phys. 121, 023106 (2017).
S. Lee, R. An, and J. A. Hunt, Nat. Nanotech. 5, 412 (2010).
E. Verneuil, M. L. Cordero, F. Gallaire, and Ch. N. Ba-roud, Langmuir 25, 5127 (2009).
P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Oxford Univ. Press, Oxford, 1995).
A. V. Zakharov and A. A. Vakulenko, J. Chem. Phys. 127, 084907 (2007).
R. S. Akopyan, R. B. Alaverdian, E. A. Santrosian, and Y. S. Chilingarian, J. Appl. Phys. 90, 3371 (2001).
J. L. Ericksen, Arch. Ration. Mech. Anal. 4, 231 (1960).
F. M. Leslie, Arch. Ration. Mech. Anal. 28, 265 (1968).
S. R. de Groot and P. Mazur, Non-Equilibrum Thermodynamics (North-Holland, Amsterdam, 1962; Mir, Moscow, 1964).
A. V. Zakharov and P. V. Maslennikov, Phys. Rev. E 96, 052705 (2017).
I. S. Berezin and N. R. Zhidkov, Computing Methods (Fizmatgiz, Moscow, 1964; Elsevier, Amsterdam, 1965).
A. A. Samarskii and E. S. Nikolaev, Numerical Methods for Grid Equations (Nauka, Moscow, 1978; Birkhäuser, Basel, 1989).
FUNDING
This work was supported by the Ministry of Education and Science of Russian Federation, projects 3.11888.2018/11.12 and 3.9585.2017/8.9.
Author information
Authors and Affiliations
Corresponding author
Additional information
Translated by A. Ivanov
APPENDIX
APPENDIX
The torque balance equation in the LC phase has the following form:
where Tel = \(\frac{{\delta {{\mathcal{W}}_{{{\text{el}}}}}}}{{\delta {\mathbf{\hat {n}}}}} \times {\mathbf{\hat {n}}}\) is elastic, Tvis = \(\frac{{\delta {{\mathcal{R}}^{{{\text{vis}}}}}}}{{\delta {{{{\mathbf{\hat {n}}}}}_{t}}}} \times {\mathbf{\hat {n}}}\) is viscous, and Ttm = \(\frac{{\delta {{\mathcal{R}}^{{{\text{tm}}}}}}}{{\delta {{{{\mathbf{\hat {n}}}}}_{t}}}} \times {\mathbf{\hat {n}}}\) is the thermomechanical contribution to the torque balance acting on the unit volume of the LC phase, respectively. Here, \({{{\mathbf{\hat {n}}}}_{t}}\) ≡ \(\frac{{d{\mathbf{\hat {n}}}}}{{dt}}\) is the material derivative of the vector \({\mathbf{\hat {n}}}\) = \({{n}_{x}}{\mathbf{\hat {i}}}\) + \({{n}_{z}}{\mathbf{\hat {k}}}\).
The linear balance momentum equation in the LC phase can be written in the form
where \(\frac{{d{v}}}{{dt}} = \frac{{\partial {v}}}{{\partial t}} + u{{{v}}_{{,x}}} + w{{{v}}_{{,z}}}\), σ = σel + σvis + σtm – \(\mathcal{P}\mathcal{E}\) is the full expression for the stress tensor (ST) consisting of the elastic σel = \( - \frac{{\partial {{\mathcal{W}}_{{{\text{el}}}}}}}{{\partial \nabla {\mathbf{\hat {n}}}}} \cdot {{(\nabla {\mathbf{\hat {n}}})}^{T}}\), viscous σvis = \(\frac{{\delta {{\mathcal{R}}^{{{\text{vis}}}}}}}{{\delta \nabla {\mathbf{v}}}}\) and thermomechanical σtm = \(\frac{{\delta {{R}^{{{\text{tm}}}}}}}{{\delta \nabla {\mathbf{v}}}}\) contributions to ST, respectively. Here, \(\mathcal{R}\) = \({{\mathcal{R}}^{{{\text{vis}}}}}\) + \({{\mathcal{R}}^{{{\text{tm}}}}}\) + \({{\mathcal{R}}^{{{\text{th}}}}}\) is the full Rayleigh dissipation function, \({{\mathcal{W}}_{{{\text{el}}}}}\) = \(\frac{1}{2}\)[K1 × (∇ · \({\mathbf{\hat {n}}}\))2 + K3(\({\mathbf{\hat {n}}}\) × ∇ × \({\mathbf{\hat {n}}}\))2] is the elastic energy density, K1 and K3 are the splay and bend elastic constants, \(\mathcal{P}\) is the hydrostatic pressure in the LC system, and \(\mathcal{E}\) is the unit tensor while
is viscous, but
and
are thermomechanical and thermal contributions to the full expression for the Rayleigh function \(\mathcal{R}\), respectively.
Here, α1–α6 are the Leslie viscosity coefficients, γ1(T) and γ2(T) are the rotational viscosity coefficients of the LC system, ξ is the thermomechanical constant, and λ|| and λ⊥ are the thermal conductivity coefficients of the LC system along and across the direction of the director \({\mathbf{\hat {n}}}\), respectively. The tensors Ds = \(\frac{1}{2}\)[∇v + (∇v)T] and Da = \(\frac{1}{2}\)[∇v – (∇v)T] are symmetric and asymmetric contributions to the strain rate tensor, M = \(\frac{1}{2}\)[∇\({\mathbf{\hat {n}}}\) + (∇\({\mathbf{\hat {n}}}\))T], and \({{\mathcal{M}}_{0}}\) = ∇ · \({\mathbf{\hat {n}}}\) is the scalar invariant of the tensor M.
The heat conduction equation describing the change in the temperature field T(x, z, t) under the effect of the heat flux q through the lower boundary of the HOLC channel has the following form:
where Q = \( - T\frac{{\delta \mathcal{R}}}{{\delta \nabla T}}\) is the heat flux in the LC system, and CP is the heat capacity coefficient, respectively.
Rights and permissions
About this article
Cite this article
Zakharov, A.V. Thermomechanically Excited Vortical Flow in a Hybrid-Oriented Nematic Channel. Phys. Solid State 61, 1136–1143 (2019). https://doi.org/10.1134/S1063783419060295
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1134/S1063783419060295