modelFlow1D2phChen

1-dimensional fluid flow model for 2-phase boiling flow (finite volumes, 2-phase, Chen correlation)
Diagram of Flow1D2phChen

Extends from Flow1D2ph (1-dimensional fluid flow model for water/steam (finite volumes, 2-phase)).

Information

This model extends Flow1D2ph by computing the distribution of the heat transfer coefficient gamma and making it available through an extended version of the wall connector.

The model can be used for one-phase or two-phase water/steam flow. The heat transfer coefficient is computed according to the following hypotheses:

  • If the fluid is subcooled liquid or superheated steam, Dittus-Boelter's correlation [1] is used.
  • If the fluid is a two-phase mixture with a steam fraction less than the (constant) critical value xCHF, the heat transfer coefficient is computed according to Chen's correlation [2].
  • If the fluid is wet steam with a steam fraction greater than the (constant) critical value xCHF, the heat transfer coefficient is computed according to Dittus-Boelter's correlation [1], by considering only the steam fraction of the mixture.

A smoothing algorithm is applied to the nodes which are in the neighbourhood of a transition boundary between non-boiling and boiling conditions, to avoid non-physical sudden changes of the nodal values of the heat transfer coefficient when the transition boundary passes through a node. The computed values of the heat transfer coefficient are thus a continuous function of the nodal enthalpies and pressures, so that it is not necessary to generate events in order to handle discontinuities

References

  1. J. C. Collier: Convective Boiling and Condensation, 2nd ed.,McGraw Hill, 1981, pp. 146.
  1. J. C. Collier: Convective Boiling and Condensation, 2nd ed.,McGraw Hill, 1981, pp. 215-220.

Revision history:

Parameters

TypeNameDefaultDescription
Realpi (from Flow1DBase)Modelica.Constants.pi
IntegerN (from Flow1DBase)2Number of nodes for thermal variables
IntegerNw (from Flow1DBase)N - 1Number of volumes on the wall interface
IntegerNt (from Flow1DBase)1Number of tubes in parallel
SI.DistanceL (from Flow1DBase)Tube length
SI.PositionH (from Flow1DBase)0Elevation of outlet over inlet
SI.AreaA (from Flow1DBase)Cross-sectional area (single tube)
SI.Lengthomega (from Flow1DBase)Perimeter of heat transfer surface (single tube)
SI.LengthDhyd (from Flow1DBase)omega/piHydraulic Diameter (single tube)
Medium.MassFlowRatewnom (from Flow1DBase)Nominal mass flowrate (total)
ThermoPower.Choices.Flow1D.FFtypesFFtype (from Flow1DBase)ThermoPower.Choices.Flow1D.FFtypes.NoFrictionFriction Factor Type
SI.PressureDifferencedpnom (from Flow1DBase)0Nominal pressure drop (friction term only!)
RealKfnom (from Flow1DBase)0Nominal hydraulic resistance coefficient (DP = Kfnom*w^2/rho)
Medium.Densityrhonom (from Flow1DBase)0Nominal inlet density
SI.PerUnitCfnom (from Flow1DBase)0Nominal Fanning friction factor
SI.PerUnite (from Flow1DBase)0Relative roughness (ratio roughness/diameter)
SI.PerUnitKfc (from Flow1DBase)1Friction factor correction coefficient
BooleanDynamicMomentum (from Flow1DBase)falseInertial phenomena accounted for
ThermoPower.Choices.Flow1D.HCtypesHydraulicCapacitance (from Flow1DBase)ThermoPower.Choices.Flow1D.HCtypes.DownstreamLocation of the hydraulic capacitance
BooleanavoidInletEnthalpyDerivative (from Flow1DBase)trueAvoid inlet enthalpy derivative
BooleanallowFlowReversal (from Flow1DBase)system.allowFlowReversal= true to allow flow reversal, false restricts to design direction
SI.PerUnitwnf (from Flow1DBase)0.02Fraction of nominal flow rate at which linear friction equals turbulent friction
SI.Accelerationg (from Flow1DBase)Modelica.Constants.g_n
SI.PerUnitdzdx (from Flow1DBase)H/LSlope
SI.Lengthl (from Flow1DBase)L/(N - 1)Length of a single volume
SI.VolumeV (from Flow1DBase)Nt*A*LTotal volume (all Nt tubes)
SI.Pressurepzero (from Flow1D2ph)10Small deltap for calculations
SI.Pressurepc (from Flow1D2ph)Medium.fluidConstants[1].criticalPressure
SI.SpecificEnthalpyhzero (from Flow1D2ph)1e-3Small value for deltah
RealxCHF0.9Steam quality corresponding to Critical Heat Flux
Initialisation
Choices.FluidPhase.FluidPhasesFluidPhaseStart (from Flow1DBase)Choices.FluidPhase.FluidPhases.LiquidFluid phase (only for initialization!)
Medium.AbsolutePressurepstart (from Flow1DBase)1e5Pressure start value
Medium.SpecificEnthalpyhstartin (from Flow1DBase)if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Liquid then 1e5 else if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Steam then 3e6 else 1e6Inlet enthalpy start value
Medium.SpecificEnthalpyhstartout (from Flow1DBase)if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Liquid then 1e5 else if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Steam then 3e6 else 1e6Outlet enthalpy start value
Medium.SpecificEnthalpy[N]hstart (from Flow1DBase)linspace(hstartin, hstartout, N)Start value of enthalpy vector (initialized by default)
Choices.Init.OptionsinitOpt (from Flow1DBase)system.initOptInitialisation option
BooleannoInitialPressure (from Flow1DBase)falseRemove initial equation on pressure

