Draft Claim: The author has indicated this claim is in a draft state and not yet ready for review.

Claim

0

Adaptive Nut Magnetic Drag is Within Acceptable Limits

Evidence

The adaptive nut is the moving element that couples to the variable-pitch screw flights to generate axial thrust. Any coupling method that produces thrust also produces some opposing force (drag / slip loss). The claim is that this drag is small enough (i) not to cap achievable screw flight pitch (lead ratio), and (ii) not to impose unacceptable waste heat or power loss.

The VPSL's drive system resembles a magnetic worm-gear. The screw flight defines a helical spatial constraint and the adaptive nut supplies a magnetic normal force to that constraint. Axial motion arises from this geometry as opposed to the mechanism of coupling to a propagating magnetic field - a mechanism used in many other classes of electromagnetic launchers.

We consider two coupling classes:

  1. Repulsive (inductive / electrodynamic) coupling, analogous to Inductrack
  2. Attractive (magnet-to-ferromagnet) coupling, with drag dominated by eddy currents / hysteresis unless mitigated

Analysis Common to Both Cases

Let's use the coordinate system of the digital twin for labeling variable. In this coordinate system, forward acceleration occurs along the positive y-axis. Upward is defined and being the direction away from the planet's surface and is along the positive z-axis. Rightward is defined as being along the positive x-axis.

In both cases the speed of the adaptive nut as a function of time is

vy=atv_y = at

The screw segments all rotated at a constant angular velocity ω\omega in radians/sec. The average speed of the contact patch of a screw flight is

vscrewflight=ω×rv_{screwflight} = \omega \times r

For context, in the VPSL reference design for launching spacecraft to Mars, vy=11,123 m/sv_y = 11,123 \space m/s at the end of the horizontal acceleration section, ω=1,060 rad/s\omega = 1,060 \space rad/s, the mid-contact-patch radius is r=0.45 mr=0.45 \space m, and the speed at the middle of the contact patch of the screw flight is vscrewflight=477 m/sv_{screwflight} = 477 \space m/s.

The relative speed at the interface between the grapplers and the screw flights, vrelv_{rel}, is

vrel=vy2+vscrewflight2v_{rel} = \sqrt{v_y^2 + v_{screwflight}^2}

The pitch of the screw flight (launcherScrewThreadPitchAtExit ) at the end of the horizontal acceleration section is

pitch=vyvscrewflight=17.1pitch = {v_y \over {v_{screwflight}} } = 17.1

The angle of the screw flight relative to the adaptive nut's direction of travel is

α=atan(1/pitch)=atan(1/17.1)=3.35°\alpha = atan(1/pitch) = atan(1/17.1) = 3.35 \degree

To incorporate the lift/drag ratio at the coupling interface into our calculations, its helpful to consider how the worm-drive system is analogous to an ice-boat travelling upwind on a close-hauled tack. The adaptive nut is analogous to the iceboat which is confined to travelling along a fixed path by the guide way in one case and by the ice and skates in the other. The screw flights are the driving medium analogous to the wind. The adaptive nut's grapplers are analogous to the iceboat's sail. In the case of an iceboat, the apparent wind direction is the direction that the wind seems to be travelling from the perspective of someone on the iceboat. This is analogous to the direction the screw flights on one side of the adaptive nut seem to be moving from the perspective of someone on the adaptive nut. Mathematically, this vector is

vapparent=vrel(i^sinα+j^cosα)\boldsymbol{v}_{\mathrm{apparent}} = -v_{rel}\left(\hat{\boldsymbol{i}}\sin\alpha + \hat{\boldsymbol{j}}\cos\alpha\right)

Where i^\hat{\boldsymbol{i}} is the unit vector in the rightward direction and j^\hat{\boldsymbol{j}} is the unit vector in the forward direction..

The analogy is imperfect in that the VPSL system employs two counter-rotating screws, which would correspond to two simultaneous “winds” moving in different directions. To preserve clarity, let's consider the interaction with a single screw.

