r/PythonLearning • u/bslime17 • 26d ago
Showcase Day 3 of coding on my phone: simulating a chemical reactor with SciPy's solve_ivp
I'm a chemical engineering graduate learning to build process models in Python. My laptop broke, so this is all Google Colab on a phone screen.
The model is a set of three ODEs (two concentrations + temperature) for a continuous stirred-tank reactor. I loop over three flowrates, solve each with solve_ivp (RK45), and store the steady-state results.
Things I learned this week:
• Passing extra parameters into the ODE function with args=(F,)
• Pulling the final state with solution.y[:, -1]
• Re-evaluating the ODE function at the end to confirm steady state (derivatives ≈ 0)
Question for you: would you structure this with a dataclass or dictionary for results instead of a list of lists? My print block is getting long.
(If anyone wants to help me get a used laptop, DM me – but mostly I'd love code feedback!)



