Hi Jack, The message that you’re seeing after fixing your first problem using Daniel’s advice is being emitted from your own code. Specifically, the set_neumann_bc() function in the class RMHD::Entities::VectorBoundaryConditions. You’ll have to inspect your code more carefully to understand what’s going wrong here.
Best, Jean-Paul > On 8. Jul 2022, at 22:41, jack urombo <juro...@hit.ac.zw> wrote: > > After following Dan's advice the following exception is now thrown > > ---------------------------------------------------- > Exception on processing: > -------------------------------------------------------- > An error occurred in line <595> of file > </home/dealii/rotmhd/source/boundary_conditions.cc> in funct > ion > void > RMHD::Entities::VectorBoundaryConditions<dim>::set_neumann_bc(dealii::types::boundary_id, > c > onst std::shared_ptr<dealii::TensorFunction<1, dim> >&, bool) [with int dim = > 3; dealii::types::boun > dary_id = unsigned int] > The violated condition was: > !this->closed() > Additional information: > The boundary conditions have already been closed > -------------------------------------------------------- > Aborting! > ---------------------------------------------------- > > > > On Friday, 8 July 2022 at 00:17:17 UTC+2 d.arnd...@gmail.com wrote: > Jack, > > step-45 https://www.dealii.org/current/doxygen/deal.II/step_45.html > <https://www.dealii.org/current/doxygen/deal.II/step_45.html> provides an > example for dealing with periodic boundary conditions. Even though it's a 2D > example, the interface is the same in 3D. > > Assuming that the arguments for your set_periodic_boundary_condition function > are similar to the ones for GridTools::collect_periodic_faces and you use a > unit cube with standard coloring, you probably want to call it with arguments > 0,1,0 and 2,3,1 respectively. > > Best, > Daniel > > On Thu, Jul 7, 2022 at 1:56 PM jack urombo <jur...@hit.ac.zw > <applewebdata://D38CD938-1D7B-4936-BC48-F6432469C554>> wrote: > I am trying to set periodic boundary conditions for a 3D problem like > > velocity->set_periodic_boundary_condition(0, 1, 2); > velocity->set_periodic_boundary_condition(2, 3, 1); > pressure->set_periodic_boundary_condition(0, 1, 2); > pressure->set_periodic_boundary_condition(2, 3, 1); > > The boundary conditions are not matching and that i triggering errors. > > How do you set the period BCs in 3D. > > The information in this message is confidential and legally privileged. It is > intended solely for the addressee(s). Access to this message by anyone else > is unauthorized. If received in error, please accept our apologies and notify > the sender immediately. You must also delete the original message from your > machine. If you are not the intended recipient, any use, disclosure, copying, > distribution or action taken in reliance of it, is prohibited and may be > unlawful. The information, attachments, opinions or advice contained in this > email are not the views or opinions of Harare Institute of Technology, its > subsidiaries or affiliates. Although this email and any attachments are > believed to be free of any virus or other defects which might affect any > computer or IT system into which they are received, no responsibility is > accepted by Harare Institute of Technology and/or its subsidiaries for any > loss or damage arising in any way from the receipt or use thereof. > > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <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+un...@googlegroups.com > <applewebdata://D38CD938-1D7B-4936-BC48-F6432469C554>. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/5eaf4f26-00dd-4e6d-841d-4a22e2328fb1n%40googlegroups.com > > <https://groups.google.com/d/msgid/dealii/5eaf4f26-00dd-4e6d-841d-4a22e2328fb1n%40googlegroups.com?utm_medium=email&utm_source=footer>. > > The information in this message is confidential and legally privileged. It is > intended solely for the addressee(s). Access to this message by anyone else > is unauthorized. If received in error, please accept our apologies and notify > the sender immediately. You must also delete the original message from your > machine. If you are not the intended recipient, any use, disclosure, copying, > distribution or action taken in reliance of it, is prohibited and may be > unlawful. The information, attachments, opinions or advice contained in this > email are not the views or opinions of Harare Institute of Technology, its > subsidiaries or affiliates. Although this email and any attachments are > believed to be free of any virus or other defects which might affect any > computer or IT system into which they are received, no responsibility is > accepted by Harare Institute of Technology and/or its subsidiaries for any > loss or damage arising in any way from the receipt or use thereof. > > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <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 > <mailto:dealii+unsubscr...@googlegroups.com>. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/f661ac5e-b073-4107-b086-bc0ee89dd9f3n%40googlegroups.com > > <https://groups.google.com/d/msgid/dealii/f661ac5e-b073-4107-b086-bc0ee89dd9f3n%40googlegroups.com?utm_medium=email&utm_source=footer>. -- 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/5DC98595-1249-4078-A88E-E56EDE80BC99%40gmail.com.
