[cryabroad]

Recently I've been looking into the k-epsilon model in OpenFOAM, and quite confused based on what is shown here https://www.openfoam.com/documentati...k-epsilon.html.

In the above link, it says the standard k-epsilon model solves the transport equations of k and epsilon as shown in the StandardkEpsilon.png figure, while OpenFOAM solves the ones shown in the OpenFOAMkEpsilon.png figure. I noticed that the standard one is also used by ANSYS Fluent.

For constant density flows without buoyancy forces, Gb in the k equation from the standard k-epsilon model vanishes, and div(u) = 0. The two forms of the k-epsilon model are exactly the same. However, for variable density flows without buoyancy forces, although Gb goes away, div(u) is in general not equal to 0. Does this suggest that Fluent will give a wrong solution but OpenFOAM will produce a right one because it explicitly takes into account effects of the variable density?

Thanks in advance.