pith. sign in

arxiv: 1907.05255 · v1 · pith:RNXMAZHCnew · submitted 2019-07-11 · 🧮 math.OC

Optimal control of district heating networks using a reduced order model

Markus Rein , Jan Mohring , Tobias Damm , Axel Klar This is my paper

Pith reviewed 2026-05-24 23:00 UTC · model grok-4.3

classification 🧮 math.OC
keywords optimal controldistrict heating networksreduced order modelsdifferential algebraic equationsenergy advectionflow direction changesfeed-in power
0
0 comments X

The pith

A reduced-order model based on system theory enables fast optimal control of district heating networks with changing flows.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper shows how to use a reduced-order model to solve optimal control problems for district heating networks much more quickly than with the full dynamics. The networks are governed by index-1 quadratic differential algebraic equations that describe how injected energy advects through the pipes. The optimization goal is to limit the maximum feed-in power without solving the full model online. The authors demonstrate that their system-theoretic reduction works on a real large-scale network where flow directions can reverse. Readers would care because faster optimization could let operators run heating systems closer to their efficiency limits in practice.

Core claim

The central claim is that the suggested reduced model decreases the computation time of the optimization significantly. The effectiveness of the presented approach is demonstrated for an existing, large-scale heating network including changes of flux directions. The model is based on a system theoretic description close to the underlying Euler equations and addresses the challenge of an index-1 quadratic in state differential algebraic equation.

What carries the argument

Reduced order model derived from a system-theoretic description of the index-1 quadratic differential-algebraic equation for energy advection.

If this is right

  • Computation time for finding optimal controls is reduced substantially.
  • The approach remains effective on large-scale networks.
  • Optimal controls can be determined even when flux directions change.
  • The method allows limiting maximal feed-in power as a product of control and state variables.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar reduced models could be developed for other types of energy distribution networks.
  • The technique may support real-time adjustments based on current network states.
  • Further work could test the reduced model on networks with different topologies or scales.

Load-bearing premise

The reduced-order model approximates the full dynamics closely enough to produce useful optimal controls when flow directions change.

What would settle it

Apply the controls found by the reduced model to a simulation of the original full-order model and check if the feed-in power stays below the limit.

Figures

Figures reproduced from arXiv: 1907.05255 by Axel Klar, Jan Mohring, Markus Rein, Tobias Damm.

