def kinematic_viscosity(density, dynamic_viscosity,
kinematic_viscosity_units=default_kinematic_viscosity_units,
density_units=default_density_units,
dynamic_viscosity_units=default_dynamic_viscosity_units):
"""
Return kinematic viscosity given density and dynamic viscosity.
"""
density = U.density_conv(density, density_units, 'kg/m**3')
dynamic_viscosity = U.dynamic_viscosity_conv(dynamic_viscosity, dynamic_viscosity_units, 'Pa s')
kinematic_viscosity = dynamic_viscosity / density
kinematic_viscosity = U.kinematic_viscosity_conv(kinematic_viscosity, 'm**2/s', kinematic_viscosity_units)
return kinematic_viscosity
def kinematic_viscosity(
density,
dynamic_viscosity,
kinematic_viscosity_units=default_kinematic_viscosity_units,
density_units=default_density_units,
dynamic_viscosity_units=default_dynamic_viscosity_units):
"""
Return kinematic viscosity given density and dynamic viscosity.
"""
density = U.density_conv(density, density_units, 'kg/m**3')
dynamic_viscosity = U.dynamic_viscosity_conv(dynamic_viscosity,
dynamic_viscosity_units,
'Pa s')
kinematic_viscosity = dynamic_viscosity / density
kinematic_viscosity = U.kinematic_viscosity_conv(
kinematic_viscosity, 'm**2/s', kinematic_viscosity_units)
return kinematic_viscosity
def alt2kinematic_viscosity(H, T='std', alt_units=default_alt_units, temp_units=default_temp_units, kinematic_viscosity_units=default_kinematic_viscosity_units):
"""
Return kinematic viscosity, given altitude and an optional temperature input.
"""
density = alt2density(H, alt_units, density_units='kg/m**3')
isa_temp = alt2temp(H, alt_units, 'K')
if T != 'std':
isa_temp = alt2temp(H, alt_units, 'K')
T = U.temp_conv(T, temp_units, 'K')
density *= isa_temp / T
else:
T = U.temp_conv(isa_temp, 'K', temp_units)
u = temp2dynamic_viscosity(T, temp_units, 'Pa s')
v = u / density
v = U.kinematic_viscosity_conv(v, 'm**2/s', kinematic_viscosity_units)
return v
def alt2kinematic_viscosity(
H,
T='std',
alt_units=default_alt_units,
temp_units=default_temp_units,
kinematic_viscosity_units=default_kinematic_viscosity_units):
"""
Return kinematic viscosity, given altitude and an optional temperature input.
"""
density = alt2density(H, alt_units, density_units='kg/m**3')
isa_temp = alt2temp(H, alt_units, 'K')
if T != 'std':
isa_temp = alt2temp(H, alt_units, 'K')
T = U.temp_conv(T, temp_units, 'K')
density *= isa_temp / T
else:
T = U.temp_conv(isa_temp, 'K', temp_units)
u = temp2dynamic_viscosity(T, temp_units, 'Pa s')
v = u / density
v = U.kinematic_viscosity_conv(v, 'm**2/s', kinematic_viscosity_units)
return v
