Hello Daniel,
please excuse the delaied answer. Thank you very much for the docker
file. I was able to pull the dockerfile and managed to get christoph.py
running within the dockerenvironment. I still have to get my head around
the whole docker environment though.
Regarding your question. You are fully correct, that you should get the
"ValueError" only with a slight modification of the Code. For the
Equations 3 and 4 you have to write.
eqn3 = E.grad[1] == DiffusionTerm(coeff=(-gamma), var=J) -
TransientTerm(coeff=vg, var=E)\
+TransientTerm(coeff= Beta2/2.0,var=Jstar)-
(k*n2*((E**2.0))+((J**2.0))*J)+(w0*tau*rho*J)\
-((sigma*rho*E)/2)-((multi(numerix.sqrt((E**2.0)+(J**2.0)))*E)/(2*((E**2)+(J**2.0))))
- alpha*E
eqn4 = J.grad[1] == DiffusionTerm(coeff=(gamma), var=E) -
TransientTerm(coeff=vg, var=J)\
-TransientTerm(coeff= Beta2/2.0,var=Estar)+
(k*n2*((E**2))+((J**2))*E)-(w0*tau*rho*E)\
-((sigma*rho*J)/2)-((multi(numerix.sqrt((E**2.0)+(J**2)))*J)/(2.0*((E**2)+(J**2.0))))
- alpha*J
All I did was subtract the TransientTerm on the lefthand side and then reinsert
it on the righthand side. Where you able to view the intermediate results by
uncommenting lines 132 and 133?
Thank you very much for the time you put into this again.
All the Best,
Christoph
_______________________________________________
fipy mailing list
fipy@nist.gov
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
