We analyse a nonlinear partial differential equation system describing the motion of a microswimmer in a nematic liquid crystal environment. For the microswimmer’s motility, the squirmer model is used in which self-propulsion enters the model through the slip velocity on the microswimmer’s surface. The liquid crystal is described using the well-established Beris–Edwards formulation. In previous computational studies, it was shown that the squirmer, regardless of its initial configuration, eventually orients itself either parallel or perpendicular to the preferred orientation dictated by the liquid crystal. Furthermore, the corresponding solution of the coupled nonlinear system converges to a steady state. In this work, we rigorously establish the existence of steady state and also the finite-time existence for the time-dependent problem in a periodic domain. Finally, we will use a two-scale asymptotic expansion to derive a homogenised model for the collective swimming of squirmers as they reach their steady-state orientation and speed.

Within the framework of the generalised Landau-de Gennes theory, we identify a Q-tensor-based energy that reduces to the four-constant Oseen–Frank energy when it is considered over orientable uniaxial nematic states. Although the commonly considered version of the Landau-de Gennes theory has an elastic contribution that is at most cubic in components of the Q-tensor and their derivatives, the alternative offered here is quartic in these variables. One clear advantage of our approach over the cubic theory is that the associated minimisation problem is well-posed for a significantly wider choice of elastic constants. In particular, this quartic energy can be used to model nematic-to-isotropic phase transitions for highly disparate elastic constants. In addition to proving well-posedness of the proposed version of the Landau-de Gennes theory, we establish a rigorous connection between this theory and its Oseen–Frank counterpart via a Г-convergence argument in the limit of vanishing nematic correlation length. We also prove strong convergence of the associated minimisers.

We investigate the defect structures around a spherical colloidal particle in a cholesteric liquid crystal using spectral method, which is specially devised to cope with the inhomogeneity of the cholesteric at infinity. We pay particular attention to the cholesteric counterparts of nematic metastable configurations. When the spherical colloidal particle imposes strong homeotropic anchoring on its surface, besides the well-known twisted Saturn ring, we find another metastable defect configuration, which corresponds to the dipole in a nematic, without outside confinement. This configuration is energetically preferable to the twisted Saturn ring when the particle size is large compared to the nematic coherence length and small compared to the cholesteric pitch. When the colloidal particle imposes strong planar anchoring, we find the cholesteric twist can result in a split of the defect core on the particle surface similar to that found in a nematic liquid crystal by lowering temperature or increasing particle size.

We propose an implicit finite-difference method to study the time evolution of the director field of a nematic liquid crystal under the influence of an electric field with weak anchoring at the boundary. The scheme allows us to study the dynamics of transitions between different director equilibrium states under varying electric field and anchoring strength. In particular, we are able to simulate the transition to excited states of odd parity, which have previously been observed in experiments, but so far only analyzed in the static case.

Here we investigate the kinematic transports of the defects in the nematic liquid crystal system by numerical experiments. The model is a shear flow case of the viscoelastic continuummodel simplified fromthe Ericksen-Leslie system. The numerical experiments are carried out by using a difference method. Based on these numerical experiments we find some interesting and important relationships between the kinematic transports and the characteristics of the flow. We present the development and interaction of the defects. These results are partly consistent with the observation from the experiments. Thus this scheme illustrates, to some extent, the kinematic effects of the defects.

We consider a class of non-quasi-convex frame indifferent energy densities that includes Ogden-type energy densities for nematic elastomers. For the corresponding geometrically linear problem, we provide an explicit minimizer of the energy functional satisfying a non-trivial boundary condition. Other attainment results, both for the nonlinear and the linearised model, are obtained by using the theory of convex integration introduced by Müller and Šverák in the context of crystalline solids.

We consider the dynamics of the director in a nematic liquid crystal when under the influence of an applied electric field. Using an energy variational approach we derive a dynamic model for the director including both dissipative and inertial forces.

A numerical scheme for the model is proposed by extending a scheme for a related variational wave equation. Numerical experiments are performed studying the realignment of the director field when applying a voltage difference over the liquid crystal cell. In particular, we study how the relative strength of dissipative versus inertial forces influence the time scales of the transition between the initial configuration and the electrostatic equilibrium state.

A continuum hydrodynamic model has been used to characterize flowing active nematics. The behavior of such a system subjected to a weak steady shear is analyzed. We explore the director structures and flow behaviors of the system in flow-aligning and flow tumbling regimes. Combining asymptotic analysis and numerical simulations, we extend previous studies to give a complete characterization of the steady states for both contractile and extensile particles in flow-aligning and flow-tumbling regimes. Another key prediction of this work is the role of the system size on the steady states of an active nematic system: if the system size is small, the velocity and the director angle files for both flow-tumbling contractile and extensile systems are similar to those of passive nematics; if the system is big, the velocity and the director angle files for flow-aligning contractile systems and tumbling extensile systems are akin to sheared passive cholesterics while they are oscillatory for flow-aligning extensile and tumbling contractile systems.

The ordered patterns formed by microphase-separated block copolymer systems demonstrate periodic symmetry, and all periodic structures belong to one of 230 space groups. Based on this fact, a strategy of estimating the initial values of self-consistent field theory to discover ordered patterns of block copolymers is developed. In particular, the initial period of the computational box is estimated by the Landau-Brazovskii model as well. By planting the strategy into the whole-space discrete method, several new metastable patterns are discovered in diblock copolymers.

In this paper, we investigate the effects of kinematic transports on the nematic liquid crystal system numerically and theoretically. The model we used is a “1+2” elastic continuum model simplified from the Ericksen-Leslie system. The numerical experiments are carried out by using a Legendre-Galerkin spectral method which can preserve the energy law in the discrete form. Based on this highly accurate numerical approach we find some interesting and important relationships between the kinematic transports and the characteristics of the flow. We make some analysis to explain these results. Several significant scaling properties are also verified by our simulations.