In this work, we propose a novel formulation for the solution of partial differential equations using finite element methods on unfitted meshes. The proposed relies discrete extension operator in aggregated method. This is robust with respect to location boundary/interface within cell. One can prove enhanced stability results, not only physical domain, but whole active mesh. However, constants ...