A New Solution for the Director Relaxation Problem in Twisted Nematic Film Based on Wavelet Analysis

. The present research considers and explains the application of the Haar wavelets as basis functions in solution of the propagation of perturbations in low-dimensional anisotropic media. The computations of the relaxation problem in the form of traveling waves have shown that the present approach possesses several advantages over regular methods.


Introduction
Not to mention that the study of such non-Newtonian fluids as liquid crystals LC has lead to the commercial production of widely used electronic devices, the theory of LC still has a few problems, whose solution presents a difficulty.The problem of a director relaxation in twisted nematic cells TNCs in the form of nonlinear waves still remains unsolved.TNCs are integral parts of liquid crystal displays, which are used in laptops and personal computers.TNC is a drop of LC placed between two parallel surfaces, separated by a spacer, and arranged in such a way that the director in the upper surface n of the liquid crystal cell is perpendicular to the director in the lower surface n − ; in this case, the director always lies in the plane parallel to the two bounding surfaces.In the absence of an external electric field, a linear director field distribution is characterized by the azimuthal angle Φ π/2 z and governed by van der Waals' forces 1, 2 .Here, Φ is the azimuthal angle formed by the director and the X-axis lying in the same plane, the Z-axis is perpendicular to the bounding surfaces, and the origin of the coordinate system is located at the lower surface of the twisted nematic cell.In response International Journal of Mathematics and Mathematical Sciences to an external electric field E Ej aligned parallel to the two bounding surfaces along the Yaxis, the twist re-orientation of the director field occurs toward the external electric field.The threshold value of the external electric field, when nonlinear relaxation mechanisms of the director field can occur, we will denote as E cr .
The use of TNC in production of LC monitors provided a strong encouraging stimulus for intensive research in this field.The present paper shows that angular moments acting on a unit element of nematic LC can generate traveling waves along the Z-axis.This problem can be easily reduced to the solution of dimensionless parabolic partial differential equation.
The goal of the paper is to present the director relaxation problem using the Haar wavelet method 3 .The earlier effort for the solution of this problem was realized in 2 .However, the considered mathematical model does not take into account the perpendicularity of the directors' orientations located in the plane of bounded plates at the initial moment t 0.
Methods of wavelet analysis adapted for the solution of partial differential equations have been actively developing in the past 15 years.The advantages of the application of wavelets for the solution of differential problems over regular numerical methods consist of linearly increasing computational costs because of the sparsity of intermediate matrices.Moreover, the quality of solutions obtained by analytically defined basis functions always exceeds classical numerical methods.

Problem Formulation
The dynamic equation based on the balance of elastic, electric, and hydrodynamic torques 2, 4 is written as follows: where Δ ε 0 ε a E 2 /2, ε 0 is the vacuum permittivity and ε a -dielectric anisotropy coefficient of nematic LC, K 2 is the twist elastic constant.
In order to study the problem of propagation of the traveling wave from the top plate to the bottom plate, it is convenient to present 2.1 with respect to dimensionless coordinates The purpose of our approach consists of the study of the generation process of a traveling wave running between bounded plates by means of wavelet basis.In the coordinate system, connected with the motion of disturbance in the media, we can introduce q * z − vτ, where v is the dimensionless velocity of wave propagation.Taking into account the rule of differentiation for complex functions, the left-hand side of 2.2 generates two terms.Thus, this equation takes the form vΦ q * τ * , q * Φ q * ,q * τ * , q * 1 2 sin 2Φ τ * , q * Φ τ * τ * , q * , q * ∈ q in , q fin , τ ≥ τ in .

2.3
In the case of strong anchoring of LC molecules with bounding surfaces, the torque balance for these surfaces leads to the following boundary and initial conditions: where H presents the Heaviside step function.The main feature of such an initial condition is that the solution strongly satisfies the upper boundary condition by adding sin which was not fulfilled in 2 .Otherwise, the initial condition means that at the initial moment, Φ τ, q * q * q * 2 ≈ 0, which is physically wrong.The location of the external disturbance q * 3 8 is set at our choice.The number 8.9999 means that the Gaussian cannot belong to two intervals of the Heaviside step function.Such form of the initial condition can be generated by a laser pulse 2 .

