Let us consider a rectangular domain with discontinuity(holes whose surface are insulated) meshed randomly with quadrilateral elements. I am trying to solve unsteady nonlinear diffusion equation in that domain. Which is the best numerical time integeration scheme/method? The results are always oscillatory for methods I am using.