Charles Danforth
Johns Hopkins University, Baltimore, MD 21218
Magnetic sails are related to the more familiar lightsails in that they both utilize ambient energy for propulsion. While lightsails use photons, magnetic sails use charged particles. They can operate in two orientations, with rather different results. A magnetic sail is designed using current technology and modeled to see its behavior and to determine net thrust. Advantages and disadvantages are discussed and some intriguing possibilities are put forward. 
1  Introduction:
The concepts of sailing a boat on a lake and riding a rocket to Mars are rarely associated. It appears, however, that intuition does not rule in this case as it is indeed quite possible to sail to Mars or indeed anywhere else.
There are several ways to do this. The first, lightsails, are a well
known and comfortably mundane entity to sciencefiction fans and cutting edge aerospace inventors. Just like a sail on a terrestrial boat uses wind to push it in another direction, a solar sail uses light pressure from the sun or a laser to push it in the opposite direction. Typically constructed of titanic sheets of very thin metallic foil spread over thousands of square kilometers, they trap the feeble momentum of photons. Though only one small test model was build to date, these objects are relatively well understood and welltreated elsewhere. Solar sails will not be dealt with here except in comparison.
The second way to sail in outer space is the magnetic sail. Instead of light pressure, the magnetic sail uses the redirected momenta of charged particles to power its flight.
The structure of the magsail is simple and in fact could be retrofitted onto almost any spacecraft. All that is required is a loop of wire. The spacecraft simply dumps this loop of wire overboard and sets up a current in it. Lenz law says that flux will attempt to remain constant. This will cause the area of the loop to increase until it forms a circle with little or no structural help required. All that is then required is to attach the loop to the spacecraft with a few shroud lines to couple the two (Fig 1.)
Assuming a monodirectional flux of charged particles, we can look at the two extremes of magsail orientation (angle of attack). The simplest case is the symmetric or axial caseso called because the particles approach along the axis of the magnetic field and scatter symmetrically. Particles entering at any other impact parameter than directly along the axis will experience a lateral acceleration thereby imparting part of their initial zmomentum to the field and hence to the spacecraft.
The other configuration, the asymmetric or normal orientation, is a little less obvious. With the dipole moment oriented in z (loop is in the xy plane) and the particles entering with V=V_{x}, they will be bent by the field into giving a y component to the velocity. Thus a component of x and y momentum is transfered to the sail. Thinking of this as a craft in a circular orbit around the sun, we now have the equivalent of outward thrust and lateral thrust, a very useful situation indeed. Thrust can be further vectored by varying the current and orientation of the loop. Indeed very complicated maneuvers are possible as related by (Zubrin and Andrews,91).
2  Modeling a Magsail
Design parameters
Since no magsail has ever been built, nor does the need for one yet exist, we need to determine the operating parameters of our sail from scratch. At present, no work exists on doubleloop or more complicated configurations of sails, and there is no reason to believe that such a structure would be more efficient than a single loop of current.
Clearly we wish to maximize the scattering cross section while minimizing mass. The effectiveness of the sail is related to the strength of magnetic field it can produce. Calling the radius of the sail R, the effective cross section scales as R^{2} while the magnetic field strength scales as R^{2}. The mass of the sail itself will also scale with R. Thus, let's choose an intermediate value for sail area, say 10 kilometers. (Note for those unfamiliar with the field: 10km is tiny compared to most usable light sail designs.)
Field strength also scales as I, so we should attempt to get maximum current through the cable. To minimize current loss, it would be a good idea to make the loop out of some sort of superconductor. While it has certain advantages in not needing a huge generating capability on board (solar panels, radiothermal generator, batteries, etc), this brings with it the inherent difficulty of critical temperature and current beyond which superconducting properties break down.
Recently discovered hightemperature superconductors will make this easier to accomplish. Ba_{2}Cu_{3}O_{7}Y, a new ceramictype superconductor has a critical current of about 10^{10} Amps/m^{2} (Zubrin&Andrews, 91) and a critical temperature of 90K (CRC). While in interstellar space, the background temperature of ~3K poses no threat, direct solar illumination in interplanetary space can quickly raise the temperatures above critical. With a thin layer of reflective insulation (doubling as protection from debris punctures), superconductivity can be maintained. The shroud lines need be nothing fancy at all as the accellerations experienced by the craft are not large.
The density of our superconductor is presumably something on the order of 5grams/cm^{3}. If we give it a diameter of 12mmabout the diameter of a good mountaineering ropeit could carry a megampere of current and still maintain superconductivity. The insulation can be assumed to be fairly light (~0.5 grams/cm^{3}) and fairly thin (~3mm). This gives a total mass for the current loop and insulation of 36,000 kg or 36 tons.
Table 1: Design Parameters  

