Чтобы просмотреть полный исходный файл области, в MATLAB® Командная строка, введите:
open([matlabroot '/toolbox/physmod/simscape/library/m/+foundation/+two_phase_fluid/two_phase_fluid.ssc'])
Сокращенное объявление области двухфазной жидкости показано ниже, с опущенными для читаемости значениями промежуточной интерполяционной таблицы.
domain two_phase_fluid % Two-Phase Fluid Domain % Copyright 2013-2018 The MathWorks, Inc. parameters p_min = {1e-3, 'MPa' }; % Minimum valid pressure p_max = {95, 'MPa' }; % Maximum valid pressure u_min = {0, 'kJ/kg'}; % Minimum valid specific internal energy u_max = {4000, 'kJ/kg'}; % Maximum valid specific internal energy p_crit = {22.0640, 'MPa' }; % Critical pressure p_atm = {0.101325, 'MPa' }; % Atmospheric pressure q_rev = {0.01, 'Pa' }; % Dynamic pressure threshold for flow reversal transition_range = {0.1, '1'}; % Range in vapor quality for smoothing of density partial derivatives % Default fluid property tables for water % Rows of the tables correspond to the normalized internal energy vector % Columns of the tables correspond to the pressure vector unorm_TLU = {[ -1 -0.958333333333333 ... 1.95833333333333 2 ], '1'}; % Normalized internal energy vector unorm_liq_TLU = {[ -1 -0.958333333333333 ... -0.0416666666666666 0 ], '1'}; % Normalized liquid internal energy vector unorm_vap_TLU = {[ 1 1.04166666666667 ... 1.95833333333333 2 ], '1'}; % Normalized vapor internal energy vector unorm_mix_TLU = {[ 0 0.0204081632653061 ... 0.979591836734694 1 ], '1'}; % Normalized mixture internal energy vector p_TLU = {[ 0.00100000000000000 ... 95.0000000000000 ], 'MPa' }; % Pressure vector v_TLU = {[ 0.00100020732544948 ... 0.00646806895821810 ], 'm^3/kg'}; % Specific volume table s_TLU = {[ -7.05045840376749e-13 ... 6.81117017439253 ], 'kJ/kg/K'}; % Specific entropy table T_TLU = {[ 273.159998291077 ... 1348.81511704704 ], 'K'}; % Temperature table nu_TLU = {[ 1.79172848157620 ... 0.361552941837863 ], 'mm^2/s'}; % Kinematic viscosity table k_TLU = {[ 0.555598831086800 ... 0.250538481070411 ], 'W/(m*K)'}; % Thermal conductivity table Pr_TLU = {[ 13.6057970740934 ... 0.671580102714643 ], '1' }; % Prandtl number table DrhoDp_liq_TLU = {[ 0.509755264394891 ... 1.80807388692342 ], 'kg/m^3/MPa'}; % Partial derivative of liquid density with pressure % at constant specific internal energy table DrhoDp_vap_TLU = {[ 7.36831747174794 ... 1.45923555583730 ], 'kg/m^3/MPa'}; % Partial derivative of vapor density with pressure % at constant specific internal energy table DrhoDp_mix_TLU = {[ 1350.39031864440 ... 1.80807388692342 ], 'kg/m^3/MPa'}; % Partial derivative of mixture density with pressure % at constant specific internal energy table DrhoDu_liq_TLU = {[ 0.0160937306771423 ... -0.350392782066966 ], 'kg^2/m^3/kJ'}; % Partial derivative of liquid density with specific internal energy % at constant pressure table DrhoDu_vap_TLU = {[ -1.79894990819782e-05 ... -0.0610611856834271 ], 'kg^2/m^3/kJ'}; % Partial derivative of vapor density with specific internal energy % at constant pressure table DrhoDu_mix_TLU = {[ -0.0137995431136761 ... -0.350392782066966 ], 'kg^2/m^3/kJ'}; % Partial derivative of mixture density with specific internal energy % at constant pressure table u_sat_liq_TLU = {[ 29.2976389515495 ... 2242.74587564859 ], 'kJ/kg'}; % Saturated liquid specific internal energy vector u_sat_vap_TLU = {[ 2384.48887142334 ... 2242.74587564859 ], 'kJ/kg'}; % Saturated vapor specific internal energy vector end variables p = { 0.1, 'MPa' }; % Pressure u = { 1500, 'kJ/kg' }; % Specific internal energy end variables(Balancing = true) mdot = { 0, 'kg/s' }; % Mass flow rate Phi = { 0, 'kW' }; % Energy flow rate end end
Объявление области содержит следующие переменные и параметры:
Через переменную p (абсолютное давление), в МПа
Через переменную mdot (массовый расход жидкости), в кг/с
Через переменные u (удельная внутренняя энергия), кДж/кг
Через переменную Phi (энергетическая скорость потока жидкости), в кВт
Параметр p_min, определяющий минимально допустимое давление
Параметр p_max, определяющий максимально допустимое давление
Параметр u_min, определяющий минимально допустимую удельную внутреннюю энергию
Параметр u_max, определяющий максимально допустимую удельную внутреннюю энергию
Параметрический p_crit, определяющий критическое давление
Параметр p_atm, определяющий атмосферное давление
Параметр q_rev, определяющий динамический порог давления для обращения потока
Он также содержит данные интерполяционной таблицы, объявленные в качестве параметров области, для следующих свойств жидкости:
Нормированный вектор внутренней энергии
Нормированный вектор внутренней энергии жидкости
Нормированный вектор внутренней энергии пара
Нормированный внутренний энергетический вектор смеси
Вектор давления
Конкретная таблица объемов
Специфическая таблица энтропии
Таблица температур
Таблица кинематической вязкости
Таблица теплопроводности
Таблица чисел Прандтля
Частная производная плотности жидкости с давлением при постоянной удельной внутренней энергетической таблице
Частная производная плотности паров с давлением при постоянной удельной внутренней энергетической таблице
Частное производное плотности смеси с давлением при постоянной удельной внутренней энергетической таблице
Частная производная плотности жидкости с удельной внутренней энергией по постоянному давлению таблице
Частная производная плотности паров с удельной внутренней энергией при таблице постоянного давления
Частное производное плотности смеси с удельной внутренней энергией за постоянное давление таблицу
Насыщенный жидкостно-удельный вектор внутренней энергии
Вектор внутренней энергии, специфичный для насыщенного пара
Эти объявления параметров распространяются на компоненты, соединенные с областью Двухфазной Жидкости, и поэтому вы можете использовать их в tablelookup
функция в уравнениях компонента. В частности, двухфазные блоки жидкости в библиотеке Foundation используют эти интерполяционные таблицы в целях интерполяции.
Для обращения к этой области в объявлениях собственного компонента используйте следующий синтаксис:
foundation.two_phase_fluid.two_phase_fluid