Picture map of Europe with pins indicating European capital cities

Open Access research with a European policy impact...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by Strathclyde researchers, including by researchers from the European Policies Research Centre (EPRC).

EPRC is a leading institute in Europe for comparative research on public policy, with a particular focus on regional development policies. Spanning 30 European countries, EPRC research programmes have a strong emphasis on applied research and knowledge exchange, including the provision of policy advice to EU institutions and national and sub-national government authorities throughout Europe.

Explore research outputs by the European Policies Research Centre...

Local stability and a renormalized Newton Method for equilibrium liquid crystal director modeling

Gartland, Jr, E.C. and Ramage, Alison (2012) Local stability and a renormalized Newton Method for equilibrium liquid crystal director modeling. Working paper. University of Strathclyde, Glasgow.

[img] PDF
9_strathrepar_2012.pdf

Download (253kB)

Abstract

We consider the nonlinear systems of equations that result from discretizations of a prototype variational model for the equilibrium director field characterizing the orientational properties of a liquid crystal material. In the presence of pointwise unit-vector constraints and coupled electric fields, the numerical solution of such equations by Lagrange-Newton methods leads to problems with a double saddle-point form, for which we have previously proposed a preconditioned nullspace method as an effective solver [A. Ramage and E. C. Gartland, Jr., submitted]. The characterization of local stability of solutions is complicated by the double saddle-point structure, and here we develop efficiently computable criteria in terms of minimum eigenvalues of certain projected Schur complements. We also propose a modified outer iteration (“Renormalized Newton Method”) in which the orientation variables are normalized onto the constraint manifold at each iterative step. This scheme takes advantage of the special structure of these problems, and we prove that it is locally quadratically convergent. The Renormalized Newton Method bears some resemblance to the Truncated Newton Method of computational micromagnetics, and we compare and contrast the two.