【NFX CFD与流体力学】：牛顿流体和非牛顿流体
midas NFX是集结构、热、流体和优化设计于一体的有限元仿真分析软件，而midas NFX CFD则是其中计算流体动力学求解器高级模块。它采用独特的有限单元法（FEM）进行编程，在midas NFX 程序中搭配结构求解器组成一体化的图形用户界面（GUI），并提供全中文界面。是一款简单易学，易上手，并且实用的 CFD 软件。
流体的粘度是其在给定速率下对变形的抵抗力的量度。两种不同粘度的液体自由落体时，相同的时间内粘度高的液体变形小，有一种“拉丝”的感觉，而粘度低的流体介质的连续性较差，很容易产生分离。
模拟具有不同黏度的液体，右面的液体比左面的液体粘度高
产生以上现象的原因是处于相对运动状态的相邻流体层之间的内部摩擦力，也称黏滞力，是流体受到剪应力变形或拉伸应力时所产生的阻力。
为了研究该阻力的性质，1687年，牛顿首先做了最简单的剪切流动实验。在平行平板之间充满粘性流体，平板间距为y，下板静止不动，上板以速度U在自己平面内等速平移。由于板上流体随平板一起运动，因此附在上板的流体速度为U，附在下板的流体速度为零。
牛顿平行板剪切流动试验
实验指出，两板之间的速度分布u(y)服从线性规律。作用在上板的力与板的面积A、板的运动速度U成正比，同间距y成反比。由此得出：
其中：τ为剪应力（Pa）， du/dy为速度梯度（s-1），μ为动力粘度（ Pa·s ）。
因此，将满足牛顿内摩擦定律的流体称为牛顿流体，即切应力和速度梯度之间是线性关系。例如：水、汽油、空气等低分子量的流体。
将不满足牛顿内摩擦定律的流体称为非牛顿流体，切应力和速度梯度之间不是线性关系的流体。例如：石油、牙膏、泥石流、等高分子聚合物的浓溶液和悬浮液。
NFX软件根据流体粘度将流体分为牛顿流体和非牛顿流体
在midas NFX软件中，牛顿流体的粘度是常数，不随速度而变化。这里说的常数不是说任何情况下都是常数，而是说在恒温下是一个常数。所以说牛顿流体的粘度可以随温度发生变化。
牛顿流体的粘度可以定义为随温度变化的函数
恒温下牛顿流体的粘度是一个常数。用通俗的话来说，就是在恒定的温度下你无论用多么大的力量去攻击水，流体的粘度都不会发生变化。
当用锤子去重击牛顿流体之蜂蜜时
非牛顿流体
然而，当我们快速的用锤子去锤击非牛顿流体的表面，神奇的现象发生了，根本锤不下！流体的粘度通过外力被彻底的改变了。
快速锤击淀粉糊的表面根本锤不下去
但是，当我们将锤子缓慢的放入非牛顿流体淀粉糊中，非牛顿流体很容易包裹住锤子。说明此时流体粘度较小。
锤子缓慢的放入淀粉糊很容易被包裹
展现出的这种“遇强则强”特性的非牛顿流体我们称之为“胀塑性流体”，是非牛顿流体中的一种。
胀塑性流体
表观粘度随切变速度的增加而增加的流体。在网上视频里看到的非牛顿流体大多是这一种。在很小的剪切应力下就可能流动，但在非常高的剪切应力或剪切速率下，粘度可能无限地增大，即存在剪切增稠行为。许多高分子的分散体系，如固含量很高的悬浮液、糊状物、涂料以及泥浆、淀粉、高分子凝胶，都属于胀塑性流体。玉米淀粉水溶液就是一种最常见也是最容易得到的胀塑性流体。
假（拟）塑性流动
与胀塑性流体相反，表观粘度随切变速度的增加而减小的流体。存在剪切变稀行为。大部分液态食品都是假塑性液体。如菜汤、酱油、浓糖水、番茄酱、苹果酱等高分子溶液、乳状液和悬浮液都属于假塑性流体。
番茄酱在受到晃动以后会变稀，更容易倒出来
当所受剪切应力超过临界剪切应力 τy 后，才能变形流动的流体，也称塑性流体。但一旦发生流动，其粘度保持不变，呈牛顿流体行为。如果超过临界剪切应力τy 后，呈现剪切变稀或剪切增稠的非牛顿行为，则称此流体为广义宾汉流体。
为什么称为塑性流体，是因为在屈服前，宾汉流体实际是弹性固体。牙膏就属于宾汉流体。所以要想让牙膏流动，必须要先用力挤出来。
彩色牙膏里面的彩色不会混合
选择非牛顿流体时，midas NFX CFD 支持幂律、CarreauYasuda 模型、交叉模型、Herschel-Bulkley 模型以及粘度曲线等 5 种模型。
每种模型对应的粘度表示式如下：
① 幂律： 
② Carreau Yasuda模型：
③ 交叉模型：
④ Herschel-Bulkley模型：
⑤ 粘度曲线：输入粘度与剪切速率的离散数据点。
其中， η为粘度（Pa·s）；
表示剪切速率；λ,a为材料特征常数；η_0为第一牛顿粘度，η_∞为第二牛顿粘度。
常见物质的粘度
空气
水
水的粘度可按下式计算：μ=0.001779／（1＋0.03368t＋0.0002210t2）
或者 μ=0.02939*EXP(507.88/(K-149.3))
| 液体 | 黏度 [Pa·s] | 黏度 [cP] |
|---|---|---|
| 血液（37 °C） | (3–4)×10-3 | 3–4 |
| 蜂蜜 | 2–10 | 2,000–10,000 |
| 糖浆 | 5–10 | 5,000–10,000 |
| 熔融态玻璃 | 10–1,000 | 10,000–1,000,000 |
| 巧克力酱 | 10–25 | 10,000–25,000 |
| 熔融巧克力 | 45–130 | 45,000–130,000 |
| 番茄酱 | 50–100 | 50,000–100,000 |
| 猪油 | ≈ 100 | ≈ 100,000 |
| 花生酱 | ≈ 250 | ≈ 250,000 |
| 酥油 | ≈ 250 | ≈ 250,000 |
我们知道，对于具有恒定流速的管道，要想维持管道内流体的流动，需要管道两端存在压力差来克服处于相对运动状态的流体层之间的摩擦。这个力与流体的粘度和剪切速率成正比。
流体的粘度度是由于相邻层间以不同的速度运动时产生的摩擦造成的。管中心处阻力最小，液层流动速度最大；管壁附近液层同时受到液体黏性阻力和管壁摩擦力作用，速度最小，在管壁上液层的移动速度为零（假定壁面为无滑移壁面时）。
考虑液体粘性，管轴附近的流动比在管壁附近流动得更快
流体分析时常用的流体是空气和水，都属于牛顿流体。感兴趣的朋友可以尝试做一个非牛顿流体的仿真，如下。限于篇幅，后期我们将会逐步分享更多案例，欢迎后台私信交流。
非牛顿流体仿真
参考文献：J. Sepulveda, A. Montillet, D. Della Valle, C.Loisel, A. Riaublanc, Deformation of gas-liquidinterfaces in a non-Newtonian fluid at high throughputs inside a microfluidicdevice and effect of an expansion on bubble breakup mechanisms, ChemicalEngineering Science (2019)
挤牙膏案例来源B站UP主：USim_Abaqus
