I am solving a nonlinear poisson equation of the form
$$-\nabla \cdot (\epsilon \nabla u) = (N - \rho*e^{(u-f)/c})$$
I am getting an unexpected results, and I'd like to know if I am defining the question properly. My question is, is the following the correct definition of the above problem?
u = Function(V)
v = TestFunction(V)
F = inner(eps*grad(u),grad(v)))*dx - (N - rho*exp((u-f)/c)) * v * dx
solve(F == 0, u, boundary_conditions, solver_parameters = \
{'newton_solver':{'relative_tolerance':1e6}})
Thanks