modelVolumeVLE_L4_Advanced

A 1D tube-shaped control volume considering one-phase and two-phase heat transfer in a straight pipe with detailed dynamic momentum and energy balance.

▶ Check in your browser

Extends from ClaRa.Basics.Icons.Volume_L4, ClaRa.Basics.Icons.ComplexityLevel (Displays the complexity level inside model icon ).

Information

For detailed model documentation please consult the html-documentation shipped with ClaRa.

Authorship and Copyright Statement for original (initial) Contribution

Author:

References:

For references please consult the html-documentation shipped with ClaRa.

Remarks:

This component was developed by ClaRa development team under the 3-clause BSD License.

Acknowledgements:

ClaRa originated from the collaborative research projects DYNCAP and DYNSTART. Both research projects were supported by the German Federal Ministry for Economic Affairs and Energy (FKZ 03ET2009 and FKZ 03ET7060).

CLA:

The author(s) have agreed to ClaRa CLA, version 1.0. See https://claralib.com/pdf/CLA.pdf

By agreeing to ClaRa CLA, version 1.0 the author has granted the ClaRa development team a permanent right to use and modify his initial contribution as well as to publish it or its modified versions under the 3-clause BSD License.

The ClaRa development team consists of the following partners:

TLK-Thermo GmbH (Braunschweig, Germany)

XRG Simulation GmbH (Hamburg, Germany).

Parameters

Type	Name	Default	Description
Real	Zeta_in	0	Inlet losses additional to wall friction
Real	Zeta_out	0	Outlet losses additional to wall friction
Units.DensityMassSpecific[geo.N_cv]	rho_nom	TILMedia.VLEFluid.Functions.density_phxi(medium, p_nom, h_nom)	Nominal density
Real	suppFreqCorr	if suppressHighFrequencyOscillations then 1 else 0
Fundamental Definitions
TILMedia.VLEFluid.Types.BaseVLEFluid	medium	simCenter.fluid1	Medium in the component
Boolean	frictionAtInlet	false	True if pressure loss shall be located between first cell and inlet
Boolean	frictionAtOutlet	false	True if pressure loss shall be located between last cell and outlet
Nominal Values
Units.Pressure[geo.N_cv]	p_nom	ones(geo.N_cv)*1e5	Nominal pressure
Units.EnthalpyMassSpecific[geo.N_cv]	h_nom	ones(geo.N_cv)*1e5	Nominal specific enthalpy for single tube
Units.MassFlowRate	m_flow_nom	100	Nominal mass flow for single tube
Units.Pressure	Delta_p_nom	1e4	Nominal pressure loss w.r.t. all parallel tubes
Initialisation
Integer	initOption	1	Type of initialisation
Units.EnthalpyMassSpecific[geo.N_cv]	h_start	ones(geo.N_cv)*800e3	Initial specific enthalpy for single tube
Units.Pressure[geo.N_cv]	p_start	ones(geo.N_cv)*1e5	Initial pressure
Units.MassFlowRate[geo.N_cv + 1]	m_flow_start	ones(geo.N_cv + 1)*100	Initial mass flow rate
Units.MassFraction[medium.nc - 1]	xi_start	zeros(medium.nc - 1)	Initial composition
Initialisation › Model Settings
Boolean	useHomotopy	simCenter.useHomotopy	True, if homotopy method is used during initialisation
Summary and Visualisation
Boolean	showExpertSummary	simCenter.showExpertSummary	True, if a summary shall be shown, else false
Boolean	showData	false	True, if a data port containing p,T,h,s,m_flow shall be shown, else false
Expert Settings
Boolean	suppressHighFrequencyOscillations	false	Suppress oscillations at frequencies greater than inverse travelling time of sound
Real	pressureAdvCalc	0	Parameter for a method of calculation advective pressure drop using:\|\| 0:= upstream velocities \| 1:= velocities in energy cells
Expert Settings › Mass Flow Stabilization
Boolean	useMeanEnthalpyAtInlet	false	Use mean enthalpy at inlet, stabilises zero flows
Boolean	useMeanEnthalpyAtOutlet	false	Use mean enthalpy at inlet, stabilises zero flows
Boolean	advectivePressureLoss	true
Boolean	limitMassChange	false	Set to true to limit time derivative of control volume mass. CAUTION: Precise short time dynamics is artificially changed! Can be useful if simulation stops in case of phase change.
Real	massChangeLimit	10	Limit abs(drhodt[I]/rho[I])=abs(der(mass[I])/mass[I])<= massChangeLimit

Connectors

Type	Name	Description
ClaRa.Basics.Interfaces.FluidPortIn	inlet	Inlet port
ClaRa.Basics.Interfaces.FluidPortOut	outlet	Outlet port
ClaRa.Basics.Interfaces.HeatPort_a[geo.N_cv]	heat

Components

Type	Name	Default	Description
String	complexity (from ComplexityLevel)	"??"
ClaRa.SimCenter	simCenter
Summary	summary
Units.EnthalpyMassSpecific[geo.N_cv]	h		Cell enthalpy
Units.EnthalpyMassSpecific	h_in
Units.EnthalpyMassSpecific	h_out
Units.Mass[geo.N_cv]	mass		Mass of fluid in cells
Real[geo.N_cv]	drhodt
Modelica.Units.SI.MassFraction[geo.N_cv,medium.nc - 1]	xi		Mass fraction
Real	Xi_flow		Mass flow rate of fraction
Modelica.Units.SI.MassFraction[medium.nc - 1]	xi_inlet		Inlet mass fraction of component
Modelica.Units.SI.MassFraction[medium.nc - 1]	xi_outlet		Outlet mass fraction of component
Units.Power[geo.N_cv + 1]	H_flow		Enthalpy flow rate at cell borders
Units.MassFlowRate[geo.N_cv + 1]	m_flow
Units.Velocity[geo.N_cv]	w		flow velocities within cells of energy model == flow velocities across cell borders of flow model
Units.Velocity	w_inlet		flow velocity at inlet
Units.Velocity	w_outlet		flow velocity at outlet
Units.Velocity[geo.N_cv + 1]	w_FM		flow velocities within cells of flow model == flow velocities across cell borders of energy model
Units.Velocity[geo.N_cv]	w_up		upstream flow velocities for calculation of advective pressure losses
TILMedia.VLEFluid.MixtureCompatible.VLEFluid_ph[geo.N_cv]	fluid
TILMedia.VLEFluid.MixtureCompatible.VLEFluid_ph	fluidInlet		if useHomotopy then homotopy(inStream(inlet.h_outflow), h_in) else h_in
TILMedia.VLEFluid.MixtureCompatible.VLEFluid_ph	fluidOutlet		(outlet.h_outflow + inStream(outlet.h_outflow))/2
PressureLoss	pressureLoss		Pressure loss model
HeatTransfer	heatTransfer		heat transfer model
MechanicalEquilibrium	mechanicalEquilibrium		Mechanical equilibrium model
Geometry	geo

Name	Description
Outline
Wall_L4
Summary
PressureLoss
HeatTransfer
Geometry
MechanicalEquilibrium

Revisions

For revisions please consult the html-documentation shipped with ClaRa.