Isothermal Liquid Domain
The isothermal liquid domain declaration is shown below.
domain isothermal_liquid
% Isothermal Liquid Domain
% Copyright 2019-2022 The MathWorks, Inc.
parameters
bulk_modulus_model = foundation.enum.bulk_modulus_model.const; % Isothermal bulk modulus model
% 1 - const
% 2 - linear
air_dissolution_model = false; % Model air dissolution
rho_L_atm = {998.21, 'kg/m^3'}; % Liquid density at atmospheric pressure (no entrained air)
beta_L_atm = {2.1791e9, 'Pa' }; % Liquid isothermal bulk modulus at atmospheric pressure (no entrained air)
beta_gain = {6, '1' }; % Isothermal bulk modulus vs. pressure increase gain
nu_atm = {1.0034e-6, 'm^2/s' }; % Kinematic viscosity at atmospheric pressure
air_fraction = {0.005, '1' }; % Volumetric fraction of entrained air in mixture at atmospheric pressure
rho_g_atm = {1.225, 'kg/m^3'}; % Gas (air) density at atmospheric condition
polytropic_index = {1.0, '1' }; % Air polytropic index
p_atm = {0.101325, 'MPa' }; % Atmospheric pressure
p_crit = {3, 'MPa' }; % Pressure at which all entrained air is dissolved
p_min = {1, 'Pa' }; % Minimum valid pressure
properties_range_check = simscape.enum.assert.action.error; % Pressure below minimum valid value
% 0 - none
% 1 - warn
% 2 - error
end
variables
p = {0.1, 'MPa'}; % Pressure
end
variables (Balancing = true)
mdot = {0, 'kg/s'}; % Mass flow rate
end
end
It contains the following variables and parameters:
Across variable p (absolute pressure), in MPa
Through variable mdot (mass flow rate), in kg/s
Enumerated parameter bulk_modulus_model, defining the bulk modulus parametrization, with two values:
0 — Bulk modulus is constant
1 — Bulk modulus is a function of pressure
Enumerated parameter air_dissolution_model, defining the entrained air parametrization, with two values:
0 — Entrained air is constant
1 — Entrained air is a function of pressure
Parameter rho_L_atm, defining the liquid density at atmospheric pressure, with zero entrained air
Parameter beta_L_atm, defining the liquid isothermal bulk modulus at atmospheric pressure, with zero entrained air
Parameter beta_gain, defining the ratio of bulk modulus to pressure increase, for when the bulk modulus is a function of pressure
Parameter nu_atm, defining the kinematic viscosity at atmospheric pressure
Parameter air_fraction, defining the volumetric fraction of air that is entrained at atmospheric pressure
Parameter rho_g_atm, defining the air density at atmospheric condition
Parameter polytropic_index, defining the air polytropic index
Parameter p_atm, defining the atmospheric pressure
Parameter p_crit, defining the pressure at which all entrained air is dissolved, for when the amount of entrained air is a function of pressure
Parameter p_min, defining the minimum valid pressure
To refer to this domain in your custom component declarations, use the following syntax:
foundation.isothermal_liquid.isothermal_liquid