第一节：固定式发烟机成烟（投影）面积的计算

将烟雾的形成分为两个主要阶段：①发动机主导的近场喷射阶段：决定烟雾的初始体积和浓度。② 大气环境主导的远场扩散阶段：决定烟雾的最终形状和面积。

1. 第一阶段：近场喷射模型 (发动机主导)

这个阶段的目标是计算从喷口出来的、已充分混合的烟雾气溶胶的初始体积通量和线源强度。

1.1 有效烟剂供应速率 (ṁ_smoke, effective)

这不仅仅是喂料器的设定值，而是考虑了混合效率后的值。
ṁ_smoke, effective = ṁ_feeder * η_mix

ṁ_feeder：喂料器设定的烟剂质量流量 (kg/s)。

η_mix： 文丘里管混合效率 (0-1之间的系数)。这是一个经验值，取决于文丘里管设计、粉体流动性等。优化设计可使 η_mix > 0.95。

1.2 初始烟云体积通量 (Q_initial)

假设发动机排气与吸入的空气/烟剂充分混合后，形成一个温度为T_jet，速度为V_jet的混合气柱。
Q_initial = (ṁ_air, core + ṁ_air, entrained + ṁ_smoke, effective) / ρ_jet

ṁ_air, core： 发动机核心空气质量流量 (kg/s)，与推力强相关。

ṁ_air, entrained： 文丘里管阶段卷吸的空气流量 (kg/s)。

ρ_jet： 混合射流的密度 (kg/m³)，可通过理想气体定律估算。

1.3 关键近场参数：初始扩散尺寸 (σ_y0, σ_z0)

高速湍流射流会有一个初始的扩散角度（θ）。在距离喷口下游一定距离（例如L = 10D， D为喷口直径）处，我们可以定义烟雾团的初始横向和垂直尺度（标准差）：
σ_y0 ≈ σ_z0 ≈ L * tan(θ/2)

θ： 湍流射流扩散角，通常在10-20度之间。

这个初始尺寸 σ_y0 和 σ_z0 是第二阶段模型的重要输入。

2. 第二阶段：远场扩散模型 (大气主导)

在这个阶段，烟雾团被看作一个在风速U下移动的、具有初始尺寸的连续体/线源。我们采用广泛应用于污染物扩散的高斯烟羽模型作为基础。

2.1 基本假设：

顺风方向（x轴）的传输由平均风速U主导。

横向（y轴）和垂直（z轴）的扩散服从高斯分布（钟形曲线）。

烟幕在地面（z=0）发生反射（即不考虑地面吸收）。

2.2 高斯烟羽模型公式：

在点(x, y, z)处的烟剂浓度C为：
C(x, y, z) = [S / (2π U σ_y σ_z)] * exp[-y² / (2σ_y²)] * { exp[-(z - H)² / (2σ_z²)] + exp[-(z + H)² / (2σ_z²)] }

2.3 模型参数与我们的系统关联：

源强 (S)： 这就是我们第一阶段计算的 ṁ_smoke, effective (kg/s)。它是烟雾的“源头”。

风速 (U)： 环境风速，例如3 m/s。

有效源高 (H)： 对于地面固定源，H ≈ 0。但如果喷流有向上的动量或浮力，H可能略大于0。

扩散参数 (σ_y(x), σ_z(x))： 这是模型的核心！ 它们表示在下风距离x处，烟羽在y和z方向的标准差。它们随时间（或随x）增大。

2.4 扩散参数 (σ) 的确定：

这是最复杂的部分，取决于大气稳定度（我们之前讨论的晴天/阴天）。常用的经验公式是Briggs公式，它给出了σ_y和σ_z随x变化的函数。

例如，在中性大气条件（最常见的假设）下，对于农村地形：

σ_y(x) = σ_y0 + 0.08 * x / √(1 + 0.0001 * x)

σ_z(x) = σ_z0 + 0.06 * x / √(1 + 0.0015 * x)

注意：我在这里加入了初始尺寸 σ_y0 和 σ_z0，以衔接我们的近场模型。标准模型通常从点源（σ0=0）开始。

3. 综合估算模型与计算流程

现在，将所有部分组合起来，用于估算可见烟幕的轮廓。

步骤1：定义可见阈值浓度 (C_threshold)
人眼或传感器能识别烟雾的最低浓度。这是一个经验值，需要通过实验确定。

步骤2：计算烟幕轮廓
对于每一个下风距离 x，我们寻找满足 C(x, y, z) = C_threshold 的 y 和 z 值。

烟幕宽度 W(x)： 在中心线（y=0, z=0）处，浓度最高。宽度可以近似为在特定高度（如z=1.5m，人眼高度）处，浓度降至阈值时的横向距离的2倍。
W(x) ≈ 2 * √[ -2 * σ_y(x)² * ln( C_threshold / C_max(x) ) ]
其中 C_max(x) = S / (π U σ_y(x) σ_z(x)) (地面中心线最大浓度简化式)。

烟幕长度 L_total：
烟幕的有效长度是当中心线最大浓度 C_max(x) 下降到 C_threshold 时的下风距离 x_max。
L_total ≈ x_max, 其中 C_max(x_max) = C_threshold

烟幕形状：
整个烟幕可以描绘为一个以x轴为中心，宽度为W(x)，高度通常为2.15σ_z(x)的近似棱形或椭球体。

步骤3：估算成烟面积 (A)
对于一个地面源，可见烟幕的投影面积可以近似为一个梯形或通过积分计算：
A ≈ ∫ W(x) dx （从 x=0 到 x=L_total）

一个更简单的近似是： A ≈ 0.5 * (W_initial + W_max) * L_total

W_initial 是近场结束时的宽度。

W_max 是烟幕最大宽度（通常在x_max附近）。

4. 总结与模型局限性

4.1 总结：

这个综合模型说明：成烟面积 = f(有效烟剂供应速率， 发动机射流初始尺度， 风速， 大气稳定度， 可见阈值)。这个模型的优势在于：物理意义清晰，将内部系统参数与外部环境联系起来。提供了一个定量的分析框架，可以通过实验数据对扩散参数和阈值进行校准，从而不断提高预测精度。

4.2 重要局限性：

经验系数：模型包含多个需要实验确定的系数（η_mix, θ, C_threshold， 大气稳定度类别）。

均匀假设：假设风和大气稳定度在空间和时间上是均匀的，现实中很少如此。

简化物理：忽略了颗粒的沉降、蒸发、化学反应等复杂过程。

计算复杂度：完整的计算需要迭代或数值积分，适合编程实现。

4.3 实践建议：

首先使用这个模型进行参数敏感性分析。例如，会发现：风速U是决定长度的最关键外部因素。有效源强S是决定浓度和长度的最关键内部因素。大气稳定度是决定扩散速率（宽度和高度）的最关键环境因素。

然后，在固定地点和典型气象条件下，通过几次实验来测量实际的烟幕尺寸，反向校准模型中的关键参数（如C_threshold），之后你就可以用这个校准后的模型相当准确地预测在其他设定（如不同推力、不同供料速率）下的成烟效果了。

（未完待续）