analyze-magnetic-levitation

Files (1)

• SKILL.md

SKILL.md (first lines)

# Analyze Magnetic Levitation

Determine whether a given magnetic system can achieve stable levitation, identify which physical mechanism enables or forbids it, calculate the conditions for force balance and stability, and verify that the levitation is stable against perturbations in all spatial degrees of freedom including tilting modes.

## When to Use

- Evaluating whether a proposed magnetic levitation design is physically viable
- Determining why a permanent magnet arrangement fails to levitate and identifying a workaround
- Analyzing superconducting levitation systems (Meissner effect, flux pinning, mixed-state trapping)
- Designing or troubleshooting active electromagnetic feedback levitation (maglev trains, magnetic bearings)
- Assessing diamagnetic levitation feasibility for a given material and field strength
- Understanding spin-stabilized magnetic levitation (Levitron) dynamics

## Inputs

- **Required**: Description of the levitated object (mass, geometry, magnetic moment or susceptibility)
- **Required**: Description of the field source (permanent magnets, electromagnets, superconducting coils, arrangement geometry)
- **Optional**: Operating environment (temperature, vacuum, vibration constraints)
- **Optional**: Desired levitation height or gap
- **Optional**: Stability requirements (stiffness, damping, bandwidth for active systems)

## Procedure

### Step 1: Characterize the System

Establish the complete physical description of the object and field source before any analysis:

1. **Object properties**: Record the mass m, geometry (sphere, disk, rod), magnetic moment mu (for permanent magnet objects), volume magnetic susceptibility chi_v (for paramagnetic, diamagnetic, or ferromagnetic materials), and electrical conductivity sigma (relevant for eddy current effects).
2. **Field source properties**: Describe the source configuration -- permanent magnet array (Halbach, dipole, quadrupole), electromagnet with coil parameters (turns, current, core material), or superconducting coil (critical current, critical field).
3. **Field geometry**: Determine the spatial profile of the magnetic field B(r). Identify the field gradient dB/dz along the levitation axis and the curvature d^2B/dz^2 that governs stability.
4. **Environmental constraints**: Note temperature range (cryogenic for superconductors), atmosphere (vacuum reduces damping), and vibration spectrum.

```markdown
## System Characterization
- **Object**: [mass, geometry, mu or chi_v, sigma]
- **Field source**: [type, configuration, key parameters]
- **Field profile**: [B(r) functional form or measured map]
- **Gradient**: [dB/dz at intended levitation point]
- **Environment**: [temperature, pressure, vibration]
```

**Got:** A complete specification of the object and field source sufficient to determine forces and stability without further assumptions.

**If fail:** If the magnetic susceptibility or moment is unknown, measure or estimate it from material data tables. Without this quantity, force calculations are impossible. For composite objects, compute an effective susceptibility from the volume-weighted average.

### Step 2: Apply Earnshaw's Theorem

Determine whether passive static levitation is possible for the given system:

1. **State Earnshaw's theorem**: In a region free of currents and time-varying fields, no static arrangement of charges or permanent magnets can produce a point of stable equilibrium for a paramagnetic or ferromagnetic body. Mathematically, the Laplacian of the magnetic potential energy satisfies nabla^2 U >= 0 (for paramagnetic/ferromagnetic), so U has no local minimum.
2. **Classify the object's response**: Determine whether the levitated object is paramagnetic (chi_v > 0), diamagnetic (chi_v < 0), ferromagnetic (chi_v >> 0, nonlinear), superconducting (perfect diamagnet, chi_v = -1), or a permanent magnet (fixed mu).
3. **Apply the theorem**:
   - For paramagnetic, ferromagnetic, or permanent magnet objects in a static field from permanent magnets or fixed currents: Earnshaw forbids stable levitation. At least one spatial direction will be unstable.
   - For diamagnetic objects: Earnshaw does NOT forbid levitation. nabla^2 U <= 0 allows a local energy minimum. Passive static levitation is permitted.
   - For superconductors: The Meissner effect provides perfect diamagnetism, and flux pinning can provide both levitation and lateral stability.
4. **Document the verdict**: State clearly whether the system is Earnshaw-forbidden or Earnshaw-permitted, and which material property determines the classification.

```markdown
## Earnshaw Analysis
- **Object magnetic classification**: [paramagnetic / diamagnetic / ferromagnetic / superconducting / permanent magnet]
- **Susceptibility**: chi_v = [value with units]
- **Earnshaw verdict**: [FORBIDDEN / PERMITTED]
- **Reasoning**: [which condition of the theorem applies or fails]
```

