Large Scale Lattice-Boltzmann Simulation of Liquid Crystals

Lead Research Organisation: University College London
Department Name: Chemistry

Abstract

In simple fluids, the atoms or molecules are disordered and can flow freely. In a crystal, they are arranged on a lattice, and cannot move past one another, creating a rigid material. Liquid crystals are in between: they are ordered in some senses and disordered in others. An example is the slimy mess you get when a bar of soap is left in a patch of water: this is a 'smectic' liquid crystal, in which the molecules pack into layers. Each layer lines up with the next in a crystalline stack, but in the plane of the layers the material is fluid. This is responsible for its 'slimy' feel. The sliminess can be quantified by measuring the material's 'rheology'. (Rheology is the science of flow behaviour.) Liquid crystals include many high tech materials used in laptop displays, flat-screen TVs, and other devices. In many of these devices, the flow of the material (for example in response to an electric field) is part of what makes the device work or not work. Many of these devices use 'nematic' liquid crystals in which rod-shaped molecules are lined up in the same direction but are not on a lattice; others involve 'cholesteric' or (potentially) 'blue phase' liquid crystals whose structure is more complex.For both scientific and technological reasons it is very important to understand properly the flow of liquid crystals in response to stresses and/or electric and magnetic fields. This is a very difficult task for two reasons. Firstly, there is the complicated, partially ordered structure to consider. Secondly, this structure is made even more complex in real materials by the presence of so-called 'defects'. These defects are of quite specific types, different in each type of liquid crystal. For nematics the defects are strange worm-like structures. (In fact, the name 'nematic' comes from the greek word for a worm.) In the simplest cases it is possible to solve using pen and paper the equations that describe the flow of pure liquid crystals, but when defects are present this is almost always impossible. The aim of the project is to develop and use methods for solving the relevant equations on very large computers. Only the biggest computers can provide the high resolution studies needed to address the problem of defects, since these are extended objects, large compared to the molecular scale. The work involves combining skill in simplifying the equations themselves (removing all inessential details from the description) with in-depth knowledge of how to make large computers solve such equations efficiently. For each type of liquid crystal, we plan to address both the way defects influence the flow behaviour, and the way a flow affects the organization of defects. This circle of influence is responsible for quite complex behaviour that is seen in the laboratory and, if understood, might be exploited in the next generation of liquid crystal technologies.

Publications

10 25 50
publication icon
Boghosian BM (2011) Unstable periodic orbits in the Lorenz attractor. in Philosophical transactions. Series A, Mathematical, physical, and engineering sciences

publication icon
Boghosian BM (2011) New variational principles for locating periodic orbits of differential equations. in Philosophical transactions. Series A, Mathematical, physical, and engineering sciences

publication icon
Saksena R (2009) Rheological response and dynamics of the amphiphilic diamond phase from kinetic lattice-Boltzmann simulations in Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences

publication icon
Saksena RS (2009) Real science at the petascale. in Philosophical transactions. Series A, Mathematical, physical, and engineering sciences

 
Description This project has allowed the successful simulation of the dynamics and rheology (flow and deformation behaviour) of mesophases (periodic fluid structures) observed in liquid crystals (locally liquid structures with long range orders) relevant to many technical and biomedical applications.
The employment of a kinetic lattice-Boltzmann method (LBM) for amphiphilic fluid mixtures (containing three phases including one surface active phase, e.g. water, oil and soap) allows for a bottom-up modelling approach from defined particle properties for the first time.
The comparison of simulations with experiment requires to reach realistic time and length scales of the response of liquid crystals to simple and periodic shear forces. Positive agreement establishes the numerical tools developed and applied as valid for the simulation of liquid crystal systems in complex geometries.
Various mesophases like the amphiphilic diamond, simple cubic, hexagonal, p- and gyroid phase have been successfully reproduced by the model. The successful comparison of the measured rheological response to application of shear forces to experimental results warrants confidence in the developed methods.
This scientific success was enabled by utilization of substantial HPC resources on the petascale (requiring in the order of a quadrillion floating point operations per second to compute). By this, building on technology developed in the RealityGrid project (GR/R67699, EP/C536452), this project allowed us to reach ground-breaking simulation scales solving computer science challenges.
In addition, beyond the original objectives the project has helped to inspire the investigation of chemically resolved nanocomposites (novell materials) including mesophases by coarse grained molecular dynamics simulations allowing for high impact contributions to international journals.
Exploitation Route The simulation application and tools developed in this project are of general interest for the modelling of liquid crystal systems for applications in technology, biomedical and material research applications. In particular, the method allows to investigate dynamics and rheology of mesophases in complex, intrinsically three-dimensional geometries not accessible by analytical, mathematical models.
The computer science developed in the scope of the project bridging the RealityGrid project (GR/R67699, EP/C536452) to the petascale provides valuable insights into large-scale computations on the pathway to the exascale.
Sectors Chemicals,Digital/Communication/Information Technologies (including Software),Manufacturing, including Industrial Biotechology,Pharmaceuticals and Medical Biotechnology

 
Description Scientific research on the petascale demonstrated by the applications of the code LB3D and data analysis tools developed in this project has enabled the inclusion of LB3D in the global Open Petascale Libraries application library of Fujitsu. Fujistsu UK funded 50% of the Ph.D. student Gary Davies subsequently working in the EPSRC Biocolloids project (EP/I034602) and high end computing consortium UKCOMES (EP/L00030X). The project furthermore informed decisions of Prof Coveney both in context of his chairing of the CPP steering panel and his work on the Strategy for the UK Research Computing Ecosystem. The latter lead to the establishment of the UK e-infrastructure leadership council which has been influential concerning a UK government investment volume of £500M into HPC resources over the period 2011-2016.
First Year Of Impact 2011
Sector Digital/Communication/Information Technologies (including Software),Government, Democracy and Justice
Impact Types Economic,Policy & public services

 
Description EPSRC Standard Research
Amount £519,283 (GBP)
Funding ID EP/I00887X/1 
Organisation Engineering and Physical Sciences Research Council (EPSRC) 
Sector Public
Country United Kingdom
Start 02/2012 
End 07/2015
 
