This repo contains the Python code and the Jupyter notebooks which were used to produce the simulation results shown in my master's thesis which I wrote in 2019 at FAU Erlangen.

More information about this research on port-Hamiltonian systems and Stirling-cycle devices can be found on my personal research website.

The code folder contains the Python files which are imported from the Jupyter notebooks.

The notebook butcher.ipynb was used to learn how to compute Butcher tableaus for Gauss-Legendre collocation methods.

The notebooks harmonic_oscillator.ipynb, kepler.ipynb and spring_pendulum.ipynb deal with conservative mechanical systems and merely serve to test the implementation of the Gauss-Legendre collocation methods.

The notebooks piston_animation_euler.ipynb, piston_animation.ipynb and piston_animation2.ipynb deal with the central example of the thesis. The first one just uses the explicit Euler scheme. The other two use Gauss-Legendre collocation methods and rely on automated code generation using SymPy.

The notebook ideal_gas_sym.ipynb is supposed to check the equations which describe an ideal gas.

The notebook piston.ipynb contains some symbolic computations concerning the modeling of the piston example which I also could have done on paper.

The notebook carnot_efficiency.ipynb simply plots the Carnot efficiency for different temperatures.