[Getdp] ... coordinates in integrals ...
Matt Koch
mattkoch at scitex.us
Sat Feb 24 02:35:58 CET 2007
Hi All,
thanks to Bernard Kubicek's groundwork (FaceConduct2D in GetDP Wiki
Examples), I have finally made some small progress in solving my thermal
interface problem. I am NOT using his exact implementation, and I was
able to make the two regions NOT be connected, but his line of thinking
helped me break a long existing deadlock.
To go farther, I have to learn about coordinates in integrals.
Basically, my interface is a small gap between two parallel disks, and
the integral involves radial coordinates from both sides of the gap,
i.e. radial coordinates of both disk surfaces facing the gap. The
function to be integrated on either surface is something as follows,
although the exact form really does not matter.
f(r1,r2,h) =(r1^2+r2^2+h^2)^2/[(r1^2+r2^2+h^2)^2 - 4*r1^2*r2^2]^1.5
Hereby, r1 is a radius on one surface, r2 is a radius on the other
surface and h is the separation of the two surfaces, i.e. the gap. The
actual integrals are as follows:
1) F1(r1,h) = Int{f(r1,r2,h)}dA2 with dA2 = 2*pi*r2*dr2
2) F2(r2,h) = Int{f(r1,r2,h)}dA1 with dA1 = 2*pi*r1*dr1
The important part is that the dependencies of F1 and F2 on r1 and r2,
respectively, must be retained, while those on r2 and r1, respectively,
are eliminated. So, one question is, how do I access r1 on one surface
and r2 on the other surface. I have read about $X (or X[]) and $XS (or
XS[]) in section 2.5 of the manual, but I do not know if the concept of
source and destination applies here, and if so, which would be source
and which would be destination.
Any clues as to how to attack these integrals would be much appreciated.
Once I have this answer, I will still need to extend it to where the
integral actually contains an unknown field variable in addition to the
above function f, but I'll try to cross that bridge when I get there.
Thanks and Regards,
Matt Koch