本研究由来自TOBB University of Economics and Technology的Selin Aradag、Airbus的Hyung Jo Kim以及Rutgers, The State University of New Jersey的Doyle D. Knight共同完成，并于2010年发表在期刊《Engineering Applications of Computational Fluid Mechanics》第4卷第4期，文章标题为“Two and Three Dimensional Simulations of Supersonic Cavity Configurations”。

本研究属于计算流体力学与空气动力学交叉领域，专注于超音速流动的数值模拟。研究的背景源于航空航天工程中的一个重要实际问题：飞行器内部武器舱（通常可简化为空腔构型）在超音速飞行时，流经空腔的流动会产生复杂的非定常流场，引发剧烈的压力波动和共振声学模态。这些振荡可能导致飞机和所携带武器（stores）的结构损伤，并阻碍武器的成功分离与释放。因此，深入理解空腔流动的物理机制，对于开发有效的流动控制方法以抑制振荡、保证安全至关重要。尽管已有大量研究，但空腔流动，尤其是开放空腔流动，仍有许多物理细节有待厘清。计算流体动力学（CFD）是可视化并分析这种高度非定常流动的有力工具。本研究设定了两个明确目标：第一，展示CFD在准确预测超音速空腔流动方面的能力；第二，评估二维模拟在捕捉空腔流动物理方面的充分性。为此，研究团队对超音速流经一个开放空腔的构型，同时进行了二维和三维的时均非定常雷诺平均纳维-斯托克斯（RANS）模拟，并将结果相互对比，同时与文献中的实验数据进行验证。

研究的详细工作流程如下。首先，研究团队确定了用于模拟的几何构型与流动条件。他们选择了Kaufman等人（1983）实验中的一个案例作为基准，以确保有可靠的实验数据用于对比。该空腔为矩形，长深比（L/D）为5.07，长宽比（L/W）为1.90，属于典型的开放空腔构型（L/D < 10）。来流马赫数为1.5，基于空腔长度的雷诺数为109万，空腔上游的边界层为湍流。其次，在计算方**面，研究使用了商业CFD求解器GASPEX（版本4.1.0）。该求解器采用有限体积法对三维可压缩时均非定常RANS方程进行空间离散，时间和空间上均具有二阶精度。湍流模型选用Wilcox的k-ω两方程模型，并采用萨瑟兰粘度定律。由于构型关于空腔中心线对称，计算域仅取物理域的一半以节省计算资源。边界条件的设置非常关键：入流边界条件通过数值求解平板在来流条件下的二维定常流动方程获得（使用程序EDDYBL），并通过三次样条插值映射到空腔计算域，确保了入流边界层厚度与实验匹配。固体壁面采用无滑移条件，出流和上边界采用一阶外推，中心线为对称条件，侧向边界（z/W=1处）也采用一阶外推。初始条件在空腔口上方设为入流条件，空腔内部则设为壁面值。非定常求解采用隐式双时间步进法，物理时间步长为1x10-6秒，对应约0.0012个无量纲罗斯特周期，以确保每个罗斯特周期内有约1000个时间步来充分解析非定常行为。数值通量格式采用Roe格式，并使用了minmod限制器。第三，研究构建了非均匀笛卡尔网格，计算域分为空腔内、外两个区域。网格在壁面附近进行了加密，确保了Δy+和Δz+小于1，满足精确解析粘性底层的准则。在入流边界层内约有24个网格点，剪切层下游面通常有25-30个网格点。总网格量对于内部区域为95x54x35，外部区域为177x62x70（x， y， z方向）。最后，在数据处理与分析方面，研究比较了二维和三维模拟得到的时均量，包括空腔底板和下游面的声压级（SPL）分布，以及特定测压点（位于下游面y/D=-0.4处）的压力历程。通过对压力历程进行离散傅里叶变换（DFT），得到了功率谱以分析主导频率。此外，研究还通过可视化手段分析了流场结构，包括中心面的流线图、数值纹影图（密度梯度等值线）、涡量图，以及垂直于流向的横截面上的流线、压力云图和纹影图，以揭示流动的三维特征。

