Hey all, Quick question that is completely throwing me for a loop. I'm working on a poisson equation for pressure and wanting to specify dirichlet boundary conditions. I used the penalty method (per libmesh documentation), and it works beautifully, aka matches the real solution. However, when I change to 2nd order polynomials, everything falls apart. The penalty method does implement the boundary conditions, but the solution is complete garbage.
Any explanation for changing (.., FIRST) to (.., SECOND) and the solution (the code runs fine) but the solution being garbage? Cheers, Caleb Phillips _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
