Двухфазная жидкая область

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