modelDiffuser

Pressure drop of a conical diffusor

Extends from Interfaces.SISOFlow_nonConstArea (SISOFlow considering dynamic effects due to different inlet/outlet cross sectional areas).

Information

This component models a diffuser. It computes the pressure loss of the fluid depending on the massflow rate or the massflow rate depending on a given pressure difference, some medium properties and the geometry of the diffuser. The component is valid for both incompressible and compressible calculation up to at least Ma 0.3 at pipe outlet and one phase medium.

The pipe bend component is using the partial model SISOFlow_nonConstArea implementing the common flow balances for components with non constant cross sectional areas. For the calculation of pressure loss the function dp_conicalDiffuserOverall_DP is used. It is based on the formulas and data provided in "Handbook of Hydraulic Resistance" by Idel'chik (1960). Be aware that there the book has been updated since then.

The function calculates and feeds back the pressure loss as well as total pressure loss coefficient and the partial pressure loss coefficientes due to diffuser enlargement and wall friction. For more information click here. The functions assumes pressure loss coefficient due to diffuser enlargement to be independent from Reynolds Number, constant wall fricition along diffuser length and uniform velocity profile at inlet.

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 pressure loss of the diffuser in the described in the figure is shown. (Currently not yet available)

[P. Jordan; HTWG Konstanz; 01/24]

HTWG Konstanz

Parameters

TypeNameDefaultDescription
StringinstanceName (from DropOfCommonsPlus)getInstanceName()Instance name
BooleanassumeConstantDensitytrue= true, if incompressibility is assumed (use '= false' for Ma > 0.3)
SI.Anglealphaif not setLength then alpha_par else 2*atan((d_2 - d_1)/2/l_par)Central divergence angle of diffuser walls
SI.Lengthl_difif not setLength then (d_2 - d_1)/2/tan(alpha/2) else l_parDiffuser length
SI.Lengthks_internalif material == ThermofluidStream.Processes.Internal.Material.concrete then 5e-3 elseif material == ThermofluidStream.Processes.Internal.Material.wood then 0.5e-3 elseif material == ThermofluidStream.Processes.Internal.Material.castIron then 0.25e-3 elseif material == ThermofluidStream.Processes.Internal.Material.galvanizedIron then 0.15e-3 elseif material == ThermofluidStream.Processes.Internal.Material.steel then 0.059e-3 elseif material == ThermofluidStream.Processes.Internal.Material.drawnPipe then 0.0015e-3 else ksPipe roughness
Layout
BooleandisplayInstanceName (from DropOfCommonsPlus)dropOfCommons.displayInstanceNames= true, if instance name is displayed
BooleandisplayParameters (from DropOfCommonsPlus)dropOfCommons.displayParameters= true, if displaying parameters is enabled
Advanced
ThermofluidStream.Utilities.Units.InertanceL (from SISOFlow_nonConstArea)dropOfCommons.LInertance of the flow
StateSelectm_flowStateSelect (from SISOFlow_nonConstArea)StateSelect.defaultState select for m_flow
Booleanclip_p_out (from SISOFlow_nonConstArea)If false, set dr_corr to zero
Medium.AbsolutePressurep_min (from SISOFlow_nonConstArea)dropOfCommons.p_minMinimal steady-state output pressure
ThermofluidStream.Utilities.Units.InertanceL_valuedropOfCommons.LInertance of pipe
BooleancomputeLtrueCompute L from r and l
Medium.Densityrho_mindropOfCommons.rho_minMinimal input density
StateSelectdpStateSelectStateSelect.defaultState select for dp
Initialization › Mass flow
InitializationMethodsinitM_flow (from SISOFlow_nonConstArea)ThermofluidStream.Utilities.Types.InitializationMethods.noneInitialization method for m_flow
SI.MassFlowRatem_flow_0 (from SISOFlow_nonConstArea)0Initial value for m_flow
ThermofluidStream.Utilities.Units.MassFlowAccelerationm_acceleraton_0 (from SISOFlow_nonConstArea)0Initial value for der(m_flow)
Geometry
SI.Lengthd_1Inlet diameter
SI.Lengthd_2Outlet diameter
BooleansetLengthtrue= true, if diffusor length l_par is given (= false, if divergence angle alpha_par is given)
SI.Lengthl_par0Diffusor length
SI.Anglealpha_par0Central divergence angle of diffuser walls
Roughness
ThermofluidStream.Processes.Internal.MaterialmaterialThermofluidStream.Processes.Internal.Material.otherMaterial of pipe
SI.LengthksPipe roughness

Connectors

TypeNameDefaultDescription
ThermofluidStream.Interfaces.Inletinlet (from SISOFlow_nonConstArea)
ThermofluidStream.Interfaces.Outletoutlet (from SISOFlow_nonConstArea)

Components

TypeNameDefaultDescription
SI.MassFlowRatem_flow (from SISOFlow_nonConstArea)inlet.m_flowmass flow through component
SI.Velocityc_in (from SISOFlow_nonConstArea)Inlet velocity
SI.Velocityc_out (from SISOFlow_nonConstArea)Outlet velocity
SI.Pressuredr_corr (from SISOFlow_nonConstArea)
SI.Pressuredp (from SISOFlow_nonConstArea)Static pressure difference dp = p_out - p_in
SI.Pressuredq (from SISOFlow_nonConstArea)Dynamic pressure difference dq = q_out - q_in
Medium.DensityrhoMean density
Medium.DynamicViscositymuMean dynamic viscosity
SI.ReynoldsNumberRed_1*abs(m_flow)/A_1/muReynolds Number
Internal.Types.PressureLossCoefficientzeta_difTotal pressure loss coefficient of diffuser
Internal.Types.PressureLossCoefficientzeta_expPressure loss coefficient due to diffuser enlargement
Internal.Types.PressureLossCoefficientzeta_frPressure loss coefficient of diffuser due to wall friction