Regarding implementation of Periodic Boundary Condition (PBC)

59 views
Skip to first unread message

Deepika Kushwah

unread,
Dec 26, 2022, 11:38:18 PM12/26/22
to dea...@googlegroups.com
Hello Everyone,

I am trying to apply PBC for an elasticity problem. I have already checked tutorial step 45..I tried to implement that  procedure (step-45) in my code, but I am not getting any difference in results with and without PBC (I am applying PBC on left(1) and right(2) in x direction see fig). I am attaching here the fig in which the left(1), right(2), bottom(3), top(4) are the boundary_indicators and also the piece of code. 
I think there is some issue with the DoFTools::make_periodicity_constraints() only, please help me to resolve this issue.


fig.png
image.png
image.png
The above results are the same in both cases. That means maybe PBC is not applied here. But what to do so that I can see the result with PBC? Should I change the displacement boundary conditions which are already applied?


Thanks & Regards,
Deepika

**************************************************************************
This e-mail is for the sole use of the intended recipient(s) and may
contain confidential and privileged information. If you are not the
intended recipient, please contact the sender by reply e-mail and destroy
all copies and the original message. Any unauthorized review, use,
disclosure, dissemination, forwarding, printing or copying of this email
is strictly prohibited and appropriate legal action will be taken.
************************************************************************************************

Wolfgang Bangerth

unread,
Dec 28, 2022, 6:27:23 PM12/28/22
to dea...@googlegroups.com

Deepika -- what you have tried already to find out what is going on? Have you
looked at the object that is filled by
DoFTools::make_periodicity_constraints() to see whether the constraints placed
in there are correct, for example? And are you applying the constraints? In
your code snippet, the 'constraint's object is filled by the called to
make_periodicity_constrains(), but then comes the closing } and the
constraints object is destroyed without you ever doing anything with it.

Best
W.

On 12/26/22 21:37, Deepika Kushwah wrote:
> *** Caution: EXTERNAL Sender ***
>
> Hello Everyone,
>
> I am trying to apply PBC for an elasticity problem. I have already checked
> tutorial step 45..I tried to implement that  procedure (step-45) in my code,
> but I am not getting any difference in results with and without PBC (I am
> applying PBC on left(1) and right(2) in x direction see fig). I am attaching
> here the fig in which the left(1), right(2), bottom(3), top(4) are the
> boundary_indicators and also the piece of code.
> I think there is some issue with the DoFTools::make_periodicity_constraints()
> <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fcurrent%2Fdoxygen%2Fdeal.II%2FnamespaceDoFTools.html%23a929249499b1e5624728d212e90a8e037&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C9173f58905d940ef0a1508dae7c42dbc%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638077128203671970%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=xGi8tbTwSTM9DXh0dmM%2Ba3W6ENzkIy64Jg91G%2FzF35o%3D&reserved=0> only, please help me to resolve this issue.
>
>
> fig.png
> image.png
> image.png
> The above results are the same in both cases. That means maybe PBC is not
> applied here. But what to do so that I can see the result with PBC? Should I
> change the displacement boundary conditions which are already applied?
>
>
> Thanks & Regards,
> Deepika
>
> **************************************************************************
> This e-mail is for the sole use of the intended recipient(s) and may
> contain confidential and privileged information. If you are not the
> intended recipient, please contact the sender by reply e-mail and destroy
> all copies and the original message. Any unauthorized review, use,
> disclosure, dissemination, forwarding, printing or copying of this email
> is strictly prohibited and appropriate legal action will be taken.
> ************************************************************************************************
>
> --
> The deal.II project is located at http://www.dealii.org/
> <https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.dealii.org%2F&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C9173f58905d940ef0a1508dae7c42dbc%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638077128203671970%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=1gEwpQ2%2BK7sNa%2FZaQCr74278zTuTg1PbStGmNuNdcIQ%3D&reserved=0>
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fforum%2Fdealii%3Fhl%3Den&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C9173f58905d940ef0a1508dae7c42dbc%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638077128203671970%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=kYXB8TRnMAU%2F%2Bxi85F%2F9MMKqfq%2BT58DuzBGTCJqL3iw%3D&reserved=0>
> ---
> 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
> <mailto:dealii+un...@googlegroups.com>.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/CAKTdrprev1ewczwCrcaVTcctH0Y0ZykyNGfjXJgP%3DnA0aY6L1Q%40mail.gmail.com <https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgroups.google.com%2Fd%2Fmsgid%2Fdealii%2FCAKTdrprev1ewczwCrcaVTcctH0Y0ZykyNGfjXJgP%253DnA0aY6L1Q%2540mail.gmail.com%3Futm_medium%3Demail%26utm_source%3Dfooter&data=05%7C01%7CWolfgang.Bangerth%40colostate.edu%7C9173f58905d940ef0a1508dae7c42dbc%7Cafb58802ff7a4bb1ab21367ff2ecfc8b%7C0%7C0%7C638077128203828204%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=yMCnQK%2FqyUXrblyVEddHLGX3wZtWfMqI5n9gik%2FqQDs%3D&reserved=0>.

--
------------------------------------------------------------------------
Wolfgang Bangerth email: bang...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/


Deepika Kushwah

unread,
Dec 29, 2022, 11:28:46 AM12/29/22
to dea...@googlegroups.com
Yes. I was doing wrong. Now I have corrected it.
Thank you for your suggestion.


Thanks 
Deepika

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://groups.google.com/d/msgid/dealii/24321403-2f36-c487-e312-23b26d689d32%40colostate.edu.
Reply all
Reply to author
Forward
0 new messages