Hi,

I have a problem on hand below:

Solve the potential flow in the domain:

40 40 40 40 40 .............40 38 38 36 36 . . . 0 0 0 . . 0 . 2 0 2 0 0 0 0 0 .............0

The streamline phi satisfies laplace equation. So how to program this using Gauss-Seidel iteraion since there's some zero points (a solid boundary) inside the domain? Thanks~

L

ooops, the domain doesn't show up. Actually it's a square domain with some solid boundary inside. For the inlet and the outlet, phi distributes linearly.

one easy strategy is to set the coefficient (in the matrix) corresponding to the solid points to a very large number (machine big). That's the diagonal term. This way you're solving for the "unknown" solid points and getting values that are nearly zero (machine zero). There are of course more sophisticated approaches, but try the latter first and see if it's good enough for you

Adrin