library(pipenostics)
# Calculate normative temperature drop based on Minenergo-325 for pipe segment
pipeline <- list(
year = 1968,
laying = "channel",
d = 700, # [mm]
len = 1000 # [m]
)
regime <- list(
temperature = c(130, 150), # [°C]
pressure = .588399, # [MPa]
flow_rate = 250 # [ton/hour]
)
pipe_loss_power <- do.call(
m325nhl,
c(pipeline, temperature = list(regime[["temperature"]]), duration = 1) # [kcal/hour]
)
temperature_drop <- dropt(
temperature = regime[["temperature"]], # [°C]
loss_power = pipe_loss_power # [kcal/hour]
) # [°C]
print(temperature_drop)
# [1] 1.366806 1.433840
Run the code above in your browser using DataLab