Modeling: Behavior
To get some idea of the behavior of a magnetic sail, a computer code was written (see Appendix) to find the trajectories of particles through the magnetic field. The particles in question are supplied by the solar wind which is composed mostly of free protons and electrons. Protons were choosen for our study because they carry about 1000 times the momentum of electrons and thus are much more efficient at transferring momentum to the spacecraft. The incident velocity was choosen to be 500 km/s, a good average velocity for particles at 1AU from the sun. In fact, the velocity varies between about 300 and 700 kilometers per second depending on solar conditions and the locations of coronal holes (Zelik and Smith).
Defining some geometry (fig 3), the vector from the origin (placed for convenience at the center of the loop) to the loop r' was written in terms of the equatorial angle phi
Table 2: Cross Section, MonteCarlo Results  
 
3  Conclusions
Clearly, the big advantage of magnetic sails is their ability to operate with "ambient" power sources and no heavy onboard propellant. The technology is relatively simple and is essentially possible today. Furthermore, it is a technology which can be retrofitted onto existing, or resonably nearfuture, spacecraft with little trouble.
One of the greatest hazards to humans in space, especially for long duration flights away from the earth's sheltering magnetosphere, is the harsh radiation and highenergy particle environment of space. While ionizing radiation is not affected by the magnetic sail, charged particles are. The sail's magnetosphere deflects any high energy charged particle from the center of the current loop. While this helps a bit, more conventional sheilding (such as tanks of water) must be used to absorb radiation and neutral particles.
One would think that the entire spacecraft would need to be shielded from the strong magnetic fields. However, due to the sail's large radius, the field strength at the center of a 10km sail is only 6x10^{11} tesla per amp of loop current. For a typical operating current of a megamp, this still produces a field of only 63 microtesla or 0.63 Gauss, comparable to the magnetic field of the Earth. This poses no hazard to humans or machinery.
While deployment of 2dimensional structures (such as light sails) is a very complicated matter involving special folding techniques, onedimensional structures (such as the loop of a magsail) can simply be wound around spools and deployed at a moment's notice. It should also be possible to vary the size of the loop ("reefing" in nautical terminology) providing another navigational degree of freedom.
There are two disadvantages inherent to magnetic sailing. The first, one which it shares with light sailing, is that the density of the momentumbearing medium falls off as r^{2}. There is a bright side to this, however. The solar wind tends to accellerate as it moves outward so, even though its density falls, the momentum per particle rises.
The other major disadvantage, one unique to magsails, is the unpredictability of the "propulsion system". Traditionally, timing is crucial when dealing with spacecraft trajectories. Orbits are calculated to the nth order years ahead of time. The magsail's thrust, and hence the mission time, depends largely on the velocity and number density of incoming protons. As has been pointed out previously in this paper, solar wind velocities and densities vary wildly depending on conditions on the surface of the sun. Fortunately, there is a way to make this navigational uncertainty much less debilitating. Since sails can vary their thrust and in effect hover, they can travel over a standard transfer orbit and then loiter at the destination until the target body arrives. Velocity matching is simply a matter of adjusting current and angle of attack. Details of these calculations do not come into this paper, but are breifly outlined in (Zubrin and Andrews 91).
While it is doubtful that magnetic sails could serve to accellerate a rocket ship at a high enough rate to make interstellar travel feasible, they could function nicely as a decellerating factor for interstellar vehicles. In this case, it is not the charged particles which enter the magnetic field but rather the magnetic field sweeping through a volume of particlefilled space. Interstellar space is filled with gas, some of it ionized and some not. However, the shock of crossing the magnetospheric boundary of the magsail (what is referred to as a Bow Shock in hydro and magnetodynamics) would serve to ionize the medium very effectively. Furthermore, since the relative velocity of the medium is so much higher than is encountered in interplanetary space, the efficiency of the magsail as brake rises. One can easily picture interstellar probes unfurling magnetic wings as they brake to investigate distant stars...
4  References
5  Technical Appendix
The core Fortran code is presented below. It follows the path of one particle through the magnetic field given a set of initial coordinates and velocities, calculates the magnetic field, and finds the total deltav of the particle. A data file is written giving the timesteps and the (x,y,z) coordinates and velocities of the proton at each step. It can be easily modified for many particles or for a MonteCarlo simulation. The reader is assumed to have a working knowledge of basic Fortran77.
program magtracer c by Charles Danforth for JHU 171304 c 42497 last modified 42997 c Variable declarations real b(3) (magnetic field real p,dp (integration angle, interval integer t,tmax,dt (time, max, interval real rr2,g,current,const (radius squared, loop rad real r(3,10000) (particle coordinate(time) real v(3,10000) (particle velocity(time) real a(3,10000) (particle accelleration(time) real radius (scalar dist from origin real dv(3),deltav (deltaV vector and scalar real Pi,mu0,q,m (physical constants (SI units) c Constant Definitions Pi=3.1415926535 mu0=4*Pi*.0000001 q=1.6E19 m=1.672E27 c Open data file open(1,file='magtracer.output') c &&& multiparticle loop start &&& c  INITIAL CONDITIONS  c positions in meters (1,2,3)(x,y,z) r(1,1)=50000 r(2,1)=(yinit)*10000 r(3,1)=0 c velocities in meters/second v(1,1)=500000 v(2,1)=0 v(3,1)=0 c  PARAMETERS  c Set number of time steps TMAX in units of dt tmax=250 dt=.001 c Set dp increment (radians) dp=.01 c Set mag field scale factors current=I(amps), loop radius=g(meters) current=1E6 g=10000 const=I*mu0*dp/(8*Pi**2) c write initial condition in the data file write(1,*),1,r(1,1),r(2,1),r(3,1),v(1,1),v(2,1),v(3,1) c ### Start moving the particle ############## do t=2,tmax radius=(r(1,t1)**2+r(2,t1)**2+r(3,t1)**2)**.5 c ^^^^^^ FIND MAGNETIC FIELD at current particle location ^^^^^^^ c Code offers three options for magnetic field. The first is the c full BiotSavart law which should be universally applicable. c Second is the magnetic dipole approximation which is very similar c to the first but takes much less time to run since it involves no c numerical integration routines. c Thirdly, for comparison, a uniform field is presented. c &&&&& BIOTSAVART LAW FOR CURRENT LOOP (radius 1) &&&&& do i=1,3 b(i)=0 end do c Numerically Integrate around loop, dl X r / r^3 do p=0,6.28,dp rr2=((r(1,t1)g*cos(p))**2+(r(2,t1)g*sin(p))**2+r(3,t1)**2) b(1)=b(1)+const*r(3,t1)*g*cos(p)/rr2**1.5 b(2)=b(2)+const*g*sin(p)*r(3,t1)/rr2**1.5 b(3)=b(3)const*g*(r(2,t1)*sin(p)r(1,t1)*cos(p)+1)/rr2**1.5 end do c &&&&& DIPOLE FIELD (dipole approx) &&&&& c >> Current Loop in xy plane, dipole in z direction << c if(radius.lt.0.00001) then c b(1)=0 c b(2)=0 c b(3)=bo c else c b(1)=bo*3*r(1,t1)*r(3,t1)/(radius)**5 c b(2)=bo*3*r(2,t1)*r(3,t1)/(radius)**5 c b(3)=bo*(2*r(3,t1)**2r(1,t1)**2r(2,t1)**2)/(radius)**5 c end if c c &&&&& CONSTANT FIELD &&&&& c >> Constant field in z direction << c b(1)=0 c b(2)=0 c b(3)=bo c Find the accelleration for this step (componentwise Lorenz Force Law) a(1,t)=q/m*(v(2,t1)*b(3)v(3,t1)*b(2)) a(2,t)=q/m*(v(3,t1)*b(1)v(1,t1)*b(3)) a(3,t)=q/m*(v(1,t1)*b(2)v(2,t1)*b(1)) c Move the particle in r and v do i=1,3 v(i,t)=v(i,t1)+a(i,t)*dt r(i,t)=r(i,t1)+v(i,t)*dt end do c Print everything to a file and screen write(1,*),t,r(1,t),r(2,t),r(3,t),v(1,t),v(2,t),v(3,t) print*,t,r(1,t),r(2,t),r(3,t),v(1,t),v(2,t),v(3,t),radius end do c Calculate deltav of particle (and hence of the field) do i=1,3 dv(i)=v(i,tmax)v(i,1) end do deltav=(dv(1)**2+dv(2)**2+dv(3)**2)**.5 print*,'DELTA V=',dv(1),dv(2),dv(3),deltav c &&& End of multiparticle loop &&& close(1) end