A discussion on “Numerical computations of resonant sloshing using the modified isoAdvector method and the buoyancy-modified turbulence closure model” [Appl. Ocean Res. (2019), 93, article no. 101829, doi:10.1016.j.apor.2019.05.014]

First and foremost, the discussors wish to compliment the authors, Li et al. [6], on their impressive results involving computational fluid dynamics (CFD) simulation of resonant sloshing. We have written this discussion merely to clarify what we feel are some potential mis-characterizations of our recent work [5] found within the paper. While we recognize that these have not been the main emphasis of the discussed paper, we raise the issues below with the simple hope of preventing their further propagation in the literature. In [6] it is stated that “Larsen and Fuhrman [5] proved that classical two-equation RANS closure models are unstable in the two-phase flow. Those closure models were originally developed for single-phase flow, and those turbulence models can be applied to the two-phase flow which is treated as a single continuum mixture. However, those turbulence models can lead to the overestimation of the turbulence level in the transition region at the [air-water] interface.” We fear that the statements above mis-characterizes the nature of our work regarding two important issues, as detailed below.