In addition to observations, numerical simulations with the Weather Research and Forecasting (WRF) model are used in this study to further investigate the role of mesoscale processes in the formation and evolution of the NLLJ cases observed during PECAN. Previous work identified the optimal horizontal grid spacing, vertical grid spacing, and planetary boundary layer scheme for reproducing Great Plains NLLJs with the WRF model. Comparisons between the observed and WRF simulated NLLJs allow for validation of the model configuration. Validated WRF model outputs then provide more detailed information about the observed NLLJ structure and evolution, as the increased spatial and temporal coverage of the numerical simulations permits a more thorough description of the events. This work gives an overview of the selected PECAN NLLJ cases and presents an initial comparison of the WRF simulations with high-resolution observations from the CLAMPS system. Analyses suggest that the simulated NLLJ typically forms earlier than the observed NLLJ. The WRF simulations produce NLLJs that are stronger than the observed counterparts early in the event, but weaker than the observed NLLJs later on. Sudden variations in the boundary-layer winds seen in the observations on individual nights are well captured in the WRF simulations. The WRF model predicts the height of the NLLJ maximum on a case-to-case basis; sometimes the height is captured quite well while other times the height is overestimated by the model. Potential temperature fields are well represented by the WRF model, as indicated by comparisons with the observed data. Finally, turbulent bursts often appear in the simulated flow fields very similarly to how they appear in the observed data.