主要内容

Butterfly Valve (MA)

R2026b

Butterfly valve in a moist air network

Since R2026b

  • Block icon

Libraries:
Simscape / Fluids / Moist Air / Valves & Orifices / Flow Control Valves

Description

The Butterfly Valve (MA) block models flow through a butterfly valve in a moist air network. A rotating disc inside the valve body controls the flow. When the disc aligns with the flow path, the valve is open. When the disc is perpendicular to the flow path, the valve is closed. The physical signal at port S controls the disc rotation angle.

Valve Parameterizations

The block behavior depends on the Valve parametrization parameter:

  • Cv flow coefficient — The flow coefficient Cv determines the block parameterization. The flow coefficient measures the ease with which the moist air flows when driven by a certain pressure differential.

  • Kv flow coefficient — The flow coefficient Kv, where Kv=0.865Cv, determines the block parameterization. The flow coefficient measures the ease with which the moist air flows when driven by a certain pressure differential.

  • Sonic conductance — The sonic conductance of the resistive element at steady state determines the block parameterization. The sonic conductance measures the ease with which the moist air flows when choked, which is a condition in which the flow velocity is at the local speed of sound. Choking occurs when the ratio between downstream and upstream pressures reaches a critical value known as the critical pressure ratio.

  • Orifice area — The size of the flow restriction determines the block parametrization.

Opening Area

The block calculates the valve opening area during simulation by using the input at port S. The opening area calculations depend on the Opening Characteristic parameter.

Area of Projected Ellipses

If you set Opening Characteristic to Area of projected ellipses, the block calculates the flow area as the area of a circle minus the obstruction from the disc. When the disc rotates, it creates an elliptical obstruction. The physical signal at port S controls the disc rotation. The valve is fully shut when the signal at port S is 0. If the value of the Disc thickness parameter is zero, the valve is fully open when the signal at port S is π/2 rad. Otherwise, the valve is fully open when the signal at port S is a value less than π/2 rad. This value depends on the value of the Disc thickness parameter.

The opening area is

A=α(1cos(φ)),

where:

  • α is the value of the Disc area parameter.

  • φ is the rotation of the disc, specified by the physical signal at port S.

The block normalizes the opening area by the valve maximum area,

Amax=2(R2cos1(1h/R)(Rh)2Rhh2),

where:

  • R is the disc radius.

  • h=R-t2, which is the width of each opening of either side of the disc.

  • t is the value of the Disc thickness parameter.

Tabulated Area

If you set Opening Characteristic to Tabulated, the block interpolates the valve opening from the Orifice area vector, Cv flow coefficient vector, Kv flow coefficient vector, or Sonic conductance vector parameters. The elements in these vectors correspond one-to-one to the elements in the Disc rotation vector parameter. The block interpolates between the data points by using linear interpolation and uses nearest extrapolation for points beyond the table boundaries.

Momentum Balance

The block equations depend on the Valve parametrization parameter. When you set Valve parametrization to Cv flow coefficient, the mass flow rate, m˙, is

m˙=CvSopenSMaxN6Y(pinpout)ρin,

where:

  • Cv is the value of the Maximum Cv flow coefficient parameter.

  • Sopen is the valve opening area.

  • SMax is the maximum valve area when the valve is fully open.

  • N6 is a constant equal to 27.3 for mass flow rate in kg/hr, pressure in bar, and density in kg/m3.

  • Y is the expansion factor.

  • pin is the inlet pressure.

  • pout is the outlet pressure.

  • ρin is the inlet density.

The expansion factor is

Y=1pinpout3pinFγxT,

where:

  • Fγ is the ratio of the isentropic exponent to 1.4.

  • xT is the value of the xT pressure differential ratio factor at choked flow parameter.

The block smoothly transitions to a linearized form of the equation when the pressure ratio, pout/pin, rises above the value of the Laminar flow pressure ratio parameter, Blam,

m˙=CvSopenSMaxN6Ylamρavgpavg(1Blam)(pinpout),

where:

Ylam=11Blam3FγxT.

When the pressure ratio, pout/pin, falls below 1FγxT, the orifice becomes choked and the block switches to the equation

m˙=23CvSopenSMaxN6FγxTpinρin.

When you set Valve parametrization to Kv flow coefficient, the block uses these same equations, but replaces Cv with Kv by using the relation Kv=0.865Cv. For more information on the mass flow equations when the Valve parametrization parameter is Kv flow coefficient or Cv flow coefficient, [2][3].

When you set Valve parametrization to Sonic conductance, the mass flow rate, m˙, is

m˙=CSopenSMaxρrefpinTrefTin[1(poutpinBcrit1Bcrit)2]m,

where:

  • C is the value of the Maximum sonic conductance parameter.

  • Bcrit is the critical pressure ratio.

  • m is the value of the Subsonic index parameter.

  • Tref is the value of the ISO reference temperature parameter.

  • ρref is the value of the ISO reference density parameter.

  • Tin is the inlet temperature.

The block smoothly transitions to a linearized form of the equation when the pressure ratio, pout/pin, rises above the value of the Laminar flow pressure ratio parameter Blam,

