Для просмотра полного исходного файла домена в командной строке MATLAB ® введите :
open([matlabroot '/toolbox/physmod/simscape/library/m/+foundation/+thermal_liquid/thermal_liquid.ssc'])
Сокращенное объявление термического жидкого домена показано ниже, а промежуточные значения таблицы поиска опущены для читаемости.
domain thermal_liquid
% Thermal Liquid Domain
% Copyright 2012-2016 The MathWorks, Inc.
parameters (Size=variable)
% Default liquid property tables for water
% Rows of the tables correspond to the temperature vector
% Columns of the tables correspond to the pressure vector
T_TLU = {[273.1600:10:373.16]', 'K' }; % Temperature vector
p_TLU = {[0.01, 0.1, 5:5:50], 'MPa'}; % Pressure vector
pT_validity_TLU = {[
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1
0 1 1 1 1 1 1 1 1 1 1 1
0 1 1 1 1 1 1 1 1 1 1 1
0 1 1 1 1 1 1 1 1 1 1 1
0 1 1 1 1 1 1 1 1 1 1 1
0 1 1 1 1 1 1 1 1 1 1 1
0 1 1 1 1 1 1 1 1 1 1 1
], '1'}; % Pressure-temperature validity matrix
rho_TLU = {[
999.8 999.8 ... 978.2 980.3
], 'kg/m^3'}; % Density table
u_TLU = {[
0.0002 0.0018 ... 407.1700 405.9800
], 'kJ/kg'}; % Specific internal energy table
nu_TLU = {[
1.7917 1.7914 ... 0.3000 0.3007
], 'mm^2/s'}; % Kinematic viscosity table
cp_TLU = {[
4.2199 4.2194 ... 4.1245 4.1157
], 'kJ/(kg*K)'}; % Specific heat at constant pressure table
k_TLU = {[
561.0400 561.0900 ... 703.3500 706.0000
], 'mW/(m*K)'}; % Thermal conductivity table
beta_TLU = {[
1.9649 1.9654 ... 2.3455 2.3788
], 'GPa'}; % Isothermal bulk modulus table
alpha_TLU = {1e-4 * [
-0.6790 -0.6760 ... 6.8590 6.8000
], '1/K'}; % Isobaric thermal expansion coefficient table
mu_TLU = {[
1.79134166000000 ... 0.294776210000000
], 'cP'}; % Dynamic viscosity table
Pr_TLU = {[
13.4736964762477 ... 1.71842839588810
], '1'}; % Prandtl number table
end
parameters
pT_region_flag = {1, '1' }; % Valid pressure-temperature region parameterization
% 0 - By minimum and maximum value
% 1 - By validity matrix
T_min = {273.16, 'K' }; % Minimum valid temperature
T_max = {373.16, 'K' }; % Maximum valid temperature
p_min = {0.01, 'MPa' }; % Minimum valid pressure
p_max = {50, 'MPa' }; % Maximum valid pressure
p_atm = {0.101325, 'MPa' }; % Atmospheric pressure
k_cv = {1.43e-4, 'kg/(m*s)'}; % Ratio of thermal conductivity to specific heat
max_aspect_ratio = {5, '1' }; % Maximum component aspect ratio (length/diameter) for thermal conduction
end
variables
p = {0.1, 'MPa'}; % Pressure
T = {300, 'K' }; % Temperature
end
variables (Balancing=true)
mdot = {0, 'kg/s'}; % Mass flow rate
Phi = {0, 'kW' }; % Energy flow rate
end
end
Он содержит следующие переменные и параметры:
Поперек переменного p (абсолютное давление), в МПа
Через переменный mdot (массовый расход), в кг/с
Поперек переменной T (температура), в кельвине
Через переменную Phi (расход энергии), в кВт
Параметр pT_region_flag, определяющий допустимую параметризацию области давление-температура, с двумя значениями:
0 - по минимальному и максимальному значению
1 - По матрице достоверности
Параметр T_min, определяющий минимальную допустимую температуру
Параметр p_min, определяющий минимальное допустимое давление
Параметр T_max, определяющий максимальную допустимую температуру
Параметр p_max, определяющий максимальное допустимое давление
Параметр p_atm, определяющий атмосферное давление
Параметр k_cv, определяющий отношение теплопроводности к удельной теплоте
Параметр max_aspect_ratio, определяющий максимальное соотношение сторон компонента (длина/диаметр) для теплопроводности
Он также содержит таблицы поиска, объявленные как параметры домена переменного размера, для следующих термодинамических свойств жидкости:
Плотность
Удельная внутренняя энергия
Кинематическая вязкость
Удельное тепло при постоянном давлении
Теплопроводность
Изотермический объемный модуль
Изобарный коэффициент теплового расширения
Динамическая вязкость
Номер Прандтля
Эти объявления параметров переменного размера распространяются на компоненты, связанные с областью тепловой жидкости, и поэтому их можно использовать в tablelookup функция в уравнениях компонентов. В частности, термические жидкие блоки в библиотеке Foundation используют эти таблицы поиска в целях интерполяции.
Для ссылки на этот домен в объявлениях пользовательских компонентов используйте следующий синтаксис:
foundation.thermal_liquid.thermal_liquid