User talk:Dmcdysan/sandbox

6/30/22 VERSION WITH SPECIFICS FOR A CURRENT LOOP - REPLACED W/ MAGNETIC MOMENT
The magnetic field at magnetopause at distance $$L$$ (m) for a magnetic dipole generated by a magnetic moment $$\mathbf m$$current of $$I_c$$ (A) in a loop of radius $$R_c$$ iswhere $\mathbf m=I_c \ R_c^2$ is the magnetic moment of a current loop. Substituting this into Equation $$ and solving for $$L$$ yields the following result The Force derived by a magnetic sail for a plasma environment is determined from MHD and kinetic equations is:  where $C_d$  is a coefficient of drag determined by numerical analysis,  $\rho u^2/2$  (Pa) is the dynamic wind pressure, and $$S=\pi L^2$$ (m2) is the effective blocking area of the plasma magnet sail with characteristic length $L$  (m) also known as the magnetopause radius  $R_{mp}$  (m). Note that this equation has exactly the same form as the drag equation in fluid dynamics. For a dipole magnetic field $C_d$ is a function of sail tilt angle (0 degrees perpendicular to flow, 90 degrees in line with flow) determined by numerical calculation. For a large current loop following the Biot–Savart law $C_d$ is a smaller value than that of a dipole and is also a function of tilt angle. Through analysis, numerical calculation, simulation and experimentation several conditions must be met before a magnetic sail can generate significant force. An important one is that in order to achieve significant thrust as determined by MHD equations the standoff distance $$L$$ must be significantly greater than the Gyroradius, also called the Larmor radius or cyclotron radius as follows:where $m_i$  (kg) is the ion mass, $$v_\perp$$ (m/s) is the velocity of ions perpendicular to the magnetic field, $$|q|$$ (C) is the elementary charge of the ion, $B_x$ (T) is the magnetic field strength at the point of reference $x$  and $C_{Li}$  is a constant that differs by source with $C_{Li}$  =1 and Wikipedia gyroradius and $C_{Li}$  =2  . At magnetopause let $$v_\perp=u_{pe}$$ and $$B_x=B_{mp}$$in the above equation, form the inequality   $r_g/L>a \ ,a>1$  using $L$  from Equation $$, and after rearrangement yields the following condition on magnetic moment $$\mathbf m$$ for MHD applicability:where the constant $$C_a = \frac{4 m_i^3}{\mid q \mid ^3 C_{SO} \mu_0^2} \approx 2.9\times10^-12$$ for CSO=1.

Note that this MHD applicability test depends upon the ratio of effective plasma wind velocity $$u=u_{pe}$$ and plasma density $$\rho=\rho_{pe}$$y for a specific plasma environment and use case. If this condition is not met then the effective sail blocking area $$A_e$$ can be much smaller than $$S=\pi L^2$$. In some cases $$A_e$$ can be determined from a kinetic model specific to the plasma environment. Dmcdysan (talk) 19:10, 1 July 2022 (UTC)

Solar wind example - Need to decide placement
Although the ram pressure of the solar wind is relatively small at 1-6 nPa at 1 AU, the force exerted on a magnetic sail is the product of this pressure, the effective area of the sail, which for a sail of radius 100 km is  m2, and a coefficient of drag ranging between 3.6 and 5. This corresponds to a force on the order of 2,000-8,000 Newtons, which is on the order of the acceleration of a small car. However, since no propellant is exhausted the specific impulse is effectively infinite and even with very small acceleration a very high velocity is potentially achievable. Dmcdysan (talk) 19:16, 1 July 2022 (UTC)

