lambdatot = lambda1 + lambda2
k = Functions.psDL((c-c_c), lambdatot, lambda2, ksquare, mu2, B, c_c)
c = c + h * k
timeps = time.time() - start
if(i % plotEveryNth == 0):
Cons = (c).sum()/(Nx*Ny)
mu, mugrad = Functions.mu((c-c_c), ksquare, kappa, a, b, B, c_c)
#Flux = - lambda1 * (mu + mugrad) - lambda2 * (mu + mugrad - mu2)
Flux = s
musum = (mu+mugrad).sum()/(Nx*Ny)
Fluxtot = Flux.sum()
jx, jy = Functions.flux(c, Nx, Ny, dx, dy)
print(Cons, musum, Fluxtot, np.amin(c), np.amax(c))
Plotting.plot(x,y,(c),0,1,'c',(i/plotEveryNth))
Plotting.Lineplot(x[int(Nx/2)],c[:][int(Nx/2)],int(Nx/2),'Lineplot',(i/plotEveryNth))
Plotting.plot(x,y,Flux,(np.amin(Flux)),(np.amax(Flux)),'Reaction Flux',(i/plotEveryNth))
Plotting.plot(x,y,((mu+mugrad)),(np.amin((mu+mugrad))),(np.amax((mu+mugrad))),'mutot',(i/plotEveryNth))
Plotting.plotvec(X, Y, jy, jx, 'DiffFlux', (i/plotEveryNth))
if (i > 0):
print(timeps, (times[i]-times[i-plotEveryNth]))
# In[ ]:
imp.reload(Plotting)
Plotting.video('c(x,y,t)','Droplets_ReactConc-Lattice2', 1)
Plotting.video('j(x,y,t)','Droplets_ReactJ-Lattice2', 3)
