Scrap heap that is only temporary

Initial Bubble Conditions

The DT gas bubble contains 3 mg of DT. Taking D–T molecular mass ≈5 g/mol (deuterium+tritium), the moles of gas are

n = (3×10^{−3}g)(mol/5g) = 6.0×10^{−4} mol.

At 750 K and 1 atm, we use the ideal-gas law (PV=nRT) to find the initial volume:

V_0=nRT_0/P_0
= (6.0×10^{−4} mol)(8.314 J/mol·K)(750 K)/(1.013×10^5 Pa)
≈ 3.69×10^{−5} m^3

Assuming a spherical bubble, the initial radius is

R_0=(3V_0/4π)^{1/3}≈(3×3.69×10^{−5}/4π)^{1/3}≈0.0207 m (2.07 cm).

Adiabatic Compression

The external pressure jumps to P_f=10,000 atm (≈ 1.013\times10^9 Pa). We assume the compression is fast and adiabatic (Q=0). For an ideal gas undergoing a (quasi-static) adiabatic process, P V^\gamma= constant, where \gamma is the heat-capacity ratio. For diatomic DT at high temperature (rotations excited, vibrations partially excited) we take \gamma\approx1.4. Thus

P_0V_0^γ=P_fV_f^γ ⟹ V_f=V_0(P_0/P_f)^{1/γ}.

Substituting P_0=1.013\times10^5 Pa, P_f=1.013\times10^9 Pa, and \gamma=1.4,

V_f ≈ 5.14×10^{−8}m^3.

The final bubble radius is R_f ≈ 2.3×10^{−3} m (2.3 mm).

Final Temperature

Adiabatic compression also raises the gas temperature. For an ideal gas, one can combine PV=nRT and P V^\gamma= const to get T_f/T_0=(V_0/V_f)^{γ−1} = (P_f/P_0)^{(γ−1)/γ}.

Using \gamma=1.4, T_0=750 K, P_f/P_0=10^4, we find

T_f=750×(10^4)^{0.4/1.4}
≈ 10,500 K.

Time to Maximum Compression

For a spherical gas bubble collapsing in an incompressible liquid under a sudden pressure jump, the Rayleigh collapse time gives a good first-order estimate:

T_c \approx 0.915 \cdot R_0 \sqrt{\frac{\rho_L}{P_\infty}}

Where:

  • T_c = collapse time
  • R_0 = initial bubble radius
  • ρ_L = liquid density
  • P_∞ = external overpressure driving the collapse

:pushpin: Example: DT Bubble in Molten FLiBe

  • R_0 = 2.07 cm
  • ρ_L = 2000 kg/m³ (approx. for FLiBe)
  • P_∞ = 10,000 atm ≈ 1.013×10^9 Pa

Plugging in:

T_c \approx 0.915 \cdot 0.0207 \cdot \sqrt{\frac{2000}{1.013 \times 10^9}} \approx 2.6 \times 10^{-5} \text{ s}

This 26 microseconds collapse time is consistent with a 50 μs pressure pulse duration. Thus the bubble should reach minimum radius before the pulse of pressure vacates the area.

Summary: The 3 mg DT bubble (initially at 750 K, 1 atm) has initial radius R_0\approx2.07 cm. Under a 10 000 atm adiabatic compression, the final radius is R_f\approx2.28 mm, the final temperature T_f\approx1.0\times10^4 K.

Here’s a pressure showdown across three very different domains—nuclear, mechanical, and ballistic. Let’s break it down:

:microscope: High-Temperature Light Water Reactor (LWR)

  • Pressure Range: ~155 bar (15.5 MPa or ~2,250 psi)
  • Why so high? This pressure keeps water in a liquid state even at temperatures around 315 °C (600 °F), allowing efficient heat transfer without boiling.
  • Design Implication: The reactor vessel and piping must withstand continuous high pressure and temperature for decades.

:articulated_lorry: Diesel Engine Firing Chamber (Max Compression)

  • Pressure Range: ~150–250 bar (2,200–3,600 psi) for heavy-duty engines
  • Peak Values: Some experimental or high-performance diesels can reach up to 300–400 bar (~4,000–5,800 psi)
  • Why so high? Compression ignites the fuel-air mixture without a spark plug. Higher pressure = more efficient combustion and torque.

:water_pistol: Bullet Shell Casing (During Firing)

  • Pressure Range: ~3,000–60,000 psi depending on caliber and load
  • Examples:
    • 9mm: ~35,000 psi standard, up to 62,000 psi with bullet setback
    • .223 Remington: ~55,000–60,000 psi typical chamber pressure
  • Why so extreme? The rapid combustion of gunpowder generates a shockwave that propels the bullet at high velocity. The casing must contain this pressure momentarily before the bullet exits.

:balance_scale: Pressure Comparison Table

