Three Dimensional Director Modelling
The computer modelling group at UCL has many years experience
in Electromagnetics, and Liquid Crystal Modelling. As a product of
this work a program has been developed, which is applicable to a wide
range of LC geometries. A CAD program is used to produce the mesh,
supporting complex shaped volumes, and curved alignment surfaces and
electrodes.
A variational approach to the Oseen-Frank free energy is taken, discretised
using finite elements on a mesh of tetrahedral elements. A variable
time step is used to accelerate the computation.
- Three elastic constants
- Flexoelectricity
- Vectorial or Tensorial Description of the LC director
- Weak anchoring
- Multiple alignment surfaces and electrodes
- Multiple dielectrics and LC domains
- Arbitrary voltage waveforms.
- Neumann or Periodic Boundary conditions
- Fast GMRES potential solver
The solution for the
director is alternated within each time step with the calculation
of the potential, and is iterated until a consistent solution is
achieved.
A simple interface provides a means to configure
the material properties, and gives a graphical indication of the
progress of the simulation.
Ion transport within the the LC can be taken into
account in this model, by algorithms produced at the University
of Gent. In addition a constant charge formulations is available,
as well as two-dimensional finite element and finite difference
implementations.
If you are interested in obtaining any of the programs
detailed here, please feel free to email me.
If you have any problems using the program please
check the troubleshooting page which details
some common mistakes.
Defining the Cell Geometry
WinLCD runs on a mesh of tetrahedral elements. A choice
of two mesh generator is available GID and Tetgen.
GID has a CAD like interface which allows complex structure to
be drawn. An academic version can be downloaded and used for free,
but there is a limit on the number of elements it can generate,
and so it is only feasible to model simple structures. The professional
version removes this limit, and is reasonably priced.
http://gid.cimne.upc.es/download/
An alternative is to use TetGen which uses input files as opposed
to a GUI which makes inputting of complex structures more difficult,
although it can take surface meshes from programs such as GID. It
produces high quality meshes and the author, Hang Si, is responsive
to questions and is very helpful. Tetgen is free, but has its own
licence restrictions which can be viewed on the TetGen website.
There is no limit on the number of elements it can produce.
The documentation and source code for TetGen can be found at the
following location:
http://tetgen.berlios.de
Creating a Planar Cell and Meshing
Complex 3D geometries can be modelled, but to demonstrate
how to use the programs we begin with a simple planar structure.
Before the physical structure is defined it is usual to introduce
the types of volumes supported the program as well as the boundary
conditions that can be applied to surfaces. A material number is
assigned to all surfaces and volumes that we define.
There are two supported types of volumes dielectric and LC volumes.
A structure may contain up to 7 different types of dielectric and
7 types of LC material!
There are more types of surfaces that can be defined. An LC device
will have alignment layers which are treated in some way, giving
a preferred direction of orientation at the surface. A surface may
be specified as an alignment surface by applying a material number
FixLC-x, where x={1,2,..}. There will be some energy penalty associated
with deflection from this preferred direction, the magnitude of
which depends on the surface treatment. The nature of the alignment
on the surface is specified in elastic.i, which will be discussed
in more detail in the next section.
Devices will also contain electrodes to reorient the LC. Any surface
can be assigned as an electrode using the material names Electrode-x.
In a real device the alignment layer is a thin insulating layer
that separates the LC from the electrodes or glass substrate. If
the alignment layer is sufficiently thin compared to the thickness
of the LC layer it is possible to fix the LC director and pin the
potential on the same surface mesh. To describe such surfaces will
dual behaviour it is necessary to be able to combine material name,
e.g. FixLC-1_Electrode-2. This method is necessary because two material
numbers can't be applied to a surface.
In GID the user only needs to know the material names discussed
above. So that GID knows these names, a problem_type is required,
which can be downloaded here.
This zip file contains a directory that should be copied into the
GID problem type directory. In Tetgen the user is required to know
the actual material number, which can be derived from the material
number using the table below, which specifies the bits of the number
that represent each property. For example if bit 128 only of the
material number is set, the corresponding name is Electrode-2
32768 |
... |
2048 |
1024 |
... |
64 |
32 |
16 |
8 |
4 |
2 |
1 |
Alignment Surface
e.g.
FixLC-1 = 2048 |
Electrode
e.g.
Electrode-1 = 64
Electrode-2 = 128 |
Dielectric(1) / Domain(0) |
Volume Number
e.g. Domain-1 = 4 |
Neumann |
Periodic |
There structure must fulfull one requirement. The program expects
the LC layer to extend in x-y plane. This assumption is in face
only used when setting the initial director configuration and will
probably be removed in a future revision. For a simple structure
this means that aliment surfaces will be normal to z.
The following form can be used to find the material numbers (requires
JavaScript):
Surface Material Number |
0 |
To create a simple structure with GID click here
For more GID tips click here. For
more advanced structures click here
To create a simple structure with Tetgen click here
To create a more complex IPS stucture in Tetgen click here
After following these instruction the mesh will have been generated.
A text file mesh.txt will be obtained containing
a list of node, tetrahedra and triangles.
Input
Files: Material Properties and Voltage Waveforms
Elastic.i
v2
Material Parameters Input File
=====================================================
z Scaling
Factor (multiplies all z mesh coordinates)
0.5
=====================================================
DOMAINS
=====================================================
-----------------------------------------------------
Domain 1
-----------------------------------------------------
k11[pN] 12.0
k22[pN] 12.0
k33[pN] 12.0
e11 0.0
e33 0.0
gamma 0.09
delta_epsilon 6.8
epsilon_perpen 5.0
pitch[um] 0.0
=====================================================
ALIGNMENT SURFACES
=====================================================
-----------------------------------------------------
FixLC 1
-----------------------------------------------------
anchoring strong (strong/weak)
pretilt 2.0
pretwist 80.0
weak anchoring parameters
we 50.0
wn 50.0
psi 0.0 (default zero)
=====================================================
DIELECTRICS
=====================================================
-----------------------------------------------------
Dielectric 1
-----------------------------------------------------
delta_epsilon 0.0
epsilon_perpen 4.0
tilt 0.0
twist 0.0
|
Isotropic elements are created by most mesh generators,
and since LC cells are commonly long and thin, it is useful
to create a thicker structure and scale it down later.
Each domain has its own elastic constants and viscosity (gamma
[Pa.s])
Domain 1 is applied to any volumes defines as Domain 1 in
GID or Tetgen. You can add more domains if you want, just
copy and paste the domain section and change the title to
Domain 2.
Alignment surfaces can either have strong or weak anchoring.
For weak anchoring there is a different energy penalty associated
with azimuthal and zenithal deformations. Hence, there are
two anchoring coefficients: we and wn, which correspond to
w_\eta and w_\zeta of [1]. These can be related to the Raphini
Papular coefficients.
The units of the anchoring coefficents are
*10^{-6} Jm^{-2}
The twist and tilt are relative to x, i.e. the twist is atan(y/x)
and the tilt is asin(z)
Anisotropic dielectrics can be entered into the program,
hence the entries for the twist and the tilt
|
VOLTDATA.i
v2
Voltage Functions Definition Input File
=====================================================
Result
File to Commence Run from: (or 'none')
none
=====================================================
VOLTAGE FUNCTIONS
define the number of shots then list the voltage
then the length of each pulse in seconds
=====================================================
-----------------------------------------------------
Electrode 1
-----------------------------------------------------
nshots 1
voltage 5.0 0.0
length 0.050 0.100
=====================================================
OTHER PARAMETERS
=====================================================
ground potential[V]
0.d0
time step[s]
5d-6
results look time[s] (start:step:end)
20.0d-6
0.001:0.001:0.150 |
If a result file name is entered the program starts
from the director configuration stored in the file, and the
initial time is set to the time stored in the result file
Define the number of time the function repeats.
Define the voltage and the length of each pulse (in seconds)
This example applies an 50ms 5V pulse, followed by a 100ms
0V pulse.
Entries for more electrodes can be added as required
Keep the time step at about this value.
Results are output at these times. MATLAB style formatting
allowed.
|
Running the Simulation
A simulation directory must contain 4 files
- elastic.i: Material properties
- voltdata.i: Voltage waveforms
- mesh.txt: Tetrahedral Mesh
- WinLCD.exe: Simulation program
The first three files should be configured as described above,
then load up WinLCD, press Play and wait! As the program runs it
will output result files at the result look times specfied in voltdata.i.
The file names will be of the form resultxxx.dat.
These files contain the director and potential values at each node
of the mesh. The program will run until the last look time is reached
or until the user closes the program. If the program is accidenttally
closed, a simulation can be restarted from an existing result file.
Just modify the 6th line of voltdata.i.
Viewing Results
The result files can be opened up in Matlab using the code available
here, or the Visualisation Program,
which can be found here.
[1] W. Zhao, C.-X. Wu, and M. Iwamoto, “Analysis of weak-anchoring
effect in nematic liquid crystals,” Phys. Rev. E, vol. 62,
no. 2, pp. 1481–1484, 2000.
This page last modified
10 November, 2006
by r.james
|