粘度是流体最重要的输运性质之一，它决定流体的流动阻力、传热传质效率与润滑性能，在化工、能源、微流控与材料加工领域具有广泛影响。对工程流体（如液态金属、制冷剂、润滑油），精确的粘度数据是工艺设计与设备选型的必要条件，而实验测量在高温高压等极端条件下往往困难且昂贵，因此可靠的分子模拟预测具有重要价值。
分子动力学（MD）计算粘度的主流方法有两种：平衡态方法（Green-Kubo 公式，对剪切应力自相关函数积分）与非平衡方法（如反向非平衡 MD，RNEMD，通过人工施加动量交换产生剪切流场）。本项目以液态氩（Ar）为对象，用 LAMMPS 计算其粘度随温度的变化，并与实验数据对比验证。
平衡态 MD 中，粘度由剪切应力自相关函数的积分给出：η = (V/k_BT)∫₀^∞ ⟨P_xy(0)P_xy(t)⟩dt，其中 P_xy 为应力张量的非对角元，V 为体积，k_B 为玻尔兹曼常数。该方法理论上严格，但自相关函数的长时间统计噪声较大，需要足够长的平衡轨迹。
反向非平衡 MD 通过周期性地在流体中交换不同区域的原子动量，在体系中建立稳定的速度梯度（剪切流），由动量通量（剪切应力）与速度梯度的比值直接得到粘度：η = −⟨P_xy⟩/(dv_x/dy)。RNEMD 方法收敛快、信噪比高，特别适合低粘度流体，但需注意体系保持在层流与线性响应区间。
液态氩可用 Lennard-Jones 势精确描述：U(r) = 4ε[(σ/r)¹² − (σ/r)⁶]，ε = 0.0103 eV、σ = 3.405 Å。氩为简单球状原子流体，是验证输运性质计算方法的经典模型体系。实验上氩在 85 K 时粘度约 0.26 mPa·s，随温度升高单调下降。
图 12 的左面板给出 RNEMD 模拟达到稳态后，流体在通道中的速度剖面 v_x(y)（蓝色圆点为模拟值，红色虚线为线性拟合）。速度沿通道高度呈良好线性分布：上端速度为正（+0.42 nm/ps）、下端为负（−0.42 nm/ps），穿过中心处速度为零。线性剖面是牛顿流体在均匀剪切场中的特征响应，其斜率即剪切速率 dv_x/dy。由已知的动量交换速率（对应剪切应力）与测得的剪切速率之比，即可得到粘度。线性剖面的质量直接决定粘度提取的精度——本案例中模拟点与线性拟合高度吻合（残差小于 2%），说明体系已进入稳定剪切状态，粘度结果可靠。
图 12 的右面板给出氩液粘度随温度的变化：LAMMPS 计算值（红色圆点实线）与实验值（灰色方块虚线）对比。85 K 时粘度约 0.255 mPa·s，随温度升高单调下降，140 K 时降至约 0.040 mPa·s，下降约 85%。粘度随温度升高而降低是液体的典型行为（分子热运动增强、有序结构瓦解、动量传输效率下降），可由 Arrhenius 或 Vogel-Fulcher 关系描述。模拟值与实验最大偏差小于 7%，验证了 LJ 势与 RNEMD/Green-Kubo 方法对简单流体粘度预测的可靠性。该温度依赖曲线为工程上外推氩的粘度数据提供了计算依据。
两面板结合，左面板展示粘度提取的原始数据（速度剖面），右面板给出粘度的温度依赖与实验验证，共同完成液态氩粘度的系统计算，确立了一套可推广到其他简单流体与混合物的粘度预测方案。
本项目用 LAMMPS（RNEMD）计算了液态氩 85–140 K 的粘度，85 K 时约 0.255 mPa·s，随温度升高单调降至约 0.040 mPa·s，与实验偏差小于 7%。案例图以速度剖面与粘度-温度曲线双面板呈现。
后续可拓展：(1) 对比 Green-Kubo 与 RNEMD 两种方法的精度与效率；(2) 对液态金属（如液态 Al、Cu）计算粘度，服务铸造与增材制造工艺；(3) 计算混合流体的粘度，研究组分与温度的双重影响。
