Modeling and Simulation of Transitional Rayleigh-Taylor Flow with Partially-Averaged Navier-Stokes Equations

TL;DR

This paper demonstrates that the partially-averaged Navier-Stokes (PANS) equations can effectively simulate transitional Rayleigh-Taylor flows, capturing complex flow features with lower computational cost than traditional methods.

Contribution

The study validates the PANS BHR-LEVM closure for modeling transitional RT flows, showing its efficiency and accuracy at reduced resolutions compared to LES and DNS.

Findings

PANS accurately predicts RT flow development.

Lower resolution PANS matches LES and DNS results.

Proper parameter selection is crucial for fidelity.

Abstract

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_\phi$ 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 large-eddy simulations and direct numerical 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 $\phi$ of the turbulent closure, $f_\phi$. 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.