m˙=CSopenSMaxρrefTrefTavg[1(BlamBcrit1Bcrit)2]m(pinpout1Blam).

When the pressure ratio, pout/pin, falls below the critical pressure ratio, Bcrit, the orifice becomes choked and the block switches to the equation

m˙=CSopenSMaxρrefpinTrefTin.

Only use the Sonic conductance setting of the Valve parameterization parameter for pneumatic applications. If you use this setting for moist air with high levels of trace gases or are modeling a fluid other than air, you may need to scale the sonic conductance by the square root of the mixture specific gravity. For more information on the mass flow equations when the Valve parametrization parameter is Sonic conductance, see [1].

When you set Valve parametrization to Orifice area based on geometry, the mass flow rate, m˙, is

m˙=CdSopen2γγ1pinρin(poutpin)2γ[1(poutpin)γ1γ1(SopenS)2(poutpin)2γ],

where:

  • Sopen is the valve opening area.

  • S is the value of the Cross-sectional area at ports A and B parameter.

  • Cd is the value of the Discharge coefficient parameter.

  • γ is the isentropic exponent.

The block smoothly transitions to a linearized form of the equation when the pressure ratio, pout/pin, rises above the value of the Laminar flow pressure ratio parameter, Blam,

m˙=CdSopen2γγ1pavg2γγρavgBlam2γ[1Blamγ1γ1(SopenS)2Blam2γ](pinγ1γpoutγ1γ1Blamγ1γ).

When the pressure ratio, pout/pin, falls below(2γ+1)γγ1, the orifice becomes choked and the block switches to the equation

m˙=CdSopen2γγ+1pinρin1(γ+12)2γ1(SopenS)2.

For more information on the mass flow equations when the Valve parametrization parameter is Orifice area based on geometry, see [4].

Mass Balance

The block conserves mass through the valve

m˙A+m˙B=0m˙wA+m˙wB=0m˙gA+m˙gB=0m˙dA+m˙dB=0

where is the mass flow rate and the subscript w denotes water vapor, the subscript g denotes trace gas, and the subscript d denotes water droplets.

Energy Balance

Because the block is an adiabatic component, no heat exchange occurs between the fluid and the wall that surrounds it. No work is done on or by the fluid as it traverses from inlet to outlet. Energy can flow only by advection through ports A and B. By the principle of conservation of energy, the sum of the port energy flows is always equal to zero,

ϕA+ϕB=0,

where ϕ is the energy flow rate into the valve through ports A or B.

Assumptions and Limitations

This block does not model supersonic flow.

Ports

Input

expand all

Physical signal port associated with the position of the valve, in rad. A value of 0 represents a fully shut valve and a value of π/2 represents a fully open valve.

Programmatic Use

Port: S

Conserving

expand all

Moist air conserving port associated with the valve inlet or outlet.

Programmatic Use

Port: A

Moist air conserving port associated with the valve inlet or outlet.

Programmatic Use

Port: B

Parameters

expand all

Method to use to calculate the mass flow rate of the valve:

  • Cv flow coefficient — The flow coefficient Cv determines the block parameterization.

  • Kv flow coefficient — The flow coefficient Kv, where Kv=0.865Cv, determines the block parameterization.

  • Sonic conductance — The sonic conductance of the resistive element at steady state determines the block parameterization.

  • Orifice area — The size of the flow restriction determines the block parametrization.

Programmatic Use

Parameter: valve_parameterization
Values: "fluids.moist_air.valves_orifices.flow_control_valves.enum.OrificeParameterization.Geometry" | "fluids.moist_air.valves_orifices.flow_control_valves.enum.OrificeParameterization.Cv" | "fluids.moist_air.valves_orifices.flow_control_valves.enum.OrificeParameterization.Kv" | "fluids.moist_air.valves_orifices.flow_control_valves.enum.OrificeParameterization.SonicConductance"

Method to use to parameterize the block by using the method of area of projected ellipses or by using tabulated data.

Programmatic Use

Parameter: opening_characteristic
Values: "fluids.moist_air.valves_orifices.flow_control_valves.enum.ButterflyValveSpec.Tabulated" | "fluids.moist_air.valves_orifices.flow_control_valves.enum.ButterflyValveSpec.ProjectedEllipse"

Flow area for a given disc rotation. Each element corresponds one-to-one with the elements in the Disc rotation vector parameter. The first element is the valve leakage and must be nonzero.

Dependencies

To enable this parameter, set Valve parameterization to Orifice area based on geometry and Opening characteristic to Tabulated.

Programmatic Use

Parameter: area_TLU

Value of the Cv flow coefficient when the restriction area available for flow is at a maximum. This parameter measures the ease with which the moist air traverses the resistive element when driven by a pressure differential.

Dependencies

To enable this parameter, set Valve parameterization to Cv flow coefficient and Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: Cv_max

Vector of Cv flow coefficients. Each coefficient corresponds to an element in the Disc rotation vector parameter. This parameter measures the ease with which the moist air traverses the resistive element when driven by a pressure differential. The size of the vector must be the same as the Disc rotation vector parameter.

Dependencies

