3D Hyperbolic Grid Generation - Help Requested...

Marcus Lobbia

I'm currently developing an automated grid generation program for complex 3D configurations using overset grids (it will be generating grids for my 3D Euler/NS compressible flow solver). Anyway, I'm having some difficulties getting the 3D field grids generated using hyperbolic methods. [For reference, I'm attempting to code the algorithm from "Handbook of Grid Generation" pp. 5.3-5.10, and "Enhancements of a Three-Dimensional Hyperbolic Grid Generation Scheme" by Chan & Steger, App. Math. and Comp. Vol. 51 pp. 181-205].

I'm using a block tridiagonal solver to integrate the linearized system of equations. However, I'm having problems with stability (especially on highly-clustered surface grids). Any suggestions or comments from anyone with experience in developing these types of programs is appreciated.

Also, if anyone can point me in the direction of some source code for setting up the block-tridiagonal system (or for the algorithm in general), I'd appreciate it - it might help me figure out some of the problems in my own code.

Thanks in advance...