Equation line formatting and numbering in Wikipedia
Template MathFormula$$r_{Li}=\frac {m_i v_\perp}{|q| B_x C_{Li}}$$r_{Li}=\frac {m_i v_\perp}{|q| B_x C_{Li}

Template NumBlk with Template " $$ "Template NumBlk with Template " " Issues: Equations with "|" or "}}" create ??

Use "" for "|", "} }" for "}}"

True first parameter as a unique name for the page searchable in Source editing mode. Use second for display number, use a unique prefix to minimize renumbering.

See $$

Equation template

$$$$

See Equation 1

See Equation 2

Insert 1 row x 3 col Table, col 1 is indent, col 2 is Math, col 3 is anchor Table formatting is rather complex. Dmcdysan (talk) 20:14, 4 July 2022 (UTC)

Ashida 2014 Thesis Issues
Ashida 2014 thesis has issues stemming from Equation (2).Dmcdysan (talk) 18:34, 11 July 2022 (UTC)

In 2011[Citation] and 2014[Citation] Ashida and others documented Particle In Cell (PIC) simulation results for a kinematic model for cases where $$r_g >> L$$ where MHD is not applicable.

Their model for magnetosphere radius $$L$$ in Equation (2) for characteristic length L had a form unlike that of related papers that make it difficult to easily compare with other results.

Equation (12) of their study included the additional force of the injected plasma jet $$F_{jet}$$ comprised of momentum and static pressure of ions and electrons and defined thrust gain as $F_{MPS}/(F_{mag}+F_{jet})$, which differs from the definition of a term by the same name in other studies. [Cite Funaki 2012, 2013]. It represents the gain of MPS over that of simply adding the magnetic sail force and the plasma injection jet force.

For the values cited in the conclusion, $$F_{MPS}/F_{mag}$$ is 7.5 in the radial orientation. Note that the attack angle defined by Ashida is the angle of the magnetic moment and not the orientation of the coil as defined by Nishida and is therefore differs by 90 degrees.

DIFFERENT $$C_d$$ From Ashida 2014 with erroneous Equation (2) see Equation (8) from the force simulation results for a magnetosphere radius $$L$$ of 4300 m, a coil radius $$R_c $$ of 75 m and $$I_c $$ selected to yield a magnetic moment $\mathbf m = I_c \pi R_c^2$ corresponding to the specified value of $$L$$ in accordance with Equation $$. The coefficient of drag $$C_d$$ determined the relative thrust with an attack angle $\alpha _t=90 $  degrees $$C_d=1.8 $$ and with  $\alpha _t=0 $   degrees $$C_d=1. $$Dmcdysan (talk) 19:36, 11 July 2022 (UTC)

MPS Edits
Edits from MPS not used or cuts. Dmcdysan (talk) 19:33, 17 July 2022 (UTC) MPS edits ??? More details on a proposed demonstrator spacecraft - Paper available on request. Some info in Funaki 2015 press.

[Ueno?] Preliminary results reported in 2019 indicated a 50% increase above the theoretical magnetosphere size. [Citation]

[Ashida 2014 with Error in Eqn 2] Note that the attack angle defined by Ashida is the angle of the magnetic moment and not the orientation of the coil as defined by Nishida and is therefore differs by 90 degrees. MODIFIED IN ANOTHER 2014 PAPER TO ALIGN.

Dmcdysan (talk) 19:34, 17 July 2022 (UTC)

Plasma magnet edits
Plasma magnet text equations not used. Dmcdysan (talk) 19:35, 17 July 2022 (UTC) NOT USED

Equation (15) of [Slough06] defines the magnetic field BRMF (T) near an antenna coil of radius RA (m) in terms of the electron skin depth $$\delta$$ (m) as follows where RA is the antenna coil radius (m). Since the RMF falls off as 1/r2. undefinedNote that for the solar wind the skin depth  d is only a few meters depending upon wRMF (rad/s). — Preceding unsigned comment added by Dmcdysan (talk • contribs) 19:40, 17 July 2022 (UTC)


 * MHD.6
 * $$F_w(f_o)=\frac {C_d}{2} \, \pi \, p_w\,

\biggl(\frac {\mathbf m}{\pi} \sqrt{\frac {\mu_0}{2 C_{SO} p_w} }\biggr)^{2/f_{o} }$$
 * MHD.1
 * $$p_w=\rho u^2 =p_B=\frac{B_{mp}^2}{C_{SO} \mu _0}$$
 * $$F_w(f_o)=\frac {C_d}{2} \, \pi \, p_w\,

\biggl(\frac {\mathbf m}{\pi} \sqrt{\frac {\mu_0}{2 C_{SO} p_w} }\biggr)^{2/f_{o} }$$
 * Note that Equation (B-1) defines the standoff magnetic field  with r (kg/m3) the plasma density and u=w-vs (m/s) the apparent wind speed and leaves selection of Rmp(m) as a design point.
 * If Rmp/R0 is held constant, then wRMF would be constant as well. If Rmp is held constant, then increasing R0 would mean that wRMF should be decreased.
 * === Cuts, Equations not used ===
 * $$P_w= u \, F_w = C_d\ \rho  \frac{u^3}{2} \pi R_{mp}^2 = \frac {C_d u \pi}{2 \mu_0 C_{SO} } R_{mp}^2 B_{mp}^2$$
 * $$= \frac {Z e}{m_i} \, \frac {C_{BR}(f_o)}{(R_c)^{f_o} } $$
 * $$\omega_{RMF}>\omega_{ci} = \frac {Z e B(0)}{m_i}= \frac {Z e}{m_i} \, \frac {C_{BR}(f_o)}{(R_0)^{f_o} } $$
 * $$\frac {C_{BR} }{(R_c)^{f_o} }$$
 * $$L \approx R_{mp} \approx \biggl(\frac {B(R_0)}{B_{mp} } \biggr)^{1/f_o} = \biggl(\frac {\sqrt{\mu_0} \, \, \mathbf m}{4\,  \pi \, u \, \sqrt{C_{SO}  \rho} }  \biggr)^{1/f_o}= \biggl(\frac {\sqrt{\mu_0} \, \, \mathbf m}{4\,  \pi \, \sqrt{C_{SO}  p_w} }  \biggr)^{1/f_o}$$  — Preceding unsigned comment added by Dmcdysan (talk • contribs) 19:39, 17 July 2022 (UTC)

Magnetic Field Model
Edits to MFM. Dmcdysan (talk) 18:48, 24 July 2022 (UTC)


 * When the magnetic field source strength $$B(R_0) $$ is specified, substituting $$B_{mp} $$ from the pressure balance analysis from Equation $$ into the above and solving for $R_{mp} $ yields the following:


 * $$L \approx R_{mp} \approx R_0 \biggl(\frac {B(R_0)}{B_{mp} } \biggr)^{1/f_o} \approx \biggl(\frac {C_0 \sqrt{\mu_0} \, \, \mathbf m}{4\,  \pi \, u \, \sqrt{C_{SO}  \rho} }  \biggr)^{1/f_o}$$
 * This is the expression for $$L$$ when $$f_o=3 $$ with $$C_{SO}=1/2 $$ for Equation (4), with $$C_{SO}=2 $$ for Equation (4), and the magnetopause distance of the Earth. This equation shows directly how the falloff rate $$f_o $$ dramatically increases the effective sail area $$S=\pi R_{mp}^2 $$ for decreasing values of $f_o $ for a given field source magnetic moment $$\mathbf m$$ and $$B_{mp}$$determined from the pressure balance equation $$. Substituting this into Equation $$ yields the plasma wind force as a function of falloff rate $$f_o $$ as follows


 * where $$p_w = \rho u^2$$ is the wind pressure from Equation $$. This is the same expression as Equation (10b) when $$f_o=3 $$ and $$C_{SO}=1/2 $$. Note that force increases as falloff rate decreases. Dmcdysan (talk) 18:49, 24 July 2022 (UTC)

Cut overview text from magnetic sail
Physical principles involved include: plasma characteristics for the Solar wind, a planetary ionosphere and the interstellar medium; interaction of magnetic fields with charged particles in a plasma; analogies with the Earth's magnetopause; and performance measures; such as, force achieved, energy requirements and the mass of the magnetic sail system. A number of proposals using this concept have been developed since 1988 starting with the Magsail proposed by Andrews and Zubrin, which although analyzed by many sources as theoretically possible has the practical disadvantage of a mass on the order of 100 tonnes. Subsequent designs strove to achieve similar benefits with reduced mass by injecting plasma, use of a rotating magnetic field, and combining the technology of a magnetic sail with that of an electric sail. This article summarizes the literature regarding performance of these designs in theoretical models, numerical analyses, simulations and laboratory experiments. Some criticisms regarding these results have been rebutted but other issues remain unresolved. Trials for several of these designs have been proposed. This article concludes with a performance comparison of these magnetic sail designs with other each other as well as other technologies for the use cases of: acceleration or deceleration in a stellar plasma wind, deceleration in the interstellar medium, deceleration in a planetary ionosphere, a comparison with electric and solar sails, and a summary of advantages and disadvantages. Dmcdysan (talk) 01:33, 1 September 2022 (UTC)

A simulation of magnetic field lines around a circular current
Deleted by Constant314 from Biot-Savart Law on 9/26/2021.

Assuming a circular current is in xy-plane with center at the origin, thus the magnetic field in yz-plane will have x-component = 0 due to symmetrical cancellation from the circular current. The current runs counterclockwise. At a point, P, in yz-plane (y, z), the contribution from a segment of the circular current at R is,$$d \mathbf{B}=\mathbf{j} d J + \mathbf{k} d K$$with,$$d J=\frac{zR\sin\theta \ d\theta}{|R^2+P^2-2yR\sin\theta|^{\frac32}}$$$$d K=\frac{(-yR\sin\theta+R^2) d\theta}{|R^2+P^2-2yR\sin\theta|^{\frac32}}$$The magnetic field B can be obtained by line integration over the circular current with $\theta = [0,\ 2\pi]$. Thus,$$\mathbf{B}=J\mathbf{j}+K\mathbf{k}$$$$J=\int_C d J=\int_0^{2\pi} \frac{zR\sin\theta\ d\theta}{|R^2+P^2-2yR\sin \theta|^{\frac32}}$$$$K=\int_C d K=\int_0^{2\pi} \frac{(-yR\sin\theta+R^2) d\theta}{|R^2+P^2-2yR\sin\theta|^{\frac32}}$$:Here is a short Python function to execute the numerical calculation:

To draw complete field lines with direction markers, the following Python function is used: 

Dmcdysan (talk) 03:14, 6 September 2022 (UTC)