`International Journal of Mathematics and Mathematical SciencesVolume 2011, Article ID 746918, 8 pageshttp://dx.doi.org/10.1155/2011/746918`
Research Article

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

1Department of Mathematics and Computer Science, University of Salerno, Via Ponte Don Melillo, 84084 Fisciano (SA), Italy
2Department of Engineering Physics, Bashkir State University, 100 Mingazheva Street, Ufa 450014, Russia

Received 4 December 2010; Accepted 24 January 2011

Copyright © 2011 Aleksey Kudreyko and Nail Migranov. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

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.

#### 1. 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 of the liquid crystal cell is perpendicular to the director in the lower surface ; 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 and governed by van der Waals' forces [1, 2]. Here, is the azimuthal angle formed by the director and the -axis lying in the same plane, the -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 to an external electric field aligned parallel to the two bounding surfaces along the -axis, 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 .

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 -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 .

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.

#### 2. Problem Formulation

The dynamic equation based on the balance of elastic, electric, and hydrodynamic torques [2, 4] is written as follows: where , is the vacuum permittivity and -dielectric anisotropy coefficient of nematic LC, 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 where , -dimensionless time and , -dimensionless coordinate [2].

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 , where 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 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 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 (which was not fulfilled in [2]). Otherwise, the initial condition means that at the initial moment, , which is physically wrong. The location of the external disturbance 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].

#### 3. 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 is defined as follows: where integer numbers , show the approximation level of a wavelet is the dilation parameter. Index can be computed as follows: , and its minimum value is 2. The maximum value of is . The scaling function corresponds to In order to proceed the solution of (2.3), we must compute [3] the following integrals: 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: Also, let us introduce matrices with its elements , , and . It is necessary to underline that the dimension of these matrices equals to .

#### 4. Solution of the Base Equation

The unknown function in (2.3) is searched within the interval as the following series: where the wavelet coefficients represent coordinates of a vector of dimension . If differs from , then coefficients 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 . It thus appears the discretization of the space variable, and the plot of can be presented by means of scaling functions and without any losses of the corresponding data.

Since the Haar wavelets are defined for , then we must present (2.3) in new variables, that is, In view of formulae (4.2), (2.3) takes the form At this point we have everything to proceed to the study of the relaxation process in a TNC governed by the external electric field . Let us integrate (3.1) with respect to on interval , and twice with respect to on . Then, we get The boundary condition on the bottom plate yields , while the one on the top plate gives . Then, we obtain The substitution of (4.5) into (4.4) and the corresponding discretization of independent variables gives the following system of equations: The initial values of , 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.

Figure 1: Propagation of the traveling wave in a TNC for the moving coordinate system. Parameters of the method: , ; parameters of the problem: , , , , .

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.

Figure 2: Propagation of the traveling wave with respect to the laboratory coordinate system versus . Parameters of the method: , ; parameters of the problem: , , , , .

The criterion for relaxation of a director was chosen in the form: . 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.

#### 5. 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 space. Another reason for our choice of basis functions is low computational costs conditioned [8] by the sparsity of matrix .

#### Acknowledgments

The work of A. A. Kudreyko is supported by the Istituto Nazionale di Alta Matematica Francesco Severi (Rome-IT) under scholarship U 2010/000139, 1 October 2010. A. A. Kudreyko is thankful to Professor Carlo Cattani (University of Salerno) for many valuable discussions.

#### References

1. P. G. de Gennes and J. Prost, The Physics of Liquid Crystals, Oxford Science, 2nd edition, 1995.
2. A. V. Zakharov and A. A. Vakulenko, “Relaxation of the director field in the form of a traveling wave in twisted nematic cells,” Physics of the Solid State, vol. 50, no. 3, pp. 573–578, 2008.
3. Ü. Lepik, “Numerical solution of evolution equations by the Haar wavelet method,” Applied Mathematics and Computation, vol. 185, no. 1, pp. 695–704, 2007.
4. A. V. Zakharov and A. A. Vakulenko, “Intrinsic torsional reorientations in a twisted nematic liquid crystal cell,” Physical Review E, vol. 72, no. 2, Article ID 021712, 7 pages, 2005.
5. C. Cattani and J. Rushchitsky, Wavelet and Wave Analysis as Applied to Materials with Micro or Nanostructure, vol. 74 of Series on Advances in Mathematics for Applied Sciences, World Scientific, Hackensack, NJ, USA, 2007.
6. C. Cattani, “Multiscale analysis of wave propagation in composite materials,” Mathematical Modelling and Analysis, vol. 8, no. 4, pp. 267–282, 2003.
7. Ü. Lepik, “Numerical solution of differential equations using Haar wavelets,” Mathematics and Computers in Simulation, vol. 68, no. 2, pp. 127–143, 2005.
8. I. Daubechies, Ten Lectures on Wavelets, vol. 61 of CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, Pa, USA, 1992.