Are there any simulations that can be obtained that use smooth particle hydrodynamics and can be configured to include different initial conditions? I wish to simulate planetary collisions and their impact. I wish to exactly replicate the following link.


  • $\begingroup$ It would be awesome if you could tell a bit more precisely, what do you really wish for. There are plenty of algorithms, inlcuding SPH and its variations, many codes written and many nice outputs produced, planetary collisions included. $\endgroup$ – Alexey Bobrick Jan 17 '14 at 18:01
  • 1
    $\begingroup$ The algorithm that I am looking for should be able to simulate collisions like the following link : youtube.com/watch?v=Fwl_JBQtH9o $\endgroup$ – Artemisia Jan 18 '14 at 2:44
  • $\begingroup$ I am trying to exactly duplicate the code used for that YouTube link. $\endgroup$ – Artemisia Jan 18 '14 at 2:45
  • $\begingroup$ Could you find Canup's SPH code? I'm a PhD Brazilian student and I'm looking for a code like you.. Thanks =) $\endgroup$ – Andre Amarante May 24 '17 at 6:22
  • $\begingroup$ @AndreAmarante, unfortunately, I was unable to find Canup's code. Other SPH codes did provide a lot of insight into how to figure it out, so I would recommend checking those out. Cheers! $\endgroup$ – Artemisia Jun 5 '17 at 23:55

I would recommend MPA Garching's Gadget code for cosmological simulations of structure formation. It's primarily gravitational, but I do believe you can include gas effects as well:

GADGET computes gravitational forces with a hierarchical tree algorithm (optionally in combination with a particle-mesh scheme for long-range gravitational forces) and represents fluids by means of smoothed particle hydrodynamics (SPH). The code can be used for studies of isolated systems, or for simulations that include the cosmological expansion of space, both with or without periodic boundary conditions. In all these types of simulations, GADGET follows the evolution of a self-gravitating collisionless N-body system, and allows gas dynamics to be optionally included. Both the force computation and the time stepping of GADGET are fully adaptive, with a dynamic range which is, in principle, unlimited.

I don't know if this is what is commonly used for tasks such as yours, but it might be a place to start looking.

  • $\begingroup$ Thanks :) it is a good headstart. I wish to exactly simulate youtube.com/watch?v=Fwl_JBQtH9o $\endgroup$ – Artemisia Jan 18 '14 at 2:46
  • 1
    $\begingroup$ There many different SPH codes around the astronomy community (I have one myself), but AFAIK none except gadget are publicly available for download (simply because these codes are very specific and easy to break if you don't know what you're doing). The problem with gadget is that it implements an out-of-date SPH algorithm (>10 years old) and was never really meant for planet collisions (I bet you would need a particular equation of state for planet collisions). $\endgroup$ – Walter Jan 18 '14 at 10:13
  • $\begingroup$ Ah yes. I plan to use the ANEOS analytical equations of state. $\endgroup$ – Artemisia Jan 18 '14 at 10:22
  • $\begingroup$ I would like a simulation with a more updated version of the SPH algorithm, no problem if I have to purchase it :) Any suggestions? $\endgroup$ – Artemisia Jan 18 '14 at 10:23
  • 1
    $\begingroup$ @Walter: there are many other good publically available SPH codes. Check any recent review on SPH, e.g. Springel, Price or Rosswog. $\endgroup$ – Alexey Bobrick Jan 19 '14 at 1:04

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.