- #1
TheCanadian
- 367
- 13
I've been trying to solve a set of coupled, nonlinear PDEs with the general form:
## \frac {\partial A}{\partial t} = aAC + bBD - cE ##
## \frac {\partial B}{\partial t} = - c(E+F) ##
## \frac {\partial C}{\partial t} = dAD - eE ##
## \frac {\partial D}{\partial t} = dBC + - eF ##
## \frac {\partial E}{\partial z} = cCD - fE ##
## \frac {\partial F}{\partial z} = cCD - fF ##
(Upper case letters stand for the solutions sought; lower case letters represent constants)
I've tried existing packages such as FiPy and tried my own basic solver using FEM, but neither has been successful aside from a few nontrivial cases. With the two independent variables (z and t) especially, this seems to complicate the solver. Any advice you have on existing libraries or numerical methods I should implement would be of great interest. (I am most familiar with Python, but am open to responses for any programming language.)
## \frac {\partial A}{\partial t} = aAC + bBD - cE ##
## \frac {\partial B}{\partial t} = - c(E+F) ##
## \frac {\partial C}{\partial t} = dAD - eE ##
## \frac {\partial D}{\partial t} = dBC + - eF ##
## \frac {\partial E}{\partial z} = cCD - fE ##
## \frac {\partial F}{\partial z} = cCD - fF ##
(Upper case letters stand for the solutions sought; lower case letters represent constants)
I've tried existing packages such as FiPy and tried my own basic solver using FEM, but neither has been successful aside from a few nontrivial cases. With the two independent variables (z and t) especially, this seems to complicate the solver. Any advice you have on existing libraries or numerical methods I should implement would be of great interest. (I am most familiar with Python, but am open to responses for any programming language.)