Parallel implementation of the algorithm describing the behavior of liquid crystals under the action of electric field

Mathematical model for description of the dynamic processes in liquid crystals under weak mechanical, temperature and electric perturbations was proposed. In this work the parallel computational algorithm for analysis of the influence of electric field on a layer of a liquid-crystal medium is represented. In the exterior part of the layer, a potential of electric field satisfies the Laplace equation, which solution is constructed by the method of straight lines. Inside the layer, the equation taking into account anisotropy of a medium is solved by means of iterative method. In the presence of electric field, the molecules of the crystal experience the action of bulk forces and couple stresses, which can be used as right-hand sides of the basic equations of mathematical model. The algorithm is implemented using the CUDA technology for computer systems with graphics accelerators. Parallel program contains modules realizing the method of straight lines with the use of a three-point sweep, the fast Fourier transform and the SLAE solution method by means of the LU-decomposition, the iterative method of solving the equation for potential.