pith. sign in

arxiv: 2304.02062 · v4 · submitted 2023-04-04 · 🧮 math.NA · cs.NA

An A Posteriori Error Estimator for Electrically Coupled Liquid Crystal Equilibrium Configurations

J.H. Adler , D. B. Emerson This is my paper

Pith reviewed 2026-05-24 08:58 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords a posteriori error estimatorliquid crystalsFrank-Oseen modeladaptive finite elementspenalty methodelectrically coupled systemsmultilevel iteration
0
0 comments X

The pith

An a posteriori error estimator provides reliable global bounds and efficient local indicators for approximations of electrically coupled liquid crystal models using a penalty constraint.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives an a posteriori error estimator for the first-order optimality conditions of the electrically and flexoelectrically coupled Frank-Oseen model. It extends prior elastic-system results by incorporating a penalty term for the unit-length constraint and proves that the estimator bounds global approximation error while efficiently marking local error for refinement. Numerical tests in a multilevel nested-iteration setting show that adaptive meshes guided by the estimator reduce computational work while producing improved physical properties in the solutions. A reader would care because the estimator makes high-accuracy simulations of electrically driven liquid crystals feasible without uniform fine grids everywhere.

Core claim

The paper establishes an a posteriori error estimator for the nonlinear optimality system of the electrically and flexoelectrically coupled Frank-Oseen model under a penalty formulation of the unit-length constraint. Theory proves that the estimator is reliable for the global error and efficient for local error, allowing its direct use in adaptive refinement. Experiments confirm that the resulting adaptive multilevel scheme improves efficiency and yields better physical behavior than uniform refinement on challenging electrically coupled problems.

What carries the argument

The a posteriori error estimator constructed from the residual of the first-order optimality conditions of the penalty-regularized coupled Frank-Oseen model.

If this is right

  • The estimator can be inserted into existing multilevel nested-iteration solvers to drive adaptive mesh refinement.
  • Local error indicators produced by the estimator focus degrees of freedom on regions of high variation while leaving smooth regions coarse.
  • The same reliability and efficiency proofs apply to both global error control and local marking strategies.
  • Numerical tests demonstrate measurable reductions in total work and improvements in computed physical quantities such as defect structure.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The estimator could be combined with goal-oriented adaptivity to control specific quantities of interest such as switching voltages rather than the full energy norm.
  • Extension to three-dimensional domains would require only verification that the penalty parameter scaling remains stable under the same residual arguments.
  • The approach may transfer to other pointwise constraints in variational models of materials, such as incompressibility in liquid crystal polymers.

Load-bearing premise

Results previously shown for purely elastic systems extend directly to the electrically coupled case once a penalty term enforces the unit-length constraint.

What would settle it

A sequence of uniformly refined meshes on which the computed true error exceeds the estimator value by a factor that grows with refinement, or an adaptive sequence in which the estimator fails to reduce the true error at the expected rate.

Figures

Figures reproduced from arXiv: 2304.02062 by D. B. Emerson, J.H. Adler.

Figure 1
Figure 1. Figure 1: (a) Fine-mesh computed solution (restricted for visualization). (b) Resulting mesh patterns after four levels of AMR overlaid on the value of n1. In order to quantify efficiency, an approximate work unit (WU) is calculated for each simulation. A WU roughly approximates the work required by any full NI hierarchy in terms of assembling and solving a single linearization step for the Hessian on the finest uni… view at source ↗
read the original abstract

This paper derives an a posteriori error estimator for the nonlinear first-order optimality conditions associated with the electrically and flexoelectrically coupled Frank-Oseen model of liquid crystals, building on previous results for elastic systems. The estimator is proposed for a penalty approach to imposing the unit-length constraint required by the model. Moreover, theory is proven establishing that the estimator provides a reliable estimate of global approximation error and an efficient measure of local error, suitable for use in adaptive refinement. Numerical experiments demonstrate significant improvements in efficiency with adaptive refinement guided by the proposed estimator in a multilevel, nested-iteration framework and superior physical properties for challenging electrically coupled systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper derives an a posteriori error estimator for the nonlinear first-order optimality conditions of the electrically and flexoelectrically coupled Frank-Oseen model of liquid crystals. It employs a penalty approach for the unit-length constraint, proves that the estimator is reliable for the global approximation error and efficient for local error indicators, and demonstrates its performance in adaptive multilevel finite-element refinement for challenging electrically coupled configurations.

