Finite-difference solution of an elliptic partial differential equation with discontinuous coefficients

An algoritm is presented for point relaxation of an elliptic partial diferential equation at a grid station where its coefficient are discontinuous. Such equations can arise from the use of skewed co-ordinate system based on discotinuous reference functions (lengths) to map a complicated physical geometry onto a simple computational domain. Efficary of the scheme is illustrated by examples involving potential flow in a sharply bent channel and in a conical diffuser.