有了上一篇的详细分析后,这里其实没有太多的新东西,主要是提出了温度的影响,加了两个力,提出mac网格,提出共轭梯度法求解。大概列一下
一、力
1、浮力与重力
烟雾会由于重力下沉,由于浮力上升,对于这两个力本文建立了一个模型,令:
α、β是正常量,z是向上的单位向量( 0 , 0 , 1 ),Tamb是空气温度。顺带一提,这里的温度场也遵循advection流动。
2、旋转力
另一个力是控制旋涡,在本文之前有一些方法是随机地添加一个无散度的速度场来制造一些细节,刚好看过一个图拿过来:
本文研究的是物理正确的细节。首先旋度 w = ∇ × u,单位化向量 η= ∇ |w|,N = η / |η|,最后