We will present a diffuse domain approach for solving PDEs in complex, involving domain. The approach extends the original complex domain into a larger and regular domain, where the original PDEs are extended and solved. The original boundary conditions are absorbed as a source term into the reformulated equations. Several numerical examples including multiphase flows, the epitaxial growth of two-dimensional material with anisotropy are shown to illustrate the approach.