I'll never admit to how long I stared at this to find the error, but make sure you're distributing in the correct place in your assembly.
Thanks to anyone who took a peek! Cheers, Aiden On Wednesday, December 6, 2023 at 4:52:02 PM UTC-5 Aiden Huffman wrote: > Hi Everyone, > > I want to thank Luca for his reponse to an earlier post, helped me figure > out a bit about inheritence in C++. > > I am having issues with a new iteration on a Cahn-Hilliard solver I wrote. > I am trying to use the block Schur decomposition illustrated in several > tutorials but I get this strange asymmetric error. I am implementing a > linearly stabilised splitting scheme from here > <https://www.cambridge.org/core/services/aop-cambridge-core/content/view/0672FBD318BBE2621A51AE0F2C9C2FE3/S1815240600003169a.pdf/div-class-title-phase-field-models-for-multi-component-fluid-flows-div.pdf> > (sec > 3.2.5): > [image: Screenshot from 2023-12-06 16-17-58.png] > where F'(phi^n) = (phi^n)^3 - 2phi^n and yields the following block > structure > [image: Screenshot from 2023-12-03 18-04-56.png] > However, I get this strange accumulating error along the diagonal which > the -1 phase drifting down and to the left. I've included a quick snapshot > of the asymmetric error at the phase transition: > [image: cahn-hilliard-sol-snap1.png] [image: cahn-hilliard-sol-snap2.png] > I thought I was incorrectly applying boundary conditions when doing the > Schur complement but trying to solve the system directly as I've done in > the current version of main.cpp resulted in the same error. > > I've also tried forgoing the finite element system completely but I am > having trouble correctly applying the boundary conditions through the > linear operator classes. So I figured I should stick to one question, on my > most successful implementation. > > Any suggestions on how to debug this, or if people have had similar issues > working with finite element systems would be awesome > > Kind Regards, > Aiden > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/298368e6-dd34-48f2-b9ae-9bd5f5de92dan%40googlegroups.com.
