a certain integral REVIEWED
DOCX · 55.0 KB
Open DOCX file
Short calculation note by Phil dated 3.4.14, arising from his two-cylinder computation of Az in the extreme skin effect limit with the King Helmholtz integral. He integrates by parts, substitutes z = e^{jθ}, and obtains a contour integral with up to three poles. He notes that Maple's residue call can evaluate it for each integer m, though the result is very complicated.
AI-written summary; may contain errors.
Extracted text (machine-read; may contain errors)
Calculation of an integral PhL 3.4.14
This integral appears in my attempt to compute Az for the two cylinder problem (in the extreme skin effect limit) using the King Helmholtz integral [ see "computation of Az for two cylinders case.doc" Section 10 ]. As written, I cannot find this integral in GR7 (which never has both sinθ and cosθ together for trig integrals), and Maple could only do it for small values of m. Here I convert it to a contour integral after doing parts integration and I show that one could in fact do the integral with some amount of energy and Maple coding! The sample Maple code for doing this is located in "Maple does contour integral.mws". I can see that the result is going to be very complicated, so this is probably not a good way to compute Az. But it is good to know that Maple can in effect do any contour integral using its "residue" call, and I just added an example to my Maple User's Guide.
The integral in question is this
I(m) = !Syntax Error, Idθ cos(mθ) ln [A + Bsinθ + Ccosθ]
where m is an integer. I can write
cos(mθ) = (1/m)∂θsin(mθ)
and do parts integration. The parts vanish since sin(mθ) vanishes and we get then
I(m) = !Syntax Error, Idθ sin(mθ) [A + Bsinθ + Ccosθ]-1 ( Bcosθ - Csinθ)
Now write everything in terms of eiθ :
I(m) = !Syntax Error, Idθ {(1/2)(ejmθ - e-jmθ)}[A + B{(1/2)(ejθ - e-jθ)}+ C{(1/2)(ejθ + e-jθ)}]-1
({ ( B{(1/2)(ejθ + e-jθ)} - C{(1/2)(ejθ - e-jθ)})
or
I(m) = (1/2) !Syntax Error, Idθ (ejmθ - e-jmθ)[2A + B(ejθ - e-jθ)+ C (ejθ + e-jθ)]-1
( ( B{ejθ + e-jθ)} - C(ejθ - e-jθ))
Now set z = ejθ so that dz = jzdθ and then we have
I(m) = -j (1/2) dz/z (zm - z-m)[2A + B(z- z-1)+ C (z + z-1)]-1 [ ( B{z+ z-1)} - C(z - z-1) ]
Multiply top and bottom by z to get
I(m) = -j (1/2) dz/z (zm - z-m)[2Az + B(z2- 1)+ C (z2 + 1)]-1 [ ( B{z2+ 1)} - C(z2 - 1) ]
This has the general form a = B-C b = B+C
I(m) = -j (1/2) dz/z (zm - z-m)[2Az + B(z2- 1)+ C (z2 + 1)]-1 [ ( B{z2+ 1)} - C(z2 - 1) ]
= -j (1/2) dz + [ same with m→-m ]
This is a completely doable contour integral showing two three possible pole locations.
Can Maple help with contour integration? Here is some text from a book on Maple
It shows how in a very simple case you can use the solve and residue calls to assist the evaluation, and it also shows a few other Maple tricks which are new to me like lhs and rhs.
I wrote some test code to do the above integral
dz
and it works just fine, though I have to do each m separately. For each m, then, I can get an analytic expression for this integral
I(m) = !Syntax Error, Idθ cos(mθ) ln [A + Bsinθ + Ccosθ]
and therefore I can get a finite sum approximation for Az(r,θ) for the two cylinders problem!