Now in the adaptive nut's frame of reference, the screw flights on one side are moving along the vapparent\boldsymbol{v}_{\mathrm{apparent}} vector, and the screw flights on its other side are moving along a mirror image vector. If repulsive coupling is used, then the two vectors make a V shape, with the V pointing aft. If attractive coupling is used, then the two vectors still make a V shape, but in this case the V points forward.

If the screws were constant pitch, then from the adaptive nut's perspective, this arrangement would be static and it would advance forward at a steady speed. But because the screws' pitch is variable,

Repulsive Coupling - Inductrack method

To derive the math that calculates the adaptive nut's magnetic drag from other known parameters, let's start with a simplified scenario. Image that there are two angled plates (stand-ins for the flights on the counterrotating twin screws) that are squeezing the adaptive nut in the manner of a watermelon seed being squeezed between two fingers. Imagine that the vehicle's forward velocity is constant, but that there is an opposing force, such as aerodynamic drag, that is resisting the forward motion of the vehicle.

To maintain speed, the angled plates must apply lateral "squeezing" forces, FLateralF_{Lateral}. Each lateral force is converted into a forward force, FThrustF_{Thrust}, by the geometry of the screw, so that

FThrust=FLateral/tan(θ)F_{Thrust} = F_{Lateral} / tan(\theta)

At the coupling interface, there is a repulsive force perpendicular to the plate which is

FThrust=FLiftcos(θ)FLift=FThrust/cos(θ)F_{Thrust} = F_{Lift} \cdot cos(\theta) \\ F_{Lift} = F_{Thrust} / cos(\theta)

There is also a magnetic drag force which is parallel to the plate which is related to the lift force by the magnetic levitation system's lift-over-drag ratio, LODLOD

FLift=FDragLODF_{Lift} = F_{Drag} \cdot LOD

