On strategies for enforcing interfacial constraints and evaluating jump conditions with the extended finite element method

We consider a problem stemming from recent models of phase transitions in stimulus-responsive hydrogels, wherein a sharp interface separates swelled and collapsed phases. Extended finite element methods that approximate the local solution with an enriched basis such that the mesh need not explicitly ‘fit’ the interface geometry are emphasized. Attention is focused on the weak enforcement of the normal configurational force balance and various options for evaluating the jump in the normal component of the solute flux at the interface. We show that as the reciprocal interfacial mobility vanishes, it plays the role of a penalty parameter enforcing a pure Dirichlet constraint, eventually triggering oscillations in the interfacial velocity. We also examine alternative formulations employing a Lagrange multiplier to enforce this constraint. It is shown that the most convenient choice of basis for the Lagrange multiplier results in oscillations in the multiplier field and a decrease in accuracy and rate of convergence in local error norms, suggesting a lack of stability in the discrete formulation. Under such conditions, neither the direct evaluation of the gradient of the approximation at the phase interface nor the interpretation of the Lagrange multiplier field provide a robust means to obtain the jump in the normal component of solute flux. Fortunately, the adaptation and use of local, domain-integral methodologies considerably improves the flux evaluations. Several example problems are presented to compare and contrast the various techniques and methods