The ability of accurately predicting transitional fluid flow behavior is important to the design of many engineering systems. From a design standpoint, Reynolds Averaged NavierStokes based CFD represents a reasonable compromise between accuracy and expense. More recent research efforts for RANS based CFD have focused on models specifically developed to predict transitional behavior.
The kklomega model is a new transitional turbulence model which is described in [1]. However in 2013 it is appeared that this paper has several errors which must be corrected. In OpenFOAM these corrections have been applied with the version 2.4.0 and later.
Since there are many studies based on the 2008 version of the kklomega model by OpenFOAM (between 2.1.0 to 2.3.1) and probably by some other CFD tools which have used the old implementation, it may be interesting to see the real difference of the two implementations as resultwise.
In the previous study a simple comparison between the kkLOmega implementations via the default test case came with the OpenFOAM. This time a new comparison have been setup against an experimental case: naca0012.
NACA Technical Note 3361 [2] which covers the naca0012 airfoil on the full range of angle off attacks, has been used for the experimental data for this comparison as smooth surface, without any trip at Re 1,8e06 as the case.
Since the mean comparison is between the OpenFOAM versions adding the experimental data may give a better insight about the situation. Two OpenFOAM versions have been used for representing the two different kklomega turbulence model implementations;
A single case file has been made and all the mesh, boundary conditions, schemes etc. were identical so the only difference was on the turbulence model implemented.
The grid has been made by the Construct2D and imported into OpenFOAM. It consists of 52.001 hexahedral cells. Yplus values are min: 0.010183  max: 1.53139  average: 0.335399; so no wall functions since the kkLOmega is a lowRe model. Here is the checkMesh results:
/**\  =========  \\ / F ield OpenFOAM: The Open Source CFD Toolbox  \\ / O peration Version: 2.4.0  \\ / A nd Web: www.OpenFOAM.org  \\/ M anipulation \**/ Build : 2.4.0f0842aea0e77 Exec : checkMesh Date : Mar 15 2016 Time : 11:35:29 Host : "x" PID : 2268 Case : /x/x/naca0012/x/kklw.0/2.4/kklw9derece nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring runtime modified files using timeStampMaster allowSystemOperations : Allowing usersupplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time
Create polyMesh for time = 0
Time = 0
Mesh stats points: 104700 internal points: 0 faces: 208353 internal faces: 103653 cells: 52001 faces per cell: 6 boundary patches: 4 point zones: 0 face zones: 0 cell zones: 0
Overall number of cells of each type: hexahedra: 52001 prisms: 0 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0
Checking topology... Boundary definition OK. Cell to face addressing OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces... Patch Faces Points Surface topology kanat 349 698 ok (nonclosed singly connected) auto1 52001 52350 ok (nonclosed singly connected) auto2 52001 52350 ok (nonclosed singly connected) inletoutlet 349 698 ok (nonclosed singly connected)
Checking geometry... Overall domain bounding box (22.0904 24.6866 0) (24.7684 24.7408 1) Mesh (nonempty, nonwedge) directions (1 1 0) Mesh (nonempty) directions (1 1 0) All edges aligned with or perpendicular to nonempty directions. Boundary openness (4.81705e19 0 6.57157e15) OK. Max cell openness = 9.11061e15 OK. Max aspect ratio = 506.186 OK. Minimum face area = 3.40903e09. Maximum face area = 1.90642. Face area magnitudes OK. Min volume = 3.40903e09. Max volume = 0.827306. Total volume = 1824.66. Cell volumes OK. Mesh nonorthogonality Max: 18.9702 average: 6.15241 Nonorthogonality check OK. Face pyramids OK. Max skewness = 0.320559 OK. Coupled point location match (average 0) OK.
Mesh OK.
End
The following results have been obtained as steadystate via the simpleFoam solver. A transient comparison would also be interesting tough but lack of time...
Figure.1) The structured Ogrid used around the naca0012 geometry for the study. The details can be seen above on the checkMesh report.
Figure.2) Computed Lift coefficients. Until about 5degrees of AoA, both models gives similar results. On the other hand at 9degrees the difference becomes gradually increasing. Overall performance of the corrected kklomega implementation reasonably good for the prestall region.
Figure.3) Computed drag coefficients. Interestingly the corrected implementation performed worse than the wrong one for this case. This may be the result of relatively coarse mesh for accurate drag prediction but it seems the old implementation works reasonably well with the same grid.
Figure.4) Computed pressure distribution on the naca0012 section at 9degrees of angle of attack. There are important differences on suction side.
Figure.5) Computed wall shear stress on the naca0012 section at 9degress of angle of attack. While pressure side is identical for both kkLOmega versions there is a big difference on the suction side which possibly leads to nonaccurate drag prediction. A finer mesh must be studied.
References1. A ThreeEquation EddyViscosity Model for ReynoldsAveraged Navier–Stokes Simulations of Transitional Flow, 2008, D. Keith Walters  Davor Cokljat 2. Aerodynamic Characteristics of NACA 0012 Airfoil Section at Angles of Attack From 0 to 180  TN 3361, 1955, Chris Critzos  Harry H. Henson  Robert W. Boswinkle