Significance. If the reliability and efficiency proofs hold for the coupled system, the estimator would provide a practical tool for adaptive simulation of electro-elastic liquid-crystal equilibria, extending prior elastic-only results to a setting with Maxwell coupling and flexoelectric torques. The numerical experiments already indicate efficiency gains and improved physical fidelity under adaptive refinement.

major comments (2)
  1. [§4] §4 (reliability analysis): the proof that the estimator remains reliable after augmentation by electric and flexoelectric terms relies on direct extension of the elastic residual bounds; the additional residuals arising from the electric potential equation and the flexoelectric torque (Eq. (2.3)–(2.5)) are not shown to be controlled by the same estimator quantities, so the global reliability constant may depend on the coupling strength.
  2. [§3.2, §4.2] §3.2 (estimator definition) and §4.2 (efficiency): the penalty term for the unit-length constraint introduces a consistency error whose interaction with the Maxwell equations is not bounded explicitly; if this term is absorbed into the elastic residuals without a separate estimate, the efficiency claim for local indicators may fail when the penalty parameter is moderate.
minor comments (2)
  1. [§2] Notation for the electric potential and its finite-element space should be introduced earlier (currently first appears in §2.3) to improve readability of the coupled weak form.
  2. [Figure 5] Figure 5 caption should state the penalty parameter value used in the adaptive runs; without it the comparison to uniform refinement is harder to interpret.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback on our manuscript. The two major comments raise valid points about the explicitness of certain bounds in the reliability and efficiency analyses. We address each below and will incorporate the necessary clarifications and additions in a revised version.

read point-by-point responses

  1. Referee: [§4] §4 (reliability analysis): the proof that the estimator remains reliable after augmentation by electric and flexoelectric terms relies on direct extension of the elastic residual bounds; the additional residuals arising from the electric potential equation and the flexoelectric torque (Eq. (2.3)–(2.5)) are not shown to be controlled by the same estimator quantities, so the global reliability constant may depend on the coupling strength.

    Authors: We agree that the current presentation of the reliability proof in §4 could be strengthened by making the control of the electric and flexoelectric residuals more explicit rather than relying on a direct extension argument. In the revised manuscript we will insert a short auxiliary lemma (between the existing elastic residual estimates and the final reliability theorem) that bounds these additional terms using the same estimator quantities, while explicitly tracking the dependence of the global constant on the coupling parameters. This does not change the validity of the result but improves clarity. revision: yes

  2. Referee: [§3.2, §4.2] §3.2 (estimator definition) and §4.2 (efficiency): the penalty term for the unit-length constraint introduces a consistency error whose interaction with the Maxwell equations is not bounded explicitly; if this term is absorbed into the elastic residuals without a separate estimate, the efficiency claim for local indicators may fail when the penalty parameter is moderate.

    Authors: The referee correctly notes that the interaction between the penalty consistency error and the Maxwell residuals is not isolated in the efficiency proof. We will revise §4.2 to include a separate estimate for this interaction term (using the same local indicator quantities) that remains valid for moderate penalty parameters. The revised efficiency statement will therefore hold uniformly with respect to the penalty parameter within the range used in the numerical experiments. revision: yes

Circularity Check

0 steps flagged

Builds on prior elastic results with independent proofs for coupled case; no reduction by construction

full rationale

The paper states it builds on previous results for elastic systems but explicitly claims and proves new theory for the electrically and flexoelectrically coupled Frank-Oseen model with penalty constraint. No quoted derivation step reduces the reliability/efficiency bounds or the estimator itself to a fitted input or self-citation chain by construction. The central reliability claim has independent mathematical content beyond the elastic citations, consistent with a minor non-load-bearing self-citation pattern.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters, axioms, or invented entities; the unit-length penalty and extension from elastic results are treated as background assumptions rather than new postulates.

