-
Notifications
You must be signed in to change notification settings - Fork 1
/
cell_parameters.m
55 lines (43 loc) · 1.85 KB
/
cell_parameters.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
function params = cell_parameters(j)
% Set the parameter values for the cell (which are not model-specific) and
% store them in a structure.
% Set the nominal capacity and rate parameters
Qnom = 2.9; % nominal capacity (Ah)
% Note that Qnom can be different from the theoretical capacity and is only
% used to determine the nominal C-rate (=Qnom).
Crate = Qnom; % nominal 1C current density (A)
% Total time
mn = 60; % one minute (s)
hr = 3600; % one hour (s)
% Define SOC limits
Sd = 0; % SOC in discharged state (non-dim.)
Sc = 1; % SOC in fully charged state (non-dim.)
S0 = Sd; % initial SOC (non-dim.)
% Temperatures
CtoK = 273.15; % Celsius to Kelvin conversion (K)
Tamb = 25+CtoK; % ambient temperature (K)
Tref = 25+CtoK; % reference temperature (K)
% Safety limits
Um = 3*Crate; % maximum current amplitude (A)
Vmax = 4.2; % upper voltage limit (V)
Vcut = 2.5; % lower voltage limit (V)
Vrng = Vmax-Vcut; % voltage range (V)
Tmax = 40; % upper temperature limit (deg. C)
Tmin = 25; % lower temperature limit (deg. C)
Trng = Tmax-Tmin; % temperature range (K)
TtoK = Tmin+CtoK; % temperature conversion factor (K)
%% Open-circuit potential
% Select either two electrode potentials or one OCV function: stored either
% as a Matlab function or as a data file using the name of the file(s)
% including the file format. For example,
% OCP_filename = {'Positive.parquet','Negative.parquet'};
% in this order, or
% OCP_filename = {'Example_OCV.parquet'};
% Set OCP filename(s)
OCP_filename = {'Raj2020_NCA.parquet','Raj2020_Graphite.parquet'};
% OCP_filename = {'LGM50_NMC811','LGM50_GraphiteSiOx'};
% OCP_filename = {'Hu2012_LiNMC'};
%% Compile all parameters into the params structure
vars = setdiff(who,{'params','vars'});
for i=1:length(vars), params.(vars{i}) = eval(vars{i}); end
end