search

EmelineBolmont / mercury-90

Automatically exported from code.google.com/p/mercury-90
3 stars 2 forks source link

implement a more accurate equation for viscosity #4

Closed GoogleCodeExporter closed 9 years ago

GoogleCodeExporter commented 9 years ago 
viscosity is for the moment used as a constant. That is, the viscosity profile 
is assumed to be as steep as (alpha - 0.5) and a common value for alpha is 0.5. 

The value at R=1 itself must be more precise. The current version is determined 
at R=1 with alpha = 0.005 and a alpha-prescription for the viscosity : 
viscosity = alpha * c_s * H (with c_s = H * Omega_p)

thus, with h=0.05, alpha = 0.005 and at R=1 we get : 
viscosity = 6e-10 in (AU^2.day-1)

Original issue reported on code.google.com by ccos...@gmail.com on 8 Aug 2011 at 1:14

GoogleCodeExporter commented 9 years ago 
viscosity has now a alpha-prescription. 

We assume that c_s = H * Omega_p. Thus, the equation become : 
viscosity = alpha * Omega_p * H**2

Original comment by ccos...@gmail.com on 9 Aug 2011 at 11:46