Commit R2023b

Author: jjha

Components/Controller/Controller.slx

@@ -0,0 +1,28 @@
+component (Propagation = blocks) compressorPower
+% This block computes theoretical compressor power.  
+% This block Computes theoretical compressor power. Carnot coefficient of performance is computed from the input evaporator and condenser temperature.
+% Ideal compressor power is obtained by dividing condenser heat Qcond with Carnot coefficient of performance.
+% The output compressor power is obtained by dividing ideal compressor with system efficiency.
+
+% Copyright 2024 The MathWorks, Inc.
+
+inputs
+    compOn = {1,'1'};
+    Tevap = {0, 'K'};
+    Qcond = {0, 'W'};
+    Tcond = {0, 'K'};
+end
+outputs
+    compressPow = {0, 'W'};
+end
+parameters
+    systemEfficiency = {0.7, '1'}; % Overall system efficiency
+end
+intermediates (ExternalAccess = none)
+    denomAdd =  {((Tcond - Tevap)==0)*eps,'K'}
+    carnotCOP = Tcond/((Tcond - Tevap) + denomAdd); % Carnot CoP for an Ideal refrigeration cycle
+end
+equations
+    compressPow == compOn*(Qcond*(Qcond>0))/(systemEfficiency*carnotCOP); % Compressor power
+end
+end

Components/Controller/Controller_HVAC.slx

@@ -0,0 +1,67 @@
+%% Parameters for BEV model Thermal system - run impirical refrigeration test harness 
+
+% Copyright 2024 The MathWorks, Inc.
+
+%% Initial conditions
+
+vehicleThermal.cabin_p_init = 0.101325; % [MPa] Initial air pressure
+vehicleThermal.cabin_RH_init = 0.4; % Initial relative humidity
+vehicleThermal.cabin_CO2_init = 4e-4; % Initial CO2 mole fraction
+vehicleThermal.cabin_T_init = 298.15; % [K] Initial Temperature
+vehicleThermal.cabin_battVoltage = 100; % [V] Battery Voltage
+vehicleThermal.cabin_battVoltage = 1; % Battery SOC
+vehicleThermal.coolant_p_init = 0.101325; % [MPa] Initial coolant pressure
+vehicleThermal.ambient = 298.15; % [K] Ambient temperature 
+vehicleThermal.CabinSpTp = 293.15 ; % [K] Cabin setpoint temperature
+vehicleThermal.AConoff = 1; % Set AC compressor on/off
+%% Vehicle Cabin
+
+vehicleThermal.cabin_duct_area = 0.04; % [m^2] Air duct cross-sectional area
+vehicleThermal.cabin_duct_volume= 0.0015; % [m^3] Aier duct volume
+vehicleThermal.cabin_perPersonCO2 = 0.01; % [g/s] Adult CO2 exhale rate 
+vehicleThermal.cabin_perPersonMoisture = 0.04; % [g/s] Adult moisture exhale rate
+vehicleThermal.cabin_perPersonHeat = 70; % Average heat tranfer from human body
+vehicleThermal.cabin_exhaleTemperature = 30; % [degC] Adult exhale air temperature
+vehicleThermal.cabin_numberPassanger = 1; % Number of onboard passangers
+vehicleThermal.cabin_airRecirculation = 0.5; % Cabin air recirculation 
+
+%% Liquid Coolant System
+
+vehicleThermal.coolant_pipe_D = 0.019; % [m] Coolant pipe diameter
+vehicleThermal.coolant_channel_D = 0.0092; % [m] Coolant jacket channels diameter
+vehicleThermal.coolant_valve_displacement = 0.0063; % [m] Max spool displacement
+vehicleThermal.coolant_valve_offset = 0.001; % [m] Orifice opening offset when spool is neutral
+vehicleThermal.coolant_valve_D_ratio_max = 0.95; % Max orifice diameter to pipe diameter ratio
+vehicleThermal.coolant_valve_D_ratio_min = 1e-3; % Leakage orifice diameter to pipe diameter ratio
+vehicleThermal.pump_displacement = 0.02; % [l/rev] Coolant pump volumetric displacement
+vehicleThermal.pump_speed_max = 1000; % [rpm] Coolant pump max shaft speed
+vehicleThermal.coolant_tank_volume = 2.5 / 2; % [l] Volume of each coolant tank
+vehicleThermal.coolant_tank_area = 0.11^2; % [m^2] Area of one side of coolant tank
+%% Refrigeration System
+vehicleThermal.refrigeration.evapTubeMass = 1; % [kg]
+
+%% Radiator
+
+vehicleThermal.radiator_L = 0.6; % [m] Overall radiator length
+vehicleThermal.radiator_W = 0.015; % [m] Overall radiator width
+vehicleThermal.radiator_H = 0.2; % [m] Overal radiator height
+vehicleThermal.radiator_N_tubes = 25; % Number of coolant tubes
+vehicleThermal.radiator_tube_H = 0.0015; % [m] Height of each coolant tube
+vehicleThermal.radiator_fin_spacing = 0.002; % Fin spacing
+vehicleThermal.radiator_wall_thickness = 1e-4; % [m] Material thickness
+vehicleThermal.radiator_wall_conductivity = 240; % [W/m/K] Material thermal conductivity
+
+vehicleThermal.radiator_gap_H = (vehicleThermal.radiator_H - vehicleThermal.radiator_N_tubes*vehicleThermal.radiator_tube_H) / (vehicleThermal.radiator_N_tubes - 1); % [m] Height between coolant tubes
+vehicleThermal.radiator_air_area_flow = (vehicleThermal.radiator_N_tubes - 1) * vehicleThermal.radiator_L * vehicleThermal.radiator_gap_H; % [m^2] Air flow cross-sectional area
+vehicleThermal.radiator_air_area_primary = 2 * (vehicleThermal.radiator_N_tubes - 1) * vehicleThermal.radiator_W * (vehicleThermal.radiator_L + vehicleThermal.radiator_gap_H); % [m^2] Primary air heat transfer surface area
+vehicleThermal.radiator_N_fins = (vehicleThermal.radiator_N_tubes - 1) * vehicleThermal.radiator_L / vehicleThermal.radiator_fin_spacing; % Total number of fins
+vehicleThermal.radiator_air_area_fins = 2 * vehicleThermal.radiator_N_fins * vehicleThermal.radiator_W * vehicleThermal.radiator_gap_H; % [m^2] Total fin surface area
+vehicleThermal.radiator_tube_Leq = 2*(vehicleThermal.radiator_H + 20*vehicleThermal.radiator_tube_H*vehicleThermal.radiator_N_tubes); % [m] Additional equivalent tube length for losses due to manifold and splits
+
+vehicleThermal.fan_area = 0.25 * 2; % [m^2] Fan flow area (2 fans)
+
+%% Workspace variables
+
+
+
+

