We have operators already written to solve in the context of embedded boundaries
Poisson equation
Helmholtz equation
Conductivity (variable coefficient Helmholtz)
Viscous Tensors.