On the Numerical Treatment of Partial Differential Equations in the Neighborhood of Isolated Singularities with Applications