Components/HVAC/HVACimpericalRef.slx Components/HVAC/HVACEmpiricalRef.slxComponents/HVAC/HVACimpericalRef.slx renamed to Components/HVAC/HVACEmpiricalRef.slx

@@ -0,0 +1,35 @@
+classdef ImpiricalRef_UnitTest_MQC < matlab.unittest.TestCase
+%% Class implementation of unit test
+%
+% These are tests to achieve the Minimum Quality Criteria (MQC).
+% MQC is achieved when all runnables (models, scripts, functions) run
+% without any errors.
+%
+% You can run this test by opening in MATLAB Editor and clicking
+% Run Tests button or Run Current Test button.
+% You can also run this test using test runner (the *_runtests.m script)
+
+
+% Copyright 2024 The MathWorks, Inc.
+
+
+%% Simple tests ... just run runnables
+
+methods (Test)
+
+function MQC(~)
+  close all
+  bdclose all
+  mdl = "HVACEmpericalRefHarnessModel";
+  load_system(mdl)
+  BatteryTestHarnessParam;
+  sim(mdl);
+  close all
+  bdclose all
+end
+
+end
+
+end  % classdef
+
+

Components/HVAC/HVACthermal.slx

@@ -0,0 +1,35 @@
+%% Run unit tests
+% This script runs unit test for battery testbench
+
+% Copyright 2024 The MathWorks, Inc.
+
+RelStr = matlabRelease().Release;
+disp("This is MATLAB " + RelStr + ".")
+
+TopFolder = fullfile(currentProject().RootFolder, "Components", "HVAC");
+
+%% Test Suite & Runner
+
+% suite_1 = matlab.unittest.TestSuite.fromFile( ...
+%   fullfile(TopFolder, "Test", "BatteryHV_UnitTest.m"));
+
+suite1 = matlab.unittest.TestSuite.fromFile( ...
+  fullfile(TopFolder, "Test", "ImpiricalRef_UnitTest_MQC.m"));
+
+suite= [suite1];
+
+runner = matlab.unittest.TestRunner.withTextOutput( ...
+  OutputDetail = matlab.unittest.Verbosity.Detailed );
+
+%% JUnit Style Test Result
+
+TestResultFile = "ImpiricalRef_TestResults_" + RelStr + ".xml";
+
+plugin = matlab.unittest.plugins.XMLPlugin.producingJUnitFormat( ...
+  fullfile(TopFolder, "Test", TestResultFile));
+
+addPlugin(runner, plugin)
+
+%% Run tests
+results = run(runner, suite);
+assertSuccess(results)