Description EU CPCSA - Combined Collaborative Project and Coordination and Support Action
Amount € 3,272,777 (EUR)
Funding ID RI-261507 
Organisation European Commission 
Department Seventh Framework Programme (FP7)
Sector Public
Country European Union (EU)
Start 10/2010 
End 09/2013
 
Description EU FP7-INFRASTRUCTURES-2010-2
Amount £345,000 (GBP)
Funding ID 261507 
Organisation European Commission 
Sector Public
Country European Union (EU)
Start 10/2010 
End 09/2013
 
Description Qatar National Research Fund, as part of the National Priorities Research Program (NPRP).
Amount $1,200,000 (USD)
Organisation Qatar Foundation 
Sector Charity/Non Profit
Country Qatar
Start 11/2010 
End 10/2013
 
Title LB3D: Lattice-Boltzmann three dimensional simulation of fluids 
Description LB3D provides functionality to simulate three-dimensional simple, binary oil, water and ternary oil, water and amphiphile fluids using the Shan-Chen model for binary fluid interactions. The boundary conditions available include periodic boundaries, body forcing, and bounce-back boundaries as Lees-Edwards shearing for simple and binary fluid mixtures. The software is written in Fortran 90 and parallelized using MPI. It supports XDR and HDF5 format for I/O and provides checkpoint and restart for long-running simulations. It has been ported to many supercomputers worldwide, where it has shown excellent scalability. Most recently it has been shown to scale linearly on up to 294,000 cores on the European Blue Gene/P system Jugene. LB3D is developed by the Centre for Computational Chemistry at University College London, University of Stuttgart, and the Technical University of Eindhoven. EPCC has contributed to the development of LB3D most recently as part of the "Biocolloids" project funded by EPSRC. 
Type Of Technology Software 
Year Produced 2015 
Open Source License? Yes  
Impact LB3D has been used to study self-assembly of cubic phases, micro-mixing, flow through porous media, fluid surface interactions and other problems in complex fluidics. In addition, we have implemented the moment propagation method in LB3D, which allows us to efficiently measure the effective diffusion, dispersion, and other transport properties of ionic species in electrokinetic flows. These developments have enabled us to reveal interesting phenomena in charged rocks, namely a non-monotonic dispersion coefficient and a crossover in the effective diffusion of charged tracers. The publication of these novel scientific insights is currently under preparation [Schiller and Coveney]. This research has benefitted considerably from the generous allocation of computing time on ARCHER through UKCOMES. Dr. Schiller will continue to develop and use LB3D at Clemson University, USA, thus broadening the international scientific impact of LB3D. Furthermore, the LB3D code is one of the application codes in the European H2020 Compat project http://www.compat-project.eu/, where it will be used to enable bleeding edge simulations for computational materials science research on emerging exascale HPC systems. 
URL http://ccpforge.cse.rl.ac.uk/gf/project/lb3d