Hi everyone
I am working in RxD to find Na and k concentration inside axon, I know that RxD equation consists many terms:
\[\frac{\partial u}{\partial t} = d \nabla^2 u + f(u, t)\]
but I would like to add a new term to the right side, I know that the first term on the left side of the equation above " d \nabla^2 u" is one term of NernstPlanck flux. I would like to add the other term of NernstPlanck flux to the RxD equation. Can anyone help me please?
Add term to RxD solution.
Moderators: hines, wwlytton, ramcdougal

 Posts: 49
 Joined: Thu Jun 07, 2018 10:57 am
Re: Add term to RxD solution.
Can anyone help me please ? is there answer to my question ?

 Posts: 207
 Joined: Fri Nov 28, 2008 3:38 pm
 Location: Yale School of Public Health
Re: Add term to RxD solution.
The bad news: The reactiondiffusion module rxd currently only directly supports reactions and diffusions. So you could add arbitrary local kinetics, but you'd have to modify the source to add the type of dynamics you describe.
The good news: besides storing parameters, all of the changes would be confined to neuron.rxd.species.section1d.Section1d._setup_diffusion_matrix which ought to be cleaned up anyways. There's a little complexity based on the use of algebraic/0volume/conservationofmass nodes at section boundaries, but basically it would be a matter of setting the right entries in the Jacobian matrix g.
The good news: besides storing parameters, all of the changes would be confined to neuron.rxd.species.section1d.Section1d._setup_diffusion_matrix which ought to be cleaned up anyways. There's a little complexity based on the use of algebraic/0volume/conservationofmass nodes at section boundaries, but basically it would be a matter of setting the right entries in the Jacobian matrix g.