Made available in DSpace on 2015-05-14T12:14:03Z (GMT). No. of bitstreams: 1 parte1.pdf: 6372308 bytes, checksum: db3e915edd1663a97d16d8935fc5becf (MD5) Previous issue date: 2012-03-02 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / In this work we study numerically the generation and dynamics of topological defects in nematic liquid crystals. Our study is based on a Ginzburg-Landau model describing the evolution of the orientational order of a liquid crystal in terms of a symmetric, traceless, second-rank tensor. This phenomenological model allows studies of nematic phases at scales ranging from few nanometers to few micrometers (mesoscopic scale). Within this framework we developed a software named LICRA (Liquid CRystal Algorithm) that combines standard finite difference algorithm for the spatial derivatives with a Runge-Kutta temporal integration to solve the relaxational equations of nematodynamics without thermal fluctuations and hydrodynamic flow. Using this software we investigate the coarsening dynamics of defects of two- and three-dimensional uniaxial nematic liquid crystals. The time dependences of the structure factor and characteristic length scale were computed. The characteristic length scale is expected to grow as a power law in time, L ∝ tα. From dimensional analysis α = 1/2 and we found α = 0, 45±0, 01 in two-dimensions and α = 0, 350±0, 003 in three-dimensions. Furthermore, in all cases Porod s law is satisfied for large values of wave number k. We also investigate, using LICRA, the coarsening dynamics of liquid crystal textures in a two-dimensional nematic under applied electric fields. We consider both positive and negative dielectric anisotropies and two different possibilities for the orientation of the electric field parallel and perpendicular to the two-dimensional lattice. We determine the effect of an applied electric field pulse on the evolution of the characteristic length scale and other properties of the liquid crystal texture network. In particular, we show that different types of defects are produced after the electric field is switched on, depending on the orientation of the electric field and the sign of the dielectric anisotropy. Finally, we present the effect of the rotation of an external electric field on the dynamics of half-integer disclination networks in two and three dimensional nematic liquid crystals with a negative dielectric anisotropy. We show that a rotation of π of the electric field around an axis of the liquid crystal plane continuously transforms all half-integer disclinations of the network into disclinations of opposite sign via twist disclinations. We also determine the evolution of the characteristic length scale, thus quantifying the impact of the external electric field on the coarsening of the defect network. / Neste trabalho estudamos numericamente a formação e dinâmica de defeitos topológicos em cristais líquidos nemáticos. Nosso estudo é baseado no modelo de Ginzburg- Landau, o qual descreve a evolução da ordem orientacional de um cristal líquido em termos de um tensor de segunda ordem simétrico e com traço nulo. Este modelo fenomenológico permite estudar a fase nemática em escalas que vão de poucos nanômetros até poucos micrômetros (escala mesoscópica). Para tal estudo numérico, desenvolvemos um programa de computador que denominamos de LICRA (Liquid CRystal Algotithm). Este programa combina o algoritmo de diferença finita para calcular derivadas espaciais com a integração temporal de Runge-Kutta para resolver a equação de relaxação da nematodinâmica, sem a presença de flutuações térmicas e fluxos hidrodinâmicos. Usando este programa de computador investigamos a dinâmica de coalescência em duas e três dimensões em um cristal líquido nemático uniaxial. Tanto o fator de estrutura quando a escala de comprimento característico foram calculadas no tempo. Espera-se que esta escala cresça como uma lei de potências do tempo, L ∝ tα, onde, a partir de uma análise dimensional, α = 1/2. Encontramos os valores de α = 0, 45 ± 0, 01 em duas dimensões e α = 0, 350 ± 0, 003 em três dimensões. Além disso, em todos os casos verificamos que a lei de Porod é satisfeita para número de ondas k de grandes valores. Utilizando LICRA, investigamos também a dinâmica de coalescência de cristais líquidos nemáticos em duas dimensões submetidos a um campo elétrico externo. Consideramos a anisotropia dielétrica positiva e negativa e duas diferentes possibilidades de orientação do campo elétrico: paralelo e perpendicular ao plano da rede bidimensional. Determinamos os efeitos de um pulso de campo elétrico na evolução da escala do comprimento característico e as alterações nas texturas dos cristais líquidos. Em particular, mostramos que os diferentes tipos de defeitos que são produzidos após o campo elétrico ser aplicado dependem da orientação do campo elétrico e do sinal da anisotropia dielétrica. Finalmente, apresentamos os efeitos da rotação de um campo elétrico externo na dinâmica de uma rede de defeitos semi-inteiros em cristais líquidos nemáticos em duas e três dimensões com anisotropia dielétrica negativa. Mostramos que, girando o campo elétrico por um ângulo π ao redor de um eixo pertencente a plano da rede, ocorre uma transformação contínua de todas as desclinações semi-inteiras da rede em desclinações com sinal oposto. Esta transformação é intermediada por desclinações do tipo torção. Além disso, determinamos a evolução da escala de comprimento característico quantificando o impacto do campo elétrico externo na dinâmica de coalescimento da rede.

