S. P. MacLachlan

J. H. Adler, D. B. Emerson, S. P. MacLachlan et al.

This paper investigates energy-minimization finite-element approaches for the computation of nematic liquid crystal equilibrium configurations. We compare the performance of these methods when the necessary unit-length constraint is enforced by either continuous Lagrange multipliers or a penalty functional. Building on previous work in [1,2], the penalty method is derived and the linearizations within the nonlinear iteration are shown to be well-posed under certain assumptions. In addition, the paper discusses the effects of tailored trust-region methods and nested iteration for both formulations. Such methods are aimed at increasing the efficiency and robustness of each algorithms' nonlinear iterations. Three representative, free-elastic, equilibrium problems are considered to examine each method's performance. The first two configurations have analytical solutions and, therefore, convergence to the true solution is considered. The third problem considers more complicated boundary conditions, relevant in ongoing research, simulating surface nano-patterning. A multigrid approach is introduced and tested for a flexoelectrically coupled model to establish scalability for highly complicated applications. The Lagrange multiplier method is found to outperform the penalty method in a number of measures, trust regions are shown to improve robustness, and nested iteration proves highly effective at reducing computational costs.

J. H. Adler, D. B. Emerson, P. E. Farrell et al.

Multiple equilibrium states arise in many physical systems, including various types of liquid crystal structures. Having the ability to reliably compute such states enables more accurate physical analysis and understanding of experimental behavior. This paper adapts and extends a deflation technique for the computation of multiple distinct solutions arising in the context of modeling equilibrium configurations of nematic and cholesteric liquid crystals. The deflation method is applied as part of an overall free-energy variational approach and is modified to fit the framework of optimization of a functional with pointwise constraints. It is shown that multigrid methods designed for the undeflated systems may be applied to efficiently solve the linear systems arising in the application of deflation. For the numerical algorithm, the deflation approach is interwoven with nested iteration, creating a dynamic and efficient method that further enables the discovery of distinct solutions. Finally, four numerical experiments are performed demonstrating the efficacy and accuracy of the algorithm in detecting important physical phenomena, including bifurcation and disclination behaviors. The final numerical experiment expands the algorithm to model cholesteric liquid crystals and illustrates the full discovery power of the deflation process.

J. H. Adler, I. Lashuk, S. P. MacLachlan et al.

This paper analyzes the discrete energy laws associated with first-order system least-squares (FOSLS) discretizations of time-dependent partial differential equations. Using the heat equation and the time-dependent Stokes' equation as examples, we discuss how accurately a FOSLS finite-element formulation adheres to the underlying energy law associated with the physical system. Using regularity arguments involving the initial condition of the system, we are able to give bounds on the convergence of the discrete energy law to its expected value (zero in the examples presented here). Numerical experiments are performed, showing that the discrete energy laws hold with order $\mathcal O\left(h^{2p}\right)$, where $h$ is the mesh spacing and $p$ is the order of the finite-element space. Thus, the energy law conformance is held with a higher order than the expected, $\mathcal{O}\left(h^p\right)$, convergence of the finite-element approximation. Finally, we introduce an abstract framework for analyzing the energy laws of general FOSLS discretizations.