Skip to main content

The Alcubierre Warp Drive Tophat Function and Open Science with Sage

I transferred yesterday's Mathematica file with the Alcubierre warp drive[2] line element and space curvature calculations to the +Sage Mathematical Software System today, (the files been added to the public repository[3]).  If you haven't used Sage before, it's a Python based software package that's similar in functionality to Mathematica.  Oh, and it' free.  I also worked a little more on understanding the theory, but frankly, I made far more progress with the software than the theory.  What follows will be a little more of the Alcubierre theory, plus, a cool Sage interactive demo of one of the Alcubierre functions[1], as well as a bit about my first experience with using Sage.

The theory is fun, but it's moving slowly.  Here's the chalk board from this morning's discussion

Alcubierre setup the derivation using something called the 3+1 formalism which means we consider space to be flat, (in this case), slices that are labelled with the advancing value of time.  The Greek letters alpha and beta above, are called the lapse, and the shift respectively.  The lapse is related to how the spacetime moves from time slice to time slice, and the shift is related to how space changes as time advances.

Understanding the Alcubierre Drive Equations: the shift and the tophat function
When you're cruising across the sky at faster than the speed of light, you don't want to drag too much other stuff with you in your space warp.  It's not energy efficient and it just looks silly!

The shift mentioned in the section above is what describes how the ship advances through space and it is a function of the velocity of the moving spaceship.  Alcubierre wanted to derive a set of equations that would advance the spaceship through space without advancing anything outside the size of the spaceship.  To make sure that the shift moves just the ship and not too much of anything else that happens to be near the ship, Alcubierre tailored the shift function to be very localized in space.  Consequently, the shift is the velocity of the ship times a very tophat looking sort of function denoted as f.  The tophat ensures that the shift variable is only non-zero and proportional to velocity near the spaceship.  In the following equation,  r is the distance away from the center of the spaceship, and R is a parameter that should be a little bit large than the ship.  Sigma is a parameter that determines how localized the velocity of the spaceship is.  You'll see that as sigma increases, the function becomes a nice little square pulse centered on 0 and about the width of the spaceship.

Now here's the really cool thing about Sage. Instead of me just telling you about the function and how it behaves with respect to sigma to become sharper and with respect to R to give you a larger warp for your spaceship, you can actually adjust sigma and R and see for yourself.  Voila...[1].  Make sure to hit the button  labeled, 'Analyze the Alcubierre tophat function', to start the demo.

My Experience with Sage
Sage rocks!  I was able to do everything I'd done with Mathematica, only better!  I was also able to publish the demo above with only a little fuss and muss.  I'll give you, I haven't learned how to properly size the demo cells yet, but I'll get there I bet.  For an example of how to publish Sage demos to the web, just check out the source of this web page.  Also, if you want to just play with Sage, please notice that you can type any Sage you like into the code box above.  If things become hopelessly intractable, just refresh the page to get back to the original starting point.

Oh, and if you asked if it can make the cool warp graph, well, of course it can!!!  Jr.'s into the color pink this month,so here's you're Alcubierre space warp in pink courtesy of Sage!

The Sage Install Experience
If you're using Windows and you want to install Sage to run locally, you're gong to swear a little bit at the gods o' Linux.  They've kind of obviously written the documentation pages.  However, like me, after you calm down a bit, you'll realize there are only really two steps.  First, download the gigantic Sage binary from .  It's an ova file, (if you don't know what that is don't worry yet).  Next, download VirtualBox and just take all the installation defaults.  Select the box that says something like 'Run VirtualBox' as the last installation step.  A nice shiny,new GUI will appear asking you which machine you want to run, (this won't make any sense at fit either necessarily).  Now, click on the .ova download in your browser.  If you're as lucky as Iwas, it will come up in the VirtualBox GUI and just start Sage in your brand-spanking new Linux virtual machine.  For a detailed explanation see: .

1.  permalink to blog with the cool demo! (You're here, this is for folks reading outside the blog)

2.  Alcubierre's original open access paper

3.  Alcubierre Drive Open Science Repository on GitHub


Anonymous said…
Love it, this short but sweet post helped me see why alcubierre chose a top hat function as he did.

Popular posts from this blog

Cool Math Tricks: Deriving the Divergence, (Del or Nabla) into New (Cylindrical) Coordinate Systems

Now available as a Kindle ebook for 99 cents ! Get a spiffy ebook, and fund more physics The following is a pretty lengthy procedure, but converting the divergence, (nabla, del) operator between coordinate systems comes up pretty often. While there are tables for converting between common coordinate systems , there seem to be fewer explanations of the procedure for deriving the conversion, so here goes! What do we actually want? To convert the Cartesian nabla to the nabla for another coordinate system, say… cylindrical coordinates. What we’ll need: 1. The Cartesian Nabla: 2. A set of equations relating the Cartesian coordinates to cylindrical coordinates: 3. A set of equations relating the Cartesian basis vectors to the basis vectors of the new coordinate system: How to do it: Use the chain rule for differentiation to convert the derivatives with respect to the Cartesian variables to derivatives with respect to the cylindrical variables. The chain

Superconductors and Friction

I spent most of yesterday travelling between Austin and College Station, but I did manage to get a little bit of research in.  I'm looking for other experiments related to Hirsch's theories we can do in conjunction to the search for H-rays.  One of the measurable predictions of Hirsch's theory is a change in the coefficient of sliding friction when a material enters it's superconducting state. All about friction [6] New to the lab book?  Scroll to the bottom for background and a summary of the experiment. Hirsch mentions [1] that a finding of reduced friction in superconductors might be evidence in support of his hole theory of superconductivity. Furthermore, the electronic layer outside the surface is likely to affect the friction properties of the superconductor, by providing a ‘lubricating layer’ on top of which another material would slide. As a matter of fact, an abrupt drop in sliding friction between a lead surface and solid nitrogen has be

Law of Cosines and the Legendre Polynomials

This is so cool!!! A few days ago I extolled the virtues of the law of cosines, taught in high schools the world over, and claimed that it turned up in all kinds of problems that you run into later in physics.  I gave one example of an electrostatics problem I was working on, but I had a nagging feeling in the back of my head that there were even cooler examples I'd forgotten about.  It turns out that there is a way cooler example of the importance of the law of cosines!  The law of cosines can be used to calculate the Legendre polynomials!!! OK, so what are the Legendre polynomials?  They turn up repeatedly in graduate physics classes.  First of all, they're used to solve electrostatics problems [4].  The most noticeable place I saw them was in quantum mechanics, where they were derived as a special case of spherical harmonics.  Spherical harmonics are used to describe the wavefunction of an electron in a hydrogen atom, and ultimately to come up with graphs of the wave fun