This model describes the mass flow rate and pressure difference relation of an orifice in the form
V̇ = C sign(Δp) |Δp|m
where V̇ is the volume flow rate, C > 0 is a flow coefficient Δ p is the pressure drop and m ∈ [0.5, 1] is a flow coefficient. The equation is regularized for |Δp| < Δpt, where Δpt is a parameter. For turbulent flow, set m=1 ⁄ 2 and for laminar flow, set m=1.
The model is used for the interzonal air flow models.
For |Δp| < Δpt, the equation is regularized so that it is twice continuously differentiable in Δp, and that it has an infinite number of continuous derivatives in m and in k.
If m is not a function of time, then a, b, c and d can be pre-computed. In this situation, use AixLib.Airflow.Multizone.BaseClasses.powerLawFixedM, which allows to compute these values outside of this function, for example as parameters of a model.
function powerLaw extends Modelica.Icons.Function; input Real C "Flow coefficient, C = V_flow/ dp^m"; input Modelica.Units.SI.PressureDifference dp(displayUnit = "Pa") "Pressure difference"; input Real m(min = 0.5, max = 1) "Flow exponent, m=0.5 for turbulent, m=1 for laminar"; input Modelica.Units.SI.PressureDifference dp_turbulent(min = 0, displayUnit = "Pa") = 0.001 "Pressure difference where regularization starts"; output Modelica.Units.SI.VolumeFlowRate V_flow "Volume flow rate"; end powerLaw;
C
for volume flow coefficient (C = V_flow/dp^m),
and k
for mass flow coefficient (k = m_flow/dp^m).
This is for consistency with
AixLib.Fluid.BaseClasses.FlowModels.