added two parameters including cost_perN and F_max to econ and simula…

…tion
symbiotic-engineering · Feb 17, 2024 · 8779b55 · 8779b55
1 parent 3828aae
commit 8779b55
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 6 deletions.
diff --git a/inputs/parameters.m b/inputs/parameters.m
@@ -50,6 +50,6 @@
             'LCOE_max', .5,...                  % maximum LCOE ($/kWh)
             'power_max', Inf,...                % maximum power (W)
             'eff_pto',   0.80,...               % PTO efficiency (-)
-            'eff_array', 0.95*0.98 );           % array availability and transmission efficiency (-)
-
+            'eff_array', 0.95*0.98,...          % array availability and transmission efficiency (-)
+            'cost_perN', 0.1 );                 % cost per Newton   
 end
diff --git a/simulation/modules/econ.m b/simulation/modules/econ.m
@@ -1,5 +1,5 @@
 
-function [LCOE, capex, opex] = econ(m_m, M, cost_m, N_WEC, P_elec, FCR, efficiency)
+function [LCOE, capex, opex] = econ(m_m, M, cost_m, N_WEC, P_elec, FCR, cost_perN, F_max, efficiency)
 
 structural_cost = m_m.* cost_m;
 devicestructure =  N_WEC * structural_cost(M);
@@ -9,7 +9,7 @@
 development     = 4553000;
 infrastructure  = 990000;
 mooring         = N_WEC * 525000;
-pto             = N_WEC * 623000;
+pto             = N_WEC * (623000 + F_max * cost_perN);
 profitmargin    = 356000;
 installation    = 5909000;
 contingency     = 1590000;

diff --git a/simulation/simulation.m b/simulation/simulation.m
@@ -48,7 +48,7 @@
                                     in.M, in.h_s, in.T_s, in.rho_w, in.g, in.sigma_y, A_c, ...
                                     A_lat_sub, r_over_t, I, in.E);
 
-LCOE = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.eff_array);
+LCOE = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.cost_perN, in.F_max, in.eff_array);
 
 %% Assemble constraints g(x) >= 0
 g = zeros(1,18);
@@ -83,7 +83,7 @@
                                             in.t_fc, in.t_fb, in.t_sr, in.t_dt,...
                                             D_d, in.D_dt, in.theta_dt, in.T_s, in.h_d, ...
                                             in.M, in.rho_m, in.rho_w, in.m_scale);
-    [~,capex,opex] = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.eff_array);
+    [~,capex,opex] = econ(m_m, in.M, in.cost_m, in.N_WEC, P_elec, in.FCR, in.cost_perN, in.F_max, in.eff_array);
     val.mass_f  = mass(1);
     val.mass_vc = mass(2);
     val.mass_rp = mass(3);