Hi Wolfgang, Thanks for the feedback! I have also not found any literature dealing with this periodic boundary conditions problem -- I doubt I will try to get to the bottom of it experimentally unless I run into technical issues solving my particular problem. With both the lifting method and interior penalty method, it appears that you need knowledge of adjacent cells in assembly; I suspect this will be difficult (or at least annoying) with periodic boundaries if it's not already implemented in deal.II.
For the triharmonic equation, this is the only thing I've been able to find <https://arxiv.org/pdf/1804.03793.pdf>. It seems quite similar to step-82. In any case, I will try the method of introducing auxiliary variables (so that it's just a bunch of coupled Laplace equations). Thanks again for the help! - Lucas On Monday, March 20, 2023 at 3:35:30 PM UTC-5 Wolfgang Bangerth wrote: > > Lucas: > > > these two tutorials, I didn't see anything about periodic boundary > > conditions, however. Two potential ways I can think to enforce periodic > > boundary conditions would be (in 1D, for notational simplicity): > > > > 1. u(0) = u(L), du/dx (0) = du/dx (L) > > 2. u(0) = u(L), d^2 u/dx^2 (0) = d^2 u/dx^2 (L) > > > > For the former case, it appears that you would have to extend the method > > from step-47 or step-82 to accommodate these conditions somehow. In the > > latter case, it seems feasible that you could write it as two coupled > > Laplace equations (supposing something like what is explained in the > > step-47 note holds here). > > > > My question is: is that right? And if so, how do the solutions gotten by > > imposing these boundary conditions differ? For an analytic solution, it > > seems like these would be equivalent (periodic boundary conditions > > assert periodicity of derivatives at all orders, yes?) so I'm not sure > > what to make of it. > > Like you, I suspect that you can do it either way to achieve > periodicity, and that you can choose the way that is most convenient for > your formulation. But I am not certain, nor do I know literature that > would have covered the question. It is, however, something that can be > experimentally determined via the method of manufactured solutions. > > I do think that it is not crazy to go with the assumption that either > approach mentioned above will work. > > > > Just as an aside, the actual problem that I'm trying to solve involves a > > triharmonic operator. I suspect the situation is similar (although of > > course the methods in the tutorials wouldn't immediately generalize), > > but if there's something qualitatively different that's obvious, I'd > > like to know about that. > > I don't think there is a good reason to believe that the approaches for > the biharmonic and triharmonic equations should be substantially different. > > I do think that it would be fun to have your triharmonic program be part > of the code gallery one day! > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > -- 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/9bb52d5a-709a-4aa1-b7d2-008981d2bf11n%40googlegroups.com.
