modelPlugFlowULg

Validation against data from Université de Liège

Extends from Modelica.Icons.Example (Icon for runnable examples).

Information

The example contains experimental data from a real district heating network.

This model compares the results with the original Modelica Standard Library pipes.

The pipes' temperatures are not initialized. Therefore, results of outflow temperature before approximately the first 10000 seconds should not be considered.

Test bench schematic

Schematic of test rig at ULg

Calibration

There are some uncertainties about the heat loss coefficient between pipe and surrounding air as well as regarding the heat conductivity of the insulation material. With the given data, the length specific thermal resistance is R = 2.164 ((m K)/W), calculated as follows:

R=((1/(2*pipe.kIns)*log((0.0603+2*pipe.dIns)/(0.0603)))+1/(5*(0.0603+2*pipe.dIns)))/Modelica.Constants.pi

U = 1/R = 0.462 W/(m K)

Parameters

TypeNameDefaultDescription
Modelica.Units.SI.MassFlowRatem_flow_nominal1Nominal mass flow rate, used for regularization near zero flow
Modelica.Units.SI.TemperatureT_start_inpipeDataULg.T_start_in + 273.15Initial temperature at pipe inlet
Modelica.Units.SI.TemperatureT_start_outpipeDataULg.T_start_out + 273.15Initial temperature at pipe outlet
Modelica.Units.SI.SpecificHeatCapacitycp_defaultMedium.specificHeatCapacityCp(state = sta_default)Heat capacity of medium
Medium.ThermodynamicStatesta_defaultMedium.setState_pTX(T = Medium.T_default, p = Medium.p_default, X = Medium.X_default)Default medium state

Components

TypeNameDefaultDescription
Fluid.Sources.MassFlowSource_TWaterCityNetworkMass flow source
Fluid.HeatExchangers.Heater_TBoilerBoiler with adjustable outlet temperature
Fluid.Sources.Boundary_pTSewer1Mass flow sink
Fluid.Sensors.TemperatureTwoPortsenTem_outTemperature sensor
Fluid.Sensors.TemperatureTwoPortsenTem_inTemperature sensor
Modelica.Blocks.Sources.CombiTimeTableDataReaderMeasurement data
Modelica.Blocks.Math.UnitConversions.From_degCToutAmbient temperature in degrees
Buildings.HeatTransfer.Sources.FixedTemperatureTBouFixed boundary condition
Modelica.Blocks.Math.UnitConversions.From_degCTinInput temperature into pipe
Data.PipeDataULg151202pipeDataULg
Modelica.Blocks.Math.GaingainGain to test variations of mass flow rate within measurement uncertainty
PlugFlowPipepipePipe
Fluid.Sensors.EnthalpyFlowRatesenEntOutOutlet enthalpy sensor
Modelica.Blocks.Math.AddheatLossSimDifference between inflowing and outflowing enthalpy streams
Modelica.Blocks.Continuous.IntegratoreneLosIntIntegrate model error
Fluid.Sensors.EnthalpyFlowRatesenEntInInlet enthalpy sensor
Modelica.Blocks.Math.AdddelTTemperature difference between in- and outlet
Modelica.Blocks.Math.ProductheaLosMeaHeat loss from measurement (mflow*cp*DeltaT)
Modelica.Blocks.Math.Gaingain3Specific heat of water
Modelica.Blocks.Math.FeedbackheaLosDiffDifference between simulated and measurement enthalpy flow difference
Modelica.Blocks.Sources.ConstantTambAmbient temperature in Kelvin

Contents

NameDescription
Medium

Revisions

  • April 29, 2026, by Michael Wetter:
    Changed configuration of table to cause the parameters to be evaluated, as this leads to more efficient code.
    This is for IBPSA, #2111.
  • March 7, 2020, by Michael Wetter:
    Replaced measured data from specification in Modelica file to external table, as this reduces the computing time.
    This is for #1289.
  • November 24, 2016 by Bram van der Heijde:
    Add pipe thickness for wall capacity calculation and expand documentation section.
  • April 2, 2016 by Bram van der Heijde:
    Change thermal conductivity and put boundary condition in K.
  • Januar 26, 2016 by Carles Ribas:
    First implementation.