flame_odefun

flame_odefun, a Python code which uses the multiple-precision arithmetic package mpmath() to approximate the solution of an ordinary differential equation (ODE) which models the growth of a ball of flame in a combustion process.

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Languages:

flame_odefun is available in a Python version.

Related Data and codes:

flame_ode, a Python code which sets up an ordinary differential equation (ODE) that models the growth of a ball of flame in a combustion process.

odefun_test, a Python code which uses odefun() from the multiple precision package mpmath() to solve some simple ordinary differential equations (ODE).

python_ode, Python codes which sets up various systems of ordinary differential equations (ODE).

reaction_twoway_odefun, a Python code which uses the multiple-precision arithmetic package mpmath() to approximate the solution of the ordinary differential equations (ODE) which model a two-way chemical reaction between species W1 and W2.

References:

1. Shirley Abelman, Kailash Patidar,
   Comparison of some recent numerical methods for initial-value problems for stiff ordinary differential equations,
   Computers and Mathematics with Applications,
   Volume 55, Number 4, 2008, pages 733-744.
2. Guy Rouleau,
   Guy on Simulink: Why do we need stiff ODE solvers?,
   "https://blogs.mathworks.com/simulink/2012/07/03/why-do-we-need-stiff-ode-solvers/?s_tid=srchtitle",
   Posted 03 July 2012.
3. Cleve Moler,
   Cleve's Corner: Stiff Differential Equations,
   MATLAB News and Notes,
   May 2003, pages 12-13.
4. Cleve Moler,
   Cleve's Corner: Ordinary Differential Equations, Stiffness,
   "https://blogs.mathworks.com/cleve/2014/06/09/ordinary-differential-equations-stiffness/?s_tid=srchtitle",
   Posted 09 June 2014.

Source Code:

• flame_odefun.py, the source code.
• flame_odefun.sh, runs all the tests.
• flame_odefun.txt, the output file.

• flame_exact.png, a plot of the exact solution, for delta = 0.01;
• flame_odefun.png, a plot of the computed solution, for delta = 0.01;