Figure 1
Figure 1. Figure 1: Topology of an existing heating network supplying a district. The flow part of the [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Optimal control problem for TC1 at −3 ◦C comparing the initial control (red) to the optimized control (green), obtained by the reduced model ROM1. Part (a) shows controls(solid lines) and the total volume flow injected at the power plant(dashed lines). Part (b) presents the feed-in power for both controls together with the mean consumption (lower dashed line) and the feed-in constraint P¯ (upper dashed lin… view at source ↗
Figure 3
Figure 3. Figure 3: Result of the optimal control problem for TC1 at [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Optimal control problem for TC2 at 3 ◦C comparing the initial control (red) to the optimized control (green), obtained by the reduced model ROM1. For a detailed explanation we refer to fig. 2. be zero as well. Hence, the absolute error in the power approximation remains small. FOM 0 FOM 1 ROM 1 FOMR runtime opt./s 2525.0 3210.0 596.2 18300.0 runtime sim./s 2455.0 3146.7 546.2 18100.0 # solves DAE 15 11 11 … view at source ↗
Figure 5
Figure 5. Figure 5: Temperature signal for TC3 at the consumer exhibiting the largest relative [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Optimal control problem for TC3 at 7.5 ◦C comparing the initial control (red) to the optimized control (green), obtained by the reduced model ROM1. For a detailed explanation we refer to fig. 2. of the DAE with respect to energy densities. Since only the thermal transport is reduced, ch is identical for both full and reduced order models. The cost for the evaluation of the ODE scales with the number of ent… view at source ↗
read the original abstract

We study the optimal control of district heating networks using a reduced order model based on a system theoretic description close to the underlying Euler equations. In the presented scenarios, the central task is to limit the maximal feed-in power occurring as a product of control and state variables. The underlying dynamics of heating networks acting as optimization constraints pose the central computational complexity, prohibiting the determination of an optimal control online. The advection of the injected energy density on the network results in an index-1, quadratic in state differential algebraic equation, challenging to reduce. The suggested reduced model decreases the computation time of the optimization significantly. The effectiveness of the presented approach is demonstrated for an existing, large-scale heating network including changes of flux directions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper develops a reduced-order model (ROM) derived via system-theoretic methods from the Euler equations for district heating networks. The ROM is applied to an optimal control problem whose objective is to limit maximal feed-in power (a bilinear product of control and state). The underlying dynamics are an index-1 quadratic DAE; the ROM is asserted to cut optimization runtime substantially while remaining effective on a real large-scale network that exhibits flux-direction reversals.

Significance. If the ROM preserves the advection dynamics and bilinear terms sufficiently well under flow reversals, the work would supply a practical route to online optimal control of district heating systems, a setting where full-order DAE optimization is currently prohibitive. The demonstration on an existing large network is a positive feature; however, the absence of quantitative fidelity metrics leaves the practical utility unverified.

major comments (3)
  1. [Abstract] Abstract: the central claim that the ROM 'decreases the computation time of the optimization significantly' and is 'effective' on a network with flux-direction changes is unsupported by any reported speedup factor, trajectory error, or comparison of the ROM-derived control against the full-order index-1 DAE; without such metrics the claim that the computed controls remain feasible for the true network cannot be assessed.
  2. [Abstract] Abstract (and implied § on model reduction): the system-theoretic reduction is stated to be 'close to the underlying Euler equations,' yet no a-posteriori error bound, residual estimate, or numerical validation of the ROM state trajectory versus the full quadratic DAE is supplied, especially across the flow-reversal events that alter the advection structure; this directly affects whether the ROM controls satisfy the feed-in power limit on the original system.
  3. [Abstract] Abstract: the demonstration on a 'large-scale heating network including changes of flux directions' is presented as evidence of effectiveness, but the manuscript supplies neither a baseline comparison (e.g., full-order optimization where feasible, or alternative reduction techniques) nor any quantitative measure of sub-optimality or constraint violation when the ROM control is applied to the unreduced model.
minor comments (2)
  1. Notation for the quadratic DAE and the bilinear feed-in power term should be introduced with explicit equation numbers in the main text rather than left implicit in the abstract.
  2. [Abstract] The abstract refers to 'the presented scenarios' without indicating how many networks or parameter regimes were tested; a brief enumeration would improve clarity.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the thorough review and constructive suggestions. We agree that the abstract would benefit from more quantitative details and will revise the manuscript accordingly. Our responses to the major comments are as follows.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that the ROM 'decreases the computation time of the optimization significantly' and is 'effective' on a network with flux-direction changes is unsupported by any reported speedup factor, trajectory error, or comparison of the ROM-derived control against the full-order index-1 DAE; without such metrics the claim that the computed controls remain feasible for the true network cannot be assessed.

    Authors: We acknowledge that the abstract does not include explicit numerical values. In the revised manuscript, we will incorporate specific speedup factors observed in our experiments. However, a direct comparison between the ROM-derived control and the full-order optimal control is not feasible, as solving the full-order optimization problem is computationally prohibitive for the large-scale network—this being the primary motivation for developing the ROM. We will clarify this limitation and provide available validation metrics. revision: partial

  2. Referee: [Abstract] Abstract (and implied § on model reduction): the system-theoretic reduction is stated to be 'close to the underlying Euler equations,' yet no a-posteriori error bound, residual estimate, or numerical validation of the ROM state trajectory versus the full quadratic DAE is supplied, especially across the flow-reversal events that alter the advection structure; this directly affects whether the ROM controls satisfy the feed-in power limit on the original system.

    Authors: We agree that explicit a-posteriori error bounds or residual estimates are not provided. We will add numerical validation comparing the ROM state trajectories to those of the full quadratic DAE, with particular attention to flow-reversal events, to better assess the approximation quality and its impact on the control performance. revision: yes

  3. Referee: [Abstract] Abstract: the demonstration on a 'large-scale heating network including changes of flux directions' is presented as evidence of effectiveness, but the manuscript supplies neither a baseline comparison (e.g., full-order optimization where feasible, or alternative reduction techniques) nor any quantitative measure of sub-optimality or constraint violation when the ROM control is applied to the unreduced model.

    Authors: We agree that quantitative measures of sub-optimality and constraint violations would strengthen the claims. In the revision, we will include such metrics by applying the ROM-derived controls to the full-order model and reporting any violations. Baseline comparisons to alternative reduction techniques will also be considered if space permits. As noted, full-order optimization itself is not feasible as a baseline. revision: partial

standing simulated objections not resolved
  • Direct comparison of ROM optimal controls to full-order optimal controls on the large-scale network due to the prohibitive computational cost of the latter.

Circularity Check

0 steps flagged

No circularity: ROM derived from system-theoretic reduction of index-1 DAE

full rationale

The derivation chain begins from the Euler equations discretized as an index-1 quadratic DAE, applies standard system-theoretic model reduction (e.g., moment matching or balanced truncation) to obtain the ROM, and then uses the ROM inside an optimal-control formulation whose objective is the feed-in power product. None of these steps is shown to be self-definitional, a fitted input renamed as prediction, or dependent on a load-bearing self-citation whose own justification collapses. The single large-scale demonstration with flow reversals is presented as empirical validation rather than a mathematical identity. The central claim therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No free parameters, axioms, or invented entities are identifiable from the abstract. The approach rests on an unspecified system-theoretic reduction of Euler equations to an index-1 quadratic DAE, but no explicit assumptions or fitted quantities are stated.

pith-pipeline@v0.9.0 · 5645 in / 1116 out tokens · 25901 ms · 2026-05-24T23:00:14.341455+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages · 1 internal anchor

  1. [1]

    District heating and cooling: Review of technology and potential enhancements

    Rezaie B, Rosen MA. District heating and cooling: Review of technology and potential enhancements. Appl. Energy 2012; 93: 2–10

    work page 2012

  2. [2]

    A Review of District Heating Systems: Modeling and Optimization

    Talebi B, Mirzaei PA, Bastani A, Haghighat F. A Review of District Heating Systems: Modeling and Optimization. Front. Built Environ. 2016; 2: 1–14

    work page 2016

  3. [3]

    Model order reduction for PDE constrained optimiza- tion

    Benner P, Sachs E, Volkwein S. Model order reduction for PDE constrained optimiza- tion. In: Springer. 2014 (pp. 303–326)

    work page 2014

  4. [4]

    Explicit Model Predictive Control for Large-Scale Systems via Model Reduction

    Hovland S, Gravdahl JT, Willcox KE. Explicit Model Predictive Control for Large-Scale Systems via Model Reduction. J. Guid. Control Dyn. 2008; 31(4): 918–926

    work page 2008

  5. [5]

    Model Order Reduction Techniques with a Posteriori Error Control for Nonlinear Robust Optimization Governed by Partial Differential Equations

    Lass O, Ulbrich S. Model Order Reduction Techniques with a Posteriori Error Control for Nonlinear Robust Optimization Governed by Partial Differential Equations. SIAM J. Sci. Comput. 2017; 39(5): S112–S139

    work page 2017

  6. [6]

    Optimal Control in Networks of Pipes and Canals

    Colombo RM, Guerra G, Herty M, Schleper V. Optimal Control in Networks of Pipes and Canals. SIAM J. Control Optim. 2009; 48(3): 2032–2050

    work page 2009

  7. [7]

    Electric transmission lines: Control and numerical discretization

    G¨ ottlich S, Herty M, Schillen P. Electric transmission lines: Control and numerical discretization. Optim. Control Appl. Meth. 2016; 37(5): 980–995

    work page 2016

  8. [8]

    Mixed integer linear models for the optimization of dynamical transport networks

    Geißler B, Kolb O, Lang J, Leugering G, Martin A, Morsi A. Mixed integer linear models for the optimization of dynamical transport networks. Math. Meth. Oper. Res. 2011; 73(3): 339–362

    work page 2011

  9. [9]

    On Structure- Preserving Model Reduction for Damped Wave Propagation in Transport Networks

    Egger H, Kugler T, Liljegren-Sailer B, Marheineke N, Mehrmann V. On Structure- Preserving Model Reduction for Damped Wave Propagation in Transport Networks. SIAM J. Sci. Comput. 2018; 40(1): A331–A365

    work page 2018

  10. [10]

    Model Order Re- duction of Differential Algebraic Equations Arising from the Simulation of Gas Trans- port Networks

    Grundel S, Jansen L, Hornung N, Clees T, Tischendorf C, Benner P. Model Order Re- duction of Differential Algebraic Equations Arising from the Simulation of Gas Trans- port Networks. In: Sch¨ ops S, Bartel A, G¨ unther M, ter Maten EJW, M¨ uller PC. , eds. Progress in Differential-Algebraic Equations Springer Berlin Heidelberg; 2014; Berlin, Heidelberg: 183–205

    work page 2014

  11. [11]

    Model Reduction for Circuit Simulation

    Benner P, Hinze M, ter Maten EJW. Model Reduction for Circuit Simulation . 74 of Lecture Notes in Electrical Engineering. Springer Netherlands . 2011. 23

    work page 2011

  12. [12]

    Model reduction for DAEs with an application to flow control

    Borggaard JT, Gugercin S. Model reduction for DAEs with an application to flow control. In: Springer. 2015 (pp. 381–396)

    work page 2015

  13. [13]

    Model and Discretization Error Adaptivity Within Stationary Gas Transport Optimization

    Mehrmann V, Schmidt M, Stolwijk JJ. Model and Discretization Error Adaptivity Within Stationary Gas Transport Optimization. Vietnam J. Math. 2018

    work page 2018

  14. [14]

    Solving power-constrained gas transportation problems using an MIP-based alternating direction method

    Geißler B, Morsi A, Schewe L, Schmidt M. Solving power-constrained gas transportation problems using an MIP-based alternating direction method. Comput. Chem. Eng. 2015; 82: 303–317

    work page 2015

  15. [15]

    Challenges in Optimal Control Problems for Gas and Fluid Flow in Networks of Pipes and Canals: From Modeling to Industrial Applications

    Hante FM, Leugering G, Martin A, Schewe L, Schmidt M. Challenges in Optimal Control Problems for Gas and Fluid Flow in Networks of Pipes and Canals: From Modeling to Industrial Applications. In: Manchanda P, Lozi R, Siddiqi AH. , eds. Industrial Mathematics and Complex Systems Springer Singapore. 2016 (pp. 77–122)

    work page 2016

  16. [16]

    MIP-based instantaneous control of mixed-integer PDE-constrained gas transport problems

    Gugat M, Leugering G, Martin A, Schmidt M, Sirvent M, Wintergerst D. MIP-based instantaneous control of mixed-integer PDE-constrained gas transport problems. Com- put. Optim. Appl. 2017; 70(1): 267–294

    work page 2017

  17. [17]

    Towards simulation based mixed-integer optimization with differential equations

    Gugat M, Leugering G, Martin A, Schmidt M, Sirvent M, Wintergerst D. Towards simulation based mixed-integer optimization with differential equations. NETWORKS 2018; 72(1): 60–83

    work page 2018

  18. [18]

    Optimal operation of large district heating networks through fast fluid-dynamic simulation.Energy 2016; 102: 586–595

    Guelpa E, Toro C, Sciacovelli A, Melli R, Sciubba E, Verda V. Optimal operation of large district heating networks through fast fluid-dynamic simulation.Energy 2016; 102: 586–595

    work page 2016

  19. [19]

    Predictive control of a complex district heating network

    Sandou G, Font S, Tebbani S, Hiret A, Mondon C. Predictive control of a complex district heating network. In: . 44. Citeseer. ; 2005: 7372

    work page 2005

  20. [20]

    Optimal Temperature Control of Large Scale District Heating Networks

    Bavi` ere R, Vall´ ee M. Optimal Temperature Control of Large Scale District Heating Networks. Energy Procedia 2018; 149: 69–78

    work page 2018

  21. [21]

    Fundamentals of momentum, heat, and mass transfer

    Welty J, Rorrer GL, Foster DG. Fundamentals of momentum, heat, and mass transfer. In: Wiley. 5 ed. 2008

    work page 2008

  22. [22]

    Numerical methods for conservation laws

    LeVeque RJ. Numerical methods for conservation laws . Lectures in mathematics- Birkh¨ auser. 2 ed. 2008. OCLC: 552210683

    work page 2008

  23. [23]

    Model order reduction of hyperbolic systems at the example of district heating networks

    Rein M, Mohring J, Damm T, Klar A. Model order reduction of hyperbolic systems at the example of district heating networks. arXiv:1903.03342 [math] 2019

    work page internal anchor Pith review Pith/arXiv arXiv 1903

  24. [24]

    Approximation of Large-Scale Dynamical Systems

    Antoulas A. Approximation of Large-Scale Dynamical Systems . Society for Industrial and Applied Mathematics . 2005

    work page 2005

  25. [25]

    Structure preserving model reduction of port- Hamiltonian systems by moment matching at infinity.Automatica 2011; 46(4): 665–672

    Polyuga RV, van der Schaft A. Structure preserving model reduction of port- Hamiltonian systems by moment matching at infinity.Automatica 2011; 46(4): 665–672. 24

    work page 2011

  26. [26]

    Gugercin S, Polyuga RV, Beattie C, Schaft v. dA. Structure-preserving tangential in- terpolation for model reduction of port-Hamiltonian systems. Automatica 2012; 48(9): 1963–1974

    work page 2012

  27. [27]

    A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems

    Benner P, Gugercin S, Willcox K. A Survey of Projection-Based Model Reduction Methods for Parametric Dynamical Systems. SIAM Rev. 2015; 57(4): 483–531

    work page 2015

  28. [28]

    Reduced basis method for finite volume approximations of parametrized linear evolution equations

    Haasdonk B, Ohlberger M. Reduced basis method for finite volume approximations of parametrized linear evolution equations. ESAIM Math. Model. Numer. Anal. 2008; 42(2): 277–302

    work page 2008

  29. [29]

    H2 Model Reduction for Large-Scale Linear Dy- namical Systems

    Gugercin S, Antoulas AC, Beattie C. H2 Model Reduction for Large-Scale Linear Dy- namical Systems. SIAM J. Matrix Anal. Appl. 2008; 30(2): 609–638

    work page 2008

  30. [30]

    Anwendung von Standardlastprofilen zur Belieferung nichtleistungsgemessener Kunden

    Bundesverband der deutschen Gas- und Wasserwirtschaft (BGW) . Anwendung von Standardlastprofilen zur Belieferung nichtleistungsgemessener Kunden. tech. rep., Bun- desverband der deutschen Gas- und Wasserwirtschaft (BGW); 2006. 25

    work page 2006