Чтобы просмотреть полный доменный исходный файл, в Командной строке 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 (давление), в Mpa
Через переменную mdot (массовая скорость потока жидкости), в kg/s
Через переменную T (температура), в кельвине
Через переменную Phi (энергетическая скорость потока жидкости), в kW
Параметр 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