modelFlow1Dim_MD

1-D fluid flow model (finite volume discretization - real fluid model) based on the exact integration of the mean volume density

Information

Model Flow1Dim_2ph describes the flow of fluid through a discretized one dimentional tube.

Enthalpy and pressure are selected as state variable

In the model the node variables are characterized by the following syntax: _node.

The model is discretized using the collocated grid method: the nodes variable are the state variables and the cell variable are deduced

.

UpWind discretization is considered i.e. h_node[j+1] = h_node[j].

Numerical options

In this tab several options are available to make the model more robust:

  • Mdotconst: assume constant mass flow rate at each node.
  • max_der: if true the density derivative is truncated during phase change
  • filter_dMdt: if true a first order filter is applied to the fast variations of the density with respect to time
  • max_drhodt: it represents the maximum value of the density derivative. It activates when using max_der is set to true
  • TT: it represents the integration time of the first order filter. It activates when filter_dMdt is set to true

The model is characterized by a SummaryClass that provide a quick access to the following variables once the model is simulated:

  • Temperature at the center of each cell
  • Temperature at each node
  • Enthalpy at each node
  • Enthalpy at the center of each cell
  • Density at the center of each cell
  • Massflow at each nodes
  • Vapor quality at each nodes
  • Pressure in the tube

The model is based on:

Casella, F. Object-oriented modelling of two-phase fluid flows by the finite volume method. Proceedings 5th Mathmod Vienna, Austria, Sep 2006, p. 68.

Parameters

TypeNameDefaultDescription
Modelica.SIunits.SpecificEnthalpyhzero1e-3Small value for deltah
RealpiModelica.Constants.pipi-greco
IntegerN10Number of cells
IntegerNt1Number of tubes in parallel
Modelica.SIunits.AreaA16.18Lateral surface of the tube: heat exchange area
Modelica.SIunits.VolumeV0.03781Volume of the tube
Modelica.SIunits.VolumeViV/NVolume of a single cell
Modelica.SIunits.AreaAiA/NLateral surface of a single cell
Modelica.SIunits.MassFlowRateMdotnom0.2588Nominal fluid flow rate
Modelica.SIunits.CoefficientOfHeatTransferU_nom100if HTtype = LiqVap : Heat transfer coefficient, liquid zone
Initialization
Modelica.SIunits.PressurepstartFluid pressure start value
Medium.TemperatureTstart_inletInlet temperature start value
Medium.TemperatureTstart_outletOutlet temperature start value
Medium.SpecificEnthalpy[N + 1]hstartlinspace(Medium.specificEnthalpy_pT(pstart, Tstart_inlet), Medium.specificEnthalpy_pT(pstart, Tstart_outlet), N + 1)Start value of enthalpy vector (initialized by default)
Numerical options
BooleanMdotconstfalseSet to yes to assume constant mass flow rate at each node (easier convergence)
Booleanmax_derfalseSet to yes to limit the density derivative during phase transitions
Booleanfilter_dMdtfalseSet to yes to filter dMdt with a first-order filter
Realmax_drhodt100Maximum value for the density derivative
Modelica.SIunits.TimeTT1Integration time of the first-order filter
Initialization › Intialization options
Booleansteadystatetrueif true, sets the derivative of h (working fluids enthalpy in each cell) to zero during Initialization

Connectors

TypeNameDefaultDescription
ThermoCycle.Interfaces.Fluid.FlangeAInFlow
ThermoCycle.Interfaces.Fluid.FlangeBOutFlow
ThermoCycle.Interfaces.HeatTransfer.ThermalPortWall_int

Components

TypeNameDefaultDescription
Medium.SaturationPropertiessat
Modelica.SIunits.SpecificEnthalpyh_lsaturated liquid enthalpy
Modelica.SIunits.SpecificEnthalpyh_vsaturated vapor enthalpy
Medium.Densityrho_l
Medium.Densityrho_v
Modelica.SIunits.DerDensityByPressuredrldpBubble point density derivative with respect to pressure
Modelica.SIunits.DerDensityByPressuredrvdpDew point density derivative with respect to pressure
Modelica.SIunits.DerEnthalpyByPressuredhldpBubble point enthalpy derivative with respect to pressure
Modelica.SIunits.DerEnthalpyByPressuredhvdpDew point enthalpy derivative with respect to pressure
Medium.ThermodynamicState[N + 1]fluidState
Medium.AbsolutePressurep
Modelica.SIunits.MassFlowRateM_dot_su
Medium.SpecificEnthalpy[N]hFluid specific enthalpy at the center cells
Modelica.SIunits.SpecificEnthalpy[N + 1]h_nodeEnthalpy state variables at each node
Medium.Temperature[N]T
Medium.Temperature[N + 1]T_nodeFluid temperature at the nodes
Modelica.SIunits.Temperature[N]T_wallInternal wall temperature
Medium.Density[N]rhoFluid cell density at center cell
Medium.Density[N + 1]rho_nodeFluid cell density at the nodes
Modelica.SIunits.DerDensityByEnthalpy[N + 1]drdh_nodeDerivative of density by enthalpy at each node
Modelica.SIunits.DerDensityByEnthalpy[N]drdh1Derivative of average density by left enthalpy
Modelica.SIunits.DerDensityByEnthalpy[N]drdh2Derivative of average density by right enthalpy
Modelica.SIunits.DerDensityByPressure[N + 1]drdp_nodeDerivative of density by pressure at each node
Modelica.SIunits.DerDensityByPressure[N]drdpDerivative of density by pressure at center cell
Real[N]dMdtTime derivative of mass in each cell between two nodes
Modelica.SIunits.HeatFlux[N]qdotheat flux in each cell between two nodes
Modelica.SIunits.MassFlowRate[N + 1]Mdot_node
Modelica.SIunits.MassFlowRate[N]Mdot
Real[N + 1]x_nodeVapor quality at each node
Modelica.SIunits.PowerQ_totTotal heat flux exchanged by the thermal port
Modelica.SIunits.MassM_totTotal mass of the fluid in the component
RealAAVariable for the calculation of the mean density derivative
RealBBVariable for the calculation of the mean density derivative
SummaryClassSummary

Contents

NameDescription
Medium
SummaryClass