def parameter_sweep(sensitivity_input):
start = time.time()
# create system
inputs = ICAES.get_default_inputs()
for variable in sensitivity_input.index:
inputs[variable] = sensitivity_input[variable]
system = ICAES(inputs=inputs)
# run single cycle and analyze
system.single_cycle()
results = system.analyze_performance()
end = time.time()
results['solve_time'] = end - start
# combine inputs and results to return in single series
single_output = pd.concat([sensitivity_input, results])
return single_output
def parameter_sweep(sweep_input):
start = time.time()
# create system
inputs = ICAES.get_default_inputs()
inputs['steps'] = sweep_input['steps'] # [kPa]
inputs['p_store_min'] = sweep_input['p_min'] # [kPa]
inputs['p_store_init'] = sweep_input['p_min'] # [kPa]
inputs['p_store_max'] = sweep_input['p_max'] # [kPa]
inputs['eta_storage'] = sweep_input['eta_storage'] # [-]
inputs['V_res'] = sweep_input['V_res'] # [m3]
inputs['phi'] = sweep_input['phi'] # porosity [-]
inputs['Slr'] = sweep_input['Slr'] # liquid residual fraction [-]
system = ICAES(inputs=inputs)
# run single cycle and analyze
system.single_cycle()
results = system.analyze_performance()
end = time.time()
results['solve_time'] = end - start
# combine inputs and results to return in single series
single_output = pd.concat([sweep_input, results])
return single_output
from caes import ICAES, plot_series
import matplotlib.pyplot as plt
# create system
inputs = ICAES.get_default_inputs()
inputs['depth'] = 5550.0
system = ICAES(inputs=inputs)
# system.debug_perf()
#
# # run single cycle and analyze
system.single_cycle()
results = system.analyze_performance()
results.to_csv('single_cycle_performance.csv')
system.data.to_csv('single_cycle_timeseries.csv')
print(results)
# # plot results
system.plot_overview()
system.plot_pressures()
system.plot_pressure_losses()
def parameter_sweep(sweep_input, debug=True):
start = time.time()
# convert inputs to model units
kW_out = sweep_input['capacity_MW'] * 1e3
kWh_out = sweep_input['capacity_MW'] * sweep_input['duration_hr'] * 1e3
if debug:
print('\nStarting sizing')
print("kW_out (desired) : " + str(round(kW_out, 3)))
print("kWh_out (desired) : " + str(round(kWh_out, 3)))
print("\nDepth (m) : " + str(sweep_input['depth_m']))
print("Thickness (m) : " + str(sweep_input['thickness_m']))
print("Porosity (-) : " + str(sweep_input['porosity']))
print("Permeability (mD): " + str(sweep_input['permeability_mD']))
# allowable calculation error
error = 1e-6
# maximum number of iterations
count_max = 200
# initial guesses
m_dot = 10.0
r_f = 10.0
# initial results
kW_out_actual = 0.0
kWh_out_actual = 0.0
count = 0
while abs(kW_out_actual - kW_out) / kW_out + abs(
kWh_out_actual - kWh_out) / kWh_out > error:
if debug:
print("\nIteration : " + str(count))
print("m_dot : " + str(round(m_dot, 3)))
print("r_f : " + str(round(r_f, 3)))
# create system
inputs = ICAES.get_default_inputs()
# user inputs
inputs['depth'] = sweep_input['depth_m'] # porosity depth [m]
inputs['h'] = sweep_input['thickness_m'] # porosity thickness [m]
inputs['phi'] = sweep_input['porosity'] # formation porosity [-]
inputs['k'] = sweep_input[
'permeability_mD'] # formation permeability [mD]
# current guess/iteration
inputs['m_dot'] = m_dot # [kg/s]
inputs['r_f'] = r_f # [m]
system = ICAES(inputs=inputs)
# run single cycle and analyze
system.single_cycle()
results = system.analyze_performance()
# extract results of interest
kW_out_actual = results['kW_out_avg']
kWh_out_actual = results['kWh_out']
# update guesses
tau = 0.5 # solver time constant
m_dot = m_dot * (1.0 + tau * (kW_out - kW_out_actual) / kW_out
) # m_dot linearly linked to kW
r_f = r_f * (1.0 + tau**2 * (kWh_out - kWh_out_actual) / kWh_out_actual
) # r_f exponentially linked to kWh
count = count + 1
if count > count_max: # sizing unsuccessful
results['errors'] = True
break
if debug:
print("MW_out_avg : " + str(round(kW_out_actual / 1e3, 3)))
print("MWh_out : " + str(round(kWh_out_actual / 1e3, 3)))
end = time.time()
results['solve_time'] = end - start
results['iterations'] = count
results['m_dot'] = m_dot
results['r_f'] = r_f
# print out RTE
print(results['RTE'])
# combine inputs and results to return in single series
single_output = pd.concat([sweep_input, results])
return single_output