Abstract
In this paper, we solve a constrained Willmore problem coupled with an electrical field using IsoGeometric Analysis (IGA) to simulate the morphological evolution of vesicles subjected to static electrical fields. The model consist of two phases, the lipid bilayer and the electrolyte. The two-phases problem is modeled using the phase-field method, a subclass of the diffusive interface models. The bending, flexoelectric and dielectric energies of the model are reformulated using the phase-field parameter. A modified Augmented-Lagrangian approach was used to satisfy the constraints while maintaining numerical stability and a relatively large time step. This approach guarantees the satisfaction of the constraints at each time step over the entire temporal domain. The results show the superiority of the isogeometric analysis in solving high-order differential operators without the need for additional intermediate equations to account for classical mesh-based methods limited continuity. On the physical side, the morphological evolution of the vesicles can be simulated accurately using IGA, even when considering the flexoelectric response of the biomembrane, which adds another layer of numerical complexity to the system. The effect of the flexoelectricity, the conductivity ratio and other aspects of the problem are studied through several 3D numerical examples.
•Isogeometric analysis has been used to solve a fourth-order phase-field PDE coupled with an electric potential second-order PDE, in 3D settings.•The free energy of the system is reformulated in terms of the phase-field variable which includes bending, flexoelectric and dielectric energies. Constraints are imposed to maintain constant surface area and volume.•A modified Augmented Lagrangian method is used to satisfy the constraints at each time step while maintaining stability and accuracy.•Vesicle temporal evolution was studied through a Willmore problem, an electromechanical problem without the flexoelectric effect, and with the flexoelectric effect. The conductivity effect was also studied. All examples were based on a staggered scheme.•The examples solved in 3D settings reflect the robustness of IGA as a numerical tool to solve high-order PDEs.