Skip to contents

Let’s consider a tiny district heating model which consists with the next 4-segment tracing path:

Figure: Regime tracing
Figure: Regime tracing

Suppose we have the next sensor readings for forward tracing:

t_fw <- 130         # [°C]
p_fw <-   0.588399  # [MPa]
g_fw <- 250         # [ton/hour]

Let’s discharges to network for each pipeline segment are somehow determined as

discharges <- seq(0, 30, 10)  # [ton/hour]
print(discharges)
#> [1]  0 10 20 30

Then the calculated regime (red squares) for forward tracing is

library(pipenostics)

regime_fw <- m325traceline(t_fw, p_fw, g_fw, discharges, forward = TRUE)
print(regime_fw)
#> $temperature
#> [1] 129.1698 128.4075 127.9377 127.3039
#> 
#> $pressure
#> [1] 0.5877977 0.5873084 0.5870757 0.5868732
#> 
#> $flow_rate
#> [1] 250 240 220 190
#> 
#> $loss
#> [1] 352.2900 351.4092 350.6004 350.1019
#> 
#> $flux
#> [1] 183.4259 182.9672 182.5461 182.2866
#> 
#> $Q
#> [1] 5072976 4469925 2524323 2940856

ℹ Read Concepts and useful notes for a deeper dive into the topic.