Programming numeric simulations of thermodynamic systems in Julia

Hi Julia Community,

I’m really new to the Julia language and this is my first post in the forum. Moreover english is not my mother tongue, so I hope you can understand my question and help me to increase my skills with Julia.

I try to use the Julia language to calculate numeric simulations of different thermodynamic systems. I have thought of a possibility to structure my script, but I’m not really sure, if it’s the way to go in Julia to reach a good performance. Here is a little example what I wanna do:

Part 1:
I want to simulate different types of thermodynamic systems with different kinds of fluids and connections between the systems. So in the first step I have to create my system. I think that structs will be a good option to structure my system:

abstract type System end
abstract type Fluid end

struct Air <: Fluid end

struct Container <: System
    fluid::Fluid
    T::Vector{Float64}
    m::Vector{Float64}    
    p::Vector{Float64}
    V::Vector{Float64}
end

struct ThermoSys
    systems::Vector{System}
    connections::Vector{Any}
end

function calc_mass(V, T)
    return 1*10^5 * V / (0.287069 *10^3 * T)
end

function create_thermo_systems()
    system1 = Container(Air(), [300], [calc_mass(0.01, 300)], [1*10^5], [0.01])
    system2 = Container(Air(), [300], [calc_mass(0.05, 300)], [1*10^5], [0.05])
    ExampleSystem = ThermoSys([system1, system2], [1, 2, 0.05])
    return ExampleSystem
end

Thats just a small example. The systems will be much more complex in the final program. In this example is a function which creates two containers of different volume and filled with air. The fluid has a temperature, a pressure and a mass. The containers are connected with each other. I decided to use vectors for the definition of temperature, mass, etc… to update the data in the simulation and use the created system as a data storage.

Part 2:
Now, I hand the created System over to the simulation module. In this part, the program iterates through all systems and the solver does a stepwise calculation of it with some differential equations. The right calculations will be dependent of the type of the system or the chosen fluid with multiple dispatch. After the calculation all data vectors will be updated with the new values.

function calc_system(system_i::System)
    T_new = calc_new_temperature_of_fluid(system_i.fluid)
    push!(system_i.T, T_new)
end

In this example is a functions, which will be called of an argument of the type system and calculates a new temperature step for the fluid in the system. The new value is added to the vector after the calculation and can be used in the next calculation as starting value with “system.T[end]”. Moreover I’m able to plot the changes of all values in every system after the calculation.

Question:
I’m really sure that my chosen way to structure the simulation is possible. But I’m struggling with the performance, if the systems will become more complex and maybe there are multiple fluids in one zone. Is it a good option to create my systems with the structs and use it as data storage in the stepwise calculation? Will it become slow in complex calculations? How can I change it to perform better? I’m really happy about help and thank you in advance.

Greetings
Eulerian

2 posts - 2 participants

Read full topic