Hassan,

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?

Yes. We have certainly done simulations with many more particle properties.


2- I want to have a particle at each vertex on the triangulation. How to create them?

Take a look at some of the particle tutorials that show how to create particles at vertex locations.


3- How to create and apply constraints?

You'll have to be more specific. Constraints on what?


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).

You'll have to have data structures that allow you to identify which particles are close to which others. The usual one is an rtree, for which deal.II has wrappers.

The only deal.II-based code I know of that does these sorts of things is Bruno Blais' 'lethe' code. You should take a look at how this is done inside lethe.

Best
 W.

--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 [email protected]
                           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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/919948b8-3046-5f53-acc7-a7f1c51511fd%40colostate.edu.

Reply via email to