Abstract: Fixed (frozen) particle configurations are used to study CFD-DEM convergence behavior in the limit of infinite grid resolution. A filtering scheme based on solving a diffusion equation, i.e., a Gaussian kernel, is used to interpolate particle information to the continuum grid, free of cell size limitations. A regression method is used for Richardson extrapolation of the discrete solutions to the grid-free solution. The regression method is shown to predict accurate extrapolation curves in the presence of the observed oscillatory convergence behavior. The extrapolation curve is used to complete the solution verification process by quantifying the numerical uncertainty present in the calculations. Several numerical and model parameters are studied including the spatial discretization method, the kernel width and the mean drag law. A manageable level of discretization error is found with relative discretization errors generally <4%.