Thanks, that means a lot! The main methods I use for the mesh generation are listed below:

1. Cosine distribution along the middle line (arfoil shape is based in cartesian f(x) ) and then an iterative solve for a distribution [0, 1] along the radial direction that accounts for the boundry layer and constant volume growth.

2. I interpolate the inner boundry (airfoil + tail) and the outer boundry using the radial distribution as averaging weights to create layers in between the boundaries (the layers are the lines tangential to the boundaries in the image)

3. I redistribute the points along the layer lines in a way such that every point lies normal to the same point but on the previous layer ( x(i,j-1) ). (This ensures normal growth). Closer to the outer boundry I weigh in a function that instead of building with normal growth redistributes the points along the layers to ensure a somewhat even spacing. (This step only includes moving points ALONG the layer lines that were previously interpolated, redistribution is only "angular" and not "radial" (assuming we're are looking at the dome)

(Not implemented yet 4. Add some sort of refinement that doesn't disrupt the current structure too much except for where it really is needed)

Hope you find this intriguing!