Thank you Profs. Prof. Bangerth, I probably wasn't clear in my question. I know how to apply BCs for an interior penalty problem weakly. But, for a non-linear problem like step 15, one initialises the solution vector and sets it to the value of the BCs at the boundary before solving for the Newton update. This is done explicitly in function 'set_boundary_values()' in step 15 with the function ' VectorTools::interpolate_boundary_values <https://www.dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#ab2562d41bb26f362043f9719a8cd9b87> '. I cannot use this function with DG elements apparently and I am wondering if there are any alternatives.
Prof. Hiester, don't you think that subtracting small and manipulating large numbers in the residual will cause some problems? I'll give it a try for sure though. Thanks again On Monday, December 5, 2022 at 3:01:32 AM UTC+2 Timo Heister wrote: > I would consider the boundary conditions part of the nonlinear > iteration (in contrast to step-15, where they are strongly enforced > and as such are only needed in step 0). Every step you can evaluate > your nonlinear residual which contains a residual in the boundary > condition and that difference needs to be applied in every step > (assuming you solve for an update in each Newton step). > > On Sat, Dec 3, 2022 at 10:56 AM Wolfgang Bangerth > <bang...@colostate.edu> wrote: > > > > On 12/3/22 07: 11, Abbas wrote: > > Something like this doesn't work > with DG. > Loosely speaking, my initial approach would be to solve a system > with just the > BC terms but I am not sure. Do I have other options? In DG > methods, > > ZjQcmQRYFpfptBannerStart > > This Message Is From an External Sender > > Use caution when opening links or attachments if you do not recognize > the sender. > > > > ZjQcmQRYFpfptBannerEnd > > > > On 12/3/22 07:11, Abbas wrote: > > > > > > Something like this doesn't work with DG. > > > Loosely speaking, my initial approach would be to solve a system with > just the > > > BC terms but I am not sure. Do I have other options? > > > > In DG methods, you impose boundary values weakly, using the same > strategy with > > which you impose continuity between cells weakly. You might want to look > at > > papers on DG methods: pretty much every single one will show the jump > terms > > corresponding to Dirichlet boundary values. > > > > Best > > W. > > > > -- > > ------------------------------------------------------------------------ > > Wolfgang Bangerth email: bang...@colostate.edu > > www: > https://urldefense.com/v3/__http://www.math.colostate.edu/*bangerth/__;fg!!PTd7Sdtyuw!Vf-IQx7HQQ5bIq8xBQu1BzCdpZiz9Q3UC6emGlg1Z50xBT-oA9UeUBtyamyADtrTV0zLEZpa0RBQ_lV6jh5B8Uo$ > > > > > > -- > > The deal.II project is located at > https://urldefense.com/v3/__http://www.dealii.org/__;!!PTd7Sdtyuw!Vf-IQx7HQQ5bIq8xBQu1BzCdpZiz9Q3UC6emGlg1Z50xBT-oA9UeUBtyamyADtrTV0zLEZpa0RBQ_lV6CDq7NmU$ > > For mailing list/forum options, see > https://urldefense.com/v3/__https://groups.google.com/d/forum/dealii?hl=en__;!!PTd7Sdtyuw!Vf-IQx7HQQ5bIq8xBQu1BzCdpZiz9Q3UC6emGlg1Z50xBT-oA9UeUBtyamyADtrTV0zLEZpa0RBQ_lV6qJBbcDU$ > > --- > > 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+un...@googlegroups.com. > > To view this discussion on the web visit > https://urldefense.com/v3/__https://groups.google.com/d/msgid/dealii/c1d413d2-3c55-0523-19fd-9f399d1781da*40colostate.edu__;JQ!!PTd7Sdtyuw!Vf-IQx7HQQ5bIq8xBQu1BzCdpZiz9Q3UC6emGlg1Z50xBT-oA9UeUBtyamyADtrTV0zLEZpa0RBQ_lV6CppyxvM$ > > . > > > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ > -- 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/e6b575d7-9bfd-4137-b30b-3e3d2a335144n%40googlegroups.com.
