modelEdgedOrifice

Pressure drop due to sharp edged orifice using Modelica.Fluid.Dissipation.PressureLoss.Orifice

Extends from ThermofluidStream.Interfaces.SISOFlow (Base Model with basic flow eqautions for SISO).

Information

This orifice model computes the pressure loss of the fluid depending on the massflow or the massflow depending on a given pressure difference, some medium properties and the geometry of the orifice.

If the resistance coefficient of a real component is known, it is recommended to calibrate the model to this value varying the Darcy friciton parameter.

This component is an adaptation of ThickEdgedOrifice by Modelica to make it compatible with ThermofludiStream library.

The model is usable for both incompressible and compressible calculation up to at least Ma 0.3 at pipe outlet and one phase medium. The best performance is achieved when using steady state or slowly changing boundary conditions. Numerical stability is best by given mass flow rate and one given pressure boundary. When using two pressure boundaries deviations due to inertia have to be accepted. The Model is not valid for hydraulic shock calculation (sudden change of pressure or mass flow rate).

The pipe bend component is using the partial model SISOFlowBend implementing the common flow balances. For the calculation of pressure loss the function dp_thickEdgedOverall_DP by Modelica is implemented and extended by an additional input which allows the manual adjustment of the Darcy friction factor. (This was removed to increase the compatibility with Modelica.Fluid but may be added again) The input records dp_curvedOverall_IN_con & dp_thickEdgedOverall_IN_var are overwritten with the input parameters defining the orifice geometry and fluid properties. For more information on the underlying pressure loss function, click here. To improve the accuracy when compressible media are used, center state fluid properties (mean dynamic viscosity & mean density) are defined and refered to in the pressure loss function.

The following figure, a resistance coefficient charts representing the pressure loss model under common conditions is shown. (Currently not yet available)

[P. Jordan; HTWG Konstanz; 10/23]

HTWG Konstanz

Parameters

TypeNameDefaultDescription
StringinstanceName (from DropOfCommonsPlus)getInstanceName()Instance name
BooleanassumeConstantMaterialPropertiestrue= true, if constant density and dynamic viscosity is assumed (use '= false' e.g. for Ma > 0.3)
Modelica.Fluid.Dissipation.PressureLoss.Orifice.dp_thickEdgedOverall_IN_conIn_conInput record constants
Layout
BooleandisplayInstanceName (from DropOfCommonsPlus)dropOfCommons.displayInstanceNames= true, if instance name is displayed
BooleandisplayParameters (from DropOfCommonsPlus)dropOfCommons.displayParameters= true, if displaying parameters is enabled
Advanced
BooleanconsiderInertance (from SISOFlow)dropOfCommons.considerInertance=true, if transient momentum (inertance) term is considered; disable only for advanced use
Utilities.Units.InertanceL (from SISOFlow)dropOfCommons.LInertance
StateSelectm_flowStateSelect (from SISOFlow)StateSelect.defaultState selection for mass flow rate
Booleanclip_p_out (from SISOFlow)= false, if dr_corr=0 (correction of inertial pressure difference)
Medium.AbsolutePressurep_min (from SISOFlow)dropOfCommons.p_minMinimum steady-state output pressure
Medium.Densityrho_mindropOfCommons.rho_minMinimal input density
Initialization › Mass flow rate
InitializationMethodsinitM_flow (from SISOFlow)ThermofluidStream.Utilities.Types.InitializationMethods.noneInitialization method for mass flow rate
SI.MassFlowRatem_flow_0 (from SISOFlow)0Initial value for mass flow rate
Utilities.Units.MassFlowAccelerationm_acceleration_0 (from SISOFlow)0Initial value for derivative of mass flow rate
Geometry
SI.Lengthd_1Perimeter of pipe
SI.Lengthd_0Perimeter of vena contraction
SI.Lengthl_0Length of vena contraction
Initialization › dp
ThermofluidStream.Utilities.Types.InitializationMethodsinitdpThermofluidStream.Utilities.Types.InitializationMethods.noneInitialization method for dp
SI.Pressuredp_00Initial value for dp
ThermofluidStream.Utilities.Units.MassFlowAccelerationdp_acceleraton_00Initial value for der(dp)
Initialization
SI.Pressuredp_smooth1Start linearisation for decreasing pressure loss
Advanced › Inertance
BooleancomputeLtrue= true, if inertance L is computed from the geometry
ThermofluidStream.Utilities.Units.InertanceL_valuedropOfCommons.LInertance

Connectors

TypeNameDefaultDescription
Inletinlet (from SISOFlow)
Outletoutlet (from SISOFlow)

Components

TypeNameDefaultDescription
SI.MassFlowRatem_flow (from SISOFlow)inlet.m_flowMass flow rate
SI.Pressuredr_corr (from SISOFlow)Correction of inertial pressure difference
SI.Pressuredp (from SISOFlow)Pressure difference
Medium.AbsolutePressurep_in (from SISOFlow)Medium.pressure(inlet.state)Inlet pressure
Medium.SpecificEnthalpyh_in (from SISOFlow)Medium.specificEnthalpy(inlet.state)Inlet specific enthalpy
Medium.MassFraction[Medium.nXi]Xi_in (from SISOFlow)Medium.massFraction(inlet.state)Inlet mass fractions
Medium.AbsolutePressurep_out (from SISOFlow)Outlet pressure
Medium.SpecificEnthalpyh_out (from SISOFlow)Outlet specific enthalpy
Medium.MassFraction[Medium.nXi]Xi_out (from SISOFlow)Outlet mass fractions
Medium.DensityrhoMean density
Medium.DynamicViscositymuMean dynamic viscosity
Modelica.Fluid.Dissipation.PressureLoss.Orifice.dp_thickEdgedOverall_IN_varIn_varInput record variables
SI.ReynoldsNumberRed_0*abs(m_flow)/A_0/muReynoldsnumber at vena contraction