TechOverflow Logo
TechOverflow

如何使用 UliEngineering 在 Python 中计算 Andrade 粘度

English Deutsch Français 简体中文

你可以使用 UliEngineering Python 库轻松计算液体的 Andrade 粘度。Andrade 方程是一个 Arrhenius 型模型，用于描述粘度对温度的依赖关系：

andrade_viscosity.py
from UliEngineering.Physics.Viscosity import andrade_viscosity, CommonLiquids
from UliEngineering.EngineerIO import *

# 计算 20°C 时水的 Andrade 粘度
T = 20 + 273.15  # 转换为开尔文
eta = andrade_viscosity(T, CommonLiquids.Water.andrade)
print(f"Water viscosity at 20°C: {format_value(eta, 'Pa·s')}")

# 计算 50°C 时乙醇的 Andrade 粘度
T = 50 + 273.15
eta = andrade_viscosity(T, CommonLiquids.Ethanol.andrade)
print(f"Ethanol viscosity at 50°C: {format_value(eta, 'Pa·s')}")

# 使用自定义 Andrade 常数
from UliEngineering.Physics.Viscosity import AndradeConstants
custom = AndradeConstants(name="Custom liquid", A=1.0e-6, B=2000.0)
eta = andrade_viscosity(300.0, custom)
print(f"Custom liquid viscosity at 300K: {format_value(eta, 'Pa·s')}")

示例输出

andrade_viscosity_output.txt
Water viscosity at 20°C: 1.00 mPa·s
Ethanol viscosity at 50°C: 1.05 µPa·s
Custom liquid viscosity at 300K: 2.35 mPa·s

Andrade 粘度图

Andrade 方程由下式给出：

$$ \eta = A \cdot \exp\left(\frac{B}{T}\right) $$

其中 $\eta$ 为动力粘度，$T$ 为开尔文绝对温度，$A$ 和 $B$ 为材料特定的常数。参数 $B$ 与粘性流动的活化能有关。

上图展示了在 0°C 到 100°C 温度范围内若干常见液体的 Andrade 粘度模型。注意粘度随温度呈指数下降，这是大多数液体的典型特征。

UliEngineering 库通过 CommonLiquids 类为水、乙醇、甲醇、甘油、橄榄油、汞、丙酮和苯等常见液体提供了预定义的 Andrade 常数。

相关文章

绘图生成脚本

plot_andrade_viscosity.py
#!/usr/bin/env python3
import matplotlib.pyplot as plt
import numpy as np
import sys
sys.path.insert(0, '/home/uli/dev/UliEngineering')

from UliEngineering.Physics.Viscosity import (
    andrade_viscosity,
    CommonLiquids,
)

# 用于绘图的摄氏温度范围
T_C = np.linspace(0, 100, 200)  # 0 到 100°C
T_K = T_C + 273.15  # 转换为开尔文

# 创建图形
plt.figure(figsize=(10, 6))

# 绘制不同常见液体的曲线
liquids_to_plot = [
    ('Water', CommonLiquids.Water.andrade, 'blue'),
    ('Ethanol', CommonLiquids.Ethanol.andrade, 'green'),
    ('Methanol', CommonLiquids.Methanol.andrade, 'red'),
    ('Glycerol', CommonLiquids.Glycerol.andrade, 'purple'),
    ('Mercury', CommonLiquids.Mercury.andrade, 'orange'),
]

for name, constants, color in liquids_to_plot:
    eta = andrade_viscosity(T_K, constants) * 1000  # 转换为 mPa·s
    plt.plot(T_C, eta, label=name, color=color, linewidth=2)

plt.xlabel('Temperature (°C)', fontsize=12)
plt.ylabel('Dynamic Viscosity (mPa·s)', fontsize=12)
plt.title('Andrade Viscosity Model for Common Liquids', fontsize=14, fontweight='bold')
plt.legend(loc='upper right', fontsize=10)
plt.grid(True, alpha=0.3)
plt.yscale('log')  # 对数刻度，因为粘度跨越多个数量级

plt.tight_layout()
plt.savefig('andrade_viscosity_plot.svg', format='svg', dpi=300)
print("Plot saved to andrade_viscosity_plot.svg")
Check out similar posts by category: Physics, Python