Flow in a circular duct: Simsons Rule

1. Oct 9, 2009

mcooper

I want to integrate log(ri)log(rj) over the unit circle where i,j=1,2,....,N and ri=distance between a point (x,y) and a singularity outside the domain. The answers will form the aij of a matrix.

The method I am using is to pattern the singuarities so that they form a circle around the domain. Therefore every singularity is a distance 2 from the centre of the domain. So does every log(ri/j)=2?

Also what are the limits of intergration? Are they -1 and 1 for dy and -sqrt(1-y2) and sqrt(1-y2) for dx?

I understand this is not specifically a Simpsons Rule question but once these confusions have been sorted I will be going on to solve the integral using Simpsons on Matlab.

Thank you

2. Oct 9, 2009

LCKurtz

Just trying to understand your question. Here's my take on it for you to clear up:

You apparently have N singularity points:

P1(x1,y1),P2(x2,y2)...PN(xN,yN)

For each aij you want to do the following integral, where P(x,y) represents a point in the unit circle and || || is the Euclidean norm?

$$a_{ij} = \int \int_C \log(||P-P_i||\cdot \log||P-P_j||\ dA$$

$$= \int \int_C \log(\sqrt{(x-x_i)^2+(y-y_i)^2})\cdot \log(\sqrt{(x-x_j)^2+(y-y_j)^2})\ dA$$

To answer your question about the limits, yes, you have them correct for a dx dy integral

3. Oct 9, 2009

mcooper

Correct. I am unsure whether I have to integrate for different points (x,y) in the unit circle (domain) or just choose one, in which case the easiest point to choose would be (0,0) and log(ri) would be 2 for every i=1,...n.

Hopefully I have managed to attach the report that I am basing my problem on. Mine is a more simple problem of the circular pipe but the calculations are similar. The question I have asked is regarding the matrix A at the top of page 3. I am trying ot find the unknowns cj.

Thanks

File size:
268.6 KB
Views:
51