Finally, the magnetic drag force can be broken down into its x and y components. (Note that this diagram's coordinate system differs from that of the digital twin).

We then get

FDrag,x=FDragcos(π/2+θ)FDrag,y=FDragsin(π/2+θ)F_{Drag,x} = F_{Drag} \cdot cos(\pi/2 + \theta) \\ F_{Drag,y} = F_{Drag} \cdot sin(\pi/2 + \theta)

Note that FDrag,xF_{Drag,x} and FAerodynamicF_{Aerodynamic} are negative.

To maintain a steady velocity, we need to satisfy

FAerodynamic+2FDrag,x+2FThrust=0F_{Aerodynamic} + 2F_{Drag,x} + 2F_{Thrust} = 0

... where the factors of two are introduced because there are two plates squeezing the adaptive nut.

Now it's possible to calculate FDragF_{Drag}, which value that will be needed to calculate how much power will be lost to magnetic and aerodynamic friction.

First we substitute in FDrag,x=FDragcos(π/2+θ)F_{Drag,x} = F_{Drag} \cdot cos(\pi/2 + \theta) and FThrust=FLiftcos(θ)F_{Thrust} = F_{Lift} \cdot cos(\theta) to get

FAerodynamic+2FDragcos(π/2+θ)+2FLiftcos(θ)=0F_{Aerodynamic} + 2F_{Drag} \cdot cos(\pi/2 + \theta) + 2F_{Lift} \cdot cos(\theta) = 0

Then we substitute in FLift=FDragLODF_{Lift} = F_{Drag} \cdot LOD to get

FAerodynamic+2FDragcos(π/2+θ)+2FDragLODcos(θ)=0F_{Aerodynamic} + 2F_{Drag} \cdot cos(\pi/2 + \theta) + 2F_{Drag} \cdot LOD \cdot cos(\theta) = 0

Finally we can factor and rearrange to solve for FDragF_{Drag}

FDrag(2cos(π/2+θ)+2LODcos(θ))=FAerodynamicF_{Drag} (2 cos(\pi/2 + \theta) + 2 LOD \cdot cos(\theta)) = -F_{Aerodynamic} FDrag=FAerodynamic2(cos(π/2+θ)+LODcos(θ))F_{Drag} = {-F_{Aerodynamic} \over 2(cos(\pi/2 + \theta) + LOD \cdot cos(\theta))}

We could also solve for FLateralF_{Lateral} to determine how much force needs to be applied, but that formula is not relevant to this particular claim.

To compute the actual drag force, we need to evolve this analysis from a steady speed case to the accelerating case. This simply requires adding in an additional term for the inertial force, FInertialF_{Inertial}, that the thrust also needs to overcome to achieve the desired acceleration. This force can be calculated with

FInertial=maF_{Inertial}=-ma

So the resulting formula is

FDrag=FAerodynamicFInertial2(cos(π/2+θ)+LODcos(θ))F_{Drag} = {-F_{Aerodynamic} - F_{Inertial} \over 2(cos(\pi/2 + \theta) + LOD \cdot cos(\theta))}

Lastly, we need to calculate the power loss associated with FDragF_{Drag}. The speed at the coupling interface is (vrelv_{rel}) the speed of the adaptive nut divided by sin(θ)sin(\theta), so the power loss is

PMagneticDrag=Fv=FDragvnutsin(θ)P_{MagneticDrag} = Fv = F_{Drag} \cdot {v_{nut} \over sin(\theta)}

Now that we have derived the needed formula, we must determine the numbers that we need to plug into it.

Inductrack lift/drag is a function of relative speed at the interface

For an Inductrack-style RL track model, Post & Ryutov give (their Eq. 12): LD=ωLR=2πvλLR\frac{L}{D} = \frac{\omega L}{R} = \frac{2\pi v}{\lambda}\frac{L}{R} where

  • λ\lambda = spatial wavelength of the magnetic field pattern seen by the loops (along the motion direction)
  • LL = effective inductance of a representative track circuit (self + mutual, per their model)
  • RR = effective resistance of that circuit
  • vv = relative speed along the wavelength direction

Mapping to the grappler–flight interface, use v=vrelv = v_{\mathrm{rel}}:

(LD)ind(vrel)=2πvrelλLR{\left(\frac{L}{D}\right)_{\mathrm{ind}}(v_{\mathrm{rel}}) = \frac{2\pi v_{\mathrm{rel}}}{\lambda}\frac{L}{R}}

Source: Post & Ryutov, “The Inductrack Approach to Magnetic Levitation,” Eq. (12).

  1. Relating “lift” to “thrust” in the worm-gear constraint

Let

  • NN = magnetic “lift-like” normal force between nut and flight (normal to the local flight surface)
  • TT = axial thrust delivered to the nut

Geometric resolution (ideal constraint, ignoring losses) gives: T=NsinαN=TsinαT = N\sin\alpha \quad\Rightarrow\quad N = \frac{T}{\sin\alpha}

Inductrack “drag” force (dissipative) is then:

Dm(vrel,vz)=N(vz)(L/D)ind(vrel)D_m(v_{\mathrm{rel}}, v_z) = \frac{N(v_z)}{(L/D)_{\mathrm{ind}}(v_{\mathrm{rel}})}

so

Dm(vrel,vz)=T/sinα(vz)2πvrelλLR=Tsinα(vz)λR2πL1vrelD_m(v_{\mathrm{rel}}, v_z) = \frac{T/\sin\alpha(v_z)}{\frac{2\pi v_{\mathrm{rel}}}{\lambda}\frac{L}{R}} = \frac{T}{\sin\alpha(v_z)}\frac{\lambda R}{2\pi L}\frac{1}{v_{\mathrm{rel}}}

Attractive Coupling

In the ideal limit (if the adaptive nut were infinitely long) the magnetic field would be static in the nut frame and static in the screw frame. In practice the adaptive nut's length is finite, so in the screw's frame (the frame observed from of any given position along the screws) the magnetic field varies as the adaptive nut passes by, first increasing, then staying steady, then decreasing.

Reviews

The following reviews are limited in scope to the validity of the claim made above, and do not imply that the reviewer has taken a position regarding any other claim or the overall feasibility of a concept that is supported by this claim.

No reviews yet.