.Annex60.Fluid.MixingVolumes.BaseClasses.PartialMixingVolume

Information

This is a partial model of an instantaneously mixed volume. It is used as the base class for all fluid volumes of the package Annex60.Fluid.MixingVolumes.

Typical use and important parameters

Set the constant sensibleOnly=true if the model that extends or instantiates this model sets mWat_flow = 0.

Set the constant simplify_mWat_flow = true to simplify the equation

  port_a.m_flow + port_b.m_flow = - mWat_flow;

to

  port_a.m_flow + port_b.m_flow = 0;

This causes an error in the mass balance of about 0.5%, but generally leads to simpler equations because the pressure drop equations are then decoupled from the mass exchange in this component.

To increase the numerical robustness of the model, the constant prescribedHeatFlowRate can be set by the user. This constant only has an effect if the model has exactly two fluid ports connected, and if it is used as a steady-state model. Use the following settings:

• Set prescribedHeatFlowRate=true if the only means of heat transfer at the heatPort is a prescribed heat flow rate that is not a function of the temperature difference between the medium and an ambient temperature. Examples include an ideal electrical heater, a pump that rejects heat into the fluid stream, or a chiller that removes heat based on a performance curve. If the heatPort is not connected, then set prescribedHeatFlowRate=true as in this case, heatPort.Q_flow=0.
• Set prescribedHeatFlowRate=false if there is heat flow at the heatPort computed as K * (T-heatPort.T), for some temperature T and some conductance K, which may itself be a function of temperature or mass flow rate.
  If there is a combination of K * (T-heatPort.T) and a prescribed heat flow rate, for example a solar collector that dissipates heat to the ambient and receives heat from the solar radiation, then set prescribedHeatFlowRate=false.

Set the parameter use_C_flow = true to enable an input connector for the trace substance flow rate.

Implementation

If the model is (i) operated in steady-state, (ii) has two fluid ports connected, and (iii) prescribedHeatFlowRate=true or allowFlowReversal=false, then the model uses Annex60.Fluid.Interfaces.StaticTwoPortConservationEquation in order to use the same energy and mass balance implementation as is used as in steady-state component models. In this situation, the functions inStream are used for the two flow directions rather than the function actualStream, which is less efficient. However, the use of inStream has the disadvantage that hOut has to be computed, in Annex60.Fluid.Interfaces.StaticTwoPortConservationEquation, using

if allowFlowReversal then
  hOut = Annex60.Utilities.Math.Functions.regStep(y1=port_b.h_outflow,
                                                  y2=port_a.h_outflow,
                                                  x=port_a.m_flow,
                                                  x_small=m_flow_small/1E3);
else
  hOut = port_b.h_outflow;
end if;

Hence, for allowFlowReversal=true, if hOut were to be used to compute the temperature that drives heat transfer such as by conduction, then the heat transfer would depend on upstream and the downstream temperatures for small mass flow rates. This can give wrong results. Consider for example a mass flow rate that is positive but very close to zero. Suppose the upstream temperature is 20ˆC, the downstream temperature is 10ˆC, and the heat port is connected through a heat conductor to a boundary condition of 20ˆC. Then, hOut = (port_b.h_outflow + port_a.h_outflow)/2 and hence the temperature heatPort.T is 15ˆC. Therefore, heat is added to the component. As the mass flow rate is by assumption very small, the fluid that leaves the component will have a very high temperature, violating the 2nd law. To avoid this situation, if prescribedHeatFlowRate=false, then the model Annex60.Fluid.Interfaces.ConservationEquation is used instead of Annex60.Fluid.Interfaces.StaticTwoPortConservationEquation.

For simple models that uses this model, see Annex60.Fluid.MixingVolumes.

Revisions

• February 19, 2016 by Filip Jorissen:
  Added outputs U, m, mXi, mC for being able to check conservation of quantities. This if or issue 247.
• January 22, 2016 by Michael Wetter:
  Updated model to use the new parameter use_mWat_flow rather than sensibleOnly.
• January 17, 2016, by Michael Wetter:
  Removed protected block masExc as this revision introduces a conditional connector for the moisture flow rate in the energy and mass balance models. This change was done to use the same modeling concept for the moisture input as is used for the trace substance input.
• December 2, 2015, by Filip Jorissen:
  Added conditional input C_flow for handling trace substance insertions.
• July 17, 2015, by Michael Wetter:
  Added constant simplify_mWat_flow to remove dependencies of the pressure drop calculation on the moisture balance.
• July 1, 2015, by Filip Jorissen:
  Set prescribedHeatFlowRate=prescribedHeatflowRate for Annex60.Fluid.Interfaces.StaticTwoPortConservationEquation. This results in equations that are solved more easily. See issue 282 for a discussion.
