Projection-Free Method for the Full Frank-Oseen Model of Liquid Crystals (2405.03145v1)
Abstract: Liquid crystals are materials that experience an intermediate phase where the material can flow like a liquid, but the molecules maintain an orientation order. The Frank-Oseen model is a continuum model of a liquid crystal. The model represents the liquid crystal orientation as a vector field and posits that the vector field minimizes some elastic energy subject to a pointwise unit length constraint, which is a nonconvex constraint. Previous numerical methods in the literature assumed restrictions on the physical constants or had regularity assumptions that ruled out point defects, which are important physical phenomena to model. We present a finite element discretization of the full Frank-Oseen model and a projection free gradient flow algorithm for the discrete problem in the spirit of Bartels (2016). We prove Gamma-convergence of the discrete to the continuous problem: weak convergence of subsequences of discrete minimizers and convergence of energies. We also prove that the gradient flow algorithm has a desirable energy decrease property. Our analysis only requires that the physical constants are positive, which presents challenges due to the additional nonlinearities from the elastic energy.
- An energy-minimization finite-element approach for the frank–oseen model of nematic liquid crystals. SIAM Journal on Numerical Analysis, 53(5):2226–2254, 2015.
- Paraview: An end-user tool for large data visualization. The visualization handbook, 717(8), 2005.
- Determination of the liquid-crystal surface elastic constant k 24. Physical review letters, 67(11):1442, 1991.
- François Alouges. A new algorithm for computing liquid crystal stable configurations: the harmonic mapping case. SIAM journal on numerical analysis, 34(5):1708–1726, 1997.
- Minimizing oseen-frank energy for nematic liquid crystals: algorithms and numerical results. In Annales de l’IHP Physique théorique, volume 66, pages 411–447, 1997.
- Approximation of fractional harmonic maps. arXiv preprint arXiv:2104.10049, 2021.
- Sören Bartels. Stability and convergence of finite-element approximation schemes for harmonic maps. SIAM journal on numerical analysis, 43(1):220–238, 2005.
- Sören Bartels. Numerical methods for nonlinear partial differential equations, volume 47. Springer, 2015.
- Sören Bartels. Projection-free approximation of geometrically constrained partial differential equations. Mathematics of Computation, 85(299):1033–1049, 2016.
- Benchmarking numerical algorithms for harmonic maps into the sphere. arXiv preprint arXiv:2209.13665, 2022.
- Stable gradient flow discretizations for simulating bilayer plate bending with isometry and obstacle constraints. IMA Journal of Numerical Analysis, 42(3):1903–1928, 2022.
- Quasi-optimal error estimates for the approximation of stable harmonic maps. arXiv preprint arXiv:2209.11985, 2022.
- ΓΓ\Gammaroman_Γ-convergent LDG method for large bending deformations of bilayer plates. arXiv preprint arXiv:2301.03151, 2023.
- A structure-preserving fem for the uniaxially constrained q-tensor model of nematic liquid crystals. Numerische Mathematik, 145(4):837–881, 2020.
- Yunmei Chen. The weak solutions to the evolution problems of harmonic maps. Mathematische Zeitschrift, 201(1):69–74, 1989.
- Minimum energy configurations for liquid crystals: Computational results. In Theory and Applications of Liquid Crystals, pages 99–121. Springer, 1987.
- Relaxation and gradient methods for molecular orientation in liquid crystals. Computer Physics Communications, 53(1-3):455–465, 1989.
- Finite element analysis of the landau–de gennes minimization problem for liquid crystals. SIAM Journal on Numerical Analysis, 35(1):336–362, 1998.
- The physics of liquid crystals. Number 83. Oxford university press, 1993.
- JL Ericksen. Nilpotent energies in liquid crystal theory. Archive for Rational Mechanics and Analysis, 10(1):189–196, 1962.
- JL Ericksen. Inequalities in liquid crystal theory. The physics of Fluids, 9(6):1205–1207, 1966.
- JL Ericksen. Liquid crystals with variable degree of orientation. Archive for Rational Mechanics and Analysis, 113(2):97–120, 1991.
- Frederick C Frank. I. liquid crystals. on the theory of liquid crystals. Discussions of the Faraday Society, 25:19–28, 1958.
- Vsevolod Fréedericksz and V Zolina. Forces causing the orientation of an anisotropic liquid. Transactions of the Faraday Society, 29(140):919–930, 1933.
- An operator-splitting method for a liquid crystal model. Computer physics communications, 152(3):242–252, 2003.
- Mathematical questions of liquid crystal theory. In Theory and applications of liquid crystals, pages 151–184. Springer, 1987.
- Existence and partial regularity of static liquid crystal configurations. Communications in mathematical physics, 105(4):547–570, 1986.
- A saddle point approach to the computation of harmonic maps. SIAM Journal on Numerical Analysis, 47(2):1500–1523, 2009.
- A newton-penalty method for a simplified liquid crystal model. Advances in Computational Mathematics, 40(1):201–244, 2014.
- Tunable dynamic topological defect pattern formation in nematic liquid crystals. Advanced Optical Materials, 8(1):1900991, 2020.
- Second variation of liquid crystal energy at x/|x|𝑥𝑥x/|x|italic_x / | italic_x |. Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 437(1900):475–487, 1992.
- Gamma-convergent projection-free finite element methods for nematic liquid crystals: The ericksen model. SIAM Journal on Numerical Analysis, 60(2):856–887, 2022.
- A finite element method for nematic liquid crystals with variable degree of orientation. SIAM Journal on Numerical Analysis, 55(3):1357–1386, 2017.
- Carl W. Oseen. The theory of liquid crystals. Transactions of the Faraday Society, 29(140):883–899, 1933.
- Solution of sparse indefinite systems of linear equations. SIAM journal on numerical analysis, 12(4):617–629, 1975.
- Martin Schadt. Liquid crystal materials and liquid crystal displays. Annual review of materials science, 27(1):305–379, 1997.
- Numerical method for the equilibrium configurations of a maier-saupe bulk potential in a q-tensor model of an anisotropic nematic liquid crystal. Journal of Computational Physics, 441:110441, 2021.
- Joachim Schöberl. C++ 11 implementation of finite elements in ngsolve. Institute for Analysis and Scientific Computing, Vienna University of Technology, 30, 2014.
- Epifanio G Virga. Variational theories for liquid crystals, volume 8. CRC Press, 1995.
- Shawn W Walker. A finite element method for the generalized ericksen model of nematic liquid crystals. ESAIM: Mathematical Modelling and Numerical Analysis, 54(4):1181–1220, 2020.
- Modelling and computation of liquid crystals. Acta Numerica, 30:765–851, 2021.
- Augmented lagrangian preconditioners for the oseen–frank model of nematic and cholesteric liquid crystals. BIT Numerical Mathematics, 61(2):607–644, 2021.