Based on the nodal discontinuous Galerkin method, a new way to solve the lattice Boltzmann equation(LBE), i.e. a nodal discontinuous Galerkinlattice Boltzmann method (NDGLBM), is presented in this study. In this method, the collision process and streaming process in LBE are split into two substeps. To implement the collision process, the multirelaxation time (MRT) model in LBM is adopted. Meanwhile, the streaming process can be converted into the advection equation, which is then solved via the the nodal discontinuous Galerkin method. Here, the space is discretized on the unstructured grids, and the time discretization is performed by using the lowstorage fourthorder, fivestage RungeKutta scheme. To validate the current NDGLBM, the simulations of liddriven cavity flow, flows over a stationary circular cylinder, two rotatingstationary cylinders and a NACA0012 airfoil with high Reynolds number are carried out. By comparing the obtained numerical results with previous data in literature, the good agreement is achieved.