Чтобы просмотреть полный доменный исходный файл, в Командной строке 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-2015 The MathWorks, Inc.
parameters
p_min = { 0.01, 'MPa' }; % Minimum valid pressure
p_max = { 10, 'MPa' }; % Maximum valid pressure
u_min = { 83, 'kJ/kg' }; % Minimum valid specific internal energy
u_max = { 3000, 'kJ/kg' }; % Maximum valid specific internal energy
p_atm = { 1, 'atm' }; % Atmospheric pressure
G_min = { 1e-4, 'kg/s' }; % Minimum thermal conductance coefficient (in terms of specific internal energy)
end
parameters (Size = variable)
% Default lookup tables as a function of pressure and normalized
% specific internal energy. Default values are given for water.
unorm_TLU = {[
-1
-0.965517241379310
...
1.965517241379310
2] , '1' }; % Normalized specific internal energy vector
p_TLU = {[
0.0100000000000000
0.0110069417125221
...
9.08517575651687
10 ]', 'MPa' }; % Pressure vector
v_TLU = {[
0.00100179322007424 0.00100179275967088 ... 0.0346304234950199 0.0314287403997878
], 'm^3/kg' }; % Specific volume table
T_TLU = {[
292.932206861359 292.932221312314 ... 745.161185789778 749.863646934846
], 'K' }; % Temperature table
nu_TLU = {[
1.00879736586632 1.00879622386086 ... 1.84426245592388 1.86263886789473
], 'J/(g*K)' }; % Specific heat at constant volume table
k_TLU = {[
0.597572631285823 0.597573251441156 ... 0.0722108618593745 0.0737522416007857
], 'W/(m*K)' }; % Thermal conductivity table
Pr_TLU = {[
7.05143974468479 7.05142233958132 ... 0.996712568607034 1.00323505646151
], '1' }; % Prandtl number table
u_liq = {[
191.795842042090
199.694279536627
...
1354.80955706624
1393.53799592228]', 'kJ/kg' }; % Saturated liquid specific internal energy vector
u_vap = {[
2437.15737300173
2439.67287156956
...
2557.44803624027
2545.19234394635]', 'kJ/kg' }; % Saturated vapor specific internal energy vector
end
variables
p = { 0.101325, 'MPa' }; % Pressure
u = { 83.905793864039, 'kJ/kg' }; % Specific internal energy
end
variables(Balancing = true)
mdot = { 0, 'kg/s' }; % Mass flow rate
Phi = { 0, 'kW' }; % Heat flow rate
end
end
Доменное объявление содержит следующие переменные и параметры:
Через переменную p (давление), в MPa
Через переменную mdot (массовая скорость потока жидкости), в kg/s
Через переменную u (определенная внутренняя энергия), в kJ/kg
Через переменную Phi (уровень теплового потока), в kW
Параметр p_min, задавая минимальное допустимое давление
Параметр p_max, задавая максимальное допустимое давление
Параметр u_min, задавая минимальную допустимую определенную внутреннюю энергию
Параметр u_max, задавая максимальную допустимую определенную внутреннюю энергию
Параметр p_atm, задавая атмосферное давление
Параметр G_min, задавая минимальный тепловой коэффициент проводимости, с точки зрения определенной внутренней энергии
Это также содержит данные об интерполяционной таблице, объявленные как доменные параметры переменного размера, для следующих жидких свойств:
Нормированный определенный внутренний энергетический вектор
Вектор давления
Определенная таблица объема
Температурная таблица
Удельная теплоемкость за постоянной таблицей объема
Таблица теплопроводности
Таблица числа Прандтля
Влажный жидкий определенный внутренний энергетический вектор
Влажный пар определенный внутренний энергетический вектор
Эти объявления параметра переменного размера распространяют к компонентам, соединенным с Двухфазной Жидкой областью, и поэтому можно использовать их в функции tablelookup в уравнениях компонента. В частности, двухфазные жидкие блоки в библиотеке Foundation используют эти интерполяционные таблицы в целях интерполяции.
Чтобы относиться к этой области в ваших объявлениях собственного компонента, используйте следующий синтаксис:
foundation.two_phase_fluid.two_phase_fluid