pjt222/analyze-magnetic-field

Magnetfeld analysieren

Berechnen the magnetic field produced by a given current distribution by characterizing die Quelle geometry, selecting the appropriate law (Biot-Savart for arbitrary geometries, Ampere's law for high-symmetry configurations), evaluating field integrals, checking limiting cases, incorporating magnetic material effects where relevant, and visualizing das Ergebnising field-line topology.

Wann verwenden

• Computing the B-field from an arbitrary current-carrying conductor (wire loop, helix, irregular path)
• Exploiting cylindrical, planar, or toroidal symmetry to apply Ampere's law directly
• Estimating far-field behavior via the magnetic dipole approximation
• Superposing fields from multiple current sources
• Analyzing magnetic materials: linear permeability, B-H curves, hysteresis, saturation

Eingaben

• Erforderlich: Current distribution specification (geometry, current magnitude and direction)
• Erforderlich: Region of interest where das Feld wird benoetigt (observation points or volume)
• Optional: Material properties (relative permeability, B-H curve data, coercivity, remanence)
• Optional: Desired accuracy level (exact integral, multipole expansion order, numerical resolution)
• Optional: Visualization requirements (2D cross-section, 3D field lines, magnitude contour map)

Vorgehensweise

Schritt 1: Characterize Current Distribution and Geometry

Fully specify die Quelle vor selecting a method:

1. Current path: Beschreiben the geometry of every current-carrying element. For line currents, specify the path as a parametric curve r'(t). For surface currents, specify the surface current density K (A/m). For volume currents, specify J (A/m^2).
2. Coordinate system: Waehlen coordinates aligned with the dominant symmetry. Cylindrical (rho, phi, z) for wires and solenoids. Spherical (r, theta, phi) for dipoles and loops at large distances. Cartesian for planar sheets.
3. Symmetry analysis: Identifizieren translational, rotational, and reflection symmetries. A symmetry of die Quelle is a symmetry of das Feld. Dokumentieren which components of B are nonzero by symmetry and which vanish.
4. Current continuity: Sicherstellen, dass the current distribution satisfies div(J) = 0 (steady state) or div(J) = -d(rho)/dt (time-varying). Inconsistent current distributions produce unphysical fields.

## Source Characterization
- **Current type**: [line I / surface K / volume J]
- **Geometry**: [parametric description]
- **Coordinate system**: [and justification]
- **Symmetries**: [translational / rotational / reflection]
- **Nonzero B-components by symmetry**: [list]
- **Current continuity**: [verified / issue noted]

Erwartet: A complete geometric description of the current distribution with coordinate system chosen, symmetries cataloged, and current continuity verified.

Bei Fehler: If the geometry is too complex for a closed-form parametric description, discretize into short straight segments (numerical Biot-Savart). If current continuity is violated, add displacement current or return charge accumulation terms vor proceeding.

Schritt 2: Auswaehlen Appropriate Law

Waehlen die Methode that matches das Problem's symmetry and complexity:

1. Ampere's law (high symmetry): Verwenden wenn the current distribution has sufficient symmetry that B kann pulled out of the line integral. Applicable cases:

   • Infinite straight wire (cylindrical symmetry) -> circular Amperian loop
   • Infinite solenoid (translational + rotational) -> rectangular Amperian loop
   • Toroid (rotational about the ring axis) -> circular Amperian loop
   • Infinite planar current sheet (translational in two directions) -> rectangular loop

2. Biot-Savart law (general): Use for arbitrary geometries where Ampere's law cannot simplify:

   • dB = (mu_0 / 4 pi) * (I dl' x r_hat) / r^2
   • For volume currents: B(r) = (mu_0 / 4 pi) * integral of (J(r') x r_hat) / r^2 dV'

3. Magnetic dipole approximation (far field): Verwenden wenn the observation point is far from die Quelle (r >> source dimension d):

   • Berechnen the magnetic dipole moment: m = I * A * n_hat (for a planar loop of area A)
   • B_dipole(r) = (mu_0 / 4 pi) * [3(m . r_hat) r_hat - m] / r^3
   • Valid when r/d > 5 for ~1% accuracy

4. Superposition: For multiple sources, compute B from each independently and sum vectorially. Linearity of Maxwell's equations guarantees this is exact.

## Method Selection
- **Primary method**: [Ampere / Biot-Savart / dipole]
- **Justification**: [symmetry argument or distance criterion]
- **Expected complexity**: [closed-form / single integral / numerical]
- **Fallback method**: [if primary fails or for cross-validation]

Erwartet: A justified choice of method with a clear statement of why the chosen law is appropriate for das Problem's symmetry level.

Bei Fehler: If Ampere's law is chosen but the symmetry is insufficient (B cannot be extracted from the integral), fall back to Biot-Savart. If die Quelle geometry is too complex for analytic Biot-Savart, discretize numerically.

Schritt 3: Set Up and Bewerten Field Integrals

Ausfuehren the calculation using die Methode selected in Step 2:

1. Ampere's law path: Fuer jede Amperian loop:

   • Parameterize the loop path and compute the line integral of B . dl
   • Berechnen the enclosed current I_enc by counting all currents threading the loop
   • Solve: contour_integral(B . dl) = mu_0 * I_enc
   • Extrahieren B from the integral using the symmetry established in Step 1

2. Biot-Savart integration: Fuer jede field point r:

   • Parameterize die Quelle: dl' = (dr'/dt) dt or express J(r') over the volume
   • Berechnen the displacement vector: r - r' and its magnitude |r - r'|
   • Bewerten the cross product: dl' x (r - r') or J x (r - r')
   • Integrieren over die Quelle (line, surface, or volume)
   • For analytic evaluation: exploit symmetry to reduce dimensionality (e.g., on-axis field of a loop involves only one integral)
   • For numerical evaluation: discretize into N segments, compute the sum, and check convergence by doubling N

3. Dipole calculation:

   • Berechnen total magnetic moment: m = (1/2) integral of (r' x J) dV' for volume currents, or m = I * A * n_hat for a planar loop
   • Anwenden the dipole field formula at each observation point
   • Schaetzen der Fehler: the next multipole (quadrupole) correction scales as (d/r)^4

4. Superposition assembly: Sum contributions from all sources at each observation point. Verfolgen components separately to preserve cancellation accuracy.

## Field Calculation
- **Integral setup**: [explicit expression]
- **Evaluation method**: [analytic / numeric with N segments]
- **Result**: B(r) = [expression with units]
- **Convergence check** (if numerical): [N vs. 2N comparison]

Erwartet: An explicit expression for B(r) at the observation points, with correct units (Tesla or Gauss) and a convergence check for numerical results.

Bei Fehler: If the integral diverges, check for a missing regularization (e.g., das Feld on the wire itself diverges for an infinitely thin wire -- use finite wire radius). If numerical results oscillate with N, the integrand has a near-singularity that requires adaptive quadrature or analytical subtraction of the singular part.

Schritt 4: Check Limiting Cases

Verifizieren das Ergebnis gegen known physics vor trusting it:

1. Far-field dipole limit: At large r, any localized current distribution should produce a field that matches the magnetic dipole formula. Berechnen B from your result in the limit r -> infinity and compare with (mu_0 / 4 pi) * [3(m . r_hat) r_hat - m] / r^3.

2. Near-field infinite-wire limit: Schliessen to a long straight segment of the conductor (distance rho << length L), das Feld should approach B = mu_0 I / (2 pi rho). Check this for the relevant portion of your geometry.

3. On-axis special cases: For loops and solenoids, the on-axis field has simple closed forms:

   • Single circular loop of radius R at distance z on axis: B_z = mu_0 I R^2 / [2 (R^2 + z^2)^(3/2)]
   • Solenoid of length L, n turns per length: B_interior = mu_0 n I (for L >> R)

4. Symmetry consistency: Sicherstellen, dass components predicted to vanish by symmetry (Step 1) are indeed zero in the computed result. A nonzero forbidden component indicates an error.

5. Dimensional analysis: Sicherstellen, dass B has units of Tesla. Every term should carry mu_0 * [current] / [length] or equivalent.

## Limiting Case Verification
| Case | Condition | Expected | Computed | Match |
|------|-----------|----------|----------|-------|
| Far-field dipole | r >> d | mu_0 m / (4 pi r^3) scaling | [result] | [Yes/No] |
| Near-field wire | rho << L | mu_0 I / (2 pi rho) | [result] | [Yes/No] |
| On-axis formula | [geometry] | [known result] | [result] | [Yes/No] |
| Symmetry zeros | [component] | 0 | [result] | [Yes/No] |
| Units | -- | Tesla | [check] | [Yes/No] |

Erwartet: All limiting cases match. The field has the correct units, symmetry, and asymptotic behavior.

Bei Fehler: A failed limit indicates an error in the integral setup or evaluation. The most common causes are: wrong sign in the cross product, missing factor of 2 or pi, incorrect limits of integration, or a coordinate system mismatch zwischen source and field point parameterizations.

Schritt 5: Incorporate Magnetic Materials and Visualize

Erweitern die Analyse to include material effects and produce field visualizations:

1. Linear magnetic materials: Replace mu_0 with mu = mu_r * mu_0 inside the material. Anwenden boundary conditions at material interfaces:

   • Normal component: B1_n = B2_n (continuous)
   • Tangential component: H1_t - H2_t = K_free (surface free current)
   • In the absence of free surface currents: H1_t = H2_t

2. Nonlinear materials (B-H curves): For ferromagnetic cores:

   • Use the material's B-H curve to relate B and H at each point
   • For design purposes, approximate with piecewise linear segments: linear region (B = mu H), knee region, and saturation region (B ungefaehr constant)
   • Account for hysteresis if the operating point cycles: remanent magnetization B_r and coercive field H_c define the loop

3. Demagnetization effects: For finite-geometry magnetic materials (e.g., short rods, spheres), the internal field is reduced by the demagnetization factor N_d: H_internal = H_applied - N_d * M.

4. Field visualization:

   • Plot field lines using the stream function or by integrating dB/ds along das Feld direction
   • Plot magnitude contours (|B| as a color map)
   • For 2D cross-sections, indicate the current direction (dots for out-of-page, crosses for into-page)
   • Sicherstellen, dass field lines form closed loops (div B = 0) -- open field lines indicate a visualization or calculation error

5. Physical intuition check: Bestaetigen, dass das Feld pattern makes qualitative sense. The field sollte strongest near the current source, should circulate around currents (right-hand rule), and should decay with distance.

## Material Effects and Visualization
- **Material model**: [vacuum / linear mu_r / nonlinear B-H / hysteretic]
- **Boundary conditions applied**: [list interfaces]
- **Visualization**: [field lines / magnitude contour / both]
- **Div B = 0 check**: [field lines close / verified numerically]

Erwartet: A complete field solution einschliesslich material effects where relevant, with a visualization that shows closed field lines consistent with div B = 0 and qualitative behavior matching physical intuition.

Bei Fehler: If field lines nicht close, the calculation has a divergence error -- recheck the integral or numerical method. If the material introduces unexpected field amplification, verify that mu_r is applied only inside the material volume and that boundary conditions are korrekt enforced at every interface.

Validierung

• [ ] Current distribution is fully specified with geometry, magnitude, and direction
• [ ] Current continuity (div J = 0 for steady state) is verified
• [ ] Coordinate system is aligned with the dominant symmetry
• [ ] Method selection (Ampere / Biot-Savart / dipole) is justified by symmetry analysis
• [ ] Field integrals are set up with correct cross products and limits
• [ ] Numerical results show convergence (N vs. 2N test)
• [ ] Far-field dipole limit is verified
• [ ] Near-field and on-axis limits match known formulas
• [ ] Forbidden symmetry components are zero
• [ ] Units are Tesla durchout
• [ ] Material boundary conditions are korrekt applied (if applicable)
• [ ] Field lines form closed loops (div B = 0)

Haeufige Stolperfallen

• Wrong cross-product direction: The Biot-Savart cross product is dl' x r_hat (source to field), not r_hat x dl'. Getting this backward flips the entire field direction. Use the right-hand rule as a quick check.
• Confusing B and H: In vacuum B = mu_0 H, but inside magnetic materials B = mu H. Ampere's law in Bezug auf H uses free current only; in Bezug auf B it includes bound (magnetization) currents. Mixing conventions produces factors-of-mu_r errors.
• Applying Ampere's law ohne sufficient symmetry: Ampere's law is always true but only useful when symmetry allows B to be extracted from the integral. If B varies along the Amperian loop, the law gives a single scalar equation for a spatially varying function -- underdetermined.
• Ignoring the finite length of "infinite" wires: Real solenoids and wires have ends. The infinite-wire or infinite-solenoid formula is valid only far from the ends (distance from end >> radius). Near the ends, use the full Biot-Savart integral or finite-solenoid corrections.
• Neglecting demagnetization in finite geometries: A magnetized sphere or short rod nicht have the same internal field as a long rod in the same applied field. The demagnetization factor can reduce the effective internal field by 30-100% abhaengig von the aspect ratio.
• Non-physical field lines: If a visualization shows field lines that begin or end in free space (not on a current source or at infinity), the calculation or plotting algorithm has an error. Magnetic field lines always form closed loops.

Verwandte Skills

• solve-electromagnetic-induction -- use the computed B-field to analyze time-varying flux and induced EMF
• formulate-maxwell-equations -- generalize to the full set of Maxwell's equations einschliesslich displacement current and wave propagation
• design-electromagnetic-device -- apply magnetic field analysis to the design of electromagnets, motors, and transformers
• formulate-quantum-problem -- quantum treatment of magnetic interactions (Zeeman effect, spin-orbit coupling)