Damped sound wave

[spossann]

• Move the factor 1/Np into the definition of the weights. This affects sph and accumulation routines; it makes the sph eval kernels easier (no need to differentiate between tesselation and monte-carlo), as all info is in the weights.

• Fix bug in mpi_sort_markers: add missing self.update_ghost_particles() to get the correct number of valid_mks.

• Improve unit test for Poisson accumulation.

• New tutorials and verification tests for SPH:

tutorials/tutorial_viscous_euler_sph.ipynb
tutorials/tutorial_velocity_diffusion_sph.ipynb
tutorials/tutorial_hagen_poiseuille_sph.ipynb
tutorials/tutorial_dam_break_sph.ipynb

Run verif tests for example with (run dam break only with 1 MPI proc):

mpirun -n 2 python src/struphy/models/tests/verification/test_verif_ViscousEulerSPH.py

• New files for sph kernels: eval_kernels_sph.py and pusher_kernels_sph.py

• New sections in doc for developerss reference: feec, pic and sph.

• Check out the dam break from the new verification test: https://datashare.mpcdf.mpg.de/s/sNYbZZes5GLXKqx

• Check out the Hagen-Poiseuille flow from the new verification test: https://datashare.mpcdf.mpg.de/s/9s9mgGLG8yzrZmj

[spossann]

Hi @aminrai2000 - please have a look at the new tutorials in this branch:

tutorials/tutorial_viscous_euler_sph.ipynb
tutorials/tutorial_velocity_diffusion_sph.ipynb
tutorials/tutorial_hagen_poiseuille_sph.ipynb
tutorials/tutorial_dam_break_sph.ipynb

You will find some pictures there that you can modify/use. Alternatively, you can run

mpirun -n 2 python src/struphy/models/tests/verification/test_verif_ViscousEulerSPH.py

to run these tests in parallel and generate these figures (change name ==" main" to switch test). We can discuss on Monday if you want - see you!