SUMMARYA key challenge while employing non-interpolatory basis functions in finite-element methods is the robust imposition of Dirichlet boundary conditions. The current work studies the weak enforcement of such conditions for B-spline basis functions, with application to both second-and fourth-order problems. This is achieved using concepts borrowed from Nitsche's method, which is a stabilized method for imposing constraints on surfaces. Conditions for the stability of the system of equations are derived for each class of problem. Stability parameters in the Nitsche weak form are then evaluated by solving a local generalized eigenvalue problem at the Dirichlet boundary. The approach is designed to work equally well when the grid used to build the splines conforms to the physical boundary of interest as well as to the more general case when it does not. Through several numerical examples, the approach is shown to yield optimal rates of convergence.

A chemo-mechanical model is used to capture the formation and evolution of microdomains on the deforming surface of giant unilamellar vesicles. The model is intended for the regime of vesicle dynamics characterized by a distinct difference in time scales between shape change and species transport. This is achieved by ensuring that shape equilibrium holds away from chemical equilibrium. Conventional descriptions are used to define the curvature and chemical contributions to the vesicle energetics. Both contributions are consistently non-dimensionalized. The phase-field framework is used to cast the coupled model in a diffuse-interface form. The resulting fourth-order nonlinear system of equations is discretized using the finite- element method with a uniform cubic spline basis, which satisfies global higher-order continuity. Two-dimensional and axisymmetric numerical examples of domain evolution coupled to vesicle shape deformation are presented. Curvature-dependent domain sorting and shape deformation dominated by line tension are also considered.