研究的主要结果如下。首先，在压力振荡特性方面，二维和三维模拟显示出显著差异。位于下游面测压点的压力历程表明，二维模拟中的压力波动幅值明显大于三维模拟。功率谱分析进一步揭示，在三维模拟中，主导频率是基频（约1066Hz），而在二维模拟中，主导频率则变为第一谐波（约2133Hz）。其次，在时均声压级分布上，空腔底板、前壁和后壁的SPL分布趋势相似，但二维模拟的结果与实验数据偏差更大，其SPL值整体更高，这再次印证了二维模拟夸大了压力振荡的强度。第三，在流场结构上，中心面的流线图显示，三维模拟中空腔下游部分主要由一个大涡占据，而二维模拟中则同时存在两到三个共存的涡结构。数值纹影图和涡量图清晰地展示了流动的反馈机制：来流边界层在空腔前缘分离，剪切层因开尔文-亥姆霍兹不稳定性产生涡旋，这些涡向下游运动并撞击后壁，产生压力波，压力波向上游传播再次影响剪切层的稳定性，形成一个自持的振荡循环。第四，也是本研究重点阐述的，是流动的三维特征。分析发现，空腔口沿展向（z方向）的压力分布并不均匀，下游压力高于上游，且从中心线向两侧迅速降低，这直接证明了流动的三维特性。在垂直于流向的横截面（如x/L=0.5平面）上，可以观察到比中心面更厚的剪切层，并且存在明显的展向流动和流向涡结构。这些涡结构对应着空腔内部的低压区。随着时间推移，剪切层在横截面内发生上下摆动和变形，高、低压区向空腔侧壁移动，伴随着涡的生成、演变和耗散，这些复杂的三维现象是二维模拟完全无法捕捉的。

基于以上结果，研究得出了明确的结论。在超音速空腔流动这一特定问题中，三维效应不可忽略。虽然二维模拟在计算速度和内存消耗上具有巨大优势，并且能够提供关于流动纵向模态的初步认识，但它不足以理解问题的真实物理本质。三维模拟揭示的展向流动简化了空腔内的流向涡结构，从而削弱了壁面的压力波动幅值，但也使得压力演化因展向流体环流的参与而更加不规则和难以预测。空腔的宽度与长度量级相当，这意味着横向声波与纵向波在流动结构中的参与度是相当的。因此，对于本研究采用的构型（L/W=1.90），进行二维模拟是不合理的，也不应使用二维模拟来评估CFD模拟空腔流动的准确性。不过，研究也指出，如果通过增加空腔宽度来减小L/W的比值，二维模拟可能仍有其参考价值。

本研究的科学价值在于，通过系统的对比模拟，定量和定性地揭示了二维模拟在预测超音速空腔流动振荡特性方面的局限性，并深入剖析了导致这种差异的三维物理机制。其应用价值在于为航空航天工程中武器舱流动的CFD模拟实践提供了重要指导：在进行此类问题的工程设计或流动控制方案评估时，应优先或必须采用三维模拟，以避免基于不准确的二维结果做出误判。研究强调了在追求计算经济性与保证物理准确性之间需要做出的权衡。

本研究的亮点在于：第一，研究设计直接而有力，通过并行执行同一构型下完全对应的二维和三维RANS模拟，并与实验数据交叉对比，使结论非常具有说服力。第二，对三维流动特征的分析非常深入和全面，不仅关注了中心面的传统流场，还通过多个横截面的瞬态可视化，生动地展示了展向剪切层发展和流向涡结构的复杂时空演化，这是理解三维效应关键所在。第三，结论明确且具有实践指导意义，明确指出对于特定几何比例的空腔，二维模拟不足以胜任，这对CFD工程应用社区是一个重要的提醒。第四，研究流程规范，对网格独立性、边界层分辨率、湍流模型、数值格式等CFD关键设置均有详细说明，确保了模拟的可信度。此外，研究还为进一步的流动控制研究（如文中提及的激光能量沉积控制）提供了可靠的基准流场和深入的理解基础。