Haar Wavelets
From all families of wavelets, we can distinguish a few easily defined functions, which represent a powerful tool in the solution of partial differential equations.Such wavelets are: Shannon wavelets 5 , harmonic wavelets 6 , and Haar wavelets.The latter wavelets are wellstudied in 3, 7, 8 and references therein.
The family of Haar wavelets for q ∈ 0; 1 is defined as follows: where integer numbers m 2 j , j 0, 1, . . ., J show the approximation level of a wavelet k 0, 1, . . ., m − 1 is the dilation parameter.Index i can be computed as follows: i m k 1, International Journal of Mathematics and Mathematical Sciences and its minimum value is 2. The maximum value of i is 2M 2 J 1 .The scaling function corresponds to i 1 1, q ∈ 0, 1 , 0, elsewhere.

3.2
In order to proceed the solution of 2.3 , we must compute 3 the following integrals:

3.3
Thus, Haar wavelets are transformed into piecewise polynomials of the second order.To get the solution of 2.3 with initial and boundary conditions 2.4 , we will employ the collocation method.The idea of this method consists of the following: differential equation 2.3 is fulfilled in equally spaced points the so-called collocation points within the space bounded by the plates.Define the location of collocation points as follows: , l 1, 2, . . ., 2M.

3.4
Also, let us introduce matrices H, P, R with its elements H i, l h i q l , P i, l p i q l , and R i, l r i q l .It is necessary to underline that the dimension of these matrices equals to 2M × 2M.
International Journal of Mathematics and Mathematical Sciences 5

Solution of the Base Equation
The unknown function in 2.3 is searched within the interval τ ∈ τ s , τ s 1 as the following series: where the wavelet coefficients a s i represent coordinates of a vector of dimension 2M.If τ differs from τ s , τ s 1 , then coefficients a s i change their values.
It is obvious to present a question about our choice of basis functions.As we know, many other bases exist, which could be suitable for the solution of this problem.The answer to this question should be originated from the limited opportunities of disclosure of the angle between two directors located in parallel planes, and separated by distance Δz Δq * .It thus appears the discretization of the space variable, and the plot of Φ τ s , q can be presented by means of scaling functions p i q and r i q without any losses of the corresponding data.

4.3
At this point we have everything to proceed to the study of the relaxation process in a TNC governed by the external electric field E Ej.Let us integrate 3.1 with respect to τ on interval τ s , τ s 1 , and twice with respect to q on 0; q .Then, we get Propagation of the traveling wave Φ τ, q in a TNC for the moving coordinate system.Parameters of the method: J 7, Δτ 0.005; parameters of the problem: v 2, q 1 0, 2q 2 10, σ 0.25, q 3 8.

4.6
The initial values of Φ 1 l , Φ 1 l , and Φ 1 l can be computed from the initial condition 2.4 and the corresponding change of variables 4.2 .
The solution of 2.3 with initial and boundary conditions 2.4 in the moving coordinate system was obtained by the Haar wavelet method and was illustrated in Figure 1.The initial disturbance of the director's field in a TNC was generated close to the top plate.7, Δτ 0.005; parameters of the problem: q 1 7 − vτ, q 2 10, v 2, σ 0.25, E/E cr 10.
The simulations have shown that this approach for the problem of relaxation of a director toward its equilibrium position allows its the computation of different dynamic modes in TNC.The relaxation of a director with respect to the laboratory coordinate system is shown in Figure 2.
The criterion for relaxation of a director was chosen in the form: ε |Φ τ s 1 − Φ τ s | ≤ 10 −4 .In other words, if the azimuthal angle between two time steps is less than ε, then the algorithm stops.It is worth emphasizing that the Haar wavelet method permits to get solutions with a smaller error.However, in this case, lower values of ε do not play any significant role and only increase the CPU load.

Conclusions
The present paper describe the phenomenon of spatial relaxation in a TNC for the case of strong anchoring of LC molecules with the bounded plates.It was suggested to employ the Haar wavelet method for getting a function, which describes the evolution of the azimuthal angle.The advantage of such an approach over regular methods is that the unknown function is searched as a superposition of analytically defined functions, which constitute a complete basis in the L 2 0; 1 space.Another reason for our choice of basis functions is low computational costs conditioned 8 by the sparsity of matrix a s i .