Re: [deal.II] How to implement an orthotropic material on differnet parts
Thanks a lot for your reply, I'm working on a PHD thesis(Rostock university, Germany) in the domain of nonlinear model order reduction. the thesis is attached to a European project in the field of microelectronics. the board I attached is a graphics board and has three heat sources, 2 memories and 1 CPU under the heat sink. I should solve the transient heat problem after that as step-26 without adaptive mesh. but I wanted to learn and validate first the static case with step-5. I also plan in the future to have a transient coupled thermo-mechanical, I didn't find any step file for such application, any recommendation? I'm sorry if my question wasn't clear, but as you see in the attached photo the heat sources doesn't influence the other board components which seemed not connected together. On Thursday, August 19, 2021 at 4:39:41 AM UTC+2 Wolfgang Bangerth wrote: > On 8/18/21 9:20 AM, ibrahim zawra wrote: > > Thanks a lot that helped, I applied symmetric<2,dim> tensor because it > is a > > heat problem. I also implemented the convection boundary condition. I > have a > > problem in the results, it seems the contact between surfaces is not > detected > > automatically, I want to model such contact.Please I will be grateful > for your > > guidance. > > You will have to explain more what you are doing, what the symptoms are, > what > you have already tried, what worked and what doesn't. We can't see from > just > one picture what might possibly be wrong. > > 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/53c911fb-91ac-4e13-9bd4-cdf559a95632n%40googlegroups.com.
Re: [deal.II] How to implement an orthotropic material on differnet parts
Thanks professor, that helped. It works now. For the record I used Boolean fragments (coherence), please find the updated result figure. [image: updated tempreture.png] On Thursday, August 19, 2021 at 6:57:43 AM UTC+2 Jean-Paul Pelteret wrote: > Hi Ibrahim, > > My guess is that the mesh has been exported without the connectivity > between the parts being enforced. So even though the physical vertices of > the mesh overlap where you have different parts, there are two or more > distinct vertices that lie on top of one another and the “intersecting” > parts of the geometry are therefore disconnected from one another. You can > normally fix this within the meshing software by using an “Imprint” > operation (it might be called something else in your specific software). > > Another approach which might work (depending on which finite elements you > are using) is that you can use the AffineConstraints to effectively tie the > two overlapping degrees of freedom from each region together. But you’d > have to manage this yourself. You can look at > > https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-do-i-get-the-degree-of-freedom-indices-at-vertices > to see how you can get the indices for DoFs with support at vertices, and > these functions might also be of use in this regard: > > https://dealii.org/developer/doxygen/deal.II/namespaceDoFTools.html#a5514e4f59ea659f63953d62ca429eaff > > https://dealii.org/developer/doxygen/deal.II/namespaceDoFTools.html#a06b3c33925c1a1f15de20deda20b4d21 > > https://dealii.org/developer/doxygen/deal.II/namespaceDoFTools.html#a8b97e816b29ecf963370a9d8b349828f > > I hope that this helps, > Jean-Paul > > On 19. Aug 2021, at 04:39, Wolfgang Bangerth > wrote: > > On 8/18/21 9:20 AM, ibrahim zawra wrote: > > Thanks a lot that helped, I applied symmetric<2,dim> tensor because it is > a heat problem. I also implemented the convection boundary condition. I > have a problem in the results, it seems the contact between surfaces is not > detected automatically, I want to model such contact.Please I will be > grateful for your guidance. > > > You will have to explain more what you are doing, what the symptoms are, > what you have already tried, what worked and what doesn't. We can't see > from just one picture what might possibly be wrong. > > 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+un...@googlegroups.com. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/d8166e4d-06e6-7da4-e258-9c936fed2ce8%40colostate.edu > . > > > -- 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/8e58c07f-2747-4efa-bb2c-f53b5a5893afn%40googlegroups.com.
Re: [deal.II] Re: MPI error while running inside a Docker container
The error appears inside write_vtu_in_parallel, which uses MPI I/O. It might be that the type of filesystem (docker shared volume) does not allow parallel IO or that the docker installation does not have MPI IO enabled. On Tue, Aug 17, 2021 at 10:15 AM Bruno Turcksin wrote: > > David, > > It depends where the build files are written. Do you see the build > files on your computer or do they exist only in the Docker container, > i.e, after you've killed Docker do the files still exist somewhere on > your computer ? > > Best, > > Bruno > > Le mar. 17 août 2021 à 09:00, David Montiel Taboada > a écrit : > > > > Hello, Bruno. > > > > Thank you, I will check. I assumed that was not the issue because when I > > compile inside docker files are created. > > > > David > > > > > > On Tue, Aug 17, 2021, 4:30 AM Bruno Turcksin > > wrote: > >> > >> David, > >> > >> Have you try to create a new file and save it inside docker? What is your > >> user name in docker? Are you root? It could be a problem with your > >> permission. > >> > >> Best, > >> > >> Bruno > >> > >> On Saturday, August 7, 2021 at 3:07:56 PM UTC-4 dmon...@umich.edu wrote: > >>> > >>> Sorry in my original email I wrote "This appears to happen if I try to > >>> run the code in serial." > >>> > >>> I meant it happens whether I try to run the code in serial or parallel > >>> (e.g., using "./main" or "mprun -n 4 main"). It makes no difference. > >>> > >>> Thank you, > >>> > >>> David > >>> > >>> On Sat, Aug 7, 2021 at 2:57 PM David Montiel Taboada > >>> wrote: > > This is how I am launching the container, by the way: > > $ docker run -ti -v > ~/DockerWorkspace/phaseField/applications:/home/dealii/phaseField/applications > prismspf/prismspf:2.2 > > > > On Sat, Aug 7, 2021 at 2:55 PM David Montiel Taboada > wrote: > > > > Hello, > > > > I am running a container of PRISMS-PF (which is built upon the > > dealii:v9.2.0-focal image). > > > > Inside the container, the code compiles successfully, but when I run > > it, I get the following error message: > > > > > > PRISMS-PF: Exception on processing: > > > > > > An error occurred in line <6909> of file > > in > > function > > void dealii::DataOutInterface > spacedim>::write_vtu_in_parallel(const string&, MPI_Comm) const [with > > int dim = 2; int spacedim = 2; std::string = > > std::__cxx11::basic_string; MPI_Comm = ompi_communicator_t*] > > The violated condition was: > > ierr == MPI_SUCCESS > > Additional information: > > deal.II encountered an error while calling an MPI function. > > The description of the error provided by MPI is "MPI_ERR_ACCESS: invalid > > access mode". > > The numerical value of the original error code is 20. > > > > > > Aborting! > > > > > > This appears to happen if I try to run the code in serial. > > > > However I only have the issue on a Linux machine (I have no problem if > > I run on a Mac or Windows machine, for instance) which leads me to > > believe it may be a conflict between the Docker container and the > > underlying OS. > > > > Has anyone encountered something similar? Any idea of what may be going > > on? > > > > Thank you! > > > > David > >> > >> -- > >> 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/368cf5ae-ac31-4883-91f2-8ef1b794a684n%40googlegroups.com. > > > > -- > > 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 a topic in the > > Google Groups "deal.II User Group" group. > > To unsubscribe from this topic, visit > > https://groups.google.com/d/topic/dealii/wkK_L6a3PPM/unsubscribe. > > To unsubscribe from this group and all its topics, send an email to > > dealii+unsubscr...@googlegroups.com. > > To view this discussion on the web visit > > https://groups.google.com/d/msgid/dealii/CAJdNbL-JNq-gGJo8F4TT4_YDmz_jYoDjQsqs4GF-UAX7SnzOjQ%40mail.gmail.com. > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/
Re: [deal.II] Exception handling for multithreaded case
Can you illustrate in a small test case what you are trying to do? Sorry, it was my mistake. Just to reassure myself, such a stack-trace is only printed for Assert() and not for AssertThrow(). I realised this while creating the attached MWE. Interesting. I had no idea that there was a difference, and why. I'm tracking this here: https://github.com/dealii/dealii/issues/12693 I think I know how to fix this, but you'll only get to see this in the next release (or with developer versions), of course. The object could be a scalar or a dealii::Tensor<1,dim,double>. I was looking for something on the lines of #pragma omp for using globalSum as a shared variable and a '+' reduction operation. You can try to use std::atomic and Tensor<1,dim,std::atomic> To me it seemed that there would be an openmp equivalent procedure for such simple reduction tasks. Even standard C++ :-) Best Wolfgang -- Wolfgang Bangerth email: bange...@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/0966fcc3-1ed3-a6c5-cb3e-9c0f1173edb1%40colostate.edu.
