modelFlow1D2phDB
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.
This simplified 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, boiling heat transfer is assumed with a heat transfer coefficient equal to gamma_b.
- 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, 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
- J. C. Collier: Convective Boiling and Condensation, 2nd ed.,McGraw Hill, 1981, pp. 146.
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Real | pi (from Flow1DBase) | Modelica.Constants.pi | |
| Integer | N (from Flow1DBase) | 2 | Number of nodes for thermal variables |
| Integer | Nw (from Flow1DBase) | N - 1 | Number of volumes on the wall interface |
| Integer | Nt (from Flow1DBase) | 1 | Number of tubes in parallel |
| SI.Distance | L (from Flow1DBase) | Tube length | |
| SI.Position | H (from Flow1DBase) | 0 | Elevation of outlet over inlet |
| SI.Area | A (from Flow1DBase) | Cross-sectional area (single tube) | |
| SI.Length | omega (from Flow1DBase) | Perimeter of heat transfer surface (single tube) | |
| SI.Length | Dhyd (from Flow1DBase) | omega/pi | Hydraulic Diameter (single tube) |
| Medium.MassFlowRate | wnom (from Flow1DBase) | Nominal mass flowrate (total) | |
| ThermoPower.Choices.Flow1D.FFtypes | FFtype (from Flow1DBase) | ThermoPower.Choices.Flow1D.FFtypes.NoFriction | Friction Factor Type |
| SI.PressureDifference | dpnom (from Flow1DBase) | 0 | Nominal pressure drop (friction term only!) |
| Real | Kfnom (from Flow1DBase) | 0 | Nominal hydraulic resistance coefficient (DP = Kfnom*w^2/rho) |
| Medium.Density | rhonom (from Flow1DBase) | 0 | Nominal inlet density |
| SI.PerUnit | Cfnom (from Flow1DBase) | 0 | Nominal Fanning friction factor |
| SI.PerUnit | e (from Flow1DBase) | 0 | Relative roughness (ratio roughness/diameter) |
| SI.PerUnit | Kfc (from Flow1DBase) | 1 | Friction factor correction coefficient |
| Boolean | DynamicMomentum (from Flow1DBase) | false | Inertial phenomena accounted for |
| ThermoPower.Choices.Flow1D.HCtypes | HydraulicCapacitance (from Flow1DBase) | ThermoPower.Choices.Flow1D.HCtypes.Downstream | Location of the hydraulic capacitance |
| Boolean | avoidInletEnthalpyDerivative (from Flow1DBase) | true | Avoid inlet enthalpy derivative |
| Boolean | allowFlowReversal (from Flow1DBase) | system.allowFlowReversal | = true to allow flow reversal, false restricts to design direction |
| SI.PerUnit | wnf (from Flow1DBase) | 0.02 | Fraction of nominal flow rate at which linear friction equals turbulent friction |
| SI.Acceleration | g (from Flow1DBase) | Modelica.Constants.g_n | |
| SI.PerUnit | dzdx (from Flow1DBase) | H/L | Slope |
| SI.Length | l (from Flow1DBase) | L/(N - 1) | Length of a single volume |
| SI.Volume | V (from Flow1DBase) | Nt*A*L | Total volume (all Nt tubes) |
| SI.Pressure | pzero (from Flow1D2ph) | 10 | Small deltap for calculations |
| SI.Pressure | pc (from Flow1D2ph) | Medium.fluidConstants[1].criticalPressure | |
| SI.SpecificEnthalpy | hzero (from Flow1D2ph) | 1e-3 | Small value for deltah |
| SI.CoefficientOfHeatTransfer | gamma_b | 20000 | Coefficient of heat transfer for boiling flow |
| Real | xCHF | 0.9 | Steam quality corresponding to the Critical Heat Flux |
| Initialisation | |||
| Choices.FluidPhase.FluidPhases | FluidPhaseStart (from Flow1DBase) | Choices.FluidPhase.FluidPhases.Liquid | Fluid phase (only for initialization!) |
| Medium.AbsolutePressure | pstart (from Flow1DBase) | 1e5 | Pressure start value |
| Medium.SpecificEnthalpy | hstartin (from Flow1DBase) | if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Liquid then 1e5 else if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Steam then 3e6 else 1e6 | Inlet enthalpy start value |
| Medium.SpecificEnthalpy | hstartout (from Flow1DBase) | if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Liquid then 1e5 else if FluidPhaseStart == Choices.FluidPhase.FluidPhases.Steam then 3e6 else 1e6 | Outlet enthalpy start value |
| Medium.SpecificEnthalpy[N] | hstart (from Flow1DBase) | linspace(hstartin, hstartout, N) | Start value of enthalpy vector (initialized by default) |
| Choices.Init.Options | initOpt (from Flow1DBase) | system.initOpt | Initialisation option |
| Boolean | noInitialPressure (from Flow1DBase) | false | Remove initial equation on pressure |
Connectors
| Type | Name | Default | Description |
|---|---|---|---|
| FlangeA | infl (from Flow1DBase) | ||
| FlangeB | outfl (from Flow1DBase) | ||
| ThermoPower.Thermal.DHT | wall (from Flow1D2ph) |
Components
| Type | Name | Default | Description |
|---|---|---|---|
| ThermoPower.System | system (from Flow1DBase) | System wide properties | |
| SI.Power | Q (from Flow1DBase) | Total heat flow through the lateral boundary (all Nt tubes) | |
| SI.Time | Tr (from Flow1DBase) | Residence time | |
| Medium.ThermodynamicState[N] | fluidState (from Flow1D2ph) | Thermodynamic state of the fluid at the nodes | |
| Medium.SaturationProperties | sat (from Flow1D2ph) | Properties of saturated fluid | |
| SI.Length | omega_hyd (from Flow1D2ph) | Wet perimeter (single tube) | |
| SI.Pressure | Dpfric (from Flow1D2ph) | Pressure drop due to friction | |
| SI.Pressure | Dpstat (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 | |
| Real | dwdt (from Flow1D2ph) | Dynamic momentum term | |
| Medium.AbsolutePressure | p (from Flow1D2ph) | Fluid pressure for property calculations | |
| SI.Pressure[N - 1] | dpf (from Flow1D2ph) | Pressure drop due to friction between two nodes | |
| SI.MassFlowRate | w (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.Temperature | Ts (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.SpecificEnthalpy | hl (from Flow1D2ph) | Saturated liquid temperature | |
| Medium.SpecificEnthalpy | hv (from Flow1D2ph) | Saturated vapour temperature | |
| Real[N] | x (from Flow1D2ph) | Steam quality | |
| Medium.Density[N] | rho (from Flow1D2ph) | Fluid density | |
| Units.LiquidDensity | rhol (from Flow1D2ph) | Saturated liquid density | |
| Units.GasDensity | rhov (from Flow1D2ph) | Saturated vapour density | |
| SI.Mass | M (from Flow1D2ph) | Fluid mass | |
| SI.Power[N - 1] | Qj (from Flow1D2ph) | ||
| Medium.ThermodynamicState | bubble | Bubble point state | |
| Medium.ThermodynamicState | dew | Dew point state |
Revisions
- 24 Mar 2005
by Francesco Casella:
FFtypes package and NoFriction option added. - 16 Dec 2004
by Francesco Casella:
Standard medium definition added. - 1 Nov 2004
by Francesco Casella:
Equations revisited. - 24 Sep 2004
by Francesco Casella:
Adapted to Modelica.Media. - 11 Feb 2004
by Francesco Casella:
Computation of the h.t.c. extended to all the thermodynamic conditions of the fluid. - 1 Oct 2003
by Francesco Casella:
First release.