项目概览
项目研究垂直叠置的双翼型如何通过相互作用改变升力、阻力和升阻比。在较小间距下，下翼型引导的流动会影响上翼型；间距扩大后，两者逐渐接近独立工作。因此，布局与攻角需要共同优化，不能把单翼型的优选攻角直接用于双翼型。
课程先开展 NACA 0012 单翼型的组内研究，我的个人案例在此基础上转向双翼型。工作包括数值方法比较、网格与湍流模型选择、流场解释和间距 / 攻角参数分析。报告中的实验对照来自参考数据，并非本人另行完成的风洞实验。
任务与设计要求
| 项目 | 条件或依据 | 类型 / 状态 |
|---|---|---|
| 研究对象 | 两个相同翼型垂直叠置，个人案例以 40 mm 间距作为初始构型。 | 课程任务 |
| 几何约束 | 优化讨论采用最大 150 mm 间距限制，结果应与该约束一起解释。 | 题设约束 |
| 评价指标 | 比较 CL、CD 和 CL/CD，并用压力、速度及流动结构解释性能变化。 | 分析任务 |
| 数值设置 | 比较一阶 / 二阶迎风，结合网格和湍流模型选择后续求解方案。 | 方法要求 |
| 收敛检查 | 报告设定最多 100 次迭代，监测残差与气动系数，不能只看最终一张云图。 | 求解设定 |
我的职责
- 我在单翼型组内研究基础上完成双翼型个人案例，整理离散方法、网格与湍流模型对照，改变攻角和间距，并分析翼型之间的气动干扰。
- 单翼型组内工作及参考实验数据构成方法选择的背景；双翼型参数分析作为个人报告展示。仿真结论限定在报告的二维几何、边界与模型条件内。
完成工作与方法
01 / 先检查离散方法对结果的影响
以 40 mm 间距、0° 攻角工况比较一阶与二阶迎风。报告中一阶方案的 CD 收敛至约 0.160，二阶约 0.143；CL 分别约 −0.025 与 −0.022。收敛曲线显示不同离散方法会改变数值扩散和迭代过程，所以必须先统一方法，再比较几何方案。
后续采用二阶迎风与稳态求解，以获得时间无关的气动系数并控制计算成本。报告通过松弛因子组织稳定性，设定残差降至约 10⁻³ 及系数变化趋于稳定作为检查条件。这里的收敛代表数值迭代满足设定判据，不自动证明网格误差或模型误差可以忽略。
02 / 比较计算域、网格与湍流模型
单翼型背景采用半圆入口与矩形出口，入口速度 45 m/s、出口表压 0 Pa，矩形区域上下边界使用滑移条件。两种域划分方法改变翼面和尾缘附近的单元方向及连续性：第一种在尾缘附近存在网格衔接问题，第二种具有更平滑的粗细过渡，但仍需关注远场的不对称性。
报告将两种网格方法与 k–ω SST、Transition SST 组合，并与参考实验气动系数对比。方法二配合 k–ω SST 在所比较工况下更接近参考数据，因此作为后续设置。这样的比较为方法选择提供依据；它与正式的多级网格独立性或不确定度量化并不是同一件事。
双翼型的翼间窄区域和近壁面需要更细的流动分辨。局部网格影响边界层、压力梯度与尾迹的描述，远场大小和边界也参与数值方案的选择。原始网格图提供了气动系数比较的几何与离散背景。
03 / 用流场说明翼型之间的相互作用
压力场提供翼面受压分布的空间信息，速度场则用于观察翼间加速、流动偏转和尾流。将这些图与 CL、CD 一起看，才能判断升阻比变化来自哪些流动特征。单独增加升力未必提高效率，因为阻力也可能同步增加。
小间距使两个翼型的流动更强烈耦合。下翼型对流动方向的改变可能对上翼型产生有利干扰，但这种关系取决于攻角和间距；不能把某一张云图的现象推广到所有布局。较大间距下相互影响减弱，系统评价逐渐接近对单翼型工作状态的评价。
04 / 在约束内选择攻角与间距
报告通过改变间距和攻角，比较气动指标与流场。当最大间距限制为 150 mm 时，约 8° 攻角在报告工况中得到较好的升阻比。这个结论同时包含几何约束和运行条件，不能只摘录一个攻角作为普适最优值。
更宽布局的优选攻角逐渐减小，并趋向约 5.5°，反映翼型间干扰减弱后向相对独立工作状态的过渡。二维参数分析用于筛选构型；后续可根据实际展长、雷诺数和流动时变特征扩展模型。
结果与要求验证
| 要求 | 采用方法 | 证据或结果 | 验证程度 |
|---|---|---|---|
| 统一可比较的数值方法 | 在相同初始工况比较两种迎风方法 | 二阶方案给出约 CD 0.143、CL −0.022，并作为后续设置。 | 仿真支持 |
| 选择网格与湍流模型 | 比较网格方法及参考实验系数 | 报告选择方法二与 k–ω SST；参考数据不等于本人实测。 | 仿真支持 |
| 解释双翼型干扰 | 将压力、速度和气动系数对应 | 观察布局变化对翼间流动及尾迹的影响。 | 仿真支持 |
| 在 150 mm 约束内优化 | 联合比较间距与攻角 | 约 8° 为报告条件下优选，较宽配置趋近约 5.5°。 | 仿真支持 |
收获与改进
让数值设置服务于可靠的几何比较
同一 40 mm 间距、0° 攻角下，一阶与二阶迎风的阻力系数约为 0.160 和 0.143。这个对照帮助我理解数值扩散对气动结果的影响，也成为后续统一采用二阶方法的依据。网格、尾缘过渡和计算域与离散格式共同影响压力及尾迹描述，因此我会围绕最影响方案选择的工况安排计算资源。下一轮优先对候选优选点及其相邻工况进行网格加密，记录升力、阻力和升阻比的变化，以数值稳定性支持几何排序，并确定值得进一步细化的参数区域。
双翼型的最佳攻角来自相互干扰的变化
报告在最大 150 mm 间距约束下选择约 8°，更宽布局的优选攻角则趋向约 5.5°。这两组结果使我把攻角视为与间距耦合的变量：翼间流动偏转会改变局部来流和压力分布，增大间距后，两个翼型逐渐接近独立工作。下一轮分析会分别积分上下翼型的升阻力，并对照表面压力和尾迹，判断总升阻比改善由哪一片翼型贡献。随后在候选区域细化攻角与间距组合，同时观察最优点附近的变化幅度，选择性能较好且对安装偏差更稳健的布局。
按流动特征安排模型的扩展顺序
稳态压力和速度场为理解翼间加速、流动偏转和尾迹提供了直接线索，也让二维模型能够高效筛选几何组合。随着攻角和间距变化，分离及尾迹相互作用可能增加气动力的时间变化，因此后续模型深度可以围绕候选方案的流动特征安排。我会先结合残差与气动力历程定位值得关注的工况，再通过瞬态计算分析平均载荷及波动，并检查时间步与统计时长。若进入实体设计，则进一步按实际展长研究三维效应，使计算工作逐步服务于稳定性和结构载荷评价。
项目来源
cfd 作业.pdf（与 cfd.pdf 内容相同，合并展示），重点依据网格、流场与布局优化章节。
UNSW 官方课程介绍 ↗