celestial (version 1.4.6)

cosorb: Orbital 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.

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.

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.

See Also

cosvol, cosmap, cosdist, cosgrow

Examples

# NOT RUN {
  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)
# }