3
$\begingroup$

(I just might solve this one fairly shortly and post an answer here if no one else does. And maybe even if someone else does.)

$$ \begin{align} & R > r > 0 \\[6pt] x & = (R + r \cos v) \cos u \\[4pt] y & = (R + r \cos v) \sin u \\[4pt] z & = r \sin v \end{align} $$ The point $(u,v)$ above is in $\mathbb{R}^2/(2\pi\mathbb{Z})^2$, which is a flat torus. The point $(x,y,z)$ in on a torus embedded in $\mathbb{R}^3$. The mapping $(u,v)\mapsto (x,y,z)$ is not conformal. Putting something else in place of $\cos v$ and $\sin v$ should make it conformal (for now I prefer to leave $\cos u$ and $\sin u$ undisturbed). What functions should those be?

  • 2
    The embedded tori are not conformally equivalent between themselves, so you cannot uniformize all of them by $\mathbb R^2/(2\pi \mathbb Z)^2$.2012-05-31
  • 0
    But apparently _some_ of them are conformally equivalent to the flat torus, if I can believe WimC's answer.2012-05-31
  • 1
    @MichaelHardy: the map below is not always bijective. In fact I think the $k$ is a covering degree.2012-06-01
  • 1
    Saying "*the* flat torus" is misleading in the context of conformal equivalence. We get all tori up to conformal equivalence by taking the quotient of $\mathbb R^2$ by $\mathbb Z+\tau\, \mathbb Z$, where $\mathrm{im}\,\tau>0$. Look up "moduli space of the torus".2012-06-01
  • 0
    I am continuing to digest WimC's answer, and I will add more to my own answer, or possibly post a third answer.2012-06-10

4 Answers 4

4

If you replace $v$ by $f(v)$ where $f$ satisfies

$$ f' = \frac{R}{r} + \cos(f) $$

then your mapping becomes locally conformal. This equation has a closed form solution, which I happily leave to you (or WA) to figure out. ;-) This closed form also suggest that it is a globally conformal map on your "flat" torus only if

$$ r = \frac{R}{\sqrt{1 + k^2}} $$