Connectors

TypeNameDefaultDescription
FlangeAinfl (from Flow1DBase)
FlangeBoutfl (from Flow1DBase)
ThermoPower.Thermal.DHTwall (from Flow1D2ph)

Components

TypeNameDefaultDescription
ThermoPower.Systemsystem (from Flow1DBase)System wide properties
SI.PowerQ (from Flow1DBase)Total heat flow through the lateral boundary (all Nt tubes)
SI.TimeTr (from Flow1DBase)Residence time
Medium.ThermodynamicState[N]fluidState (from Flow1D2ph)Thermodynamic state of the fluid at the nodes
Medium.SaturationPropertiessat (from Flow1D2ph)Properties of saturated fluid
SI.Lengthomega_hyd (from Flow1D2ph)Wet perimeter (single tube)
SI.PressureDpfric (from Flow1D2ph)Pressure drop due to friction
SI.PressureDpstat (from Flow1D2ph)Pressure drop due to static head
Real[N - 1]Kf (from Flow1D2ph)Friction coefficient
Real[N - 1]Kfl (from Flow1D2ph)Linear friction coefficient
Real[N - 1]Cf (from Flow1D2ph)Fanning friction factor
Realdwdt (from Flow1D2ph)Dynamic momentum term
Medium.AbsolutePressurep (from Flow1D2ph)Fluid pressure for property calculations
SI.Pressure[N - 1]dpf (from Flow1D2ph)Pressure drop due to friction between two nodes
SI.MassFlowRatew (from Flow1D2ph)Mass flowrate (single tube)
SI.MassFlowRate[N - 1]wbar (from Flow1D2ph)
SI.Velocity[N]u (from Flow1D2ph)Fluid velocity
Medium.Temperature[N]T (from Flow1D2ph)Fluid temperature
Medium.TemperatureTs (from Flow1D2ph)Saturated water temperature
Medium.SpecificEnthalpy[N]h (from Flow1D2ph)Fluid specific enthalpy
Medium.SpecificEnthalpy[N - 1]htilde (from Flow1D2ph)Enthalpy state variables
Medium.SpecificEnthalpyhl (from Flow1D2ph)Saturated liquid temperature
Medium.SpecificEnthalpyhv (from Flow1D2ph)Saturated vapour temperature
Real[N]x (from Flow1D2ph)Steam quality
Medium.Density[N]rho (from Flow1D2ph)Fluid density
Units.LiquidDensityrhol (from Flow1D2ph)Saturated liquid density
Units.GasDensityrhov (from Flow1D2ph)Saturated vapour density
SI.MassM (from Flow1D2ph)Fluid mass
SI.Power[N - 1]Qj (from Flow1D2ph)
Medium.ThermodynamicStatebubbleBubble point state
Medium.ThermodynamicStatedewDew point state

Revisions