[发明专利]水场景目标红外尾迹特性的流体模拟方法

权利要求书

1.一种水场景目标红外尾迹特性的流体模拟方法，其特征在于包括以下步骤：

1)对场景中的模型进行体素化产生固体粒子；

2)在存在液体的区域添加液体粒子；

3)对每一时间帧，利用SPH对流体的力学特性进行数值模拟；

所述步骤3)为：

(3.1)对于每一个液体粒子L_i，根据SPH方程定义，利用所有光滑核半径内的液体粒子L_j或固体粒子S_j计算液体的密度ρ_i：

ρ_i＝ρ(r_i)＝315m/(64πh⁹)∑_j(h²-|r_i-r_j|²)²

其中r_i是L_i在空间中的位置，m是L_i的质量，h为光滑核半径，r_j为粒子L_j或S_j的空间位置；

(3.2)对于每一个粒子L_i，根据理想气体状态方程计算液体的压强：

p_iM＝ρ_iRT_i

其中ρ_i，T_i，p_i分别为L_i的密度、温度、压强，M为平均摩尔质量，R为理想气体常数；

(3.3)根据SPH公式，利用所有光滑核半径内的粒子L_j或S_j计算每一个粒子L_i的加速度a_i：

a_i＝a(r_i)＝g+m*45/(πh⁶)∑_j((p_i+p_j)/(2ρ_iρ_j)*(h-r)²*(r_i-r_j)/r)+

mμ*45/(πh⁶)∑_j(u_j-u_i)/(ρ_iρ_j)*(h-r)

其中，r＝|r_i-r_j|，ρ_j是粒子L_j或S_j的密度，u_j为粒子L_j或S_j的速度，a_i，u_i分别为粒子L_i的加速度和速度；

(3.4)利用加速度与速度，根据牛顿第二定律计算模拟每个粒子L_i或S_i的运动轨迹，重复步骤(3.1)～(3.4)即可计算得到每一时刻每个粒子的力学特性；

4)对每一时间帧，利用瞬态热传导方程对流体的热力学特性进行数值模拟及红外特征图形绘制；

所述步骤4)为：

(4.1)对于每个粒子L_i，根据笛卡尔坐标系下的热传导方程，利用狄利克雷边界条件计算以下方程得到温度变化：

其中W为光滑粒子核函数W(x)，W_ij＝W(|r_i-r_j|)，q，T为积分点的热流密度与温度；

(4.2)根据计算得到的温度场T的分布，利用Planck公式计算以下方程得到其红外辐射量E：

其中ε₀为材料的表面发射率；C₁为第一辐射常数，其值为3.742×10^-16W·m²；C₂为第二辐射常数，其值为1.4388×10^-2m·K，λ₁和λ₂为红外探测仪的探测波段的两段波长；

(4.3)若场景中最高和最低辐射温度分别为T_max和T_min，由上式算得的相应的辐射强度分别为E_max和E_min，则对于辐射强度为E_i的表面，其对应的计算灰度值为：

利用计算机图形学中的颜色均衡方法，绘制得到红外特征图形效果。