Main Content

Pipe (G)

Rigid conduit for gas flow

  • Pipe (G) block

Libraries:
Simscape / Foundation Library / Gas / Elements

Description

The Pipe (G) block models pipe flow dynamics in a gas network. The block accounts for viscous friction losses and convective heat transfer with the pipe wall. The pipe contains a constant volume of gas. The pressure and temperature evolve based on the compressibility and thermal capacity of this gas volume. Choking occurs when the outlet reaches the sonic condition.

Caution

Gas flow through this block can choke. If a Mass Flow Rate Source (G) block or a Controlled Mass Flow Rate Source (G) block connected to the Pipe (G) block specifies a greater mass flow rate than the possible choked mass flow rate, you get a simulation error. For more information, see Choked Flow.

Mass Balance

Mass conservation relates the mass flow rates to the dynamics of the pressure and temperature of the internal node representing the gas volume:

MpdpIdt+MTdTIdt=m˙A+m˙B

where:

  • Mp is the partial derivative of the mass of the gas volume with respect to pressure at constant temperature and volume.

  • MT is the partial derivative of the mass of the gas volume with respect to temperature at constant pressure and volume.

  • pI is the pressure of the gas volume.

  • TI is the temperature of the gas volume.

  • t is time.

  • m˙A and m˙B are mass flow rates at ports A and B, respectively. Flow rate associated with a port is positive when it flows into the block.

Energy Balance

Energy conservation relates the energy and heat flow rates to the dynamics of the pressure and temperature of the internal node representing the gas volume:

UpdpIdt+UTdTIdt=ΦA+ΦB+QH

where:

  • Up is the partial derivative of the internal energy of the gas volume with respect to pressure at constant temperature and volume.

  • UT is the partial derivative of the internal energy of the gas volume with respect to temperature at constant pressure and volume.

  • ΦA and ΦB are energy flow rates at ports A and B, respectively.

  • QH is heat flow rate at port H.

Partial Derivatives for Perfect and Semiperfect Gas Models

The partial derivatives of the mass M and the internal energy U of the gas volume, with respect to pressure and temperature at constant volume, depend on the gas property model. For perfect and semiperfect gas models, the equations are:

Mp=VρIpIMT=VρITIUp=V(hIZRTI1)UT=VρI(cpIhITI)

where:

  • ρI is the density of the gas volume.

  • V is the volume of gas.

  • hI is the specific enthalpy of the gas volume.

  • Z is the compressibility factor.

  • R is the specific gas constant.

  • cpI is the specific heat at constant pressure of the gas volume.

Partial Derivatives for Real Gas Model

For real gas model, the partial derivatives of the mass M and the internal energy U of the gas volume, with respect to pressure and temperature at constant volume, are:

Mp=VρIβIMT=VρIαIUp=V(ρIhIβITIαI)UT=VρI(cpIhIαI)

where:

  • β is the isothermal bulk modulus of the gas volume.

  • α is the isobaric thermal expansion coefficient of the gas volume.

Momentum Balance

The momentum balance for each half of the pipe models the pressure drop due to momentum flux and viscous friction:

pApI=(m˙AS)2(1ρI1ρA)+ΔpAIpBpI=(m˙BS)2(1ρI1ρB)+ΔpBI

where:

  • p is the gas pressure at port A, port B, or internal node I, as indicated by the subscript.

  • ρ is the density at port A, port B, or internal node I, as indicated by the subscript.

  • S is the cross-sectional area of the pipe.

  • ΔpAI and ΔpBI are pressure losses due to viscous friction.

The heat exchanged with the pipe wall through port H is added to the energy of gas volume represented by the internal node via the energy conservation equation (see Energy Balance). Therefore, the momentum balances for each half of the pipe, between port A and the internal node and between port B and the internal node, are assumed to be adiabatic processes. The adiabatic relations are:

hA+12(m˙AρAS)2=hI+12(m˙AρIS)2hB+12(m˙BρBS)2=hI+12(m˙BρIS)2

where h is the specific enthalpy at port A, port B, or internal node I, as indicated by the subscript.

