# Simulating evolution of comet

I have to present a working 3D simulation of an arbitrary comet in a close proximity to the Sun or to the Jupiter. In specific, I want to simulate comets D/1993 F2 and C/2011 W3, using the theoretical model I have proposed in my recent reseach. My deadline is 1st April 2014.

As I see it now - I make my comet of tons of small particles, define the physics of these particles, define the external physics and just run the simulation for some period of time. I need to obtain both select data in numerical form, and a 3D visualisation of the whole process.

For now I do not have to do that in Minkowski space, just regular 3D space, albeit in spherical coordinates. Physics wise - thermodynamics, gravity, mechanical properties of materials, solar wind.

I've got:

• Chunk of free time.

• A Windows/Linux laptop. Might have an access to OS X workstation.

• Solid knowledge of mathematics behind these physical phenomena (PDEs, integrals, series, complex functions).

• A little bit of experience in coding - Turbo Pascal 7, C++, Delphi (order from most to less experience).

• Budget of roughly 500 USD - for software, literature etc. I am qualified for student discounts.

• I can not use existing solutions, as I have to define the mathematical model for physics according to my research.

I do not know where to start, and from what do I begin.

Edit according to the input from Ben. I have a developed mathematical model of a system comet-star, or comet-planet (it is simplified comet-sun case). I am looking on how I can make a 3D simulation of this model, with gathering specific data and the visual representation of the processes happening.

I hope this helps. If it is still to broad to be answered, feel free to specify more details.

• Welcome to SciComp.SE! Your question is broad, and it would help if you could be more clear about what exactly you're looking for. Do you already have a model that you're now looking to solve numerically? Are you still developing the model? – Ben Jun 8 '13 at 22:51
• Edited the main post according to your input, Ben. I hope this helps and I have gotten all the points correct and understandable. – user4515 Jun 8 '13 at 23:17
• To clarify: you have a set of differential equations that, when integrated, reproduce the elliptic/parabolic/hyperbolic orbit of the comet? – J. M. Jun 9 '13 at 0:58
• Basically, I have 3 sets of equations that describe external phenomena, 3 sets of equations that describe phenomena inside the comet and a set of equations that describe the orbital motion of the comet. What I want to simulate - comets destruction (and determining where it will be, and which of external phenomena will be main cause for that). Comets destruction is governed by equilibrium between external and internal phenomena. All these equations are in spherical coordinates, and 5 sets out of seven are PDEs that have to be integated during the calculation. – user4515 Jun 9 '13 at 6:54
• For visualisation you could use glowscript or VPython. On the Glowscript page there is a gravitational/collisions simulation for hard spheres. – boyfarrell Jun 9 '13 at 10:03