Source author record

S. P. MacLachlan

S. P. MacLachlan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

4works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

4 published item(s)

preprint2016arXiv

A Deflation Technique for Detecting Multiple Liquid Crystal Equilibrium States

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.

preprint2014arXiv

An Energy-Minimization Finite-Element Approach for the Frank-Oseen Model of Nematic Liquid Crystals: Continuum and Discrete Analysis

This paper outlines an energy-minimization finite-element approach to the computational modeling of equilibrium configurations for nematic liquid crystals under free elastic effects. The method targets minimization of the system free energy based on the Frank-Oseen free-energy model. Solutions to the intermediate discretized free elastic linearizations are shown to exist generally and are unique under certain assumptions. This requires proving continuity, coercivity, and weak coercivity for the accompanying appropriate bilinear forms within a mixed finite-element framework. Error analysis demonstrates that the method constitutes a convergent scheme. Numerical experiments are performed for problems with a range of physical parameters as well as simple and patterned boundary conditions. The resulting algorithm accurately handles heterogeneous constant coefficients and effectively resolves configurations resulting from complicated boundary conditions relevant in ongoing research.

preprint2014arXiv

Constrained Optimization for Liquid Crystal Equilibria: Extended Results

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.

preprint2014arXiv

Energy Minimization for Liquid Crystal Equilibrium with Electric and Flexoelectric Effects

This paper outlines an energy-minimization finite-element approach to the modeling of equilibrium configurations for nematic liquid crystals in the presence of internal and external electric fields. The method targets minimization of system free energy based on the electrically and flexoelectrically augmented Frank-Oseen free energy models. The Hessian, resulting from the linearization of the first-order optimality conditions, is shown to be invertible for both models when discretized by a mixed finite-element method under certain assumptions. This implies that the intermediate discrete linearizations are well-posed. A coupled multigrid solver with Vanka-type relaxation is proposed and numerically vetted for approximation of the solution to the linear systems arising in the linearizations. Two electric model numerical experiments are performed with the proposed iterative solver. The first compares the algorithm's solution of a classical Freedericksz transition problem to the known analytical solution and demonstrates the convergence of the algorithm to the true solution. The second experiment targets a problem with more complicated boundary conditions, simulating a nano-patterned surface. In addition, numerical simulations incorporating these nano-patterned boundaries for a flexoelectric model are run with the iterative solver. These simulations verify expected physical behavior predicted by a perturbation model. The algorithm accurately handles heterogeneous coefficients and efficiently resolves configurations resulting from classical and complicated boundary conditions relevant in ongoing research.