From: David K. <dknez@MIT.EDU> - 2010-03-27 21:59:23
|
I've been using subdomain-only variables for a conjugate heat transfer problem and I wanted to introduce some mesh refinement "by hand" by marking some specific elements for refinement. Doing this, though, appears to corrupt the system assembly (as evidenced by non-convergence of Newton iterations). The issue appears to be independent of which elements I mark for refinement. Also, if I do a global uniform refinement then everything works fine. So I wonder if there might be a problem with the hanging node constraints for subdomain-only variables? I'd be happy to prepare a test case, but thought I'd run this by the list first to see if anyone has any thoughts/suggestions. David |
From: David K. <dknez@MIT.EDU> - 2010-03-27 22:36:30
|
I may have been a bit hasty with that last message, looks like the issue could be with my user code so please disregard for now. David Knezevic wrote: > I've been using subdomain-only variables for a conjugate heat transfer > problem and I wanted to introduce some mesh refinement "by hand" by > marking some specific elements for refinement. Doing this, though, > appears to corrupt the system assembly (as evidenced by non-convergence > of Newton iterations). The issue appears to be independent of which > elements I mark for refinement. Also, if I do a global uniform > refinement then everything works fine. So I wonder if there might be a > problem with the hanging node constraints for subdomain-only variables? > > I'd be happy to prepare a test case, but thought I'd run this by the > list first to see if anyone has any thoughts/suggestions. > > > David > > ------------------------------------------------------------------------------ > Download Intel® Parallel Studio Eval > Try the new software tools for yourself. Speed compiling, find bugs > proactively, and fine-tune applications for parallel performance. > See why Intel Parallel Studio got high marks during beta. > http://p.sf.net/sfu/intel-sw-dev > _______________________________________________ > Libmesh-devel mailing list > Lib...@li... > https://lists.sourceforge.net/lists/listinfo/libmesh-devel |
From: David K. <dknez@MIT.EDU> - 2010-03-28 00:03:49
|
Yes, it was a user code issue :) The subdomain-only variable functionality is working very nicely. - David David Knezevic wrote: > I may have been a bit hasty with that last message, looks like the issue > could be with my user code so please disregard for now. > > > > David Knezevic wrote: >> I've been using subdomain-only variables for a conjugate heat transfer >> problem and I wanted to introduce some mesh refinement "by hand" by >> marking some specific elements for refinement. Doing this, though, >> appears to corrupt the system assembly (as evidenced by >> non-convergence of Newton iterations). The issue appears to be >> independent of which elements I mark for refinement. Also, if I do a >> global uniform refinement then everything works fine. So I wonder if >> there might be a problem with the hanging node constraints for >> subdomain-only variables? >> >> I'd be happy to prepare a test case, but thought I'd run this by the >> list first to see if anyone has any thoughts/suggestions. >> >> >> David >> >> ------------------------------------------------------------------------------ >> >> Download Intel® Parallel Studio Eval >> Try the new software tools for yourself. Speed compiling, find bugs >> proactively, and fine-tune applications for parallel performance. >> See why Intel Parallel Studio got high marks during beta. >> http://p.sf.net/sfu/intel-sw-dev >> _______________________________________________ >> Libmesh-devel mailing list >> Lib...@li... >> https://lists.sourceforge.net/lists/listinfo/libmesh-devel > > |