<div dir="ltr">Increasing number of snes iterations, I get convergence.<div><br></div><div>So it is a problem of initial guess being too far from the solution of the nonlinear equation.</div><div><br></div><div>Solution can be seen here</div><div><br></div><div><a href="https://github.com/cpraveen/fenics/blob/master/1d/cosmic_ray/cosmic_ray.ipynb">https://github.com/cpraveen/fenics/blob/master/1d/cosmic_ray/cosmic_ray.ipynb</a><br></div><div><br></div><div>Green curve is solution after two time steps.</div><div><br></div><div>It took about 100 snes iterations in first time step and about 50 in second time step.</div><div><br></div><div>I use exact Jacobian and direct LU solve. </div><div><br></div><div>Thanks</div><div>praveen</div></div>