stability of scheme due to skewness/aspect ratio
hi,
i would like to know how the skewness or aspect ratio of a structured grid can affect the stability of a FVM scheme. For instance, near the boundary, the cells are usually "elongated" rectangles ie they've high aspect ratio. will this affect the stability if their ratio are too high? So is the best type of cells almost squarish ie non skew and aspect ratio of 1? thanks 
Re: stability of scheme due to skewness/aspect rat
Dear Zonexo,
The effect of skewness or aspect ratio is definitely an issue in FVM, especially in viscous flow problems. The skewed cells near the boundary, are from a physical consideration, viz. the variations in the normal direction are more important than the streamwise direction in case of boundary layers. However, standard viscous discretisation procedures suffer from stability problems on such grids. This stability loss is directly related to the loss in positivity of the viscous flux discretisation. It can be shown for instance that the quadratic reconstruction is less robust on finer meshes. It is true from this perspective of stability that a nonskew cell with unit AR is better, bu this would end up in lot more cells than actually desire and more computational effort. The trick is not to bypass the numerical issue, but tackle it, as the skewed cells are a direct consequence of the boundary layer theory. The need is to construct a robust, consistent viscous flux discretisation procedure, which behaves reasonably on stretched meshes. However, to expect a procedure to work on any mesh is not possible, even such procedures may fail on highly random and skewed grids. Hope this helps Regards, Ganesh 
Re: stability of scheme due to skewness/aspect rat
You can overcome the stability issue by treating the viscous terms in an implicit manner, e.g. CrankNicholson sheme (you should be able to find papers on this). As a result, the viscous timestepping criterion (VonNeumann) can be sidestepped.

Re: stability of scheme due to skewness/aspect rat
oh thanks for all your enlightenment!

All times are GMT 4. The time now is 16:51. 