solve Anisotropy equation
4 months ago by
by changing the bilinear form to
a = inner(sigma*grad(u), grad(v))*dx
However, note that \(\sigma\) should be positive definite for the problem to be stable. (Physically: "heat needs to flow from hot to cold".) The suggested matrix has a negative eigenvalue, and will produce a very strange-looking discrete solution.
Please login to add an answer/comment or follow this question.