R: Orbital functions

cosorb {celestial} | R Documentation

## 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.

### 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)

[Package celestial version 1.4.1 Index]