Dear Hassan, There is no way to do that directly in Deal.II Lethe adds two classes to achieve this using a coarse and then a fine contact detection. https://github.com/lethe-cfd/lethe/blob/master/include/dem/particle_particle_broad_search.h https://github.com/lethe-cfd/lethe/blob/master/include/dem/particle_particle_fine_search.h
If you are interested, we could eventually try to port them to deal.II, but these are very lethe specific at the moments. Otherwise, you can use these ideas to generate your contact list. These will work if you use the background mesh as your contact detection grid (which is what we do, it's quite efficient and robust). On Saturday, February 12, 2022 at 3:59:16 a.m. UTC-5 hassan...@gmail.com wrote: > Dear Wolfgang and Bruno > > thank you for your answers. I will study the 'lethe' code. > I want to apply constraints on linear momentum for a set of particles. > Actually, I want to build an SPH code from scratch. We are working with a > total Lagrangian formulation, so position of particles is constant. > Additionally. we are working with just a layer of neighbor particles in > each direction. I did't understand how to find neighbor particles, I was > looking for a way to access neighbor particles, something like the one that > is used for the cell in step-68 "particle->get_surrounding_cell()". > > Thanks > Hassan > > On Saturday, February 12, 2022 at 6:27:09 AM UTC+3:30 blais...@gmail.com > wrote: > >> If I may suggest, at that point, I would ask you what you are trying to >> specifically achieve? >> If you want to build an SPH code from scratch, this is definitely >> possible to do using deal.II particle classes (esp. with the exchange and >> update ghosts functionalities), but this will require a lot of work. In >> that case, it might also be a good idea to look at an SPH code designed >> from scratch for this (e.g. DualSPHPhysics is a good one) >> In our case, we built our DEM code using deal.II because >> A) DEM is simpler than SPH (shorter scale interaction) >> B) We wanted to couple our DEM to our CFD and deal.II was the perfect >> vehicle for that. >> >> It is definitely doable to build from scratch, but it will be a massive >> endeavour to make it fast. >> >> >> >> On Thursday, February 10, 2022 at 3:34:13 a.m. UTC-5 hassan...@gmail.com >> wrote: >> >>> Dear dealII users >>> >>> I want to solve a series of conservation law using smooth particles >>> hydrodynamics (SPH). I am looking for a parallel solution. could you >>> please help me to setup my code correctly. >>> >>> 1- The number of unknowns is 11 (linear momentum 3 components and >>> deformation gradient 9 components). In addition to that I have to save >>> stress tensor, smoothing length, pressure and some other state variables at >>> each particles. So the number of variable that I want to store is huge. >>> Should I store them as 'property' at each particle? >>> >>> 2- I want to have a particle at each vertex on the triangulation. How to >>> create them? >>> >>> 3- How to create and apply constraints? >>> >>> 4- Finally, (the most important question), in a solution using the SPH >>> method computation is done on each particle (Target particle) according to >>> data belong to the particles are around the target particle (Neighbor >>> particles) . How to access to the data that belong to the neighbor >>> particle. For simplicity at this moment we can consider that the neighbor >>> particles are just the particles connected to the target particle on the >>> adjacent cell (one layer). >>> >>> Beta regards, >>> Hassan >>> >> -- 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/0c9e683f-a186-47df-aee8-486611818ebbn%40googlegroups.com.