The pressure losses due to viscous friction, ΔpAI and ΔpBI, depend on the flow regime. The Reynolds numbers for each half of the pipe are defined as:

ReA=|m˙A|DhSμIReB=|m˙B|DhSμI

where:

  • Dh is the hydraulic diameter of the pipe.

  • μI is the dynamic viscosity at internal node.

If the Reynolds number is less than the Laminar flow upper Reynolds number limit parameter value, then the flow is in the laminar flow regime. If the Reynolds number is greater than the Turbulent flow lower Reynolds number limit parameter value, then the flow is in the turbulent flow regime.

In the laminar flow regime, the pressure losses due to viscous friction are:

ΔpAIlam=fshapem˙AμI2ρIDh2SL+Leqv2ΔpBIlam=fshapem˙BμI2ρIDh2SL+Leqv2

where:

  • fshape is the Shape factor for laminar flow viscous friction parameter value.

  • Leqv is the Aggregate equivalent length of local resistances parameter value.

In the turbulent flow regime, the pressure losses due to viscous friction are:

ΔpAItur=fDarcyAm˙A|m˙A|2ρIDhS2L+Leqv2ΔpBItur=fDarcyBm˙B|m˙B|2ρIDhS2L+Leqv2

where fDarcy is the Darcy friction factor at port A or B, as indicated by the subscript.

The Darcy friction factors are computed from the Haaland correlation:

fDarcyA=[1.8log10(6.9ReA+(εrough3.7Dh)1.11)]2fDarcyB=[1.8log10(6.9ReB+(εrough3.7Dh)1.11)]2

where εrough is the Internal surface absolute roughness parameter value.

When the Reynolds number is between the Laminar flow upper Reynolds number limit and the Turbulent flow lower Reynolds number limit parameter values, the flow is in transition between laminar flow and turbulent flow. The pressure losses due to viscous friction during the transition region follow a smooth connection between those in the laminar flow regime and those in the turbulent flow regime.

Convective Heat Transfer

The convective heat transfer equation between the pipe wall and the internal gas volume is:

QH=Qconv+kISsurfDh(THTI)

Ssurf is the pipe surface area, Ssurf = 4SL/Dh. Assuming an exponential temperature distribution along the pipe, the convective heat transfer is

Qconv=|m˙avg|cpavg(THTin)(1exp(hcoeffSsurf|m˙avg|cpavg))

where:

  • Tin is the inlet temperature depending on flow direction.

  • m˙avg=(m˙Am˙B)/2 is the average mass flow rate from port A to port B.

  • cpavg is the specific heat evaluated at the average temperature.

The heat transfer coefficient, hcoeff, depends on the Nusselt number:

hcoeff=NukavgDh

where kavg is the thermal conductivity evaluated at the average temperature. The Nusselt number depends on the flow regime. The Nusselt number in the laminar flow regime is constant and equal to the value of the Nusselt number for laminar flow heat transfer parameter. The Nusselt number in the turbulent flow regime is computed from the Gnielinski correlation:

Nutur=fDarcy8(Reavg1000)Pravg1+12.7fDarcy8(Pravg2/31)

where Pravg is the Prandtl number evaluated at the average temperature. The average Reynolds number is

Reavg=|m˙avg|DhSμavg

where μavg is the dynamic viscosity evaluated at the average temperature. When the average Reynolds number is between the Laminar flow upper Reynolds number limit and the Turbulent flow lower Reynolds number limit parameter values, the Nusselt number follows a smooth transition between the laminar and turbulent Nusselt number values.

Choked Flow

The choked mass flow rates out of the pipe at ports A and B are:

m˙Achoked=ρAaASm˙Bchoked=ρBaBS

where aA and aB is the speed of sound at ports A and B, respectively.

The unchoked pressure at port A or B is the value of the corresponding Across variable at that port:

pAunchoked=A.ppBunchoked=B.p

The choked pressures at ports A and B are obtained by substituting the choked mass flow rates into the momentum balance equations for the pipe:

