Nicolas Ali Libre1,2, Arezoo Emdadi2, Edward J. Kansa3,4, Mohammad Rahimian2, Mohammad Shekarchi2
CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.1, pp. 61-80, 2008, DOI:10.3970/cmes.2008.024.061
Abstract The numerical solution of partial differential equations (PDEs) with Neumann boundary conditions (BCs) resulted from strong form collocation scheme are typically much poorer in accuracy compared to those with pure Dirichlet BCs. In this paper, we show numerically that the reason of the reduced accuracy is that Neumann BC requires the approximation of the spatial derivatives at Neumann boundaries which are significantly less accurate than approximation of main function. Therefore, we utilize boundary treatment schemes that based upon increasing the accuracy of spatial derivatives at boundaries. Increased accuracy of the spatial derivative approximation can be achieved by h-refinement reducing the… More >