×

Convergence of a fully discrete finite element method for a degenerate parabolic system modelling nematic liquid crystals with variable degree of orientation. (English) Zbl 1097.35082

The main topic of paper is studying a degenerate parabolic system which models the evolution of nematic liquid crystal with variable degree of orientation of the type: \[ \partial_t s=k_1 \text{div}([1+\delta| \nabla s| ^{p-2}]\nabla s)- k_2 | \nabla \varphi| ^2 s -W^\prime (s) \]
\[ s^2\partial_t \varphi =k_3\, \text{div}(s^2 \nabla \varphi) \] with appropriate initial condition and a Dirichlet boundary condition on some bounded maximum three dimensional space domain \(\Omega\) and some time interval \((0,T)\). Here the unknown scalar function \(s\) represents averaged local orientation of molecules compared to unit vector \(\mathbf{n} (x,t)\) called director and only planar director configurations \(\mathbf{n} = (\cos \varphi, \sin \varphi)\), where \(\varphi\) is the associated angle is considered. \(k_i, i=1,2,3>0\), \(\delta\geq0\) are given real numbers and \(W\) is a smooth double well potential. It is the modification of previous models investigated by Calderer, Golovaty, Lin and Liu with \( \delta=0\) and is a special case of Erickesen’s general continuum model. A regularized system with regularization parameter in second equation is proposed and well posedness of this system is proved. Moreover the existence of a weak solution and a certain dissipative energy law is proved. Then passing to the limit for regularization parameter, the same is proved also for the weak solution of the original problem. Practical fully discrete finite element method is then derived to obtain numerical solutions of the proposed problem, and a convergence of an approximation solution to the continuous one is presented in some functional spaces. Several numerical experiments for two space dimension domain are included. On their results it is shown the formation, annihilation and evolution of line singularities in such models.

MSC:

35K65 Degenerate parabolic equations
35K55 Nonlinear parabolic equations
35Q35 PDEs in connection with fluid mechanics
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs
65M60 Finite element, Rayleigh-Ritz and Galerkin methods for initial value and initial-boundary value problems involving PDEs
76A15 Liquid crystals
35K50 Systems of parabolic equations, boundary value problems (MSC2000)
PDFBibTeX XMLCite
Full Text: DOI Numdam EuDML

References:

[1] R.A. Adams , Sobolev Spaces . Academic Press, New York ( 1975 ). MR 450957 | Zbl 0314.46030 · Zbl 0314.46030
[2] M.C. Calderer , D. Golovaty , F.-H. Lin and C. Liu , Time evolution of nematic liquid crystals with variable degree of orientation . SIAM J. Math. Anal. 33 ( 2002 ) 1033 - 1047 . Zbl 1003.35108 · Zbl 1003.35108 · doi:10.1137/S0036141099362086
[3] C.M. Elliott and S. Larsson , A finite element model for the time-dependent joule heating problem . Math. Comp. 64 ( 1995 ) 1433 - 1453 . Zbl 0846.65047 · Zbl 0846.65047 · doi:10.2307/2153363
[4] J.L. Ericksen , Liquid crystals with variable degree of orientation . Arch. Ration. Mech. Anal. 113 ( 1991 ) 97 - 120 . Zbl 0729.76008 · Zbl 0729.76008 · doi:10.1007/BF00380413
[5] R. Glowinski , Numerical Methods for Nonlinear Variational Problems . Springer-Verlag, Berlin ( 1984 ). MR 737005 | Zbl 0536.65054 · Zbl 0536.65054
[6] N.G. Meyers , An \(L^{\mathrm p}\) estimate for the gradient of solutions of second order elliptic divergence equations . Ann. Scuola Norm. Sup. Pisa 17 ( 1963 ) 189 - 206 . Numdam | Zbl 0127.31904 · Zbl 0127.31904
[7] X. Xu , Existence for a model arising from the in situ vitrification process . J. Math. Anal. Appl. 271 ( 2002 ) 333 - 342 . Zbl 1011.35135 · Zbl 1011.35135 · doi:10.1016/S0022-247X(02)00112-9
[8] E. Zeidler , Nonlinear Functional Analysis and Its Applications , Vol. II/B. Springer, New York ( 1990 ). Zbl 0684.47029 · Zbl 0684.47029
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.