plasmoid 0.3.0

Fusion and self-organising plasma physics, computed: Coulomb barriers and Gamow peaks, Bosch–Hale and p–¹¹B reactivities, Lawson, screening and muon catalysis; Z-pinches and the plasma focus; plasmoids, Taylor relaxation, FRCs and reconnection; aneutronic power balance and direct energy conversion; each proposal against its physical limit.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30

//! From a fusion shot's energy to the reactions and fuel behind it: NIF's 7 April 2025 record and
//! JET's final 69 MJ pulse ([LLNL](https://annual.llnl.gov/fy-2025/national-ignition-facility-2025);
//! [UK Government, 2024](https://www.gov.uk/government/news/jets-final-tritium-experiments-yield-new-fusion-energy-record)).
//!
//! ```text
//! cargo run -p plasmoid --example fusion_yields
//! ```

use metrology::{C, KEV};
use plasmoid::aneutronic::Q_DT;
use plasmoid::fusion::{M_D, M_T};

fn main() {
    let per_reaction = Q_DT * 1e3 * KEV; // 17.59 MeV, J
    // NIF N250406: 2.08 MJ of laser light in, 8.6 MJ of fusion out.
    let (laser, yield_nif) = (2.08e6, 8.6e6);
    let reactions = yield_nif / per_reaction;
    println!("NIF: target gain {:.2}, {reactions:.2e} D–T reactions (one neutron each)", yield_nif / laser);
    assert!((yield_nif / laser - 4.13).abs() < 0.01 && (reactions / 3.05e18 - 1.0).abs() < 0.01);

    // JET, 2023: 69.26 MJ in 5.2 s from 0.21 mg of fuel.
    let burned = 69.26e6 / per_reaction * (M_D + M_T) * KEV / (C * C); // kg of deuterium and tritium
    println!("JET: {:.1} MW for 5.2 s burned {:.3} mg of D–T (published 0.21 mg)", 69.26 / 5.2, burned * 1e6);
    assert!((burned * 1e6 - 0.205).abs() < 0.002);
}

#[test]
fn yields_match_the_published_shots() {
    main();
}