modelPartialConductionElement

Partial model of quasi-stationary mass and heat transfer

Extends from Interfaces.SISOFlow (Base Model with basic flow eqautions for SISO).

Information

This model is an element with a fixed volume (fig. 1). The mass in the volume is assumed quasi-stationary (statically computed with volume and density), and the fore massflow is coupled to the rear massflow. Because of this the ConductionElement cannot be used as a loop breaker. The advantage is that multiple ConductionElements can be put behind each other without worrying about oscillations or fast eigenvalues between their masses. The ConductionElement implements equations for conservation of mass and energy for the fluid mass contained within it.

Initialization can be done by initial temperature, initial enthalpy or by using the inlet state.

The ConductionElement makes different assumptions:

  • Quasistationary mass:
    m_flow_rear = - m_flow_fore & M = rho * V (this assumption violates the conservation of mass for changing densities, since the mass in the element can change although inflow and outflow are the same)
    der(H) = der(M*h) = M*der(h) (This assumption violates the conservation of energy for changing densities, since then the mass M of fluid in the element is no longer constant)
  • Neglection of der(p) in the energy equation
    V*der(p) = 0 (this assumption violates the conservation of energy for changing pressures. For a noticeable difference in the testcase the der(p) must be in the order of 1e5 Pa/s).
    This assumption can be turned off by setting neglectPressureChanges=false (true by default) in the Advanced tab. This option requires the fore and rear input pressures to be smooth.
  • Due to stability reasons the component exhibits a different behavior for negative massflows (see fig. 2). For negative massflows, the ingoing and outgoing massflows get decoupled from the fluid remaining in the component. The fluid exits with the same state as it enters, and the heatport is connected to the fluid left in the volume (which is decoupled to the mass flowing).

Due to these assumptions minor violations in the global energy conservation can occur. With the flag enforce_global_energy_conservation in the "Advanced" tab is set true (Default: false), long-term energy storage in the ConductionElement is sacrificed to hold global energy conservation.

fig. 1: positive massflow model

fig. 2: negative massflow model

Parameters

TypeNameDefaultDescription
StringinstanceName (from DropOfCommonsPlus)getInstanceName()Instance name
SI.VolumeV0.001Volume
Layout
BooleandisplayInstanceName (from DropOfCommonsPlus)dropOfCommons.displayInstanceNames= true, if instance name is displayed
BooleandisplayParameters (from DropOfCommonsPlus)dropOfCommons.displayParameters= true, if displaying parameters is enabled
Advanced
BooleanconsiderInertance (from SISOFlow)dropOfCommons.considerInertance=true, if transient momentum (inertance) term is considered; disable only for advanced use
Utilities.Units.InertanceL (from SISOFlow)dropOfCommons.LInertance
StateSelectm_flowStateSelect (from SISOFlow)StateSelect.defaultState selection for mass flow rate
Booleanclip_p_out (from SISOFlow)= false, if dr_corr=0 (correction of inertial pressure difference)
Medium.AbsolutePressurep_min (from SISOFlow)dropOfCommons.p_minMinimum steady-state output pressure
Medium.Densityrho_mindropOfCommons.rho_minMinimal density
BooleanneglectPressureChangestrue=true, if pressure changes are neglected
SI.MassFlowRatem_flow_assert-dropOfCommons.m_flow_regAssertion threshold for negative massflows
Initialization › Mass flow rate
InitializationMethodsinitM_flow (from SISOFlow)ThermofluidStream.Utilities.Types.InitializationMethods.noneInitialization method for mass flow rate
SI.MassFlowRatem_flow_0 (from SISOFlow)0Initial value for mass flow rate
Utilities.Units.MassFlowAccelerationm_acceleration_0 (from SISOFlow)0Initial value for derivative of mass flow rate
Initialization › Specific enthalpy
Internal.InitializationMethodsCondElementinitThermofluidStream.Processes.Internal.InitializationMethodsCondElement.inletInitialization for specific enthalpy
Medium.TemperatureT_0Medium.T_defaultInitial Temperature
Medium.SpecificEnthalpyh_0Medium.h_defaultInitial specific enthalpy
Advanced › Global energy conservation
Booleanenforce_global_energy_conservationfalse= true, if global conservation of energy is enforced
Advanced › global energy conservation
SI.TimeT_e100Time constant for global conservation of energy

Connectors

TypeNameDefaultDescription
Inletinlet (from SISOFlow)
Outletoutlet (from SISOFlow)
Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_aheatPort

Components

TypeNameDefaultDescription
SI.MassFlowRatem_flow (from SISOFlow)inlet.m_flowMass flow rate
SI.Pressuredr_corr (from SISOFlow)Correction of inertial pressure difference
SI.Pressuredp (from SISOFlow)Pressure difference
Medium.AbsolutePressurep_in (from SISOFlow)Medium.pressure(inlet.state)Inlet pressure
Medium.SpecificEnthalpyh_in (from SISOFlow)Medium.specificEnthalpy(inlet.state)Inlet specific enthalpy
Medium.MassFraction[Medium.nXi]Xi_in (from SISOFlow)Medium.massFraction(inlet.state)Inlet mass fractions
Medium.AbsolutePressurep_out (from SISOFlow)Outlet pressure
Medium.SpecificEnthalpyh_out (from SISOFlow)Outlet specific enthalpy
Medium.MassFraction[Medium.nXi]Xi_out (from SISOFlow)Outlet mass fractions
Medium.SpecificEnthalpyhVolume? specific enthalpy
Medium.ThermodynamicStatestateMedium.setState_phX(p_in, h, Xi_in)Volume thermodynamic state
Medium.TemperatureTMedium.temperature(state)Volume temperature
SI.ThermalConductancekThermal conductance
SI.EnergydeltaE_systemEnergy difference between m_flow*(h_in-h_out) and Q_flow
SI.MassMMass (of the volume)