PUH-BrianMartell puh_throat_impedance_formula_simulation.py-Updated New Py Code

puh_throat_impedance_formula_simulation.py

import numpy as np
import matplotlib.pyplot as plt

# PUH v25: Throat Impedance Formula Sim — Z vs Frequency (Sharp Dip Zero)
f = np.linspace(0.5, 1.5, 500)  # Normalized frequency
f_res = 1  # Resonance
eta_density = 100  # Rebound density factor
Z_standard = np.ones_like(f)  # Standard constant Z_0
Z_throat = 1 / (1 + (f / f_res)**2 * eta_density)  # PUH sharp dip

plt.figure(figsize=(10,6))
plt.plot(f, Z_standard, label='Standard Vacuum Constant Z_0', color='red', lw=2)
plt.plot(f, Z_throat, label='PUH Throat Sharp Resonance Dip Z\to 0', color='cyan', lw=2)
plt.axvline(f_res, color='gold', ls='--', label='Resonance f_{res}')
plt.xlabel('Normalized Frequency f / f_{res}'); plt.ylabel('Impedance Z (arb.)')
plt.title('PUH v25: Throat Impedance Formula Sim')
plt.legend(); plt.grid(alpha=0.3)
plt.tight_layout()
plt.savefig('puh_throat_impedance_formula_simulation.png', dpi=300)
plt.show()

print("Standard constant, PUH sharp dip zero at resonance — tachyonic zero-drag.")