On 03/11/17 20:15, Francis Poulin wrote:
If yes, then I think I'm in good shape. One problem is that the boundary layer is very narrow for my choice of parameters so it's hard to know whether the boundary value is really zero. I don't suppose there is an easy way to find out the maximum of the boundary value in the numerical solution?
If you are imposing strong (Dirichlet) conditions, then the solution on the boundary nodes is, by definition, whatever you set it to. If you have some other means of enforcing boundary conditions and just want to know that the values are on the boundary, you could do: V = FunctionSpace(...) selector = DirichletBC(V, 0, "on_boundary") f = Function(V) # solve into f boundary_values = f.dat.data_ro[selector.nodes, ...] Lawrence