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
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
