d_wire = d * 10**-6,
emissivity = 0.3, l_wire=l_wire*10**-2,
beam_shape="Flat", l_beam = l_beam* 10**-2,
T_cracker = T_cracker,T_atoms = T_cracker,
phi_beam= (A_beam/ 10**-4) * 10**phi_exp, # Normalized to cm**2
T_base=wire_no_beam.record_dict["T_distribution"][-1]
)
wire.simulate(n_steps=n_steps, record_steps=record_steps,
time_step=time_step)
run_name = "lw_{}_phi_{}_Tc_{}".format(l_wire,phi_exp,
T_cracker)
os.makedirs(plot_dir + "signal/", exist_ok=True)
os.makedirs(plot_dir + "R_over_t/", exist_ok=True)
wire.plot_signal(plot_dir + "signal/{}".format(run_name))
wire.plot_R_over_t(plot_dir + "R_over_t/{}".format(run_name))
os.makedirs(plot_dir + "heat_flow/", exist_ok=True)
wire.plot_heat_flow(plot_dir + "heat_flow/{}".format(run_name))
wire.save(results_dir + "{}".format(run_name))
time_after = time()
run_time = time_after - time_before
print("finished run: " + run_name + "time required: "
+ "{0:.2f} minutes".format(run_time/60))
print("total time elapsed: {0:.2f} minutes".format(
(time() - start_time)/60.0))
l_wire_list = [ #5,
2.7
] # in cm
exp_list = [14, 15, 16, 17]
T_cracker_list = [2400, 2200, 2000, 1800, 1000, 500, 300, 0]
f_arr_full = np.zeros(
(len(l_wire_list), len(exp_list), len(T_cracker_list), 6))
for n_lw, l_wire in enumerate(l_wire_list):
for n_phi, phi_exp in enumerate(exp_list):
for n_T, T_cracker in enumerate(T_cracker_list):
run_name = "lw_{}_phi_{}_Tc_{}".format(l_wire, phi_exp, T_cracker)
wire = Wire()
wire = wire.load(top_dir + "results\\" + run_name)
wire.plot_heat_flow(top_dir + "plots\\" +
"heat_flow/log_{}".format(run_name),
log_y=True)
i = wire.n_wire_elements // 2
elem = [
wire.f_el(i),
wire.f_conduction(i),
wire.f_rad(i),
wire.f_beam(i),
wire.f_beam_gas(i),
wire.f_bb(i)
]
f_arr_full[n_lw, n_phi, n_T] = elem
# Plot scaling of Heat flow with temperature with various
if True:
for n_lw, l_wire in enumerate(l_wire_list):
T_atoms=T_cracker,
pressure=pressure,
p_laser=p_laser,
phi_beam=(A_beam / 10**-4) *
10**phi_exp, # Normalized to cm**2
T_base=wire_no_beam.record_dict["T_distribution"][-1])
wire.simulate(n_steps=n_steps,
record_steps=record_steps,
time_step=time_step)
run_name = "lb_{}".format(l_beam)
os.makedirs(plot_dir + "signal/", exist_ok=True)
os.makedirs(plot_dir + "R_over_t/", exist_ok=True)
wire.plot_signal(plot_dir + "signal/{}".format(run_name))
wire.plot_R_over_t(plot_dir + "R_over_t/{}".format(run_name))
os.makedirs(plot_dir + "heat_flow/", exist_ok=True)
wire.plot_heat_flow(plot_dir + "heat_flow/{}".format(run_name))
wire.plot_heat_flow(plot_dir +
"heat_flow/log_{}".format(run_name),
log_y=True)
wire.save(results_dir + "{}".format(run_name))
time_after = time()
run_time = time_after - time_before
print("finished run: " + run_name + "time required: " +
"{0:.2f} minutes".format(run_time / 60))
print("total time elapsed: {0:.2f} minutes".format(
(time() - start_time) / 60.0))