• June 9, 2015 by Michael Wetter:
  Set start value for heatPort.T and changed type of T to Medium.Temperature rather than Modelica.SIunits.Temperature to avoid an error because of conflicting start values if Buildings.Fluid.Chillers.Carnot_y is translated using pedantic mode in Dymola 2016. This is for #426.
• June 5, 2015, by Michael Wetter:
  Moved assignment of dynBal.U.start from instance dynBal to the actual model implementation. This is required for a pedantic model check in Dymola 2016. It addresses issue 266.
• May 6, 2015, by Michael Wetter:
  Improved documentation and changed the test

   final parameter Boolean useSteadyStateTwoPort=(nPorts == 2) and
   prescribedHeatFlowRate and ...
  

  to

   final parameter Boolean useSteadyStateTwoPort=(nPorts == 2) and
   (prescribedHeatFlowRate or (not allowFlowReversal)) and ...
  

  The reason is that if there is no flow reversal, then Annex60.Fluid.Interfaces.StaticTwoPortConservationEquation computes hOut = port_b.h_outflow;, and hence it is correct to use hOut to compute temperature-driven heat flow, such as by conduction or convection. See also the model documentation.
  This is for issue #412.
• February 5, 2015, by Michael Wetter:
  Changed initalize_p from a parameter to a constant. This is only required in finite volume models of heat exchangers (to avoid consistent but redundant initial conditions) and hence it should be set as a constant.
• October 29, 2014, by Michael Wetter:
  Made assignment of mFactor final, and changed computation of density to use default medium states as are also used to compute the specific heat capacity.
• October 21, 2014, by Filip Jorissen:
  Added parameter mFactor to increase the thermal capacity.
• July 3, 2014, by Michael Wetter:
  Added parameter initialize_p. This is required to enable the coil models to initialize the pressure in the first volume, but not in the downstream volumes. Otherwise, the initial equations will be overdetermined, but consistent. This change was done to avoid a long information message that appears when translating models.
• May 29, 2014, by Michael Wetter:
  Removed undesirable annotation Evaluate=true.
• February 11, 2014 by Michael Wetter:
  Removed Q_flow and added QSen_flow. This was done to clarify what is sensible and total heat flow rate as part of the correction of issue #197.
• October 8, 2013 by Michael Wetter:
  Removed propagation of show_V_flow to instance steBal as it has no longer this parameter.
• September 13, 2013 by Michael Wetter:
  Renamed rho_nominal to rho_start because this quantity is computed using start values and not nominal values.
• April 18, 2013 by Michael Wetter:
  Removed the check of multiple connections to the same element of a fluid port, as this check required the use of the deprecated cardinality function.
• February 7, 2012 by Michael Wetter:
  Revised base classes for conservation equations in Annex60.Fluid.Interfaces.
• September 17, 2011 by Michael Wetter:
  Removed instance medium as this is already used in dynBal. Removing the base properties led to 30% faster computing time for a solar thermal system that contains many fluid volumes.
• September 13, 2011 by Michael Wetter:
  Changed in declaration of medium the parameter assignment preferredMediumStates=true to preferredMediumStates= not (energyDynamics == Modelica.Fluid.Types.Dynamics.SteadyState). Otherwise, for a steady-state model, Dymola 2012 may differentiate the model to obtain T as a state. See ticket Dynasim #13596.
• July 26, 2011 by Michael Wetter:
  Revised model to use new declarations from Annex60.Fluid.Interfaces.LumpedVolumeDeclarations.
• July 14, 2011 by Michael Wetter:
  Added start values for mass and internal energy of dynamic balance model.
• May 25, 2011 by Michael Wetter:
  
  • Changed implementation of balance equation. The new implementation uses a different model if exactly two fluid ports are connected, and in addition, the model is used as a steady-state component. For this model configuration, the same balance equations are used as were used for steady-state component models, i.e., instead of actualStream(...), the inStream(...) formulation is used. This changed required the introduction of a new parameter m_flow_nominal which is used for smoothing in the steady-state balance equations of the model with two fluid ports. This implementation also simplifies the implementation of Annex60.Fluid.MixingVolumes.BaseClasses.PartialMixingVolumeWaterPort, which now uses the same equations as this model.
  • Another revision was the removal of the parameter use_HeatTransfer as there is no noticeable overhead in always having the heatPort connector present.
• July 30, 2010 by Michael Wetter:
  Added nominal value for mC to avoid wrong trajectory when concentration is around 1E-7. See also https://trac.modelica.org/Modelica/ticket/393.
• February 7, 2010 by Michael Wetter:
  Simplified model and its base classes by removing the port data and the vessel area. Eliminated the base class PartialLumpedVessel.
• October 12, 2009 by Michael Wetter:
  Changed base class to Annex60.Fluid.MixingVolumes.BaseClasses.ClosedVolume.