NumericalMethods

Polymerization equation in 3D

Polymerization.py

Explicit Euler simulation of a exothermic self-supporting polymerization reaction. Executed on GPU using tensorflow.

Solvers for partial differential equations

LinearDifferentialEqSolvers.ipynb

Implemented:

• Explicit Euler
• Implicit Euler
• Runge-Kutta
• Predictor-Corrector

A sigle trajectory with different methods.

Experiments with different values of R

All trajectories tend to a stable point in zero. R = 1.

All trajectories tend toward one of the two stable points. R = 5.

Chaos (the trajectory changes drastically if the starting point changes slightly).

Restoring from the chaos. Trajectories tend to a single loop.