An official website of the United States government
Abstract A local discontinuous Galerkin (LDG) method for approximating large deformations of prestrained plates is introduced and tested on several insightful numerical examples in Bonito et al. (2022, LDG approximation of large deformations of prestrained plates. J. Comput. Phys., 448, 110719). This paper presents a numerical analysis of this LDG method, focusing on the free boundary case. The problem consists of minimizing a fourth-order bending energy subject to a nonlinear and nonconvex metric constraint. The energy is discretized using LDG and a discrete gradient flow is used for computing discrete minimizers. We first show $$\varGamma $$-convergence of the discrete energy to the continuous one. Then we prove that the discrete gradient flow decreases the energy at each step and computes discrete minimizers with control of the metric constraint defect. We also present a numerical scheme for initialization of the gradient flow and discuss the conditional stability of it.
more »
« less
Optimal Control of the Landau–de Gennes Model of Nematic Liquid Crystals
We present an analysis and numerical study of an optimal control problem for the Landau-de Gennes (LdG) model of nematic liquid crystals (LCs), which is a crucial component in modern technology. They exhibit long range orientational order in their nematic phase, which is represented by a tensor-valued (spatial) order parameter $Q = Q(x)$. Equilibrium LC states correspond to $$Q$$ functions that (locally) minimize an LdG energy functional. Thus, we consider an $L^2$-gradient flow of the LdG energy that allows for finding local minimizers and leads to a semi-linear parabolic PDE, for which we develop an optimal control framework. We then derive several a priori estimates for the forward problem, including continuity in space-time, that allow us to prove existence of optimal boundary and external ``force'' controls and to derive optimality conditions through the use of an adjoint equation. Next, we present a simple finite element scheme for the LdG model and a straightforward optimization algorithm. We illustrate optimization of LC states through numerical experiments in two and three dimensions that seek to place LC defects (where $Q(x) = 0$) in desired locations, which is desirable in applications.
more »
« less
- Award ID(s):
- 2111474
- PAR ID:
- 10524442
- Publisher / Repository:
- SIAM
- Date Published:
- Journal Name:
- SIAM Journal on Control and Optimization
- Volume:
- 61
- Issue:
- 4
- ISSN:
- 0363-0129
- Page Range / eLocation ID:
- 2546 to 2570
- Subject(s) / Keyword(s):
- nematic liquid crystals, defects, finite element method, adjoint equation
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This paper discusses modeling and numerical issues in the simulation of the Landau--de Gennes (LdG) model of nematic liquid crystals (LCs) with cholesteric effects. We propose a fully-implicit, (weighted) $L^2$ gradient flow for computing energy minimizers of the LdG model, and note a time-step restriction for the flow to be energy decreasing. Furthermore, we give a mesh size restriction, for finite element discretizations, that is critical to avoid spurious numerical artifacts in discrete minimizers, particularly when simulating cholesteric LCs that exhibit ``twist.'' Furthermore, we perform a computational exploration of the model and present several numerical simulations in 3-D, on both slab geometries and spherical shells, with our finite element method. The simulations are consistent with experiments, illustrate the richness of the cholesteric model, and demonstrate the importance of the mesh size restriction.more » « less
-
The paper is devoted to the study of a new class of optimal control problems governed by discontinuous constrained differential inclusions of the sweeping type involving the duration of the dynamic process into optimization. We develop a novel version of the method of discrete approximations of its own qualitative and numerical values with establishing its well-posedness and strong convergence to optimal solutions of the controlled sweeping process. Using advanced tools of first-order and second-order variational analysis and generalized differentiation allows us to derive new necessary conditions for optimal solutions of the discrete-time problems and then, by passing to the limit in the discretization procedure, for designated local minimizers in the original problem of sweeping optimal control. The obtained results are illustrated by a numerical examplemore » « less
-
null (Ed.)In this paper, we propose a local discontinuous Galerkin (LDG) method for nonlinear and possibly degenerate parabolic stochastic partial differential equations, which is a high-order numerical scheme. It extends the discontinuous Galerkin (DG) method for purely hyperbolic equations to parabolic equations and shares with the DG method its advantage and flexibility. We prove the L 2 -stability of the numerical scheme for fully nonlinear equations. Optimal error estimates ( O ( h (k+1) )) for smooth solutions of semi-linear stochastic equations is shown if polynomials of degree k are used. We use an explicit derivative-free order 1.5 time discretization scheme to solve the matrix-valued stochastic ordinary differential equations derived from the spatial discretization. Numerical examples are given to display the performance of the LDG method.more » « less
-
We introduce a characterization of disclination lines in three dimensional nematic liquid crystals as a tensor quantity related to the so called rotation vector around the line. This quantity is expressed in terms of the nematic tensor order parameter Q , and shown to decompose as a dyad involving the tangent vector to the disclination line and the rotation vector. Further, we derive a kinematic law for the velocity of disclination lines by connecting this tensor to a topological charge density as in the Halperin-Mazenko description of defects in vector models. Using this framework, analytical predictions for the velocity of interacting line disclinations and of self-annihilating disclination loops are given and confirmed through numerical computation.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1137/22M1506158
