Тепловая жидкая область

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