cosorb: Orbital functions

In celestial: Collection of Common Astronomical Conversion Routines and Functions

Description

A variety of obital analysis functions. These are useful for setting up initial conditions for merging systems etc.

Usage

  cosorbVisViva(M=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)

Arguments

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

Details

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.

Value

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.

Author(s)

Aaron Robotham, Chris Power

See Also

cosvol, cosmap, cosdist, cosgrow

Examples

  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)

celestial documentation built on May 2, 2019, 3:05 p.m.