The partially-averaged Navier-Stokes (PANS) equations are used to predict the variable-density Rayleigh-Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS BHR-LEVM closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters f φ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs. cost) of the PANS model predicting the spatio-temporal development of the RT flow. Results comparable to largeeddy simulations are obtained, at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier-Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity φ of the turbulent closure, f φ . These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations' fidelity.