This section summarizes all turbulence models and the numerical implementation of these models, which feature in the current ISNaS package. For this, several issues are considered. Our goal is to develop an efficient and unconditionally positive scheme for turbulence equations, so that the non-negativity of the turbulence quantities (e.g., k and ) is guaranteed.