Hi Jon,
I've been continuing to work on my problem. I've implemented the divergence
of the current density mathematically into the positive charge and negative
charge density equations. Now, I'm encountering issues with the boundary
conditions.. the charges should be flowing towards the center of
When taking the divergence of current density: divergence.(-mu_p *
grad.potential(x,t) * Pion(x,t)) would the divergence not also be applied
to the Pion term, since it is also a function of x. I essentially applied
the product rule to obtain the `Pion * potential.faceGrad.divergence`
part, which
Justin -
A couple of things:
- Charge Density is not Pion + Nion, it's Pion - Nion
- What are the terms `Pion * potential.faceGrad.divergence` in Pion.equation
and `Nion * potential.faceGrad.divergence` in Nion.equation? I don't believe
they should be there.
- Your equations are really not
