cosorb | R Documentation |
A variety of obital analysis functions. These are useful for setting up initial conditions for merging systems etc.
cosorbVisViva(mass=1e12, Rad=162.635, SemiMajRad=162.635, Munit=1, Lunit=1e3, Vunit=1)
cosorbFreeFall(M1=1e12, M2=1, Rad=162.635, Munit=1, Lunit=1e3, Vunit=1, Tunit=1e9)
cosorbRocheRad(M1=1e12, M2=1e10, Size=35.03865, Rfac=2.44)
cosorbRocheSize(M1=1e12, M2=1e10, Rad=396.8294, Rfac=2.44)
potential_part(part, eval = NULL, mass = 1, soft = 1, Munit = 1, Lunit = 1000,
Vunit = 1)
kinetic_part(vel, mass=NULL, Vunit = 1)
mass |
Mass in units of 'Munit'. |
M1 |
Mass of primary body in units of 'Munit'. |
M2 |
Mass of secondary body in units of 'Munit'. |
Rad |
Separation between bodies in units of 'Lunit' (for cosorbRocheSize this is in arbitrary units). |
part |
Cartesian particle positions [x,y,z] in units of 'Lunit'. |
eval |
Cartesian evalution positions [x,y,z] in units of 'Lunit'. If this is not provided then the potentials are evaluated at the positions of part. |
vel |
Cartesian particle velocities [vx,vy,vz] in units of 'Vunit'. |
soft |
Numeric softening in units of 'Lunit'. |
SemiMajRad |
The semi major radius of the orbit (a > 0 for ellipses, a = Rad for circles, 1/a = 0 for parabolas, and a < 0 for hyperbolas). |
Size |
The size radius of the secondary object. Inside of this radius the object is bound to the secondary, outside of this radius the object is stripped by the primary. |
Rfac |
The Roche factor. Approximately taken to be 2.44, but in reality it varies depending on the shape of the potentials etc. |
Munit |
Base mass unit in multiples of Msun. |
Lunit |
Base length unit in multiples of parsecs. |
Vunit |
Base velocity unit in multiples of km/s. |
Tunit |
Base time unit in multiples of years. |
These functions allow for various analytic conversions between the 3 major properties related to virial radius: the mass, velocity dispresion and size. The default properties calculate properties for 1e12 Msun halos and assume masses in Msun, velocities in km/s and distances in Kpc.
cosorbVisViva function gives the required velocity in units of Vunit to create the specified orbit.
cosorbFreeFall function gives the free fall time to static initial velocity separated bodies.
cosorbRocheRad function gives the orbital radius at which the secondary will become stripped within a specified bound radius.
cosorbRocheSize function gives the limiting bound radius of the secondary for a given system.
potential_part functions gives the total gravitational potential either at the positions of particles provided (part) or at the requested evaluation positions (eval).
kinetic_part functions gives the total kinetic energy of particles provided (vel)
Aaron Robotham, Chris Power
cosvol
, cosmap
, cosdist
, cosgrow
cosorbVisViva(M=1e15, Rad=1, Lunit=1e6)
cosorbFreeFall(M1=1e15, M2=1, Rad=1, Lunit=1e6)
cosorbRocheRad(M1=1e12, M2=1e12, Size=162.635, Rfac=2.44)
cosorbRocheSize(M1=1e12, M2=1e12, Rad=396.8294, Rfac=2.44)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.