Interactive GUI based program that generates the overall species balance, system of ODEs needed for the solve_ivp and odeint method, and calculates the Jacobian both symbolically and numerically. The resulting code can easily be copied and pasted as is to be integrated with the aforementioned SciPy functions.
Project description
Chemical System Generator
Generator for symbolic functions for both full and simple systems
Symbol Generator
chemsys.py is a Gui based chemical reaction system generator that generates the right hand side (RHS) of a chemical system to be used with the solve_ivp and odeint method, the jacobian, and the overall equations for each reactions and each individual chemical species. It will generate Latex formatted equations for the individual chemical reactions, overall mass balances for each species, and both symbolic and numerical text files for the RHS and jacobian.
Installation
pip install chemsys
1. Chemical reactions
This program will generate the Latex forms of each individual reaction and then save them as both a pdf and svg. The string forms of the equations will be saved in a text (txt) file. An example of some inital reactions are given below.
2. Overall reactions for each chemical species
This program will generate the Latex forms of each individual reaction and then save them as both a pdf and svg. The string forms of the equations will be saved in a text (txt) file. An example of some inital reactions are given below.
3. Right Hand Side
The right hand side (RHS) of the system of equations will be generated both symbolically and with initial values substituted into their respective places. This is the required system for scipy's solve_ivp and odeint method. This function requires a callable in the form of fun(t, y):
4. Jacobian
The Jacobian matrix is symbolically generated in order to improve the accuracy and speed of the solvers used in the aforementioned method.
5. Use
5.1 Generate all necessary lists and values as follows.
chemnamesl, rxnsvl, Initreactions, Eqlist, indvdf, ffpath, kk, eaf, RR = guivar.fullgui()
5.2 Calculate the jacobian and all other desired functions as follows.
C, KKS, EAS, reacts, prods, equations, slat, dlat, chem, chemD, chemw, rhs, rhsf, jac, jacnumpy, Jacmath, JacSimple, lm, latexmatrix, jacsy, jacnumpysy, jacmathsy, jacsimplesy, lmsy, latexmatrixsy = symbolgen.fullgen(chemnamesl, rxnsvl, Initreactions, Eqlist, indvdf, ffpath, kk, kJtoJ(eaf), RR)
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.