On Fri, 27 Jan 2017, Manav Bhatia wrote:
The subdomain approach sounds promising. I will look into this and
the SCALAR variable.
SCALAR you want to avoid for locally supported variables in fully
implicit solves. We generate dense rows and columns for each SCALAR
dof, which is not what you want for efficiency when you know a priori
that nearly all that coupling isn't really there.
Another question: the application in mind is topology optimization
with XFEM, where voids can form/disappear in a continuum during the
course of optimization. I don’t want to modify the mesh during this
procedure. So, if a certain portion of the domain becomes a void,
what would be the best approach to make the dofs inactive in that
void?
Can I simply iterate on the node/elem dofs and flag them as
active/inactive so that I don’t have to worry about
solving/specifying Dirichlet conditions for them?
I'd add a user constraint equation for each dof.
---
Roy
------------------------------------------------------------------------------
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