System Typical Pressure Peak Pressure Notes
Light Water Reactor (LWR) ~2,250 psi ~2,250 psi Constant pressure to prevent boiling
Diesel Engine (Heavy Duty) ~2,500–3,600 psi Up to ~4,000 psi Varies by engine size and boost
Bullet Shell (e.g. 9mm/.223) ~35,000–60,000 psi Up to ~62,000 psi Instantaneous spike during firing6

:brain: Fun Insight

The pressure inside a bullet casing during firing can be more than 25 times greater than inside a nuclear reactor’s pressure vessel. And while diesel engines flirt with reactor-level pressures, they still fall short of the ballistic blast zone.

Here’s a first-principles pressure model for a spherically converging compressional wave in molten FLiBe approaching a deuterium-tritium (DT) fuel bubble.

:high_voltage: Quick Estimation Formula

We start by modeling a strong spherical shock wave propagating inward through a molten salt medium (FLiBe), converging on a center containing a gasious bubble of DT fuel. A simple scaling relation gives pressure as a function of radius r:

P(r) \approx P_0 \left( \frac{R_0}{r} \right)^{n}

  • P(r) = pressure at distance r from the bubble center
  • P_0 = reference pressure at radius R_0 (starting point of shock)
  • n = shock steepness exponent, typically 2≤n≤3.52 depending on wave strength and nonlinearity
  • This comes from energy conservation and the geometric convergence of a shock wave in 3D space

:microscope: Physical Assumptions

To make this tractable:

  • Idealized spherical symmetry: No azimuthal variations or instabilities
  • Strong shock approximation: Neglecting viscosity and thermal diffusion in molten FLiBe
  • Adiabatic compression: No heat loss during convergence
  • Constant properties in FLiBe: Density ρ, sound speed c, and Grüneisen parameter remain constant
  • Neglecting bubble recoil or deformation during compression (until late stages)

These allow us to invoke the Guderley self-similar solution for converging shocks, which gives:

P(r) \propto \left( \frac{R_0}{r} \right)^{2\gamma/(\gamma + 1)}

For molten FLiBe, assume:

  • Adiabatic index γ≈1.25 (approximation for ionized heavy fluorides)
  • This gives n≈1.11

:warning: Limitations and Where It Breaks Down

Here’s where the realism hits a wall:

  • Material nonlinearity: FLiBe’s compressibility varies under high pressure and temperature, so assuming constant γ oversimplifies
  • Shock preheating and ionization: The wave itself alters FLiBe’s local state, changing its response
  • Radiation transport: Near the DT bubble, bremsstrahlung and line radiation modify energy transport—these are absent in the model
  • Instabilities at interface: Rayleigh–Taylor and Richtmyer–Meshkov modes disrupt the bubble edge, breaking spherical symmetry
  • Phase changes: FLiBe might partially vaporize or ionize, adding opacity and feedback mechanisms not captured here
  • Bubble motion and rebound: Late-stage dynamics of the DT bubble aren’t modeled—this oversimplifies ignition prediction

:hammer_and_wrench: Improving the Model

  • Introduce equations of state (EOS) for FLiBe under compression
  • Couple hydrodynamic simulations to radiative transfer models
  • Account for real geometries and target imperfections
  • Include shock viscosity and conduction losses
  • Use multiphysics simulations like ALE or Lagrangian hydrocodes

compare the adiabatic index γ (gamma) of molten FLiBe and liquid water, then express the difference in terms of the shock steepness exponent nn from your spherical wave model.

:magnifying_glass_tilted_left: Step 1: Define the Adiabatic Index γ

The adiabatic index is defined as:

\gamma = \frac{C_p}{C_v}

Where:

  • C_p = specific heat at constant pressure
  • C_v = specific heat at constant volume

It governs how pressure responds to volume changes in adiabatic processes and directly influences shock wave behavior.

:thermometer: Step 2: Estimate γ for Each Fluid

:large_blue_diamond: Molten FLiBe

  • FLiBe is a molten salt mixture (LiF-BeF₂), highly ionic and dense.
  • Experimental and simulation data suggest:
    • C_p ≈ 2.4 J/g·K
    • C_v ≈ 1.9 J/g·K (estimated from EOS and compressibility)
  • So:

\gamma_{\text{FLiBe}} \approx \frac{2.4}{1.9} \approx 1.26

:large_blue_diamond: Liquid Water

  • Water is a polar molecular liquid with strong hydrogen bonding.
  • At room temperature:
    • C_p \approx 4.18 \, \text{J/g·K}
    • C_v \approx 3.1 \, \text{J/g·K} (approximate, varies with pressure)
  • So:

\gamma_{\text{H}_2\text{O}} \approx \frac{4.18}{3.1} \approx 1.35

:chart_increasing: Step 3: Convert to Shock Steepness Exponent n

From the Guderley-type scaling for spherical shocks:

n = \frac{2\gamma}{\gamma + 1}

Let’s compute:

