Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Numerical integration on a triangle

  1. May 22, 2008 #1
    Hi.

    To solve an integral equation I have to numerically integrate the scalar Green's function on a triangle having vertices [tex](x,y) = (0,0), (1,0), (1,1)[/tex]:


    [tex]I = \int_0^1\int_0^x \frac{e^{-i 2\pi R}}{4\pi R} dy dx[/tex]

    where

    [tex]R=\sqrt{(x-x')^2+(y-y')^2}[/tex]

    is the distance between a fixed point outside the triangle and the integration point inside the triangle. Does anyone have any experience in this area?

    To evaluate this integral one can use precomputed weights and abcissae for triangular domains, and approximate the integral by a product rule, i.e.

    [tex]I=\sum_{n=1}^N w_i f(x_i,y_i)[/tex]

    where [tex]f(x,y)[/tex] is the Green's function. But this will only be accurate when [tex]R[/tex] is sufficiently large. When [tex]R[/tex] is small, the integrand grows increasingly singular.

    So I could use the product rule only when [tex]R[/tex] is larger than some constant, but I wonder if this is the smartest way of doing it. I guess this is a very typical problem encountered when solving any kind of real-world integral equation, but I have not been able to find any info about this problem.
     
  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Can you help with the solution or looking for help too?
Draft saved Draft deleted



Similar Discussions: Numerical integration on a triangle
  1. Numerical integration (Replies: 1)

  2. Numerical Integration (Replies: 10)

  3. Numerical Integration (Replies: 1)

  4. Numerical integration (Replies: 2)

Loading...