non-physical production of k, k-w turbulence model

David Garrido

Hi everybody:

Summary: My question is: How to avoid a non-physical growing of k and wiggles of all variables near the wake cut in a structured C-Mesh when using a 2D NavierStokes code with k-w turbulence model?

I am implementing a k-w turbulence model for a 2D Navier Stokes code (structured mesh, finite differencing scheme Beam&Warming+ADI, C-mesh around the airfoil). Near the wake cut boundary (where I just compute the mean of upper and lower values of variables explicitly after each step) there is a great production of k, and fluctuations of rho, rho*u,...Those wiggles don't cause troubles when using Balwin-Lomax or Spalart-Allmaras models, but the two equation model seems not to work properly due to them. Please, did you find this problem before? If you know the way to avoid it, or if you could help me with anything useful, THANKS!