[Libmesh-users] Dirichlet boundary condition setting variable to zero?

[Renato P. <re...@gm...>]

Hi

I am struggling with a Dirichlet boundary condition setting the
boundary to zero, instead of the value I impose (3.e7). Previously, I
was using the penalty method, and it was working fine.

Is there any known corner condition where the DirichletBoundary would
drive the variable to zero? Any idea?

This is the code:
      std::set<boundary_id_type> boundary_ids;
      std::vector<unsigned int> variables;
      boundary_ids.insert(BOUNDARY_RESERVOIR);
      variables.push_back(_pvar);
      DirichletBoundary dirichlet_bc(boundary_ids, variables,
ConstFunction<>(3.e7));
      _geomec_sys.get_dof_map().add_dirichlet_boundary(dirichlet_bc);


Note: I read some posts on subdomain issues with BCs. My elements are
in "subdomain=2", but as of this test, I removed all other subdomains,
without success. Should this matter at all?

Thanks
Renato