No description available.

In this thesis, we have tried to explain the regular and irregular(chaotic) dynamics of worm like micellar solutions on applying shear, through a detailed study of the equation of motion of a nematic order parameter tensor coupled to a hydrodynamic velocity field. We have assumed spatial variations only along one direction i.e. the gradient direction(1D model). The resulting phase diagram shows various interesting steady states or phases such as spatiotemporal chaos, temporal and spatiotemporal periodicities, and alignment of the director axis along the imposed flow field. The coupling of the orientational degrees of freedom of the order parameter with the hydrodynamic flow field holds the key to the appearance of dynamic shear bands in the system. We have solved numerically a set of coupled nonlinear equations to obtain the order parameter stress developed in the system; the magnitude of the order parameter tensor, the biaxiality parameter and the orientation of the director axis of the nemato gens under shear have also been studied in detail. To study the phase diagram obtained by time integration of the equation of motion mathematically, a stability analysis of the fixed point of motion for various parameter values has been performed so that the location of the chaotic-to-aligned phase boundary is verified. Also in the periodic region of the phase diagram, the stability of limit cycles is tested by analysing the fixed point of the corresponding Poincare map. Stability analysis of the periodic orbits leads to the observation that in the parameter space, there are regions of phase coexistence where chaotic or spatiotemporally intermittent behaviour coexists with periodic behaviour. When corrections in the imposed velocity field due to the order parameter stress were taken into account and the order parameter response was looked into at several points in the parameter space, the modified equations of motion were found to reproduce the earlier behaviour in all the different regimes if the value of a dimensionless viscosity parameter is taken to be such that the bare viscous stress overrides the order parameter stress. The phase boundaries are however different from the ones seen in the earlier model. However, for a choice of the viscosity parameter such that the order parameter stress and the bare viscous stress are comparable, we see two distinctly different attractors: a banded, periodic one that is common to both α1equalto 0, and not equal to 0 and a banded chaotic one for α1not equal to 0. Here, α1is a parameter that governs the nonlinearity in the stretching of the order parameter tensor along the direction of the applied shear. Quantitative analysis of the various chaotic attractors throws up not only positive Lyapunov exponents but also that the banded chaos is a “flip-flop” kind of chaos where the switching between two long-lived states of high and lows hear stress is chaotic, where as the behaviour in either of the two states is periodic, with either a single, isolated frequency or a bunch of harmonics. Also, the spatial correlation of the shear stress in the chaotic attractors is of much larger range than the temporal correlation, the latter being almost delta-function-like. On increasing the temperature of the system till it is above the isotropic–nematic transition temperature in the absence of shear, we find that under shear, similar attractors as those in the nematic case are observed, both for passive advection and for the full 1D hydrodynamics. This is an encouraging result since actual experiments are performed at a temperature for which the system is in the isotropic phase in the absence of shear. Thus for the 1D system, the parameter space has been explored quite extensively. Considering spatial variations only along the gradient axis of the system under shear is not enough since experiments have observed interesting behaviour in the vorticity plane in which Taylor velocity rolls were noted. Hence taking the system to 2D was necessary. Our numerical study of the 2D system under shear is incomplete because we came across computational difficulties. However, on shorter time scales we have seen a two-banded state with an oscillating interface and Taylor velocity rolls as well. The methodology used for the 2D study can also be used to reproduce the 1D results by the simple step of taking initial condition with no variation in the vorticity direction. This automatically ensures that no variation in the vorticity direction ever builds up because the equations of motion ensure that these variations in the system do not grow by themselves unless fed in at the start. Using this method, we were able to reproduce all the attractors found in the 1D calculation. Thus the 1D attractors have been observed using two different methods of calculation. Further work on the full 2D numerics needs to be done because we believe that spatiotemporally complex steady-state attractor s exist in the 2D system also for appropriate values of the parameters.

No description available.

No description available.

No description available.