[TomWaits]

I am coding up k-omega-SST on a finite difference grid (such that u,p,omega,k are defined at the grid points and not at the cell centers). The near wall boundary condition for omega is



However y=0 for a finite difference grid at the wall! And my solution is blowing up at the boundaries if I take y to be the distance from the wall to the nearest grid point in my grid and use that omega at the wall. What is the appropriate way to treat omega in the near wall case for grids which are not cell centered?


Additionally, should I expect to be very large near the wall if my yplus is small? As y->0, omega->large. But in the log layer as y->0, k->0 so omega remains bounded.


Many thanks for any help,
Tom Waits



Many thanks for any help,
Tom Waits