Чтобы просмотреть полный исходный файл области, в 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