This paper presents a direct numerical simulation (DNS) of incompressible turbulent channel flow at a Reynolds number of $Re_{\tau} = 5186$. The simulation, conducted using a Fourier-Galerkin method in the streamwise and spanwise directions and a B-spline collocation method in the wall-normal direction, resolves the flow with high spatial resolution. The results exhibit several characteristics of high Reynolds number wall-bounded turbulence, including a logarithmic variation in the mean velocity profile with a von Kármán constant of $\kappa = 0.384 \pm 0.004$, a logarithmic dependence of the variance of the spanwise velocity component, and a distinct separation of scales between large outer-layer structures and small inner-layer structures. The one-dimensional spectrum of the streamwise velocity fluctuation shows a $k^{-1}$ dependence over a short range in $k$, consistent with previous experimental observations. The bi-modal structure in the premultiplied spectra, with local peaks at wavenumbers on either side of the $k^{-1}$ range, is also observed. The simulation provides valuable data for studying high Reynolds number wall-bounded turbulence and is a reliable source of information on this phenomenon.This paper presents a direct numerical simulation (DNS) of incompressible turbulent channel flow at a Reynolds number of $Re_{\tau} = 5186$. The simulation, conducted using a Fourier-Galerkin method in the streamwise and spanwise directions and a B-spline collocation method in the wall-normal direction, resolves the flow with high spatial resolution. The results exhibit several characteristics of high Reynolds number wall-bounded turbulence, including a logarithmic variation in the mean velocity profile with a von Kármán constant of $\kappa = 0.384 \pm 0.004$, a logarithmic dependence of the variance of the spanwise velocity component, and a distinct separation of scales between large outer-layer structures and small inner-layer structures. The one-dimensional spectrum of the streamwise velocity fluctuation shows a $k^{-1}$ dependence over a short range in $k$, consistent with previous experimental observations. The bi-modal structure in the premultiplied spectra, with local peaks at wavenumbers on either side of the $k^{-1}$ range, is also observed. The simulation provides valuable data for studying high Reynolds number wall-bounded turbulence and is a reliable source of information on this phenomenon.