modelPipeFlowVLE_L4_Simple

Model for simulating a pipe flow with a VLE-fluid

Extends from ClaRa.Components.VolumesValvesFittings.Pipes.PipeFlowVLE_L4_Simple (A 1D tube-shaped control volume considering one-phase and two-phase heat transfer in a straight pipe with static momentum balance and simple energy balance.), TransiEnt.Basics.Icons.PipeFlow_L4_Simple.

Information

Purpose of model

A simple model of a pipe flow with a VLE-fluid (Vapor liquid equilibrium) based on the three balance equations of energy, momentum and mass. The mass and the energy balance are dynamic, while the momentum balance is considered static. The balance equations are discretized with a finite volume method. The model uses n+2 fluid models from the TIL media to compute the fluid states at the inlet, outlet, and inside the control volumes of the pipe. Moreover, it uses a replaceable pressure loss model and a replaceable heat transfer model. Therefore, heat and pressure losses are considered. Heat losses are transferred through a heat port if the heat port is connected. However, the heat port doesn’t have to be connected for the model to function. The pipe model can be connected to other models via two fluid ports, which pass the pressure, the specific enthalpy, the mass fraction and the mass flow rate as variables. The model allows reverse flow. In the model are instances of iCom and summary, which are not relevant for the physical understanding of the model.

Level of detail, physical effects considered, and physical insight

L4 (defined in the CodingConventions)

  • Heat losses
  • Pressure losses
  • Thermal inertia
  • Variable fluid properties
  • Mass storage due to density changes

Limits of validity

The accuracy of the dynamic effects of the model depends on the number of control volumes. The thermal inertia of the pipe wall is neglected.

Interfaces

Governing Equations

The equations have been simplified for easier description of the model.

V*rho*(dh/dt)=m_in*h_in+m_out*h_out+Q_flow_out-V*(dp/dt) (for N_cv control volumes)

drho/dt=m_in+m_out (for N_cv control volumes)

p_in-p_out=delta_p(for N_cv control volumes)

drho/dt=a*(dh/dt)+b*(dp/dt)+c*(dxi/dt)
 (The density change is calculated with the TIL media fluid model. Last term is alway zero as the mass fraction doesn't change)

T and rho are determined with the fluid models

References

This model is extended from the ClaRa library. Look at the documentation of the ClaRal library for more information.

Version History

Model was documented by Jan Westphal january 2025 (j.westphal@tuhh.de)

Parameters

TypeNameDefaultDescription
Units.DensityMassSpecific[geo.N_cv]rho_nom (from VolumeVLE_L4)TILMedia.VLEFluid.MixtureCompatible.Functions.density_phxi(medium, p_nom, h_nom)Nominal density
Fundamental Definitions
TILMedia.VLEFluid.Types.BaseVLEFluidmedium (from VolumeVLE_L4)simCenter.fluid1Medium in the component
BooleanfrictionAtInlet (from VolumeVLE_L4)falseTrue if pressure loss between first cell and inlet shall be considered
BooleanfrictionAtOutlet (from VolumeVLE_L4)falseTrue if pressure loss between last cell and outlet shall be considered
Nominal Values
Units.Pressure[geo.N_cv]p_nom (from VolumeVLE_L4)ones(geo.N_cv)*1e5Nominal pressure
Units.EnthalpyMassSpecific[geo.N_cv]h_nom (from VolumeVLE_L4)ones(geo.N_cv)*1e5Nominal specific enthalpy for single tube
Units.MassFlowRatem_flow_nom (from VolumeVLE_L4)100Nominal mass flow w.r.t. all parallel tubes
Units.PressureDifferenceDelta_p_nom (from VolumeVLE_L4)1e4Nominal pressure loss w.r.t. all parallel tubes
Initialisation
IntegerinitOption (from VolumeVLE_L4)0Type of initialisation
Units.EnthalpyMassSpecific[geo.N_cv]h_start (from VolumeVLE_L4)ones(geo.N_cv)*800e3Initial specific enthalpy for single tube
Units.Pressure[geo.N_cv]p_start (from VolumeVLE_L4)ones(geo.N_cv)*1e5Initial pressure
Units.MassFraction[medium.nc - 1]xi_start (from VolumeVLE_L4)zeros(medium.nc - 1)Initial composition
Initialisation › Model Settings
BooleanuseHomotopy (from VolumeVLE_L4)simCenter.useHomotopytrue, if homotopy method is used during initialisation
Summary and Visualisation
BooleanshowExpertSummary (from VolumeVLE_L4)simCenter.showExpertSummaryTrue, if an extended summary shall be shown, else false
BooleanshowData (from VolumeVLE_L4)falseTrue, if a data port containing p,T,h,s,m_flow shall be shown, else false
BooleancontributeToCycleSummary (from PipeFlowVLE_L4_Simple)simCenter.contributeToCycleSummaryTrue if component shall contribute to automatic efficiency calculation
BooleanheatFlowIsLoss (from PipeFlowVLE_L4_Simple)trueTrue if negative heat flow is a loss (not a process product)
Geometry
Basics.Units.Lengthlength (from PipeFlowVLE_L4_Simple)1Length of the pipe (one pass)
Basics.Units.Lengthdiameter_i (from PipeFlowVLE_L4_Simple)0.1Inner diameter of the pipe
Basics.Units.Lengthz_in (from PipeFlowVLE_L4_Simple)0.1Height of inlet above ground
Basics.Units.Lengthz_out (from PipeFlowVLE_L4_Simple)0.1Height of outlet above ground
IntegerN_tubes (from PipeFlowVLE_L4_Simple)1Number Of parallel pipes
IntegerN_passes (from PipeFlowVLE_L4_Simple)1Number of passes of the tubes
Integerorientation (from PipeFlowVLE_L4_Simple)0Main orientation of tube bundle (N_passes>1)
Discretisation
IntegerN_cv (from PipeFlowVLE_L4_Simple)3Number of finite volumes
Basics.Units.Length[N_cv]Delta_x (from PipeFlowVLE_L4_Simple)ClaRa.Basics.Functions.GenerateGrid({0}, length*N_passes, N_cv)Discretisation scheme

Connectors

TypeNameDefaultDescription
ClaRa.Basics.Interfaces.FluidPortIninlet (from VolumeVLE_L4)Inlet port
ClaRa.Basics.Interfaces.FluidPortOutoutlet (from VolumeVLE_L4)Outlet port
ClaRa.Basics.Interfaces.HeatPort_a[geo.N_cv]heat (from VolumeVLE_L4)
Basics.Interfaces.EyeOuteye (from PipeFlowVLE_L4_Simple)

Components

TypeNameDefaultDescription
Stringcomplexity (from ComplexityLevel)"??"
ClaRa.SimCentersimCenter (from VolumeVLE_L4)
Summarysummary (from VolumeVLE_L4)
PressureLosspressureLoss (from VolumeVLE_L4)Pressure loss model
HeatTransferheatTransfer (from VolumeVLE_L4)heat transfer model
Geometrygeo (from VolumeVLE_L4)
MechanicalEquilibriummechanicalEquilibrium (from VolumeVLE_L4)Mechanical equilibrium model
ClaRa.Basics.Interfaces.Connected2SimCenterconnected2SimCenter (from PipeFlowVLE_L4_Simple)