pith-pipeline@v0.9.0 · 5630 in / 1086 out tokens · 18501 ms · 2026-05-24T08:58:35.751306+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [1]

    Adler, J.H., Atherton, T.J., Benson, T.R., Emerson, D.B., MacLachlan, S.P.: Energy minimization for liquid crystal equilibrium with electric and flexoelectric effects. SIAM J. Sci. Comput.37(5), S157–S176 (2015) 10 J.H. Adler and D. B. Emerson

    work page 2015

  2. [2]

    Adler, J.H., Emerson, D.B., MacLachlan, S.P., Manteuffel, T.A.: Constrained optimization for liquid crystal equilibria. SIAM J. Sci. Comput.38(1), B50–B76 (2016)

    work page 2016

  3. [3]

    Babuska, I., Rheinboldt, W.C.: A posteriori error estimates for the finite element method. Int. J. Numer. Meth. Engng (12), 1597–1615 (1978)

    work page 1978

  4. [4]

    ACM Trans

    Bangerth, W., Hartmann, R., Kanschat, G.: deal.II – a general purpose object oriented finite element library. ACM Trans. Math. Softw.33(4), 24/1–24/27 (2007)

    work page 2007

  5. [5]

    Bank, R.E., Welfert, D.B.: A posteriori error estimators for elliptic partial differential equations. Math. Comp. 44, 283–301 (1985)

    work page 1985

  6. [6]

    RAIRO Anal

    Clément, P.: Approximation by finite element functions using local regularization. RAIRO Anal. Numér. 2, 77–84 (1975)

    work page 1975

  7. [7]

    Davis, T.A., Gartland-Jr., E.C.: Finite element analysis of the Landau-de Gennes minimization problem for liquid crystals. SIAM J. Numer. Anal.35(1), 336–362 (1998)

    work page 1998

  8. [8]

    Döfler, W.: A convergent adaptive algorithm for Poisson’s equation. SIAM J. Numer. Anal.33, 1106– 1124 (1996)

    work page 1996

  9. [9]

    Emerson, D.B.: A posteriori error estimates for the Frank-Oseen model of liquid crystals. J. Coupled Syst. Multiscale Dyn.5(2), 95–110 (2017)

    work page 2017

  10. [10]

    Lagerwall, J.P.F., Scalia, G.: A new era for liquid crystal research: Applications of liquid crystals in soft matter, nano-, bio- and microtechnology. Curr. Appl. Phys.12(6), 1387–1412 (2012)

    work page 2012

  11. [11]

    Meyer, R.B.: Piezoelectric effects in liquid crystals. Phys. Rev. Lett.22(18), 918–921 (1969)

    work page 1969

  12. [12]

    Onsager, L.: The effects of shape on the interaction of colloidal particles. Ann. NY Acad. Sci.51, 627–659 (1949)

    work page 1949

  13. [13]

    Rojas-Gómez, O.A., Romero-Enrique, J.M., Silvestre, N.M., da Gama, M.M.T.: Pattern-induced an- choring transitions in nematic liquid crystals. J. Phys. Condens. Matter29(6) (2017)

    work page 2017

  14. [14]

    Computing64, 323–338 (2000)

    Starke, G.: Gauss-Newton multilevel methods for least-squares finite element computations of variably saturated subsurface flow. Computing64, 323–338 (2000)

    work page 2000

  15. [15]

    Finite element discretizations of elliptic equations

    Verfürth, R.: A posteriori error estimates for nonlinear problems. Finite element discretizations of elliptic equations. Math. Comp.62(206), 445–475 (1994)

    work page 1994

  16. [16]

    Wiley and Teubner (1996)

    Verfürth, R.: A Review of A Posteriori Error Estimation and Adaptive Mesh-Refinement Techniques. Wiley and Teubner (1996)

    work page 1996

  17. [17]

    Wan, Y., Zhao, D.: On the controllable soft-templating approach to mesoporous silicates. Chem. Rev. 107(7), 2821–2860 (2007)

    work page 2007