Расчёт массового расхода при течении жидкости через трубы различных форм сечений

В данном примере демонстрируется применение блока Laminar Leakage библиотеки Изотермической жидкости.

Данный блок моделирует ламинарное течение через сечения различной геометрии в сети изотермической жидкости.

В примере рассматриваются такие сечения как: круг, прямоугольник, эллипс и равносторонний треугольник.

Площади поперечных сечений труб равны между собой.

Схема модели:

Определение функции для загрузки и запуска модели:

function start_model_engee()
    try
        engee.close("npn_transistor", force=true) # закрытие модели 
        catch err # в случае, если нет модели, которую нужно закрыть и engee.close() не выполняется, то будет выполнена её загрузка после catch
            m = engee.load("$(@__DIR__)/laminar_leakages.engee") # загрузка модели
        end;

    try
        engee.run(m, verbose=true) # запуск модели
        catch err # в случае, если модель не загружена и engee.run() не выполняется, то будут выполнены две нижние строки после catch
            m = engee.load("$(@__DIR__)/laminar_leakages.engee") # загрузка модели
            engee.run(m, verbose=true) # запуск модели
        end
end

start_model_engee (generic function with 1 method)

Загрузка, запуск модели и запись результатов

Запуск симуляции

try
    start_model_engee() # запуск симуляции с помощью специальной функции, реализованной выше
    catch err
    end;

Building...
Progress 0%
Progress 33%
Progress 66%
Progress 100%
Progress 100%

Выделение из переменной simout данных о массовых расходах:

sleep(5)
result = simout;
res = collect(result)

4-element Vector{WorkspaceArray}:
 WorkspaceArray("laminar_leakages/Треугольное сечение")
 WorkspaceArray("laminar_leakages/Прямоугольное сечение")
 WorkspaceArray("laminar_leakages/Круглое сечение")
 WorkspaceArray("laminar_leakages/Эллиптическое сечение")

Запись результатов работы модели в отдельные переменные:

circular_flowrate = collect(res[3])
rectangular_flowrate = collect(res[2])
elliptical_flowrate = collect(res[4])
triangular_flowrate = collect(res[1]);

Визуализация результатов моделирования

using Plots
plot(circular_flowrate[:,1], circular_flowrate[:,2], label="Круглое сечение", lw = 3, title="Массовый расход через сечения", xlabel="Время, с", ylabel="Массовый расход, кг/с", legend=:bottomright)
plot!(rectangular_flowrate[:,1], rectangular_flowrate[:,2], label="Прямоугольное сечение", lw = 3)
plot!(elliptical_flowrate[:,1], elliptical_flowrate[:,2], label="Эллиптическое сечение", lw = 3)
plot!(triangular_flowrate[:,1], triangular_flowrate[:,2], label="Треугольное сечение", lw = 3)

Вывод: