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