modelSteamPipe

Detailed thermal advection model with thermal expansion effects using IF97 water vapour

Information

This models shows the mass, energy, and momentum balance equations for 1D flow of steam in a pipe, using the finite volume method. The pressure loss is assumed to be linear with the flow rate for simplicity. The inertial term and the kinetic term are neglected in the momentum balance equations, hence the pressure wave dynamics is not represented. The pipe is adiabatic, with zero energy storage in the walls. The industry-standard IF97 model is used to compute the steam properties.

The boundary conditions at the inlet are prescribed mass flow w_in_pipe and specific enthalpy h_in_pipe. The boundary condition at the outlet is constant pressure. At time = 0 the inlet flow is zero; at time = 1 the inlet flow is changed to 2 kg/s and at time = 10 the inlet specific enthalpy is raised by 5000 J/kg.

Parameters

TypeNameDefaultDescription
IntegerN10Number of nodes
SI.VolumeVtot1Total volume
Realw_nom2Nominal mass flow rate
Realdp_nom1e5Nominal total pressure loss at full flow rate
Realp_nom10e5Nominal pressure at outlet
Realh_nom3.1e6Nominal specific enthalpy
SI.Pressurep_startp_nomInitial value of pressure states
SI.SpecificEnthalpyh_starth_nomInitial value of enthalpy states
Realkfw_nom/(dp_nom/N)Friction coefficient
SI.VolumeVVtot/NTotal volume
SI.Densityrho_avMedium.density_ph(p_nom + dp_nom/2, h_nom)Average density
SI.TimetauVtot*rho_av/w_nomEstimated transport delay at nominal conditions

Components

TypeNameDefaultDescription
SI.MassFlowRatew_in_pipe0 + (if time > 1 then w_nom else 0)Inlet mass flow rate
SI.SpecificEnthalpyh_in_pipeh_nom + (if time > 10 then 5000 else 0)Inlet specific enthalpy
SI.Mass[N]MFluid mass within each volume
SI.Energy[N]EFluid internal energy within each volume
SI.MassFlowRate[N]w_inInlet flow rate of each volume
SI.MassFlowRate[N]w_outOutlet flow rate of each volume
SI.Pressure[N]pPressure states
SI.SpecificEnthalpy[N]hSpecific enthalpy states
SI.SpecificEnthalpy[N]h_inInlet specific enthalpy of each volume
SI.SpecificEnthalpy[N]h_outOutlet specific enthalpy of each volume
SI.Density[N]rhoVolume density
Medium.ThermodynamicState[N]stateVolume thermodynamic state
SI.Pressurep_in_pipeInlet pressure
SI.Pressurep_out_pipeOutlet pressure
SI.MassFlowRatew_out_pipeInlet flow rate
SI.SpecificEnthalpyh_out_pipe Outlet flow rate

Contents

NameDescription
Medium