Model Representation and Result Quantities

How the heat network is simulated: thermo-hydraulic model, time integration, result quantities and output files

Overview

The dynamic simulation represents the heat network as a thermo-hydraulic model: a quasi-steady hydraulics (mass fluxes and pressures) coupled with transient energy balances of the fluid volumes. The time integration is performed with an adaptive step size. This page describes the model representation compactly and lists the result quantities and output files as a reference.

Model Representation

Hydraulics

In each integration step, the nonlinear equation system of mass conservation at the nodes and the pressure-loss relations of all elements is solved with a Newton method. The reference pressure is prescribed at the pressure-maintenance node; from the geodetic heights of the nodes, absolute pressures additionally result (quantities with the suffix Absolute). Controlled valves and pumps enter the equation system directly with their current control variables.

Thermics

Each pipe is subdivided into fluid volumes along its length; the maximum element length Δx\Delta x is the Length of the pipe discretization from the network settings. The number of volumes per pipe is

n=max ⁣(1, LΔx)n = \max\!\left(1,\ \left\lfloor \tfrac{L}{\Delta x} \right\rfloor\right)

with the pipe length LL. For each volume, an energy balance of heat transport with the mass flux and heat exchange with the surroundings (e.g. ground) is solved. A smaller discretization length increases the accuracy of the temperature fronts in the network, but enlarges the equation system and the computation time.

Controllers and Ground Coupling

Controllers (e.g. for valves and pumps) are computed within the network model; their current manipulated and error variables are available as result quantities (see Controllers). Pipes with ground heat exchange and ground heat sources with a coupled ground model are represented via separate ground models generated on export, which are coupled to the network model by co-simulation with the configured coupling time step (see Ground Model).

Time Integration

The default integrator is CVODE, an implicit method with an adaptive step size: the time step is automatically increased or decreased based on a local error estimate (relative tolerance 10⁻⁵ by default, maximum time step 1 h). In numerically demanding situations, the steps can temporarily become very small; the outputs are produced independently of this on the output grid (default: hourly).

Result Quantities

With the option Generate default network-related outputs (see Defining Outputs), hourly outputs are generated for all network elements. Which quantities an element provides depends on the component type:

QuantityUnitMeaning
FluidMassFluxkg/sMass flux through the element
FluidVolumeFlowm³/sVolume flow through the element
FluidVelocitym/sFlow velocity (pipes)
InletNodeTemperature / OutletNodeTemperature°CTemperature at the inlet or outlet of the element
TemperatureDifferenceKOutlet minus inlet temperature
InletNodePressure / OutletNodePressurePaPressure at the inlet or outlet node
InletNodePressureAbsolute / OutletNodePressureAbsolutePaPressure including the geodetic component
PressureDifferencePaPressure difference across the element
PressureDifferencePerLengthPa/mPressure gradient per meter of pipe length
PressureLossFittingsInlet / PressureLossFittingsOutletPaPressure loss due to fittings at the inlet/outlet
FlowElementHeatLossWHeat flux from the element to the surroundings; negative if the element injects heat into the network
HeatingPowerWHeating power supplied to the fluid (heat exchanger)
PressureHeadPaHead of the pump
ElectricalPowerWElectrical power consumption (pumps, heat pumps)
PumpEfficiencyCurrent pump efficiency
COPCoefficient of performance of the heat pump
CondenserHeatFlux / EvaporatorHeatFluxWHeat flux at the condenser or evaporator of the heat pump
RequiredBuildingHeatFluxkWHeat power requested by the building (transfer station)
HeatDeficitAbsolutekWUncovered heat demand of the transfer station
HeatDeficitRelativeRatio of heat deficit to requested heat flux
OutletTemperatureSecondary / InletTemperatureSecondary°CSecondary-side temperatures of the transfer station
StoredHeatingEnergykWhCurrently stored heat quantity of the transfer station
ControllerResultValue / ControllerErrorValueManipulated variable and control deviation of controllers

Output Files

Alongside the project file, the simulation creates a folder <ProjectName>/:

  • results/ – result files
  • log/screenlog.txt – solver log

For each result quantity, a file <NetworkName>.<Quantity>.btf (binary format, default) or .tsv (text format) is written, e.g. Network.FluidMassFlux.btf. The format is set by the option Use binary format on the Outputs page. Averaged or integrated outputs receive the suffix -mean or -integral; for an output grid other than hourly, its name is appended.

The TSV files are tab-separated text files: the first column contains the time, followed by one column per network element with the element name, ID and unit in the column header, e.g. SupplyPipe(ID=1612).FluidMassFlux [kg/s].

You evaluate the results directly in the result area of VICUS Districts (false color in the 3D model, line charts) or open them in PostProc 2 – the TSV files can also be further processed in any spreadsheet or scripting tool.

Stay up to date

New features, tutorials and updates delivered to your inbox.