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