R/cosorb.R

In asgr/celestial: Collection of Common Astronomical Conversion Routines and Functions

cosorbVisViva=function(mass=1e12,Rad=162.635,SemiMajRad=162.635,Munit=1,Lunit=1e3,Vunit=1){
  G=6.67384e-11
  msol_to_kg=1.98892e30
  pc_to_m=3.08568e16
  kms_to_ms=1e3
  g = G*msol_to_kg/(pc_to_m*kms_to_ms^2)
  g = g*Munit/(Lunit*Vunit^2)
  return(sqrt(g*mass*(2/Rad-1/SemiMajRad)))
}

cosorbFreeFall=function(M1=1e12,M2=1,Rad=162.635,Munit=1,Lunit=1e3,Vunit=1,Tunit=1e9){
  G=6.67384e-11
  msol_to_kg=1.98892e30
  pc_to_m=3.08568e16
  kms_to_ms=1e3
  g = G*msol_to_kg/(pc_to_m*kms_to_ms^2)
  g = g*Munit/(Lunit*Vunit^2)
  time=(pi/sqrt(g*(M1+M2)))*(Rad/2)^(3/2)*(Lunit/1e3)
  time=time*1e9/(Tunit*0.9778139)
  return(time)
}

cosorbRocheRad=function(M1=1e12,M2=1e10,Size=35.03865,Rfac=2.44){
  return(Rfac*Size*(M1/M2)^(1/3))
}

cosorbRocheSize=function(M1=1e12,M2=1e10,Rad=396.8294,Rfac=2.44){
  return((Rad/Rfac)*(M2/M1)^(1/3))
}

kinetic_part=function(vel, mass=NULL, Vunit = 1){
  
  if(length(mass) > 1){
    com = c(sum(vel[,1] * mass) / sum(mass),
            sum(vel[,2] * mass) / sum(mass),
            sum(vel[,3] * mass) / sum(mass)
            )
    
    vel[,1] = vel[,1] - com[1]
    vel[,2] = vel[,2] - com[2]
    vel[,3] = vel[,3] - com[3]
  }
  
  return((Vunit^2) * (vel[,1]^2+vel[,2]^2+vel[,3]^2)/2)
}