CEM and enviro modules as used for umerc by rebeccamccabe · Pull Request #96 · symbiotic-engineering/MDOcean

rebeccamccabe · 2025-06-12T04:57:17Z

No description provided.

rebeccamccabe · 2025-06-12T05:25:42Z

dev/dynamics/force_sat/bode_fit.m

+% load data from MEEM (hardcoded, will replace later)
+dat = readmatrix('mdocean/mdocean/inputs/validation/MEEM_validation/damping.csv');
+w = dat(:,1);
+mag = dat(:,2);


the data you are currently testing with is actually (k, B) rather than (w, mag).
try the following for better test data:

p = parameters(); b = var_bounds(); X = [b.X_noms; 1]; [~, ~, ~, ~, val] = simulation(X,p); w = val.w; F = val.gamma_f_over_rho_g * p.rho_w * p.g * p.Hs / (2*sqrt(2)) ; X = val.X_f; mag_matrix = X ./ F; mag = mag_matrix(1,:); % all freqs at lowest wave height

rebeccamccabe · 2025-06-12T05:26:22Z

dev/environment/component.m

This file is empty, please check and push again

codecov · 2025-06-12T08:28:45Z

❌ 40 Tests Failed:

Tests completed	Failed	Passed	Skipped
83	40	43	0

View the top 3 failed test(s) by shortest run time

test::allFiguresRun(which_figs=Fig_24_ParetoFrontWithLCOEContours,which_tabs=none24)

Stack Traces | 0.000623s run time

Error occurred in test/allFiguresRun(which_figs=Fig_24_ParetoFrontWithLCOEContours,which_tabs=none24) and it did not run to completion.
    ---------
    Error ID:
    ---------
    'MATLAB:badsubscript'
    --------------
    Error Details:
    --------------
    Index exceeds the number of array elements. Index must not exceed 1.
    
    Error in pareto_search>get_seeds_epsilon_constraint (line 123)
        P_var_min = objs_opt(2);
    
    Error in pareto_search (line 18)
        [X0,fvals,probs] = get_seeds_epsilon_constraint(p,b,num_DVs,objFcn2,turnLCOEtoPower);
    
    Error in all_figures (line 216)
            pareto_search(filename_uuid);
    
    Error in test/runAllFigsTabs (line 123)
                     f_runtime_run,t_runtime_run] = all_figures( run_figs, run_tabs, testCase.uuid );

test::allFiguresRun(which_figs=Fig_36_DampingPlatePlateAspectRatio,which_tabs=none36)

Stack Traces | 0.000624s run time

Error occurred in test/allFiguresRun(which_figs=Fig_36_DampingPlatePlateAspectRatio,which_tabs=none36) and it did not run to completion.
    ---------
    Error ID:
    ---------
    'MATLAB:UndefinedFunction'
    --------------
    Error Details:
    --------------
    Unrecognized function or variable 'BoedoPrantilAnnularPlate'.
    
    Error in all_figures (line 289)
            BoedoPrantilAnnularPlate()
    
    Error in test/runAllFigsTabs (line 123)
                     f_runtime_run,t_runtime_run] = all_figures( run_figs, run_tabs, testCase.uuid );

test::allFiguresRun(which_figs=Fig_23_DesignAndObjectiveHeuristics,which_tabs=none23)

Stack Traces | 0.000625s run time

Error occurred in test/allFiguresRun(which_figs=Fig_23_DesignAndObjectiveHeuristics,which_tabs=none23) and it did not run to completion.
    ---------
    Error ID:
    ---------
    'MATLAB:badsubscript'
    --------------
    Error Details:
    --------------
    Index exceeds the number of array elements. Index must not exceed 1.
    
    Error in pareto_search>get_seeds_epsilon_constraint (line 123)
        P_var_min = objs_opt(2);
    
    Error in pareto_search (line 18)
        [X0,fvals,probs] = get_seeds_epsilon_constraint(p,b,num_DVs,objFcn2,turnLCOEtoPower);
    
    Error in all_figures (line 216)
            pareto_search(filename_uuid);
    
    Error in test/runAllFigsTabs (line 123)
                     f_runtime_run,t_runtime_run] = all_figures( run_figs, run_tabs, testCase.uuid );

To view more test analytics, go to the Test Analytics Dashboard
_{📋 Got 3 mins? Take this short survey to help us improve Test Analytics.}

…ax=inf in max CW check

…if statement

…ring/MDOcean into anthony-enviro

…kup table

…ring/MDOcean into anthony-enviro

…cean into anthony-enviro

rebeccamccabe · 2025-09-02T05:32:38Z

combined2.m

move to appropriate foldere

rebeccamccabe · 2025-09-02T05:34:28Z

mdocean/inputs/scenario_outputs.csv

rebeccamccabe · 2025-09-02T05:34:37Z

mdocean/inputs/scenario_outputs_2.csv

rebeccamccabe · 2025-09-02T05:35:00Z

mdocean/inputs/scenario_outputs_3.csv

rename to remove _3

rebeccamccabe · 2025-09-02T05:35:30Z

mdocean/inputs/var_bounds.m

 % calibrations of nominal values
-b.F_max_nom = find_nominal_inputs(b, parameters(mode));
+
+% FIXME: reminder to uncomment after debugging


don't merge this

rebeccamccabe · 2025-09-02T05:40:18Z

mdocean/simulation/modules/grid_CEM.m

+    nP = numel(pVals);
+
+    % interpolated values
+    %varsAll   = Tsub.Properties.VariableNames;


remove commented

rebeccamccabe · 2025-09-02T05:40:46Z

mdocean/simulation/modules/grid_CEM.m

