Код: Выделить всё
import pint
pq = pint.UnitRegistry()
density_water = (1 / 1.002)* 1e3 * pq.kg / pq.m**3
dynamic_viscocity_water = 959e-6 * pq.N * pq.s / pq.m**2
tube_inner_diameter = 2e-3 * pq.m
Prandtl_number = 6.62 # Dimensionless
mean_velocity_water = 1 * pq.m / pq.s
# Reynolds number:
Reynolds_number = (density_water * mean_velocity_water * tube_inner_diameter) /dynamic_viscocity_water
print(np.mean(Reynolds_number.to_base_units()))
# Dittus-Boelter equation:
Nusselt_number = 0.0265 * (Reynolds_number) * (Prandtl_number) # For cooling
print(np.mean(Nusselt_number.to_base_units()))
Однако это дает мне размерность (метр ** 2.22045e-16 / второй ** 2.22045e-16) для Nusselt_number, где он тоже должен был быть безразмерным.
Код: Выделить всё
Nusselt_number = 0.0265 * (Reynolds_number**0.8) * (Prandtl_number**0.3) # For cooling
print(np.mean(Nusselt_number.to_base_units()))
Подробнее здесь: https://stackoverflow.com/questions/790 ... nits-while