跳转至

Infinite Energy is closer than you think


online slides version

Energy

\[ E = mc^2 \]
  • fire : \(10^{-8}\%\)

\(200\text k\) years ago

pick wood burn

smoke, CO

  • fossil fuel: \(10^{-7}\%\)

\(2 \text k\) years ago

dead body from plants, dinosaurs and even our ancestors

greenhouse gases

  • electricity: \(200\) years ago

currently a energy transmission method

  • fission: \(0.09\%\)

\(<100\) years,

heavy elements are rare in the universe (U, Pu)

environment impact (Discharge of radioactive water of the Fukushima Daiichi Nuclear Power Plant, start from this year, and continue for 30 years)

  • fusion: \(0.7\%\)

still doing research,

but sun has already doing fusion for \(5\) billion years

  • blackhole: \(33\%\)

  • antimatter: \(100\%\)

If we compress the \(13.8\) billion years into one day. The solar system formed at noon. Human appear in the last half minute. Modern Human Civilization appear in the last millisecond. --Cixin Liu

Fusion

average binding energy

Lawson criterion

\(T\): Temperature

\(n\) plasma density

\(\tau_E\) confinement time $$ P\propto Tn\tau_E $$

fusion \(Tn\tau_E(m^{-3}KeVs)\) energy(\(MeV\))
\(\ce{D + T -> n + ^4He}\) \(2.9\times 10^{21}\) \(17.589\)
\(\ce{D + ^3He -> p + ^4He}\) \(5.1\times 10^{22}\) \(18.353\)
\(\begin{aligned}&\ce{5D -> 2^4He + 2n + p}\\&\ce{n + p^{hot} -> D}\end{aligned}\) \(1.1\times 10^{23}\) \(43.2\)
\(\begin{aligned}&\ce{3D -> p + n + ^4He}\\&\ce{n + p^{hot}-> D}\end{aligned}\) \(1.3\times 10^{23}\) \(21.6\)
\(\ce{D + D -> n + ^3He}\) \(5\times10^{23}\) \(3.269\)
\(\ce{D + D -> p + T}\) \(5\times 10^{23}\) \(4.032\)
\(\ce{2T -> ^4He + 2n}\) \(2\times 10^{23}\) \(11.3\)
\(\ce{p + B_{11}-> 3^4He}\) \(4\times10^{24}\) \(8.682\)
\(\ce{p + ^6Li -> ^4He + ^3He}\) \(1.25\times10^{25}\) \(4\)

https://www.zhihu.com/question/398163463/answer/2596028611

Device

fusion reactor

Tokamak

A tokamak is a device which uses a powerful magnetic field to confine plasma in the shape of a torus

The electric current in Solenoid coils should rise linearly

tokamak

magenetic field

lamor procession

Toroidal Coils

Torus shape

Toroidal Coils

Central Solennoid Coils

linear increasing current, super conducting

toroidal coils origin

toroidal coils twist

Poloidal Coils

control the section area of the plasma

poloidal coils

Divertor

remove heavier ions, transfer fusion energy

divertor

divertor from iter

https://www.iter.org/mach/Divertor

Fusion Power

\[ E= Pt \]

JET produce short energy of \(Q=65\%\)

fusion power

Development

triple product vs temperature scatterplot

https://www.fusionenergybase.com/article/measuring-progress-in-fusion-energy-the-triple-products

Stellarators

\[ Q\rightarrow \infty \]

but complex engineering problem

Stellarators

Stellarators engineering view

Z-pinch

\[ -\nabla(p+\frac{B^2}{8\pi})+\frac{1}{4\pi}(B\cdot \nabla)B = 0 \]
\[ B = (0, B_\theta(r),B_z(r)) \]
\[ \begin{matrix} B_\theta = 0 & \theta~pinch & p(r) = \frac{-B_z^2}{8\pi}+P_0\\ B_z = 0 & z~pinch & p(r) = \frac{J_z\pi}{c^2}(a^2-r^2) \end{matrix} \]

unstable

z-pinch

Sausage Instability

z-pinch sausage instability

Kink Instability

z-pinch kink instability

https://www.fusionenergybase.com/concept/z-pinch

https://sites.uw.edu/zpinchlab/z-pinch-attributes/

\(\theta\) pinch

Lenz's Law, current of plasma is reverse of the applied current in coils

Field Reverse Configuration(FRC)

self stable torus plasma

theta-pinch

General Fusion

general_fusion

https://generalfusion.com/

Helion

\(\ce{D + He_3 -> p + He_4}\)

https://www.helionenergy.com

https://indico.cern.ch/event/776181/contributions/3376242/attachments/1856269/3049031/The_Field-Reversed_Configuration_FRC_Plasma_as_v2.pdf

Plasma

\(u_i\) : velocity of particle \(i\)

\(q_i\) : electric charge of particle \(i\)

\(m_i\) : mass of particle \(i\)

\(\rho\) : mass density \(\rho = \sum_i m_i\)

\(\gamma\) : adiabatic index

current density $$ J = \sum_i q_iu_i $$ center of mass velocity $$ v = \frac{1}{\rho}\sum_i m_i u_i $$ equation of continuity $$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho v) = 0 $$ equation of state $$ \frac{d}{dt}\left(\frac{p }{\rho^\gamma}\right) = 0 $$ equation of motion $$ \rho\left(\frac{\partial }{\partial t}+v\cdot \nabla\right)v = J\times B-\nabla p $$ $$ \frac{\partial \vec{B}}{\partial t} = \nabla \times (\vec{v} \times \vec{B} - \eta \nabla \times \vec{B}) $$

\[ \rho \frac{\partial \vec{v}}{\partial t} = - \nabla p + \frac{1}{\mu_0} (\nabla \times \vec{B}) \times \vec{B} + \vec{J} \times \vec{B} + \rho \vec{g} \]

MHD Equilibrium

$$ \nabla p = j\times B $$ - \(p\) kinematic pressure - \(j\) current density - \(B\) magnetic field

Maxwell

\[ \mu_0 J = \nabla\times B \]

Grad-Shafranov Equation

Cartesian coordinates $$ \nabla^2 A = -\mu_0\frac{d}{dA}(p+\frac{B_z^2}{2\mu_0}) $$

  • \(B = \left(\frac{\partial A(x,y)}{\partial y},-\frac{\partial A(x,y)}{\partial x}, B_z(x,y)\right)= \nabla A\times \hat z + B_z\hat z\)
  • \(A\) : vector potential
  • \(p\) : pressure in the plasma