摘要：本文通过对真实火焰物理特性的分析，建立了火焰动态燃烧的数学模型，并在此数学模型基础之上借助于DirectDraw技术对图形显示的加速，在VC++ 6.0下对火焰作了效果非常逼真的计算机动态仿真。
关键词：火焰；DirectDraw；计算机仿真
引言
计算机仿真技术的基本原理都是一样的，神秘复杂的核爆同水波、火焰、烟雾等非常平常的自然现象在仿真处理过程中并没有什么太大的区别。都是经历了从实体对象到物理特性的总结，再由此建立数学模型并在数学模型基础之上提出仿真算法，最后通过计算机将其动态仿真出来等一系列步骤。本文以火焰作为仿真对象，通过对热源、热扩散以及对流等特性的分析对其建立了数学模型及仿真算法，为了能充分发挥计算机对图形的硬件加速，使用DirectDraw技术对仿真结果显示进行了加速，使之能逼真、流畅地对火焰的燃烧过程实行动态模拟。
简单近似模型设计
虽然火焰在自然界是一种极普通的自然现象，但根据流体力学的相关知识，火焰可以表达为一个相当复杂的三维动态流体系统。如要在计算机中对这样一个复杂的流体系统做出精确的仿真将需要有相当庞大复杂的数学模型为基础，而且运算量将非常巨大，在现有的微型计算机中几乎很难保证其动态实时性，这也就失去了仿真的意义。因此，在仿真时应用尽可能简单的模型来实现尽可能逼真的效果。
从物理角度分析，要产生火焰，首先要有火源，其次为了产生"焰"的效果，需要以火源为中心向上、向四周扩散，而且由于在扩散过程中逐渐远离火源，温度会逐渐下降，表现在视觉上就是火焰的冷却变暗。另外，由于火焰的高温使周围空气受热膨胀比重下降，因此会有空气的对流出现，这将把火焰向上"吹"起，使火焰向四周扩散的距离要远小于向上扩散的距离。基于以上几点认识，可以采取对应的仿真措施：对火源的设置可以用一幅二值位图来标识，非火源以低亮度像素填充、火源点则设以高亮度像素，通过对位图像素值的判别可以断定当前点是否为火源。
While(TRUE){
for(y=1;y
n2 = buffer1(x-1, y)
n3 = buffer1(x, y+1)
n4 = buffer1(x, y-1)
p = ((n1+n2+n3+n4+p) /5)； //四临像素均值
p = p-c； //同一固定冷却衰减值相减
if(p<0)
p=0
buffer2(x,y-1)=p
}
}
copy buffer2 to the screen ; //显示下一帧
copy buffer2 to buffer1； //更新Buffer1
}
根据上述近似模型可对火焰进行一定程度上的仿真，但由于没有引入随机分布火焰往往看上去相当单调规则，而且火焰总呈线性上升，冷却速度也严格地保持统一速率。要消除以上问题，可通过引入随机非均匀因素来解决。一种途径是随机布置各点冷却值使火焰冷却过程非均匀化。但由于火焰的模拟过程是实时进行的，为确保动态模拟过程中能顺畅进行，最好用预先创建的冷却位图（见右图）来代替。一般采用在屏幕上随机撒布几千个亮度不同的点并对其应用平滑处理等方法对冷却位图加以填充。通过冷却图中获取的数值来代替原来固定的冷却衰减值效果要好的多，此时的冷却过程改进为Pixel(x,y)=Pixel(x,y)-Coolingmap(x,y) 这样的衰减结果将使火焰的冷却衰减效果更加真实：
pp = coolMap + coolMapWidth*2;
p1 = lightBuf1+imgWidth*2;
p2 = p1 - imgWidth;
p3 = p1 - 1;
p4 = p1 + 1;
p5 = p1 + imgWidth;
for(i=0;i
c1=(unsigned char)(((UINT)*p1+(UINT)*p2+(UINT)*p3+(UINT)*p4+(UINT)*p5)/5);
c2 = *pp;
if(c1>c2)
c1 -= c2;
*p = c1;
pp++,p++,p1++,p2++,p3++,p4++,p5++; //内存指针修正
}
topX = (imgWidth - seedMapWidth)/2; //seedMapWidth种子位图宽度
topY = (imgHeight - seedMapHeight)/2; //seedMapHeight种子位图高度
p = lightBuf1 + (topY+2)*imgWidth + topX; //p, unsigned char型指针
ps = seedMap + seedMapWidth*2; //ps, unsigned char型指针
for(j=0;j<(seedMapHeight-4);j++) {
p1 = p; //p1, unsigned char型指针
for(k=0;k
if(rand() < t)
*p1 = 255;
}
p1++,ps++; //指针修正
}
p += imgWidth; //指针修正
}
[1] [2]