chandrasekhar_dynamical_friction¶
-
class
eridanus_pbh.potential.dynamical_friction.chandrasekhar_dynamical_friction(GMs, rhm, dens, rmin=quantity<0.0008 kpc>, time=quantity<0 Myr>)[source]¶ Bases:
eridanus_pbh.potential.dynamical_friction.DynamicalFrictionClassWrapper for to_amuse(potential.ChandrasekharDynamicalFrictionForce).
Instead of using position and velocity of individual star in get_gravity_at_point the code use self.position and self.velocity which would be set externally as the cluster’s centre of mass.
chandrasekhar_dynamical_friction.
- Parameters
- GMs: float or astropy quantity
- rhm: float
half-mass radius
- dens: float
density
- rmin: amuse Quantity
minimum r turns off dynamical friction when r < rmin
- time: amuse Quantity
the time, not used
- .. todo::
have GMs, rhm, and dens be @properties read from the CDF
Methods Summary
circular_velocity(self, r)circular_velocity.
enclosed_mass(self, r)enclosed_mass.
evolve_model(self, t_end)evolve_model.
get_gravity_at_point(self, eps, x, y, z)get_gravity_at_point.
get_potential_at_point(self, eps, x, y, z)get_potential_at_point.
mass_density(self, eps, x, y, z)mass_density.
reset_parameters(self, \*[, GMs, rhm, dens])reset_parameters.
set_rv(self, x, y, z, vx, vy, vz)set_rv.
stop(self)stop.
Methods Documentation
-
circular_velocity(self, r)[source]¶ circular_velocity.
- Parameters
- r: float
- Returns
- circular_velocity: float
-
get_gravity_at_point(self, eps, x, y, z)[source]¶ get_gravity_at_point.
- Parameters
- eps: float
- x, y, z: float
- Returns
- ax, ay, az: float
-
get_potential_at_point(self, eps, x, y, z)[source]¶ get_potential_at_point.
- Parameters
- eps: float
- x, y, z: float
- Returns
- potential: float
-
mass_density(self, eps, x, y, z)[source]¶ mass_density.
- Parameters
- eps: float
- x, y, z: float
- Returns
- mass_density: float
-
reset_parameters(self, \*, GMs=None, rhm=None, dens=None)[source]¶ reset_parameters.
- Parameters
- GMs: float
- rhm: float
- dens: float