From 461484b6fa6d811718cf2eb076c2fc86e9bb15ae Mon Sep 17 00:00:00 2001 From: Madison Dietrich Date: Fri, 14 Jul 2023 18:32:05 -0400 Subject: [PATCH 1/5] fixed opex calculation by doing power equation --- mdocean/simulation/modules/econ.m | 1 + 1 file changed, 1 insertion(+) diff --git a/mdocean/simulation/modules/econ.m b/mdocean/simulation/modules/econ.m index 578c7aa6..39e68374 100644 --- a/mdocean/simulation/modules/econ.m +++ b/mdocean/simulation/modules/econ.m @@ -26,6 +26,7 @@ opex = operations + postinstall + shoreoperations + replacement ... + consumables + insurance; +opex=1193447.03 * N_WEC^.443321491; hr_per_yr = 8766; P_avg = N_WEC * P_elec * efficiency; From 12d0934e24ff6444035ccbde2b40647e65b09955 Mon Sep 17 00:00:00 2001 From: Madison Dietrich Date: Tue, 18 Jul 2023 16:59:31 -0400 Subject: [PATCH 2/5] Fixed syntax of opex equation and added comments --- mdocean/simulation/modules/econ.m | 13 ++++--------- 1 file changed, 4 insertions(+), 9 deletions(-) diff --git a/mdocean/simulation/modules/econ.m b/mdocean/simulation/modules/econ.m index 578c7aa6..9623833e 100644 --- a/mdocean/simulation/modules/econ.m +++ b/mdocean/simulation/modules/econ.m @@ -17,15 +17,10 @@ capex = development + infrastructure + mooring + devicestructure + pto ... + profitmargin + installation + contingency; -operations = N_WEC * 27000; -postinstall = 710000; -shoreoperations = 142000; -replacement = N_WEC * 54000; -consumables = N_WEC * 8000; -insurance = (.8 + .2*N_WEC) * 227000; - -opex = operations + postinstall + shoreoperations + replacement ... - + consumables + insurance; +% opex includes cost of operation, postinstall, replacement, consumables, and +% insurance +% opex equation was generated by Casio online power regression calculator +opex = 1193000 * N_WEC^0.4433; hr_per_yr = 8766; P_avg = N_WEC * P_elec * efficiency; From e23d9573ab2fffa832052aae1797c7eb12293d4c Mon Sep 17 00:00:00 2001 From: Madison Dietrich Date: Wed, 19 Jul 2023 12:12:39 -0400 Subject: [PATCH 3/5] fixed opex equation syntax and added comments --- mdocean/simulation/modules/econ.m | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/mdocean/simulation/modules/econ.m b/mdocean/simulation/modules/econ.m index 9623833e..cf7898a8 100644 --- a/mdocean/simulation/modules/econ.m +++ b/mdocean/simulation/modules/econ.m @@ -17,7 +17,7 @@ capex = development + infrastructure + mooring + devicestructure + pto ... + profitmargin + installation + contingency; -% opex includes cost of operation, postinstall, replacement, consumables, and +% opex reflects cost of operation, postinstall, replacement, consumables, and % insurance % opex equation was generated by Casio online power regression calculator opex = 1193000 * N_WEC^0.4433; From 7d4dd948ad67df671aca66d026e197787eed4b36 Mon Sep 17 00:00:00 2001 From: rebeccamccabe Date: Wed, 8 Nov 2023 16:18:37 -0500 Subject: [PATCH 4/5] divide power by 8 --- mdocean/simulation/modules/dynamics/dynamics.m | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/mdocean/simulation/modules/dynamics/dynamics.m b/mdocean/simulation/modules/dynamics/dynamics.m index d4431c36..599a3ee0 100644 --- a/mdocean/simulation/modules/dynamics/dynamics.m +++ b/mdocean/simulation/modules/dynamics/dynamics.m @@ -39,7 +39,7 @@ X_unsat = get_response(w,m,b,k,Fd); % confirm unsaturated response doesn't exceed maximum capture width - P_unsat = 1/2 * in.B_p * w.^2 .* X_unsat.^2; + P_unsat = 1/8 * 1/2 * in.B_p * w.^2 .* X_unsat.^2; F_ptrain_over_x = sqrt( (in.B_p * w).^2 + (K_p).^2 ); F_ptrain_unsat = F_ptrain_over_x .* X_unsat; @@ -54,7 +54,7 @@ X_sat = get_response(w,m,b_sat,k_sat,Fd); % calculate power - P_matrix = 1/2 * (mult * in.B_p) .* w.^2 .* X_sat.^2; + P_matrix = 1/8 * 1/2 * (mult * in.B_p) .* w.^2 .* X_sat.^2; X_max = max(X_sat,[],'all'); h_s_extra = (in.h_s - in.T_s - (in.h_f - in.T_f) - X_max) / in.h_s; % extra height on spar after accommodating float displacement From 6d6d0b8bc2c1a7916bf41d828841f8d7608830ff Mon Sep 17 00:00:00 2001 From: rebeccamccabe Date: Thu, 9 Nov 2023 22:24:37 -0500 Subject: [PATCH 5/5] pareto front plot changes for maha princeton talk --- .../multiobjective/pareto_curve_heuristics.m | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/mdocean/optimization/multiobjective/pareto_curve_heuristics.m b/mdocean/optimization/multiobjective/pareto_curve_heuristics.m index 4e0d7ac6..7c02558c 100644 --- a/mdocean/optimization/multiobjective/pareto_curve_heuristics.m +++ b/mdocean/optimization/multiobjective/pareto_curve_heuristics.m @@ -70,7 +70,7 @@ function pareto_curve_heuristics() P_var_solar = 125; % balanced design - [~,idx_balanced] = min(abs(Pvar-40)); + [~,idx_balanced] = min(abs(Pvar-100)); LCOE_balanced = LCOE(idx_balanced); P_var_balanced = Pvar(idx_balanced); @@ -126,32 +126,32 @@ function pareto_curve_heuristics() if showSingleObj text(LCOE(idx_best_LCOE)+.03,Pvar(idx_best_LCOE),'Cheapest','FontSize',sz) text(LCOE(idx_best_Pvar)+.03,Pvar(idx_best_Pvar)-3,'Least Variable','FontSize',sz) - text(LCOE_balanced+.02,P_var_balanced+5,'Balanced Design','FontSize',sz) + text(LCOE_balanced-.15,P_var_balanced+5,'Balanced Design','FontSize',sz) end if showImages mini_plot_size = [.2 .22]; % small corner pictures of best geometries % upper left - axes('Position',[.28 .7 mini_plot_size]) + axes('Position',[.28 .6 mini_plot_size]) box on visualize_geometry(x_best_LCOE,p,true); set(gca,'XTickLabel',[],'YTickLabel',[]) % lower right - axes('Position',[.52 .15 mini_plot_size]) + axes('Position',[.51 .23 mini_plot_size]) box on visualize_geometry(x_best_Pvar,p,true); set(gca,'XTickLabel',[],'YTickLabel',[]) % balanced - axes('Position',[.22 .25 mini_plot_size]) + axes('Position',[.10 .28 mini_plot_size]) box on visualize_geometry(x_balanced,p,true); set(gca,'XTickLabel',[],'YTickLabel',[]) % RM3 - axes('Position',[.7 .5 mini_plot_size]) + axes('Position',[.7 .53 mini_plot_size]) box on visualize_geometry(x_nom,p,true); set(gca,'XTickLabel',[],'YTickLabel',[])