非牛顿流体:倒番茄酱的困惑

2013年 4月 11日

如果你喜欢吃番茄酱,那么很有可能你已经体验过我们称之为“要么全有要么全无的番茄酱困境”。想象一下,你伸手去拿一个新的番茄酱玻璃瓶,取下瓶盖,把瓶子倒过来时,然后——“什么都没有发生”。直觉上,我们会摇晃或轻拍瓶子,然后——突然你的食物就完全被涂上了番茄酱(当然,除非你的反应真的很好)。在这篇博客文章中,我们将通过仿真模拟瓶子中的番茄酱的流量来演示倾倒番茄酱时的“全有或全无”效果。

我的番茄酱在哪里?

番茄酱除了是一种特别受欢迎的调味品以外,也是一种非牛顿流体。这是它会出现上述奇怪现象的关键。作为一种非牛顿流体,番茄酱的黏度——以及其流动能力——会随着流体的速度梯度或剪切速率变化。当第一次把番茄酱瓶子倒过来的时候,由于番茄酱的初始黏度很高,很可能你只能得到非常小的流量,如果有的话。好在番茄酱还有一个有趣的剪切稀化特性。当在重力之外施加外力时,增加的剪切力将导致番茄酱的黏度降低,更容易流动。这就是为什么我们有时必须轻拍或摇晃瓶子才能享受番茄酱加薯条。

模拟瓶子中的非牛顿流体流动

为了满足我们的好奇心,以及我们的胃口,并进一步完善我们的番茄酱倾倒技术,我们使用COMSOL Multiphysics 软件和其附加的计算流体力学模块对这个过程进行模拟。在模拟中,我们做了几个简化的假设,包括瓶子的无限延伸(即番茄酱一直都有)和单相层流,忽略了瓶内滞留的空气。所有速度均以米/秒(m/s)为单位。

我们规定了瓶子的几何形状和层流的边界,包括无滑移壁面、压力入口(代表未建模番茄酱部分的重量)和压力出口(外部大气),如下图所示。

2D geometry of a ketchup bottle
用于模拟的二维几何图形。显示了应用的边界条件。

我们还添加了体积力特性,以包括重力和沿着 y 轴摇晃瓶子 0.2 秒所产生的力。

Volume force
为了给求解器提供一致的初始条件,重力产生的体积力通过阶跃函数从 0 上升到其最大值(过渡区域为 0.05 秒)。摇动动作是以 0.2 秒宽的矩形函数模拟的。在此期间,体积力加倍。

倾倒番茄酱

我们基于Koocheki 等人的工作,我们先做了一个仿真为之后模拟番茄酱的流动做准备。为了模拟番茄酱的流动,我们使用了非牛顿幂律来表示动力黏度。番茄酱配方通过在幂律方程中指定稠度系数和流动特性指数来模拟。下面的图像显示了添加1%黄原胶的番茄酱配方中较高的初始黏度和非牛顿流动性。

A plot of the dynamic viscosity.
Viscosity measured in COMSOL Multiphysics.

:p1(见上面的模型设置)为单位测量的黏度是时间的函数,时间间隔 [00.1]秒。右图:p1(见上面的模型设置)为单位测量的黏度是时间的函数,时间间隔为[0.12]秒。番茄酱的剪切稀化特性很明显,尤其是在瓶子摇晃的[11.2]秒间隔内。黏度降低,从而允许更多的番茄酱流出瓶子。

Plot showing viscosity measured along line 1.
Graph depicting the average outlet velocity.

左:在三个不同时间沿线 1(见上面的模型设置)测量的黏度。右图:平均出口速度与时间的函数关系。可以看出,番茄酱在重力的作用下慢慢开始流动,然后在摇晃瓶子时容易增加。

An animation showing pouring ketchup.

Simulation plot depicting the velocity magnitude of pouring ketchup at 1.1 seconds.

左:速度幅度动画从 0.05 到 2 秒。右:1.1 秒时的速度大小。

“番茄酱困境”

等待流速的变化对你来说可能很熟悉。一旦流速开始改变,就会快速变化,导致番茄酱泛滥,覆盖在你的食物上。流速顺其自然地增加导致了剪切力增加,进而降低黏度,导致更高的流速。这就是为什么番茄酱最初似乎是静止的,但后来突然又像水一样流动的原因。

番茄酱从瓶子中流出这个特殊的问题比一开始看起来要复杂得多。有许多变量会影响我们体验一个“绝对的” 番茄酱惊喜。如果你有兴趣进一步研究这种效应,建议你使用 COMSOL Multiphysics 和计算流体力学模块进行模拟实验。例如,你可以研究诸如倾斜角度、瓶子里的空气、甚至瓶子的设计等变量是如何影响番茄酱的倾倒体验的。


评论 (0)

正在加载...
浏览 COMSOL 博客