+function [delta_CEM_CO2, CEM_wec_capacity, delta_CEM_grid_cost,margin_to_viability] = CEM_lookup_table(zeta, omega_n, capacity_cost, power_lim_frac, B_p, location, params)
+
+
+    %if zeta == 0.05 && omega_n == 0.4 && strcmp(location,'ISONE')


remove commented

rebeccamccabe · 2025-09-02T05:40:51Z

mdocean/simulation/modules/grid_CEM.m

+    delta_CEM_CO2 = CEM_CO2_no_wave - CEM_CO2;
+
+    %{
+    isValidLookupLocation = true;% placeholder


remove commented

rebeccamccabe · 2025-09-02T05:41:31Z

mdocean/simulation/simulation.m

+                                                in.location, in);
+
+
+% old code that doesn't normalize with the number of WECs


remove commented

rebeccamccabe · 2025-09-02T05:42:24Z

pubs/journal-paper/sections/model.tex.rej

don't merge

anthnlong and others added 3 commits June 11, 2025 15:45

initial bode fit with dummy data

99ea5d7

remove DS_Store and add to gitignore

83e1135

fix fit_result typo

5a26f0e

rebeccamccabe commented Jun 12, 2025

View reviewed changes

Base automatically changed from force-saturation-fix to main June 16, 2025 21:41

anthnlong and others added 24 commits June 18, 2025 00:43

use matlab variables instead of hardocded path

05a7163

Merge branch 'main' into anthony-enviro

3a9e239

rewrote environment optimization

ae0d435

fix dimension error in F calculation

7bb4672

add A_hull output and fix bottom plate area mistake

38074b7

inf fix warning in run_single by using p.use_sat=false instead of F_m…

5e8bd65

…ax=inf in max CW check

script for plotting and fitting dummy cem data

4e24354

move devcontainer to work for all pubs and add image editor

385f828

add X phase to val and create bode sysid plots

10e3eed

xdsm with added grid (CEM) and enviro modules and fixed padding

d0659ed

paper updates: add data flow plot, bode plot, and some words

7852ea6

manually sync overleaf .bib since patch was rejected

1dd3825

Sync pubs/UMERC-2025-grid-value with Overleaf project

b7d0985

Merge remote-tracking branch 'origin/decider-pub' into anthony-enviro

87b5bbf

move gitignore and add xdsm pdf

293cd3d

spacing, title, other cosmetic tweaks to bode plot

34d55a1

write all methods sections except LCA

a7b5c77

Sync pubs/UMERC-2025-grid-value with Overleaf project

4098a24

Sync pubs/UMERC-2025-grid-value with Overleaf project

12af81b

Sync pubs/UMERC-2025-grid-value with Overleaf project

5a1855b

write LCA section and switch PG diagram to not have beta

6e6706d

manual merge of .rej file for RAO equation

cc234bd

Sync pubs/UMERC-2025-grid-value with Overleaf project

0625059

Sync pubs/UMERC-2025-grid-value with Overleaf project

1c31e35

anthnlong and others added 27 commits August 7, 2025 16:54

comment

22cd16d

fix indentation, var names, and add error handling to non viable CEM …

0525007

…if statement

Merge branch 'anthony-enviro' of https://github.com/symbiotic-enginee…

4032353

…ring/MDOcean into anthony-enviro

remove dummy function

0894170

Merge branch 'anthony-enviro' of https://github.com/symbiotic-enginee…

034ab14

…ring/MDOcean into anthony-enviro

comment out plots in fit_from_vars

e3a320f

comment out the nan-filling in places where JPD=0

bea25a8

fit_from_vars fix

11c3d5a

fully get rid of plots and prints in grid_CEM and fit_from_vars

13e4f94

Merge branch 'anthony-enviro' of https://github.com/symbiotic-enginee…

ade29ac

…ring/MDOcean into anthony-enviro

moved CSV loading to parameters.m

85b19e9

Merge branch 'anthony-enviro' of https://github.com/symbiotic-enginee…

70d1d08

…ring/MDOcean into anthony-enviro

Updated params to feed CEM data instead of directly calling

c9a10e9

new scenario outputs csv (partially filled with nans)

28355c8

use scatteredInterpolant to better handle missing values from CEM loo…

14cff52

…kup table

update grid_CEM

5c04e1d

Merge branch 'anthony-enviro' of https://github.com/symbiotic-enginee…

8ea2425

…ring/MDOcean into anthony-enviro

next round of CEM results

6c125f4

optimize grid cost instead of old objectives

e530162

Merge branch 'anthony-enviro' of github.com:symbiotic-engineering/MDO…

61d930b

…cean into anthony-enviro

change behavior of CEM interp for nan and unviable designs

8c16ffb

next scenario outputs

c7f1c2f

add power_lim to CEM interpolation

c73fa5d

switch objectives when unviable

2dfa042

extrap zeta and use delta CEM cost if viable. gradient optim runs

da387ee

comment out ssh auth sock since not working

f17a749

Merge branch 'main' into anthony-enviro

9e10a5f

rebeccamccabe commented Sep 2, 2025

View reviewed changes

rebeccamccabe changed the title ~~Anthony enviro~~ CEM and enviro modules as used for umerc Sep 2, 2025

add plot for presentation levelization fig

3a4beb0

		function [delta_CEM_CO2, CEM_wec_capacity, delta_CEM_grid_cost,margin_to_viability] = CEM_lookup_table(zeta, omega_n, capacity_cost, power_lim_frac, B_p, location, params)


		%if zeta == 0.05 && omega_n == 0.4 && strcmp(location,'ISONE')

		in.location, in);


		% old code that doesn't normalize with the number of WECs

Conversation

rebeccamccabe commented Jun 12, 2025

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

codecov bot commented Jun 12, 2025 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

❌ 40 Tests Failed:

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

2 participants

codecov bot commented Jun 12, 2025 •

edited

Loading