modelPartialDrum2Phase2Volume

partial two phase drum model

Information


// Info below is obsolete. Kept for reference until it can be updated.

Spray stream is assumed to be saturated liquid.


This model describes the cylindrical drum of a drum boiler, without assuming thermodynamic equilibrium between the liquid and vapour holdups. Connectors are provided for surgePort inlet, steamPort outlet, downcomer outlet, riser inlet, and blowdown outlet.

The model is based on dynamic mass and energy balance equations of the liquid volume and vapour volume inside the drum. Mass and energy tranfer between the two phases is provided by bulk condensation and surface condensation of the vapour phase, and by bulk boiling of the liquid phase. Additional energy transfer can take place at the surface if the steamPort is superheated.

The riser flowrate is separated before entering the drum, at the vapour pressure. The (saturated) liquid fraction goes into the liquid volume; the (wet) vapour fraction goes into the vapour volume, vith a steamPort quality depending on the liquid/vapour density ratio and on the avr parameter.

The enthalpy of the liquid going to the downcomer is computed by assuming that a fraction of the total mass flowrate (afd) comes directly from the surgePort inlet. The pressure at the downcomer connector is equal to the vapour pressure plus the liquid head.

The metal wall dynamics is taken into account, assuming uniform temperature. Heat transfer takes place between the metal wall and the liquid phase, vapour phase, and external atmosphere, the corresponding heat transfer coefficients being gl, gv, and gext.

The drum level is referenced to the centreline.

The start values of drum pressure, liquid specific enthalpy, vapour specific enthalpy, and metal wall temperature can be specified by setting the parameters p_start, h_liquid_start, h_vapor_start, Tmstart

Modelling options

The following options are available to specify the orientation of the cylindrical drum:

  • DrumOrientation = 0: horizontal axis.
  • DrumOrientation = 1: vertical axis.

Parameters

TypeNameDefaultDescription
SI.TemperatureT_liquid_startMedium.temperature_ph(p_start, h_liquid_start)
SI.TemperatureT_vapor_startMedium.temperature_ph(p_start, h_vapor_start)
Geometry
SI.VolumeV_totalTotal volume (liquid + vapor)
Assumptions
BooleanallowFlowReversalsystem.allowFlowReversal= true to allow flow reversal, false restrics to design direction
Assumptions › Dynamics
Types.DynamicsenergyDynamicssystem.energyDynamicsFormulation of energy balance
Types.DynamicsmassDynamicssystem.massDynamicsFormulation of mass balance
DynamicssubstanceDynamicsmassDynamicsFormulation of substance balances
Initialization
TRANSFORM.Units.NonDimVfrac_liquid_start0.5Initial fraction of volume in the liquid phase
SI.Pressurep_startsystem.p_startPressure start value
SI.SpecificEnthalpyh_liquid_startMedium.bubbleEnthalpy(Medium.setSat_p(p_start))Liquid specific enthalpy start value
SI.SpecificEnthalpyh_vapor_startMedium.dewEnthalpy(Medium.setSat_p(p_start))Vapour specific enthalpy start value
Initialization › Start Value: Mass Fractions
Medium.MassFraction[Medium.nX]X_startMedium.X_defaultMass fractions m_i/m

Components

TypeNameDefaultDescription
Modelica.Fluid.SystemsystemSystem properties
SI.PressurepSurface pressure
SI.VolumeV_vaporVolume occupied by vapour
SI.VolumeV_liquidVolume occupied by liquid
TRANSFORM.Units.NonDimVfrac_liquidFractional of total volume occupied by liquid
Medium.ThermodynamicStatestate_liquidmedium_liquid.stateLiquid state
SI.Densityrho_liquidmedium_liquid.dLiquid density
SI.SpecificEnthalpyh_liquidLiquid specific enthalpy
SI.TemperatureT_liquidLiquid temperature
Medium.ThermodynamicStatestate_vapormedium_vapor.stateVapor state
SI.Densityrho_vapormedium_vapor.dVapor density
SI.SpecificEnthalpyh_vaporVapor specific enthalpy
SI.TemperatureT_vaporVapor temperature
Medium.SaturationPropertiessatMedium.setSat_p(p)State vector to compute saturation properties
SI.SpecificEnthalpyh_gsatMedium.dewEnthalpy(sat)Saturated vapor specific enthalpy
SI.SpecificEnthalpyh_fsatMedium.bubbleEnthalpy(sat)Saturated liquid specific enthalpy
ThermoPower.AbsoluteTemperatureTsatsat.TsatSaturation temperature
SI.EnergyU_liquidInternal energy of liquid
SI.Massm_liquidMass of liquid
SI.MassFlowRatemb_flow_liquidMass flows across liquid boundaries
SI.EnthalpyFlowRateHb_flow_liquidEnthalpy flow across liquid boundaries or energy source/sink
SI.HeatFlowRateQb_flow_liquidHeat flow across liquid boundaries or energy source/sink
SI.PowerWb_flow_liquidWork flow across liquid boundaries or source term
SI.EnergyU_vaporInternal energy of vapor
SI.Massm_vaporMass of vapor
SI.MassFlowRatemb_flow_vaporMass flows across vapor boundaries
SI.EnthalpyFlowRateHb_flow_vaporEnthalpy flow across vapor boundaries or energy source/sink
SI.HeatFlowRateQb_flow_vaporHeat flow across vapor boundaries or energy source/sink
SI.PowerWb_flow_vaporWork flow across vapor boundaries or source term
Medium.BasePropertiesmedium_vapor
Medium.BasePropertiesmedium_liquid
SI.MassmXi_liquidSubstance mass
Medium.MassFlowRatembXi_liquid_flowIndependent mass flow rates, source or sink
SI.MassmXi_vaporSubstance mass
Medium.MassFlowRatembXi_vapor_flowIndependent mass flow rates, source or sink

Contents

NameDescription
Medium

Revisions

<htm_liquid>
<ul>
<li><i>30 May 2005</i>
    by <a href="mailto:francesco.casella@polimi.it">Francesco Casella</a>:<br>
       Initialisation support added.</li>
<li><i>16 Dec 2004</i>
    by <a href="mailto:francesco.casella@polimi.it">Francesco Casella</a>:<br>
       Standard medium definition added.</li>
<li><i>5 Jul 2004</i>
    by <a href="mailto:francesco.casella@polimi.it">Francesco Casella</a>:<br>
       Adapted to Modelica.Media.</li>
<li><i>1 Feb 2004</i>
    by <a href="mailto:francesco.casella@polimi.it">Francesco Casella</a>:<br>
       Improved equations for drum geometry.</li>
<li><i>1 Oct 2003</i>
    by <a href="mailto:francesco.casella@polimi.it">Francesco Casella</a>:<br>
       First release.</li>
</ul>