翼型 — CAE术语解説

翼型是指翼沿着与进行方向平行的平面切割时的断面形状。它是二维形状。例如，NACA 0012是对称翼，最大厚度是弦长（从前缘到后缘的直线距离）的12%。弦长是表征翼型最基本的尺寸，是确定翼型特征的基准。

很好的问题。如果知道了二维翼型的特性（升力系数、阻力系数、俯仰力矩系数），就能推定三维翼的特性。特别是，对于展弦比大的（细长）翼，可以近似为二维翼型在翼根到翼尖沿跨向叠加的形式。支配方程的基础是二维、不可压缩、非粘性条件下的Kutta-Joukowski定理。

主要参数是迎角（攻击角）、雷诺数、马赫数以及翼型形状本身。例如，典型的NACA 4412翼型，在低雷诺数（Re=10^6）、迎角0度时的升力系数

基本假设是二维、不可压缩、定常流。计算域要在翼型前后各留出足够的长度。例如，从前缘到入口边界20倍弦长，从后缘到出口边界30倍弦长，上下边界各20倍弦长。求解器采用基于压力的定常求解器（如SIMPLE法），湍流模型从相对稳定的SST k-ω模型开始。

为了分辨边界层，需要在翼表面贴棱柱层（充气层）。第一层厚度由无量纲壁距

是重要的监视量。但仅靠数值稳定还不够，要确认残差至少下降3～4个数量级。还要目视检查翼表面压力分布、后流涡结构是否物理上合理。在Ansys Fluent中设置"Force Monitors"，绘制升力系数和阻力系数的时间序列，确认周期性变动消失、趋向定常值。

首先要明确目的。是"想知道最大升力系数"还是"想评估设计迎角下的效率（升阻比）"，这决定了解析范围。之后按以下步骤：1) 获取或生成翼型坐标数据（.dat文件等）。2) 用简易面板法代码（如XFOIL）粗略了解特性。3) CFD前处理：用粗网格测试计算域和边界条件。4) 进行网格无关性研究。5) 本体计算和结果验证。

系统地改变全局网格大小（背景网格要素尺寸）和翼周边界层网格（层数和增长率）。例如，将要素数改为10万、50万、200万，观察目标量（如迎角5度的升力系数）的变化是否控制在1%以内。同时检验壁面剪应力分布、后缘附近压力梯度在细化网格时是否不再变化。这虽然有成本，但是结果可信度的依据。

即使没有实验数据，也有几个基准可用。首先，与美国NASA兰利研究中心公开的NACA 0012翼型广泛的实验和计算数据库（如"NPARC Alliance Validation Archive"）对比。其次，用不同的CFD软件（如OpenFOAM和商用软件）在相同条件下计算交叉检验。再次，与面板法（非粘性）结果对比，定性判断粘性影响（阻力、失速特性）是否正确计算。

大体分为简易工具、高级CFD、专用优化工具。简易工具的代表是XFOIL。它结合了面板法和边界层方程，既考虑非粘性也考虑粘性，笔记本电脑上几秒就能出结果。但在逆压梯度大、流动分离的条件下精度下降。

本质计算精度没有大差别，但工作流程和成本差异很大。Ansys Fluent（特别是Workbench环境）界面功能完整，可以从网格生成到后处理一体化进行。OpenFOAM是开源免费的，但主要以命令行和配置文件方式操作。研究开发要定制求解器的选OpenFOAM，设计业务要稳定重复解析的选Fluent。两者都配置了航空领域标准的湍流模型如S-A模型和SST k-ω模型。

有。例如Dassault Systèmes的Isight或ESTECO的modeFRONTIER等过程集成和设计优化（PIDO）工具。这些工具可以组合参数化翼型生成工具（如PARSEC或CST参数化）、评估工具（如XFOIL或OpenFOAM）和优化算法（如遗传算法），自动搜索在指定升力系数下阻力系数最小的翼型形状。

首先怀疑边界条件。流入速度、密度、迎角的设置值正确吗？其次看网格，特别是后缘。后缘钝的（有厚度）翼型如果后缘附近网格太粗，会过度高估分离，导致升力降低。还有湍流模型的选择也影响。如果不考虑转捩（从层流到湍流的过渡），假设整个翼面都是湍流，摩擦阻力会过大，升阻比恶化。SST k-

定常解析中振荡可能是物理上非定常现象（如卡门涡街）发生。首先，如果能改时间步，试试用时间推进法做伪非定常计算。还不行的话，检查求解器设置。在Ansys Fluent中把压力-速度耦合从"SIMPLE"改为"Coupled"可能改善收敛。还可以降低松弛因子（从默认的1改为0.5等），但计算时间会增加。

这是数值发散的前兆。最常见的原因是网格品质差。特别检查翼远处流出边界附近的网格，是否有极端长宽比（细长）的单元。还要确认求解器离散化方案不是一阶精度。对流项应用二阶以上的迎风差分方案（如QUICK、MUSCL等）。初始条件太差的话也会出现，建议用一致的流入条件初始化所有变量后开始计算。

这说明XFOIL的求解法（粘性/非粘性相互作用的迭代计算）发散了。多数情况下是那个迎角时翼上表面发生大规模分离，面板法前提崩溃。对策有两个。一是用"GDES"命令稍微修改翼型形状（如稍微抬起后缘）使流动重新附着。二是放弃那个迎角的解析，换用精度更高的CFD求解器。不要忘记XFOIL只是初步设计工具。