To enable this parameter, set Valve parameterization to Cv flow coefficient and Opening characteristic to Tabulated.

Programmatic Use

Parameter: Cv_TLU

Disc rotation for a given flow coefficient or area. The elements in this vector must correspond one-to-one with the elements in either the Cv flow coefficient vector, Kv flow coefficient vector, Sonic conductance vector, or Orifice area vector parameters. The elements in this vector must increase monotonically from left to right.

Dependencies

To enable this parameter, set Opening characteristic to Tabulated.

Programmatic Use

Parameter: rotation_Cv_TLU

Parameter: rotation_Kv_TLU

Parameter: rotation_C_TLU

Parameter: rotation_area_TLU

Value of the Kv flow coefficient when the restriction area available for the flow is at a maximum. This parameter measures the ease with which the moist air traverses the resistive element when driven by a pressure differential.

Dependencies

To enable this parameter, set Valve parameterization to Kv flow coefficient and Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: Kv_max

Vector of Kv flow coefficients. Each coefficient corresponds to an element in the Disc rotation vector parameter. This parameter measures the ease with which the moist air traverses the resistive element when driven by a pressure differential. The size of the vector must be the same as the Disc rotation vector parameter.

Dependencies

To enable this parameter, set Valve parameterization to Kv flow coefficient and Opening characteristic to Tabulated.

Programmatic Use

Parameter: Kv_TLU

Ratio between the inlet pressure, pin, and the outlet pressure, pout, defined as (pinpout)/pin where choking first occurs. If you do not have this value, look it up in table 2 in ISA-75.01.01 [3]. Otherwise, use the default value of 0.7.

Dependencies

To enable this parameter, set Valve parameterization to Cv flow coefficient or Kv flow coefficient.

Programmatic Use

Parameter: xT_Cv

Parameter: xT_Kv

Value of the sonic conductance when the cross-sectional area available for flow is at a maximum.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance and Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: C_max

Pressure ratio at which flow first begins to choke and the flow velocity reaches its maximum, given by the local speed of sound. The pressure ratio is the outlet pressure divided by the inlet pressure.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance and Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: B_crit_C

Vector of sonic conductances. Each conductance corresponds to an element in the Disc rotation vector parameter. The size of the vector must be the same as the Disc rotation vector parameter.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance and Opening characteristic to Tabulated.

Programmatic Use

Parameter: C_TLU

Vector of critical pressure ratios at which the flow first chokes. Each critical pressure ratio corresponds to an element in the Disc rotation vector parameter. The critical pressure ratio is the fraction of downstream-to-upstream pressures at which the flow velocity reaches the local speed of sound. The size of the vector must be the same as the Disc rotation vector parameter.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance and Opening characteristic to Tabulated.

Programmatic Use

Parameter: B_crit_TLU

Empirical value to use to calculate the mass flow rate in the subsonic flow regime.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance.

Programmatic Use

Parameter: m

Temperature at standard reference atmosphere, defined as 293.15 K in ISO 8778.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance.

Programmatic Use

Parameter: T_ref_C

Density at standard reference atmosphere, defined as 1.185 kg/m3 in ISO 8778.

Dependencies

To enable this parameter, set Valve parameterization to Sonic conductance.

Programmatic Use

Parameter: rho_ref_C

Surface area of one flat side of the disc.

Dependencies

To enable this parameter, set Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: disc_area

Thickness of valve disc.

Dependencies

To enable this parameter, set Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: disc_thickness

Cross-sectional area of the valve ports.

Programmatic Use

Parameter: area

Correction factor that accounts for discharge losses in theoretical flows.

Dependencies

To enable this parameter, set Valve parameterization to Orifice area based on geometry.

Programmatic Use

Parameter: Cd

Ratio of the flow rate of the orifice when it is closed to when it is open.

Dependencies

To enable this parameter, set Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: leakage_fraction

Continuous smoothing factor that introduces a layer of gradual change to the flow response when the valve is in near-open or near-closed positions. Set this parameter to a nonzero value less than one to increase the stability of your simulation in these regimes.

Dependencies

To enable this parameter, set Opening characteristic to Area of projected ellipses.

Programmatic Use

Parameter: smoothing_factor

Pressure ratio at which the flow transitions between the laminar and turbulent flow regimes. The pressure ratio is the outlet pressure divided by inlet pressure. Typical values range from 0.995 to 0.999.

Programmatic Use

Parameter: B_lam

References

[1] ISO 6358-3. "Pneumatic fluid power – Determination of flow-rate characteristics of components using compressible fluids – Part 3: Method for calculating steady-state flow rate characteristics of systems". 2014.

[2] IEC 60534-2-3. "Industrial-process control valves – Part 2-3: Flow capacity – Test procedures". 2015.

[3] ANSI/ISA-75.01.01. "Industrial-Process Control Valves – Part 2-1: Flow capacity – Sizing equations for fluid flow underinstalled conditions". 2012.

[4] P. Beater. Pneumatic Drives. Springer-Verlag Berlin Heidelberg. 2007.

Extended Capabilities

expand all

C/C++ Code Generation
Generate C and C++ code using Simulink® Coder™.

Version History

Introduced in R2026b