A renormalized Newton method for liquid crystal director modeling

Eugene C. Gartland, Jr, Alison Ramage

Research output: Contribution to journalArticle

4 Citations (Scopus)

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 linear systems 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., SIAM J. Sci. Comput., 35 (2013), pp. B226–B247]. Here we propose and analyze 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. This brings to light some anomalies of the Truncated Newton Method.
LanguageEnglish
Pages251-278
Number of pages28
JournalSIAM Journal on Numerical Analysis
Volume53
Issue number1
Early online date22 Jan 2015
DOIs
Publication statusPublished - 2015

Fingerprint

Newton-Raphson method
Truncated Newton Method
Newton Methods
Liquid Crystal
Liquid crystals
Modeling
Micromagnetics
Lagrange Method
Variational Model
Nonlinear Systems of Equations
Unit vector
Saddlepoint
Anomaly
Electric Field
Discretization
Linear Systems
Numerical Solution
Prototype
Iteration
Linear systems

Keywords

  • reduced hessian method
  • saddle-point problems
  • unit-vector constraints
  • director models
  • liquid crystals
  • renormalized newton method

Cite this

@article{d21f452c56694aa79220089c13343d88,
title = "A renormalized Newton method for liquid crystal director modeling",
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 linear systems 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., SIAM J. Sci. Comput., 35 (2013), pp. B226–B247]. Here we propose and analyze 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. This brings to light some anomalies of the Truncated Newton Method.",
keywords = "reduced hessian method, saddle-point problems, unit-vector constraints, director models, liquid crystals, renormalized newton method",
author = "{Gartland, Jr}, {Eugene C.} and Alison Ramage",
note = "{\circledC} 2015, Society for Industrial and Applied Mathematics",
year = "2015",
doi = "10.1137/130942917",
language = "English",
volume = "53",
pages = "251--278",
journal = "SIAM Journal on Numerical Analysis",
issn = "0036-1429",
number = "1",

}

A renormalized Newton method for liquid crystal director modeling. / Gartland, Jr, Eugene C.; Ramage, Alison.

In: SIAM Journal on Numerical Analysis, Vol. 53, No. 1, 2015, p. 251-278.

Research output: Contribution to journalArticle

TY - JOUR

T1 - A renormalized Newton method for liquid crystal director modeling

AU - Gartland, Jr, Eugene C.

AU - Ramage, Alison

N1 - © 2015, Society for Industrial and Applied Mathematics

PY - 2015

Y1 - 2015

N2 - 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 linear systems 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., SIAM J. Sci. Comput., 35 (2013), pp. B226–B247]. Here we propose and analyze 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. This brings to light some anomalies of the Truncated Newton Method.

AB - 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 linear systems 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., SIAM J. Sci. Comput., 35 (2013), pp. B226–B247]. Here we propose and analyze 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. This brings to light some anomalies of the Truncated Newton Method.

KW - reduced hessian method

KW - saddle-point problems

KW - unit-vector constraints

KW - director models

KW - liquid crystals

KW - renormalized newton method

UR - http://epubs.siam.org/loi/sjnaam

U2 - 10.1137/130942917

DO - 10.1137/130942917

M3 - Article

VL - 53

SP - 251

EP - 278

JO - SIAM Journal on Numerical Analysis

T2 - SIAM Journal on Numerical Analysis

JF - SIAM Journal on Numerical Analysis

SN - 0036-1429

IS - 1

ER -