pAchoked=pI+(m˙AchokedS)2(1ρI1ρA)+ΔpAIchokedpBchoked=pI+(m˙BchokedS)2(1ρI1ρB)+ΔpBIchoked

ΔpAIchoked and ΔpBIchoked are the pressure losses due to viscous friction, assuming that the choking has occurred. They are computed similar to ΔpAI and ΔpBI, with the mass flow rates at ports A and B replaced by the choked mass flow rate values.

Depending on whether choking has occurred, the block assigns either the choked or unchoked pressure value as the actual pressure at the port. Choking can occur at the pipe outlet, but not at the pipe inlet. Therefore, if pAunchokedpI, then port A is an inlet and pA = pAunchoked. If pAunchoked < pI, then port A is an outlet and

pA={pAunchoked,if pAunchokedpAchokedpAchoked,if pAunchoked<pAchoked 

Similarly, if pBunchokedpI, then port B is an inlet and pB = pBunchoked. If pBunchoked < pI, then port B is an outlet and

pB={pBunchoked,if pBunchokedpBchokedpBchoked,if pBunchoked<pBchoked 

Variables

To set the priority and initial target values for the block variables prior to simulation, use the Initial Targets section in the block dialog box or Property Inspector. For more information, see Set Priority and Initial Target for Block Variables and Initial Conditions for Blocks with Finite Gas Volume.

Nominal values provide a way to specify the expected magnitude of a variable in a model. Using system scaling based on nominal values increases the simulation robustness. Nominal values can come from different sources, one of which is the Nominal Values section in the block dialog box or Property Inspector. For more information, see Modify Nominal Values for a Block Variable.

Assumptions and Limitations

  • The pipe wall is perfectly rigid.

  • The flow is fully developed. Friction losses and heat transfer do not include entrance effects.

  • The effect of gravity is negligible.

  • Fluid inertia is negligible.

  • This block does not model supersonic flow.

Examples

Ports

Conserving

expand all

Gas conserving port associated with the inlet or outlet of the pipe. This block has no intrinsic directionality.

Gas conserving port associated with the inlet or outlet of the pipe. This block has no intrinsic directionality.

Thermal conserving port associated with the temperature of the pipe wall. This temperature may differ from the temperature of the gas volume.

Parameters

expand all

Geometry

The length of the pipe along the direction of flow.

The internal area of the pipe normal to the direction of the flow.

Diameter of an equivalent cylindrical pipe with the same cross-sectional area.

Friction and Heat Transfer

The combined length of all local resistances present in the pipe. Local resistances include bends, fittings, armatures, and pipe inlets and outlets. The effect of the local resistances is to increase the effective length of the pipe segment. This length is added to the geometrical pipe length only for friction calculations. The gas volume depends only on the pipe geometrical length, defined by the Pipe length parameter.

Average depth of all surface defects on the internal surface of the pipe, which affects the pressure loss in the turbulent flow regime.

The Reynolds number above which flow begins to transition from laminar to turbulent. This number equals the maximum Reynolds number corresponding to fully developed laminar flow.

The Reynolds number below which flow begins to transition from turbulent to laminar. This number equals to the minimum Reynolds number corresponding to fully developed turbulent flow.

Dimensionless factor that encodes the effect of pipe cross-sectional geometry on the viscous friction losses in the laminar flow regime. Typical values are 64 for a circular cross section, 57 for a square cross section, 62 for a rectangular cross section with an aspect ratio of 2, and 96 for a thin annular cross section [1].

Ratio of convective to conductive heat transfer in the laminar flow regime. Its value depends on the pipe cross-sectional geometry and pipe wall thermal boundary conditions, such as constant temperature or constant heat flux. Typical value is 3.66, for a circular cross section with constant wall temperature [2].

References

[1] White, F. M., Fluid Mechanics. 7th Ed, Section 6.8. McGraw-Hill, 2011.

[2] Cengel, Y. A., Heat and Mass Transfer – A Practical Approach. 3rd Ed, Section 8.5. McGraw-Hill, 2007.

Extended Capabilities

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

Version History

Introduced in R2016b

expand all