Fluid γ n= \frac{2\gamma}{\gamma + 1}
FLiBe 1.26 \frac{2 \cdot 1.26}{1.26 + 1} \approx 1.12
Liquid Water 1.35 \frac{2 \cdot 1.35}{1.35 + 1} \approx 1.15

So liquid water has a slightly steeper shock exponent n than molten FLiBe, meaning pressure increases more sharply as the wave converges.

:warning: Caveats and Assumptions

  • These values are approximations based on available data and idealized behavior.
  • Real systems involve:
    • Phase transitions (especially water near boiling)
    • Ionization and radiation effects
    • Nonlinear EOS behavior at high pressures

:stopwatch: How Precise Does the Timing Need to Be?

For a multi-shock compression scheme to work effectively in bubble implosion or ignition scenarios, timing precision needs to be on the order of:

  • Microseconds for shock spacing (e.g. 10–50 μs between pulses)
  • Sub-microsecond jitter tolerance to avoid destructive interference or mistimed convergence
  • Nanosecond-level rise time if you’re shaping sharp pressure fronts

This is especially true if you’re trying to synchronize shocks with bubble collapse phases or exploit constructive interference near minimum radius.

:gear: How Accurate Are Piezoelectric Stack Actuators Today?

Modern piezo stack actuators can deliver:

  • Rise times: as fast as 10–50 μs, depending on voltage and load
  • Pulse widths: tunable down to 10–20 μs with specialized drivers
  • Force modulation: up to 70 kN, with accelerations exceeding 10,000 g
  • Repetition rates: up to kHz burst mode, with synchronization across multiple stacks

High-voltage pulser systems (e.g. 1000 V, 400 A peak) can trigger stacks with μs-level precision, and multi-stack arrays can be synchronized to deliver staged pulses.

So yes—today’s piezo systems can absolutely deliver multi-shock sequences with the timing and force profiles needed for staged compression.

:rocket: Can We Build a Multi-Shock Pulse System Today?

Yes, but with caveats:

  • You’d need custom electronics: high-voltage pulsers with programmable delay lines and μs-level triggering
  • Mechanical coupling must be tight: stacks must transmit pressure into FLiBe or similar medium without damping or delay
  • Thermal management is critical: repeated pulses can cause heating and drift
  • Synchronization across stacks must be precise—think FPGA-level control or laser-triggered timing

Systems like PIA impulse generators and piezocomposite shakers already support this kind of operation in material testing and shock generation.

:brain: TL;DR

  • Timing precision: ~10 μs spacing, <1 μs jitter
  • Piezo accuracy: μs-level pulse delivery, kHz burst rates, high force
  • Feasibility: Yes—with custom drivers, synchronized stacks, and careful mechanical design

Want to sketch out a three-pulse ignition sequence with realistic actuator specs? Or riff on how to explain this to skeptics in your forum: “Turns out, you can knock three times—with piezo fists and microsecond swagger.”

Let’s walk through a first-principles estimate of the speed of sound in molten 2LiF–BeF₂ at 750 K and 1 atm, using thermodynamic fundamentals and atomic data. This won’t match high-fidelity EOS fits, but it’ll give a transparent baseline you can tweak or stress-test.

:abacus: Step 1: Atomic Weights and Molar Mass

From online sources2:

Compound Formula Molar Mass (g/mol)
LiF Li + F 25.94
BeF₂ Be + 2F 47.01
2LiF + BeF₂ 2×25.94 + 47.01 = 98.89 g/mol

:ice: Step 2: Estimate Density at 750 K

Experimental data for FLiBe (66–34 mol%) gives densities around 1.94–2.0 g/cm³ at 750 K. Since 2LiF–BeF₂ is slightly more LiF-rich, we’ll estimate:

  • ρ ≈ 1.95 g/cm³ = 1950 kg/m³

:triangular_ruler: Step 3: Estimate Bulk Modulus

The speed of sound c in a fluid is:

c= \sqrt{\frac{K}{\rho}}

Where:

  • K = bulk modulus
  • ρ = density

For molten salts, bulk modulus is rarely tabulated directly. But we can estimate it from compressibility β_T:

K= \frac{1}{\beta_T}

Typical isothermal compressibility for molten FLiBe is:

  • β_T \approx 4.5 \times 10^{-10} \, \text{Pa}^{-1}

So:

  • K \approx 2.22 \times 10^9 \, \text{Pa}

:rocket: Step 4: Calculate Speed of Sound

Plugging into the formula:

c= \sqrt{\frac{2.22 \times 10^9}{1950}} \approx \sqrt{1.138 \times 10^6} \approx \boxed{3373 \, \text{m/s}}

:white_check_mark: Final Result

Estimated speed of sound in molten 2LiF–BeF₂ at 750 K and 1 atm: ~3370 m/s
For Americans, that translates into 1/8 inch per microsecond

This aligns well with the experimentally reported ~3400 m/s for eutectic FLiBe2. The slight deviation reflects the LiF-rich composition and assumptions in compressibility.