I am trying to solve the following problem using fipy but I am a bit overwhelmed. Hoping someone can point me in the right direction.
C and q are the dependent variables and I am solving in 2 dimensions, time (t) and length (z).
Is there any example I can follow? Where should I be looking at in the documentation? Or is there a some other library better suited for this problem?
Thanks in advance.
examples.diffusion.mesh1D is the starting point for anybody using FiPy. It runs through the basic concepts of how to set up a problem and it covers the transient and diffusive behavior of your first equation.
examples.convection.robin illustrates a static convection-diffusion-source problem with a Robin boundary condition. This covers all the terms on the right-hand side of your first equation and your first boundary condition.
The two equations can either be "swept" in succession until they converge, or they can be coupled. examples.diffusion.coupled illustrates both approaches for a pair of diffusion equations.
If these aren't enough to get you started, please come back with specific questions about where you're having trouble.