for some positive integer $k$. (This because $2 \pi$ must be a period of $\cos(f)$ .)

  • 0
    I think your way of thinking about this and mine might converge soon. To be continued......2012-06-03
  • 0
    As the solution to this differential equation I'm getting $f(v) = 2\arctan\left(\sqrt{\frac{R+r}{R-r}} \tan\left(v \frac{\sqrt{R^2 - r^2}}{2r} \right) \right)$. I think I expected $f(v)$ to go from $0$ to $\pi$ as $v$ goes from $0$ to $\pi$ (so both $0$ and $\pi$ are fixed points) but that $f(\pi/2)>\pi/2$. That doesn't seem to be how this behaves.2012-06-03
  • 0
    Oh: "WA" means Wolfram Alpha. It finally registered in my brain. I found the solution from scratch by hand. It may be the only time I've had to find an antiderivative of that form for reasons other than that it was an exercise.2012-06-03
  • 0
    @MichaelHardy: Yes, Wolfram Alpha, sorry. And actually, what I fed it was $g(v) = \cos f(v)$ which satisfies $$ g' = -\sqrt{1-g^2} \left( \frac{R}{r} + g \right) $$ with $g(0) = 1$.2012-06-04
  • 0
    OK, looking at your differential equation satisfied by $g(v) = \cos f(v)$, I still come to the conclusion that $f(v) = 2\arctan\left(\sqrt{\dfrac{R+r}{R-r}} \tan\left(v \dfrac{\sqrt{R^2 - r^2}}{2r} \right) \right)$, and also that $\cos f(v) = \dfrac{(R-r)-(R+r)\tan^2(\cdots)}{(R-r)-(R+r)\tan^2(\cdots)}$, where "$\cdots$" is the same expression whose tangent was taken above, and $\sin f(v) = \dfrac{2\sqrt{R^2-r^2}\tan(\cdots)}{(R-r)+(R+r)\tan^2(\cdots)}$. So my next question is whether you are right that satisfaction of that differential equation entails that what you get is conformal.......2012-06-07
  • 0
    @MichaelHardy: For $\tfrac{R}{r} = \sqrt{k^2+1}$ I then get $$\sin(f) = \frac{k \sin(k v)}{\sqrt{k^2+1} - \cos(k v)}$$ $$\cos(f) = \frac{\sqrt{k^2+1} \cos(k v) - 1}{\sqrt{k^2 + 1} - \cos(k v)}$$2012-06-07
  • 0
    I figured out those last two identities just after posting my last comment. Then other things demanded my attention and I haven't thought about $k$ yet.2012-06-07
  • 0
    I had examined the mapping $e^{i\theta}\mapsto \dfrac{ae^{i\theta} - b}{a - be^{i\theta}}$. That leaves $\pm 1$ pointwise fixed and leaves the circle invariant. What you've got here seems reminiscent of that, but if I'm not mistaken, it's not the same thing, and the derivative of that mapping with respect to $\theta$ didn't seem to be proportional to an affine function of the real part.2012-06-08
  • 0
    @MichaelHardy I checked for $k=1$ and in fact it *is* equal to a Mobius transform that sends $i$ to $(-1+i)/\sqrt{2}$, but only on the circle (where it matters), not on the entire disc.2012-06-08
  • 0
    Suppose $f(z)=(-3z+1)/(z-3)$ and $z=e^{i\theta}$. I was wondering whether $|(d/d\theta)f(z)|$ is proporitional to an affine function of the real part of $z$, and seemed to be needed. Notice that this function leaves $\pm1$ pointwise fixed and leaves the unit circle invariant. Now $|(d/d\theta)f(z)|$ $=|(8/(z-3)^2)||(d/d\theta)e^{i\theta}|$ $=8/|z-3|^2$ $=8/((x-3)^2+y^2)$ where $x$ and $y$ are the real and imaginary parts. Using the identity $x^2+y^2=1$, this reduces to $4/(5-3x)$. This is not proportional to an affine function of $x$. Now notice that....2012-06-08
  • 0
    ....$f(i) = (-3/5)+i(4/5)$. The tangent line to the circle at that point meets the real axis at $-5/3$, and I had thought therefore that _maybe_ the thing above would be proportional to $x-(-5/3)$. When I saw $4/(-3(x+5/3))$, I suspected I had made some clumsy mistake in algebra. But it doesn't seem so. Was I mistaken in thinking I needed a function for which $|(d/d\theta)f(z)|$ is proportional to the distance from the center of the torus to the projection of $f(z)$ onto the real axis?2012-06-08
  • 0
    Suppose $x = (R + r k(v)) \cos u$, $y = (R + r \ell(v)) \sin u$, $z = r \sin v$, where $(k(v),\ell(v))$ wraps around the unit circle as $v$ goes from $0$ to $2\pi$, and everything is smooth. Would I be right in thinking that this mapping from the flat torus to the embedded torus is conformal iff $|(\partial/\partial v)(x,y,z)|=|(\partial/\partial u)(x,y,z)|$ at every point?2012-06-08
  • 0
    OK, I am indeed guilty of a clumsy mistake, although not in algebra. I hoped to reach the conclusion that certain rates of change are proportional to certain affine functions of real parts. BUT I needed the real parts of values of $f(z)$, NOT of values of $z$.2012-06-10
  • 0
    In the case where $R=5/3$ and $r=1$, your differential equation says $f'(v) = \frac 5 3 + \cos f(v)$. But I seem to be getting $\frac 4 3 f'(v) = \frac 5 3 + \cos f(v)$. The solution satisfies $\cos f(v) = (5\cos v - 3)/(5 - 3\cos v)$ and $\sin f(v) = (4\sin v)/(5\cos v- 3)$. Is it possible you missed the possibility that the coefficient of $f'(v)$ and that of $\cos f(v)$ might differ from each other?2012-06-24
  • 0
    @MichaelHardy The directional derivatives in $u$ and $v$ directions are orthogonal. Therefore it suffices that their lengths be equal. For $v$ that length is $r \cdot f'$ and for $u$ it is $R + r\cdot \cos(f)$. So my equation looks fine it seems.2012-06-25
1

This will not be a complete answer, but it's a different point of view from other stuff posted here. Recall the conventional paramatrization of the torus:

\begin{align} & R > r > 0 \\[6pt] x & = (R + r \cos v) \cos u \\[4pt] y & = (R + r \cos v) \sin u \\[4pt] z & = r \sin v \end{align}

As a mapping from the flat torus $\mathbb{R}^2/(2\pi\mathbb{Z})^2)$ to the embedded torus, this is not conformal. I was suggesting replacing it by

\begin{align} & R > r > 0 \\[6pt] x & = (R + r k(v)) \cos u \\[4pt] y & = (R + r \ell(v)) \sin u \\[4pt] z & = r \sin v \end{align}

where $v\mapsto(k(v),\ell(v))$ is some other parametrization of the unit circle.

Let's call $u$ the longitude and $v$ the latitude. Then the meridians of longitude are closer together on the inside of the embedded torus than on the outside, but the parallels of latitude are still just as far apart for equal values of $\Delta v$. That tells us the mapping is not conformal. So the idea would be to make $(x,y,z)$ move only as fast with respect to $v$ as it moves with respect to $u$. That rate would be proportional to the distance from the center of the embedded torus.

