Twisting a Surface, or: a Doughnut in Outer Space

Do you know this delicious German pastry?

(Photo by Guido Draheim).

It’s called “ Spritzkuchen” or “Spritzring”. I thought of it when I was toying with TikZ and pgfplots. It’s similar to a twisted torus. Today I dealt with it again, when writing on the French TeX question and answer site TeXnique.fr, with the help of my French TeX friends Patrick Bideault and Denis Bitouzé.

Before I come to code and explanation, here are some of the steps and images in the construction process, using pgfplots.

Let’s get a representative cross section. In polar coordinates, the sine function sin(x) gives a circle, sin(3x) are three leaves, we add a bit radius (1.25) as piece in the middle. This gives us the function sin(3x) + 1.25 in polar coordinates:

We embed it in the space, such as in the xy plane with z=0 as (x,y,z)(t) = ( cos(t)(sin(3t)+1.25), sin(t)(sin(3t) + 1.25), 0 ):

Or differently turned:

We can move it straight through the space for example by drawing in the xz plane and let y run linearly: (x,y,z)(t) = ( cos(t)(sin(3t)+1.25), t, sin(t)(sin(3t)+1.25) ). In space it looks this way:

But we want to turn it. Moving it in a circle in space, we use a torus map, such as:

x(t,s) = (2+cos(t))cos(s+pi/2)
y(t,s) = (2+cos(t))sin(s+pi/2)
z(t,s) = sin(t)

We connect it with our original function:

x(t,s) = (6+(sin(3t)+1.25)cos(t))cos(s)
y(t,s) = (6+(sin(3t)+1.25)cos(t))sin(s)
z(t,s) = (sin(3t)+1.25)sin(t)

Here is a half section, moved half a turn:

And fully rotated:

And now to the twist. We can twist it by adding a multiple oft t resp. y to the argument of the function. So we achieve a rotation with growing y value:

Not that we did the math and the TeX, let’s frame it. Music: from the Soundtrack of “Independence Day”, as today is the 4th of July and I hang around in the US, in Miami, Florida

 

Ok, code or it didn’t happen.

All the code and the explanation can be read here – it developed over time:

• German: Drehtransformation mit pgfplots, original post on TeXwelt.de
• English: Rotation transformation of a parametrized plot, with answer by cmhughes on TeX.SE
• French: Représenter un vissage à l’aide de pgfplots, final version on TeXnique.fr