烟的视觉模拟外文文献翻译中文译文.docVIP

烟的视觉模拟 Visual Simulation of Smoke Ronald Fedkiw Stanford University Jos Stam Aliaswavefront Henrik Wann Jensen Stanford University 翻译：YeHaike (/yehaike) 摘要 废话略(请看原文)，正式翻译请看第2节。 1 介绍 废话略(请看原文)，正式翻译请看第2节。 1.1 前人工作 废话略(请看原文)，正式翻译请看第2节。 1.2 我们的模型 废话略(请看原文)，正式翻译请看第2节。 2 流体流动方程组 从一开始，我们假设我们的气体能作为非粘性、不可压、常量密度(constant density)流体进行模拟。在气体中粘性的影响是微不足道的，尤其在粗网格上，粗网格上的数值耗散比起物理粘性和分子扩散来说更占主要地位。当烟的速度远低于声速时，可压缩性的影响也微不足道，这样不可压缩的假设就很大程度上地简化了数值计算方法。因此，烟的速度模型方程，记为，由不可压缩欧拉方程(Euler equations)组[14]给出 这两个方程说明速度应该让质量守恒(conserve mass，方程1)和动量守恒(conserve momentum，方程2)。物理量p是气体的压强(pressure)，f表示外部力(external forces)。我们也任意设置流体的常量密度(constant desity)为1. 就像[7，6，17]中的那样，我们使用两步来解这些方程。首先我们通过求解方程2(经过一个时间步?t，不包含压力项)计算一个中间速度域u*: 在这一步之后，我们使用一个投影方法(projection method)[3]来强制这个域u*为不可压。这等价于从下面的泊松方程(Poisson equation)计算压强(pressure)： 在一个法向量为n的边界点处，使用纯Neumann(诺伊曼)边界条件，也就是，。(注意利用Dirichlet边界条件也很简单，其中压强直接被指定，而不是指定它的法向导数(normal derivative))。中间速度域u*之后通过把它减去压强的梯度(gradient of pressure)被修改成不可压缩的： 我们也需要用于温度T和烟的密度ρ的进化(evolution)的方程组。我们假设这两个标量(scalar quantities)简单地被沿着烟的速度移动(平流运输)。 密度和温度都影响流体的速度。由于重力，重的烟会向下掉，然而由于浮力，热气体会向上升。我们使用一个简单的模型来说明这种效果，该模型定义了一个外部力，该力直接与密度和温度成比例： 其中z=(0,0,1)点在向上垂直方向。Tamb是空气的环境温度，ɑ和β是两个具有适当单位的正符号常量，这样方程8才有物理意义。注意当ρ=0和T=Tamb时，这个力等于0。 方程2，6，7都包含平流运算符(advecton operator) 。如[17]中那样，我们使用一个半拉格朗日法(a semi-Lagrangian method)[19]来求解该项。我们使用一个迭代求值器(iterative solver)来求解用于计算压强的泊松方程(Possion equation，方程4)。我们在第4节介绍这些求值器怎样也能处理物体浸入流体(bodies immersed in the fluid)的问题。 3 漩涡限制 通常烟和空气混合物包含的速度域具有很大空间偏差，并伴随着相当大量的各种规模的旋转和狂暴的结构。非物理的数值耗散减弱了这些有趣的流动特征，我们的新方法的目的就是在粗网格上把这些漩涡加回来。把它们加回来的方式之一就是创建一个随机或假随机的小规模流动域扰乱，使用试探式的(启发式的)或基于物理的模拟。例如，可以使用Kolmogorov谱产生一个divergence free速度域并把这个加到计算出的流动域来代替丢失的小规模结构(了解Kolmogorov谱的CG应用，看[18])。虽然这给流动提供了小规模细节，但是这不会把小规模细节放在物理上正确的位置，这个正确位置就是在流动域中小规模细节丢失的地方。反而，这些细节以一种杂乱的形式被添加，这让烟变得“活跃”，以一种非物理的形式翻滚和卷曲。烟的逼真运动的关键在于使烟看起来像是被动的自然现象而不是相对的由活生生的人制造出来的烟。 我们的方法寻找流动域中小规模特征应该产生的位置，然后在粗的CG(计算机图形学)网格上，以一种基于物理的方式把小规模特征加进这些位置，促进了烟的被动滚动，这样让烟看起来有真实地混乱。具有无限的计算能力，任何符合的数值方法都可以用于获得可接受的结果，只需简单地增加网格点的数量直到想要的限制行为被观察到。但是，实际上，计算资源有限，网