Hello, I need to change my boundary conditions during the calculation, more precisely I need to remove or disable a dirichlet bc, which was added through DofMap.add_dirichlet_boundary(..). As far as I can see there is no .remove_dirichlet_boundary(..). Is there some special value like NAN that my FunctionBase<> subclass can return to indicate that no constraints should be applied? I have another issue with DirichletBoundaries when using together with PeriodicBoundary conditions: Supposed there is a cube with periodic bc on the left and right side and dirichlet bc on the top side, the nodes on the edge between the left and upper side as well as the nodes on the edge between right and upper side will always have periodic bc. As far as I found out this is because the PR-BC were applied when add_periodic_bc is called and the dirichlet-BC were applied when reiniting the equation-systems. So I found no way to get all nodes on the upper side to be constrained with dirichlet-bc in this scenario, is there any way to do this?
regards Robert ------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
