I successfully implemented parallel::fullydistributed::Triangulation in my code. The way it works right now is it requires an initial triangulation to generate the TriangulationDescription::Description< dim, dim > using create_description_from_triangulation. It means that I should first generate the triangulation. This still imposes the memory constraint on the process as you already generated the large mesh for each processor. I tried it and still require the same memory.
Is there any way that while we use GridGenerator, it can directly generate the mesh and distribute it over nodes using p::f::T? Any recommendation is highly appreciated. Thanks Reza On Thursday, April 22, 2021 at 12:55:47 PM UTC-4 Wolfgang Bangerth wrote: > On 4/22/21 10:52 AM, Mohammadreza ya'ghoobi wrote: > > > > I'm trying to solve a simpler problem of uniaxial response of a plate > with > > hole using FE_Nothing as you recommended. I re-read steps 10, 27, and 46 > to > > refresh my memory. I think I have the logic but I appreciate if you > check some > > of the steps I have a little concern about. The two main step I took are: > > 1) I assign the cells inside the hole as FE_Nothing and the rest as > regular FE. > > 2) During the stiffness matrix and residual force formation, I only go > over > > the elements in regular FE part and exclude the FE_Nothing element from > the > > calculation. > > > > My only concern are the cells on the interface. Does performing assembly > and > > residual formation just on the FE cells (and exclude FE_Nothing cells) > > guarantee the traction free surface of the cells on the surface of the > hole? > > Or I need to take extra action to guarantee the traction free surface of > the hole? > > That depends on how you create the FE_Nothing element. Its constructor > takes a > flag that determines what the solution should look like on the other side: > Should it be zero at the interface, or be unconstrained? You want it to be > unconstrained -- it's a boundary at which you intend to do nothing special > (which will then correspond to zero-traction Neumann-type boundary > condition). > > 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/9510882b-8243-413b-ac08-d6a9dbedcc58n%40googlegroups.com.
