[cfdprogrammer]

Dear users,

I am developing a simple 2D algebraic mesh generator initially based on the TFI technique.

Given a convex quadrilateral of horizontal sides AB and CD, and vertical sides BC and DA, I have a strange result when the number of nodes on the AB is different with respect to that on BC.

If you look at the ref. images (1) and (2) in the link below, you will notice in (2) how the right-bottom corner moves down from its original position, giving an unwanted crooked lower edge when the number of subdivision on BC is 10 against the 20 on AB.

The result of Fig. (1) is correct and is obtained when AB and BC have equal number of elements.

Am I not understanding fully the concept of TFI or is there actually a bug in my code?

IMAGES HERE

Thank you very much in advance,
All the best
cfd