skip to main content

Title: Minimizers for the Cahn--Hilliard energy functional under strong anchoring conditions
We study analytically and numerically the minimizers for the Cahn-Hilliard energy functional with a symmetric quartic double-well potential and under a strong anchoring condition(i.e., the Dirichlet condition) on the boundary of an underlying bounded domain. We show a bifurcation phenomenon determined by the boundary value and a parameter that describes the thickness of a transition layer separating two phases of an underlying system of binary mixtures. For the case that the boundary value is exactly the average of the two pure phases, if the bifurcation parameter is larger than or equal to a critical value, then the minimizer is unique and is exactly the homogeneous state. Otherwise, there are exactly two symmetric minimizers. The critical bifurcation value is inversely proportional to the first eigenvalue of the negative Laplace operator with the zero Dirichlet boundary condition. For a boundary value that is larger (or smaller) than that of the average of the two pure phases, the symmetry is broken and there is only one minimizer. We also obtain the bounds and morphological properties of the minimizers under additional assumptions on the domain.Our analysis utilizes the notion of the Nehari manifold and connects it to the eigenvalue problem for the negative Laplacian more » with the homogeneous boundary condition. We numerically minimize the functional E by solving the gradient-flow equation of E, i.e., the Allen-Cahn equation, with the designated boundary conditions, and with random initial values. We present our numerical simulations and discuss them in the context of our analytical results. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
SIAM journal on applied mathematics
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    In materials that undergo martensitic phase transformation, macroscopic loading often leads to the creation and/or rearrangement of elastic domains. This paper considers an example involving a single-crystal slab made from two martensite variants. When the slab is made to bend, the two variants form a characteristic microstructure that we like to call “twinning with variable volume fraction.” Two 1996 papers by Chopra et al. explored this example using bars made from InTl, providing considerable detail about the microstructures they observed. Here we offer an energy-minimization-based model that is motivated by their account. It uses geometrically linear elasticity, and treats the phase boundaries as sharp interfaces. For simplicity, rather than model the experimental forces and boundary conditions exactly, we consider certain Dirichlet or Neumann boundary conditions whose effect is to require bending. This leads to certain nonlinear (and nonconvex) variational problems that represent the minimization of elastic plus surface energy (and the work done by the load, in the case of a Neumann boundary condition). Our results identify how the minimum value of each variational problem scales with respect to the surface energy density. The results are established by proving upper and lower bounds that scale the same way. Themore »upper bounds are ansatz-based, providing full details about some (nearly) optimal microstructures. The lower bounds are ansatz-free, so they explain why no other arrangement of the two phases could be significantly better.

    « less
  2. The primal variational formulation of the fourth-order Cahn-Hilliard equation requires C1-continuous finite element discretizations, e.g., in the context of isogeometric analysis. In this paper, we explore the variational imposition of essential boundary conditions that arise from the thermodynamic derivation of the Cahn-Hilliard equation in primal variables. Our formulation is based on the symmetric variant of Nitsche's method, does not introduce additional degrees of freedom and is shown to be variationally consistent. In contrast to strong enforcement, the new boundary condition formulation can be naturally applied to any mapped isogeometric parametrization of any polynomial degree. In addition, it preserves full accuracy, including higher-order rates of convergence, which we illustrate for boundary-fitted discretizations of several benchmark tests in one, two and three dimensions. Unfitted Cartesian B-spline meshes constitute an effective alternative to boundary-fitted isogeometric parametrizations for constructing C1-continuous discretizations, in particular for complex geometries. We combine our variational boundary condition formulation with unfitted Cartesian B-spline meshes and the finite cell method to simulate chemical phase segregation in a composite electrode. This example, involving coupling of chemical fields with mechanical stresses on complex domains and coupling of different materials across complex interfaces, demonstrates the flexibility of variational boundary conditions in the context ofmore »higher-order unfitted isogeometric discretizations.« less
  3. Abstract Variable-order space-fractional diffusion equations provide very competitive modeling capabilities of challenging phenomena, including anomalously superdiffusive transport of solutes in heterogeneous porous media, long-range spatial interactions and other applications, as well as eliminating the nonphysical boundary layers of the solutions to their constant-order analogues.In this paper, we prove the uniqueness of determining the variable fractional order of the homogeneous Dirichlet boundary-value problem of the one-sided linear variable-order space-fractional diffusion equation with some observed values of the unknown solutions near the boundary of the spatial domain.We base on the analysis to develop a spectral-Galerkin Levenberg–Marquardt method and a finite difference Levenberg–Marquardt method to numerically invert the variable order.We carry out numerical experiments to investigate the numerical performance of these methods.
  4. In this work we present a systematic review of novel and interesting behaviour we have observed in a simplified model of a MEMS oscillator. The model is third order and nonlinear, and we expressit as a single ODE for a displacement variable. We find that a single oscillator exhibits limitcycles whose amplitude is well approximated by perturbation methods. Two coupled identicaloscillators have in-phase and out-of-phase modes as well as more complicated motions.Bothof the simple modes are stable in some regions of the parameter space while the bifurcationstructure is quite complex in other regions. This structure is symmetric; the symmetry is brokenby the introduction of detuning between the two oscillators. Numerical integration of the fullsystem is used to check all bifurcation computations. Each individual oscillator is based on a MEMS structure which moves within a laser-driven interference pattern. As the structure vibrates, it changes the interference gap, causing the quantity of absorbed light to change, producing a feedback loop between the motion and the absorbed light and resulting in a limit cycle oscillation. A simplified model of this MEMS oscillator, omitting parametric feedback and structural damping, is investigated using Lindstedt's perturbation method. Conditions are derived on the parameters of the modelmore »for a limit cycle to exist. The original model of the MEMS oscillator consists of two equations: a second order ODE which describes the physical motion of a microbeam, and a first order ODE which describes the heat conduction due to the laser. Starting with these equations, we derive a single governing ODE which is of third order and which leads to the definition of a linear operator called the MEMS operator. The addition of nonlinear terms in the model is shown to produce limit cycle behavior. The differential equations of motion of the system of two coupled oscillators are numerically integrated for varying values of the coupling parameter. It is shown that the in-phase mode loses stability as the coupling parameter is reduced below a certain value, and is replaced by two new periodic motions which are born in a pitchfork bifurcation. Then as this parameter is further reduced, the form of the bifurcating periodic motions grows more complex, with yet additional bifurcations occurring. This sequence of bifurcations leads to a situation in which the only periodic motion is a stable out-of-phase mode. The complexity of the resulting sequence of bifurcations is illustrated through a series of diagrams based on numerical integration.« less
  5. Strong coupling between geomechanical deformation and multiphase fluid flow appears in a variety of geoscience applications. A common discretization strategy for these problems is a continuous Galerkin finite element scheme for the momentum balance equation and a finite volume scheme for the mass balance equations. When applied within a fully implicit solution strategy, however, this discretization is not intrinsically stable. In the limit of small time steps or low permeabilities, spurious oscillations in the piecewise-constant pressure field, i.e., checkerboarding, may be observed. Further, eigenvalues associated with the spurious modes will control the conditioning of the matrices and can dramatically degrade the convergence rate of iterative linear solvers. Here, we propose a stabilization technique in which the mass balance equations are supplemented with stabilizing flux terms on a macroelement basis. The additional stabilization terms are dependent on a stabilization parameter. We identify an optimal value for this parameter using an analysis of the eigenvalue distribution of the macroelement Schur complement matrix. The resulting method is simple to implement and preserves the underlying sparsity pattern of the original discretization. Another appealing feature of the method is that mass is exactly conserved on macroelements, despite the addition of artificial fluxes. The efficacy ofmore »the proposed technique is demonstrated with several numerical examples.« less