Entrapped air in porous media can significantly affect water flow but simulations of air entrapment are still challenging. We developed a pore-network model using quasi-static algorithms to simulate air entrapment during spontaneous wetting and subsequent drainage processes. The model, implemented in OpenPNM, was tailored to replicate an experiment conducted on a medium-sized unconsolidated sand sample. We started building the model with three types of relatively small networks formed by 54,000 pore bodies which we used to calibrate basic network topological parameters by fitting the model to the water retention curve and the saturated hydraulic conductivity of the sand sample. Using these parameters, along with X-ray image data (µCT), a larger network formed by over 250,000 pore bodies was introduced in the form of stacked sub-networks where topological parameters were scaled along the z-axis. We investigated the impact of two different contact angles on air entrapment. For a contact angle of 0, the model showed good agreement with the experimental data, accurately predicting the amount of entrapped air and the saturated hydraulic conductivity. On the contrary, for a contact angle of π/4, the model provided reasonable accuracy for saturated hydraulic conductivity but overestimated the amount of entrapped air. Overall, this approach demonstrated that a reasonable match between simulated and experimental data can be achieved with minimal computational costs.