【文献阅读】2004 Goktekin NS方程加弹性 A Method for Animating Viscoelastic Fluids

概述

题目：A Method for Animating Viscoelastic Fluids

一种模拟粘弹性流体的方法

发表时间：2004年

期刊/会议：SIGGRAPH

作者：Tolga G. Goktekin Adam W. Bargteil James F. O’Brien 加州大学伯克利分校

摘要：本文描述了一种为粘弹性流体（如粘液、液体肥皂、布丁、牙膏或粘土）的行为设置动画的技术，这些流体同时具有流体和固体的特性。该技术基于之前的欧拉方法，通过在基本Navier-Stokes方程中加入额外的弹性项，为具有自由表面的不可压缩流体设置动画。弹性项是通过在整个流体中积分和平流应变率来计算的。从弹性阻力到粘性流动的转变由冯·米塞斯的屈服条件控制，随后的行为由准线性塑性模型控制。

图片：

笔记：

作者认为，流体和固体的唯一不同点在于是否存在弹性力（使材料回复原始形状）。

流体和固体的相同点在于：（1）都抵抗体积变化（2）固体阻尼力和流体粘性力都可以表示为应变率的函数。

粘弹性体是介于流体和固体之间的材料。为该材料逐步施加外力，它先表现出弹性，然后开始塑性流动。分界点在于所承受的应力大小，这里称之为弹性极限。当弹性极限设定得很高时，材料表现为固体。当弹性极限设定得很低时，材料表现为流体。介于两者之间得，如粘液、肥皂液、牙膏、粘土，被称为粘弹性流体或弹塑性固体。(note: 从这里看出来，粘弹性是从流体角度说的，弹塑性是从固体角度说得，两者其实是一种东西)

采用von Mises屈服准则确定流体/固体的转换点。

Related works:

[Fung, 1965], [Han and Reddy, 1999], and [Bird et al.,1987]这三篇文献阐述了粘弹性流体材料得不同描述方法。

本文方法主要借鉴[Enright et al., 2004]。采用欧拉网格模拟自由表面。该方法由两个部分组成。一是直角网格，用于定义流体状态。二是level set函数。

Method:

控制方程（NS方程）如

关于此方程：

1. 关键在于右端第四项，即弹性力项。
2. 该项中，${\mu }_e$是弹性模量，给定一个常数。
3. 该项中，关键在于$ϵ$，即应变（弹性应变）。应变是没法直接计算的，这也是本文难点和创新点之一。

（Notes: 没标注的应变默认为弹性应变）

应变无法直接计算是因为物体时刻处于流动状态，没法像固体那样追踪形变。所以本文提出了一种特殊的处理方法，即用积分应变率的方式计算应变。

处理方法大致如下：先计算一个总应变，总应变由总应变率的积分得到（应变率即应变的导数）。该积分的起始点是0时刻的总应变率（设定为0），如图所示。

而总应变率是可以由速度的梯度（结果为二阶张量）计算的，如下图。

只要再计算出塑性应变，弹性应变就可以由总应变减去塑性应变得到。

塑性应变同样可由塑性应变率积分得到。而塑性应变率的计算则要利用von Mises准则。该准则说明，只要当弹性应变达到一定程度（定义为屈服点）的时候，就会发生塑性流动。塑性应变率是与超出屈服点的部分成正比的。具体从数学公式上来讲：

先定义一个临时变量，称之为应变偏置。它是弹性应变的函数。它的物理意义是弹性应变去掉膨胀的部分（因为膨胀收缩这种体积变化是不被我们允许的，我们要做的是不可压缩流体）。

然后采用von Mises屈服准则：

这个公式的物理意义是：假如应变没有超过屈服点，塑性应变率为0；假如应变超过屈服点，塑性应变率与超出屈服点的应变成正比。这点可以从max()这个函数看出。其中**参数$\gamma$就是屈服点，控制了发生流动的难易程度，越大越难流动。前面的系数$\alpha$**控制了塑性流动的快慢，越大越快。通过调节这两个参数，就实现了粘弹性流体的模拟。

我们把上述公式综合起来，得到一个公式：

弹性应变是我们要求解的未知量。观察发现，右边的应变偏置${ϵ}^{\prime }$是弹性应变的函数。左边则是弹性应变率，即弹性应变的导数。所以这是一个微分方程，或者说积分方程。

（假如只用一个公式来概括附加弹性力的NS方程中的弹性力怎么算，只要上面这一个公式就够了）

额外附注：

1. 原文是采用欧拉法的，所以还需要用半拉格朗日法进行平流，弹性应变也是需要被平流的。
2. 运算速度几乎与没添加弹性项的模拟一样。
3. 由于适用交错网格，应变张量的不同分量，和速度一样，是存储在不同面的。
4. 用level set进行表面重建。level set的网格精度为计算网格的二倍。
5. 作者提到该方法也可以模拟血液，但许多血液的特性是没法表现的，这些特性与血液的微观结构有关。其他聚合物悬浮液同理。