Чтобы просмотреть полный доменный исходный файл, в Командной строке MATLAB®, введите:
open([matlabroot '/toolbox/physmod/simscape/library/m/+foundation/+gas/gas.ssc'])
Сокращенное газовое доменное объявление показывают ниже с промежуточными значениями интерполяционной таблицы, не использованными для удобочитаемости.
domain gas
% Gas Domain
% Copyright 2016 The MathWorks, Inc.
parameters
gas_spec = {1, '1'}; % Gas specification
% 1 - Perfect
% 2 - Semiperfect
% 3 - Real
% Perfect gas properties
R = {0.287, 'kJ/(kg*K)'}; % Specific gas constant
Z = {1, '1' }; % Compressibility factor
T_ref = {293.15, 'K' }; % Reference temperature for gas properties
h_ref = {420, 'kJ/kg' }; % Specific enthalpy at reference temperature
cp_ref = {1, 'kJ/(kg*K)'}; % Specific heat at constant pressure
cv_ref = {0.713, 'kJ/(kg*K)'}; % Specific heat at constant volume
mu_ref = {18, 'uPa*s' }; % Dynamic viscosity
k_ref = {26, 'mW/(m*K)' }; % Thermal conductivity
Pr_ref = {0.692307692307692, '1' }; % Prandtl number
% Semiperfect gas properties
T_TLU1 = {[150:10:200, 250:50:1000, 1500, 2000]', 'K'}; % Temperature vector
log_T_TLU1 = {[
5.01063529409626
5.07517381523383
...
7.60090245954208
], '1'}; % Log temperature vector
h_TLU1 = {[
275.264783730547
285.377054177734
...
2377.14064127409
], 'kJ/kg'}; % Specific enthalpy vector
cp_TLU1 = {[
1.01211492398124
1.01042105529234
...
1.24628356718428
], 'kJ/(kg*K)'}; % Specific heat at constant pressure vector
cv_TLU1 = {[
0.725174216111164
0.723480347422265
...
0.959342859314206
], 'kJ/(kg*K)'}; % Specific heat at constant volume vector
mu_TLU1 = {[
10.3766056544352
10.9908682444892
...
68.0682900809450
], 'uPa*s'}; % Dynamic viscosity vector
k_TLU1 = {[
14.1517155766309
15.0474512994325
...
114.486299090693
], 'mW/(m*K)'}; % Thermal conductivity vector
Pr_TLU1 = {[
0.742123270231960
0.738025627675206
...
0.740982912785154
], '1'}; % Prandtl number vector
a_TLU1 = {[
245.095563145758
253.217606015000
...
863.440849227825
], 'm/s'}; % Speed of sound vector
int_dh_T_TLU1 = {[
0
0.0652630980004620
0.126478959779276
...
2.79681971660776
], 'kJ/(kg*K)'}; % integral of dh/T vector
% Real gas properties
% Default gas property tables for dry air
% Rows of the tables correspond to the temperature vector
% Columns of the tables correspond to the pressure vector
T_TLU2 = {[150:10:200, 250:50:1000, 1500, 2000]', 'K' }; % Temperature vector
p_TLU2 = {[0.01:0.01:0.1, 0.12, 0.15, 0.2, 0.5, 1, 2, 5, 10]', 'MPa'}; % Pressure vector
log_T_TLU2 = {[
5.01063529409626
5.07517381523383
...
7.60090245954208
], '1'}; % Log temperature vector
log_p_TLU2 = {[
9.21034037197618
9.90348755253613
...
16.1180956509583
], '1'}; % Log pressure vector
log_rho_TLU2 = {[
-1.45933859209149 -0.765580954956293 ... 2.84006136461620
], '1'}; % Log density table
s_TLU2 = {[
3.85666832168988 3.65733557342939 ... 4.66584072487367
], 'kJ/(kg*K)'}; % Specific entropy table
h_TLU2 = {[
276.007989595737 275.926922934925 ... 2386.79535914098
], 'kJ/kg'}; % Specific enthalpy table
cp_TLU2 = {[
1.00320557010184 1.00416915257750 ... 1.24767439351222
], 'kJ/(kg*K)'}; % Specific heat at constant pressure table
cv_TLU2 = {[
0.715425577953031 0.715655648411093 ... 0.960303115685940
], 'kJ/(kg*K)'}; % Specific heat at constant volume table
mu_TLU2 = {[
10.3604759816291 10.3621937105615 ... 68.3249440282350
], 'uPa*s'}; % Dynamic viscosity table
k_TLU2 = {[
14.0896194596466 14.0962928994967 ... 114.905858092359
], 'mW/(m*K)'}; % Thermal conductivity table
Pr_TLU2 = {[
0.737684026417089 0.738165370950116 ... 0.741888050943969
], '1'}; % Prandtl number table
a_TLU2 = {[
245.567929192228 245.496359667264 ... 878.939999571000
], 'm/s'}; % Speed of sound table
log_drho_dp_TLU2 = {[
-10.6690690699104 -10.6678475863467 ... -13.2956456388403
], '1'}; % Log derivative of density with respect to pressure table
log_drho_dT_TLU2 = {[
-6.46809263806814 -5.77245144865022 ... -4.77782516923660
], '1'}; % Log derivative of density with respect to temperature table
drhou_dp_TLU2 = {[
5.41617782089664 5.42024592837099 ... 3.03095417965253
], '1'}; % Derivative of internal energy per unit volume with respect to pressure table
drhou_dT_TLU2 = {[
-0.195280173069178 -0.391714814336739 ... 1.27305147462835
], 'kJ/(m^3*K)'}; % Derivative of internal energy per unit volume with respect to temperature table
pT_region_flag = {1, '1'}; % Valid pressure-temperature region parameterization
% 1 - Range of gas property tables
% 2 - Specified minimum and maximum values
% 3 - Validity matrix
pT_validity_TLU2 = {ones(24, 18), '1'}; % Pressure-temperature validity matrix
T_min = {150, 'K' }; % Minimum valid temperature
T_max = {2000, 'K' }; % Maximum valid temperature
p_min = {0.01, 'MPa'}; % Minimum valid pressure
p_max = {10, 'MPa'}; % Maximum valid pressure
p_atm = {0.101325, 'MPa'}; % Atmospheric pressure
T_unit = {1, 'K' }; % Unit for log temperature
p_unit = {1, 'Pa' }; % Unit for log pressure
rho_unit = {1, 'kg/m^3' }; % Unit for log density
drho_dp_unit = {1, 'kg/(m^3*Pa)'}; % Unit for log derivative of density with respect to pressure
drho_dT_unit = {1, 'kg/(m^3*K)' }; % Unit for log derivative of density with respect to temperature
log_ZR = {5.65948221575962, '1'}; % Log of compressibility factor times specific gas constant
k_cp = {26e-06, 'kg/(m*s)'}; % Ratio of thermal conductivity to specific heat at typical operating conditions
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 (температура), в K
Через переменную Phi (энергетическая скорость потока жидкости), в kW
Параметр T_min, задавая минимальную допустимую температуру
Параметр T_max, задавая максимальную допустимую температуру
Параметр p_min, задавая минимальное допустимое давление
Параметр p_max, задавая максимальное допустимое давление
Параметр p_atm, задавая атмосферное давление
Параметр gas_spec обеспечивает выбор трех газовых моделей:
1 — Совершенный (значение по умолчанию)
2 — Полусовершенный
3 действительный
В библиотеке Foundation Gas блок Gas Properties (G) служит источником для доменных значений параметров, включая выбор газовой модели. Для получения дополнительной информации о распространении доменных параметров смотрите Работу с Доменными Параметрами.
Доменное объявление также содержит наборы параметров, которые задают газовые свойства для каждой газовой модели.
Свойства для полусовершенного и действительного газа в форме данных об интерполяционной таблице. Эти объявления параметра распространяют к компонентам, соединенным с Газовой областью, и поэтому можно использовать их в tablelookup функция в уравнениях компонента.
Чтобы относиться к этой области в ваших объявлениях собственного компонента, используйте следующий синтаксис:
foundation.gas.gas