# Coil Length (m) = 3.000000e-02
# Projectile length (m) = 3.000000e-02
# Effective gap length (m) = 2.500000e-03
# Effective core radius (m) = 4.000000e-03
# Firing capacitance (F) = 1.000000e-03
# Number of turns in coil = 3.2e+02
# Low-field core relative permeability = 5.000000e+03
# Saturation flux density of core (Wb/m^2) = 1.000000e+00
# Saturation magnetisation of core (A/m) = 7.957747e+05
# Coil loss resistance (Ohms) = 6.570000e-01
# Mass of armature (kg) = 6.000000e-03
# Simple RLC (SCR-like) drive.
# Initial capacitor voltage = 3.000000e+02
# Initial coil current (A) = 0.000000e+00
# Initial armature speed (m/sec) = 0.000000e+00
# Initial armature position (m) = 0.000000e+00
# Time step (sec) = 1.000000e-06
# Print out every 10 time steps.
# Tolerance for the non-linear solution (Newton-Raphson) = 1.000000e-08

# Maximum force on projectile (N) = 1.124648e+02
# Minimum force on projectile (N) = -4.859833e-01
# Average force on projectile (N) = 2.552838e+01
# Maximum projectile speed (m/s) = 1.758019e+01
# Projectile exit speed (m/s) = 1.751668e+01
# Minimum projectile speed (m/s) = 5.775584e-07
# Maximum current magnitude (A)= 3.026064e+02
# Maximum capacitor voltage (V)= 2.999992e+02
# Minimum capacitor voltage (V)= -7.553270e+00
# Maximum core flux density magnitude (T)= 3.504001e+00
# Time at which pulse is switched off (sec) = 3.349520e-317
# CPU seconds used in computation. = 0.110000

• 
  1. With open coils, you can play games with the gap length to simulate the effects of the open coil. This is a bit of a fudge, though. This model specifically counts on the confinement of the field by the shield.
  
  2. The model breaks down for thick, multilayer windings (like yours). Errors arise because the radially directed magnetic field becomes significant along the length of the core. EDIT: Furthermore: many of the flux lines do not couple into the armature /EDIT. The model assumes current is concentrated in a shell surrounding the core. The closer we come to this ideal, the better it works.
  
  3. The specific solution will, of course, depend strongly on the magnetisation model. If you used a different saturation model, the results will certainly be different.
  
  4. It's bloody difficult to get exactly repeatable experimental results, since even slight misplacements of the armature can vary results considerably!
  

Bravo! Glad to hear that you got my code to work on your computer..

Great work, Mero.. Sorry about the delay...

Cheers.