**Got:** A definitive classification of whether the proposed levitation is Earnshaw-forbidden or Earnshaw-permitted, with the specific physical reasoning documented.

**If fail:** If the object has mixed magnetic character (e.g., a ferromagnetic core with a diamagnetic shell), analyze each component separately. The overall stability depends on the net energy landscape, which may require numerical field computation.

### Step 3: Identify Circumvention Mechanism

If Earnshaw's theorem forbids passive static levitation, identify which of the four standard circumvention mechanisms applies:

1. **Diamagnetic levitation**: The levitated object itself is diamagnetic (chi_v < 0). Examples: pyrolytic graphite over NdFeB magnets, water droplets and frogs in 16 T Bitter magnets. Requires strong field gradients; the condition is (chi_v / mu_0) * B * (dB/dz) >= rho * g, where rho is density.

2. **Superconducting levitation**: The object is a type-I or type-II superconductor below T_c.
   - **Meissner levitation**: Complete flux expulsion provides a repulsive force. Stable but has limited load capacity and requires the superconductor to remain in the Meissner state (B < B_c1).
   - **Flux pinning** (type-II superconductors): Magnetic flux vortices are pinned at defect sites in the material. This provides both vertical levitation force and lateral restoring force, allowing the superconductor to be suspended below or above the magnet. The object is locked in 3D position relative to the field source.

3. **Active electromagnetic feedback**: Sensors measure the object's position, and a controller adjusts electromagnet currents to maintain the equilibrium. Examples: EMS maglev trains (Transrapid), active magnetic bearings. Requires power supply, sensors, and a control system with bandwidth exceeding the mechanical resonance frequency.

4. **Spin-stabilized levitation**: A spinning permanent magnet (Levitron) achieves a gyroscopic stabilization of the tilting mode that Earnshaw's theorem otherwise makes unstable. The spin must exceed a critical frequency omega_c for gyroscopic stiffness to overcome the magnetic torque. The object must also remain within a narrow mass window.

```markdown
## Circumvention Mechanism
- **Mechanism**: [diamagnetic / superconducting (Meissner or flux pinning) / active feedback / spin-stabilized]
- **Physical basis**: [why this mechanism evades Earnshaw's theorem]
- **Key requirements**: [material property, field strength, temperature, spin rate, or control bandwidth]
- **Limitations**: [load capacity, power consumption, cryogenics, mass window]
```

**Got:** Identification of the specific mechanism with its physical basis clearly explained, including quantitative requirements for the mechanism to function.

**If fail:** If the system does not clearly fit any of the four mechanisms, check for hybrid approaches (e.g., permanent magnets for the primary force with eddy current damping for stability, or diamagnetic stabilization of a paramagnetic system). Also consider whether the system uses electrodynamic levitation (moving conductors in a magnetic field), which is a distinct mechanism based on Lenz's law.

### Step 4: Calculate Levitation Conditions

Compute the force balance and quantitative conditions for stable levitation:

1. **Vertical force balance**: The magnetic force must equal gravity.
   - For a magnetic dipole in a field gradient: F_z = mu * (dB/dz) = m * g.
   - For a diamagnetic object: F_z = (chi_v * V / mu_0) * B * (dB/dz) = m * g.
   - For a superconductor (image method): Model the superconductor as a mirror and compute the repulsion between the magnet and its image.
   - For active feedback: F_z = k_coil * I(t), where I(t) is the feedback-controlled current.

2. **Solve for levitation height**: The force balance equation F_z(z) = m * g determines the equilibrium height z_0. For analytic field profiles, solve algebraically. For measured or numerically computed fields, solve graphically or numerically.

3. **Restoring force gradient (stiffness)**: Compute k_z = -dF_z/dz evaluated at z_0. For stable levitation, k_z > 0 (force decreases with increasing height). The natural frequency of vertical oscillation is omega_z = sqrt(k_z / m).

4. **Lateral stiffness**: Compute the restoring force gradient in the horizontal plane, k_x = -dF_x/dx. For Earnshaw-permitted systems (diamagnetic, superconducting), this should be positive. For feedback systems, it depends on the sensor-actuator geometry.

5. **Load capacity**: Determine the maximum mass that can be levitated by finding the field gradient at which the equilibrium becomes marginally stable (k_z -> 0 at the maximum displacement).

```markdown
## Levitation Conditions
- **Force balance equation**: [F_z(z) = m*g, explicit form]
- **Equilibrium height**: z_0 = [value]
- **Vertical stiffness**: k_z = [value, units N/m]
- **Vertical natural frequency**: omega_z = [value, units rad/s]

From pjt222/agent-almanac, MIT.

More from pjt222/agent-almanac