So look at the circle \begin{align} x & = R + r \cos v \\ z & = r \sin v \end{align} i.e. $(x-R)^2 + z^2 = r^2$. Draw any line through $(x,z)=(0,0)$ that intersects that circle twice: once on the "inside", i.e. closer to $(x,z)=(0,0)$ and thus closer to the center of the torus, and once on the "outside". As that line moves, so that both of those two points of intersection move, what is the ratio of their two rates of motion along the arc? Here's an easy result I proved: That ratio of rates of motion is the same as the ratio of their $x$-coordinates! I rashly conjectured this based only on the fact that it's obviously true in the case where the aforementioned line through the origin happens to be the $x$-axis in the $xz$-plane. There didn't seem to be any obvious geometric reason why it should happen, and yet one can verify that it does.

What we would like, then, is that points $v$ and corresponding points $\pi-v$ would be mapped by $$ v\mapsto \begin{cases} x = R + r k(v) \\ z = r \ell(v) \end{cases} $$ to two points on the circle such that if you draw the line through them it goes straight through the origin in the $xz$-plane.

A corollary is that the point $v=\pi/2$ would get mapped to the point on the circle at which the tangent line goes through the origin.

Although this tells us exactly the mapping $(k(v),\ell(v))\mapsto(k(\pi-v),\ell(\pi-v))$, it stops short of identifying $v\mapsto(k(v),\ell(v))$, except at four special points ($(\pi/2)\mathbb{Z}\mod 2\pi\mathbb{Z}$). Hence this is not a complete answer.

0

Another direction my thinking about this took, is one that WimC says is in some sense equivalent to his own answer. He approaches the problem by asking what differential equations should be satisfied by the mapping, and that seems sensible since the property of conformality is a differential equation---it's about how fast things are changing at each point.

The image point $(x,y,z)$ should move just as fast as $v$ changes, as it does as $u$ changes. But the rate of change of $(x,y,z)$ with respect to $u$ is proportional to the distance from $(x,y,z)$ to the central axis $x=y=0$. So look at just the $xz$-plane and look at the point moving around the circle of radius $r$ as $v$ changes. It should change fast at $x=R+r$ and slowly at $x=R-r$. Instead of $v$, think of the angle $\theta$, the "latitude". Latitudes $0$ and $\pi$ should map to latitudes $0$ and $\pi$ while points in between move closer to the point where $x=R-r$. So we have a mapping that fixes two points on the circle while leaving the circle invariant. But we've all seen functions like that before, namely $$ f(z) = \frac{az+b}{bz+a}. $$ This leaves $\pm1$ pointwise fixed and leaves the circle $|z|=1$ invariant. I looked at the case $a=3$ and $b=1$. It maps $i$ to $(-3/5)+(4/5)i$. The tangent line to the circle at that point meets the real axis at $-5/3$.

And lo and behold: the rates of change of $|f(z)|$ are proportional to the differences between the real parts of $f(z)$ and $-5/3$.

0

I was looking at this question because of the tokamak fields and the use of Poisson's equations to solve for the electric and magnetic fields inside a region. By mapping the boundaries of the torus to a rectangular grid, the classical second derivative of the potential of the charge can be found analytically and the boundary conditions become simpler. Hypothetically, one shape can be mapped to another. This is seen in the mapping of the stereographic projection, the semicircle to the rectangular plane, and it is proposed now that the torus can be mapped to a rectangular plane. If this is so, imagine a practice problem (E and M, Jackson 2.14) and subsequently revolving the two dimensional cylindrical shape to create a surface of revolution, a torus with four charged plates. Each plate has a voltage and then this creates a potential landscape. The curvature of the plate quadrupoles shields stray fields, but I imagine the stable or unstable equilibria within the Poisson equations solution (using the Lorentz force EOM) might indicate the preference of concave curved plates over convex plates. This problem seems to be confusing. In essence, what is the best way to solve for the Poisson equations fields inside a torus with hemispherical bosses located on the inner surface. These bosses create ELMs in real world tokamaks. I guess the questions I have are about super positioning Poisson's equations solutions to boss problems, cylinder-toroid quadrupole problems and similar ones to generate the actual field shape inside the tokamak.

To summarize: I create a torus, I map a torus to a grid (x,y,z), I place boundary conditions for Poisson's equation on my grid, I solve for the field in the grid, I map to the toroid. Then, I perturb the boundaries with hemispheres and superimpose the new fields onto the grid and subsequently onto the torus. This generates my field. Then I do the perturbation method to generate neutron magnetic moment contributions, and the poloidal field.

I worked on one of the Jackson problems in grad school and was going to do the toroid in the next few weeks. I have some code in matlab for what I am going to turn into the Poisson equation but I have yet to write the PDE ...

 a=5; %radius of circle  c=15; % radius of torus  [u,v]=meshgrid(0:10:360);% this generates the entire torus, but I need to generate a quadrupole like section  x=(c+a*cosd(v)).*cosd(u);   y=(c+a*cosd(v)).*sind(u);  z=a*sind(v);  % This where I plug in the second derivative 

My brain is full.

Steve