Recommend parameter ranges to achieve 25 kW target power.

target_power = 25000  # W

print("=" * 70)
print("  Parameter Recommendations to Achieve 25 kW")
print("=" * 70)

# Current tested ranges
f_max_tested = 199.5e3  # Hz
E_p_max_tested = 30e-3  # J
power_max_tested = f_max_tested * E_p_max_tested / 1000  # kW

print(f"\nCurrent tested ranges achieved: {power_max_tested:.2f} kW")
print(f"Target power: {target_power / 1000:.0f} kW")
print(f"Multiplier needed: {target_power / (f_max_tested * E_p_max_tested):.2f}x")

print("\n" + "-" * 70)
print("  Option 1: Increase Frequency (keep E_p = 30 mJ)")
print("-" * 70)
f_needed = target_power / 30e-3
print(f"Required frequency: {f_needed / 1e3:.1f} kHz")
print(f"Current max: {f_max_tested / 1e3:.1f} kHz")
print(f"Increase needed: {f_needed / f_max_tested:.2f}x")

print("\n" + "-" * 70)
print("  Option 2: Increase Energy per Pulse (keep f = 200 kHz)")
print("-" * 70)
E_p_needed = target_power / 200e3
print(f"Required E_p: {E_p_needed * 1e3:.1f} mJ")
print(f"Current max: {E_p_max_tested * 1e3:.1f} mJ")
print(f"Increase needed: {E_p_needed / E_p_max_tested:.2f}x")

print("\n" + "-" * 70)
print("  Option 3: Balanced Increase")
print("-" * 70)
# Double both parameters
f_balanced = 400e3  # 400 kHz
E_p_balanced = 60e-3  # 60 mJ
power_balanced = f_balanced * E_p_balanced
print(f"f = {f_balanced / 1e3:.0f} kHz, E_p = {E_p_balanced * 1e3:.0f} mJ")
print(f"Power = {power_balanced / 1e3:.1f} kW")
print()
# Alternative: smaller increase in both
f_balanced2 = 300e3  # 300 kHz
E_p_balanced2 = 83e-3  # 83 mJ
power_balanced2 = f_balanced2 * E_p_balanced2
print(f"f = {f_balanced2 / 1e3:.0f} kHz, E_p = {E_p_balanced2 * 1e3:.0f} mJ")
print(f"Power = {power_balanced2 / 1e3:.1f} kW")

print("\n" + "=" * 70)
print("  RECOMMENDATIONS")
print("=" * 70)
print("\n1. **Immediate Action**: Update parameter sweep ranges")
print("   - Frequency: 10 - 900 kHz (vs current 10-200 kHz)")
print("   - Energy per pulse: 1 - 130 mJ (vs current 1-30 mJ)")
print("\n2. **Conservative Approach** (2x increase in each):")
print("   - Frequency: 10 - 400 kHz")
print("   - Energy per pulse: 1 - 60 mJ")
print("   - Expected max power: ~24 kW (close to target)")
print("\n3. **Verify other constraints remain satisfied**")
print("   - All top candidates had positive margins on:")
print("     * Breakdown voltage")
print("     * Coverage fraction")
print("     * Energy density")
print("   - These should remain valid at higher f and E_p")
print()