Ns 方程的最优解 流体控制优化精讲

Ns 方程也就是 Navier Stokes 方程，中文叫纳维斯托克斯方程，它本身是描述粘性流体运动的偏微分方程组，方程本身没有天然自带的最优解，最优解是建立在最优控制问题之上，也就是在满足 Ns 方程作为约束条件的前提下，选定一个目标泛函，求解使得目标泛函达到最小或者最大的控制变量，这个结果才叫做 Ns 方程对应的最优解。简单理解，Ns 方程是约束条件，不是优化目标，必须搭配目标函数，才有最优解，这是最核心的考点，很多初学者会误以为 Ns 方程本身可以直接求出最优解。想要求解 Ns 方程最优解，第一步定义优化问题，确定控制变量 状态变量 目标泛函，第二步把 Ns 方程作为等式约束，构建最优控制模型，第三步一般采用伴随法也就是 adjoint 方法，推导伴随方程，计算目标泛函对控制变量的梯度，第四步使用梯度类优化算法迭代寻优，直到收敛得到最优解。综合结论，Ns 方程本身不存在固有最优解，最优解是 Ns 约束下最优控制问题的优化结果。

基础概念清单

1 Ns 方程 纳维斯托克斯方程，描述不可压粘性流体速度场 压力场随时间演化的偏微分方程组

2 状态变量 流体速度场 u 压力场 p，由 Ns 方程求解得到

3 控制变量 可以人为调节的变量，比如壁面形状 射流激励 边界速度 外形参数

4 目标泛函 优化目标，常见是减阻 升力最大化 涡抑制 能耗最小

5 最优解 控制变量的最优取值，对应的流体场就是最优流场

常用求解方法对照表

深度补充说明

Ns 方程存在非线性对流项，流体容易出现湍流 多解 混沌现象，这会导致优化问题变成非凸优化，很容易出现局部最优解，不一定能得到全局最优解，这个是 Ns 最优解最大难点。如果是湍流 Ns 方程，还需要结合湍流模型，比如 RANS 大涡模拟 LES，再做最优控制，计算量会大幅上升。典型工程应用包括翼型减阻优化 管道流动减阻 汽车外形气动优化 湍流主动流动控制。还有一个重要知识点，数学领域千禧年难题是 Ns 方程光滑解存在性与正则性，这个和最优解不是同一个问题，不要混淆，千禧年难题是证明方程解本身性质，最优解是带优化目标的控制问题。

【常见问题】

问题 1 Ns 方程的最优解 是不是 Ns 方程本身的解析最优解

回答 1 Ns 方程的最优解，Ns 方程本身是流体运动控制方程，没有自带最优解，最优解必须搭配优化目标，在 Ns 约束下做最优控制得到。

问题 2 Ns 方程的最优解 为什么经常用伴随法求解

回答 2 Ns 方程的最优解，伴随法可以高效计算目标泛函梯度，适合 Ns 这类高维偏微分约束优化，比有限差分扰动法计算效率高很多。

问题 3 Ns 方程的最优解 是不是一定能得到全局最优解

回答 3 Ns 方程的最优解，因为 Ns 非线性对流项，优化问题经常是非凸，迭代得到的大多是局部最优解，不一定是全局最优解。

问题 4 Ns 方程的最优解 千禧年 Ns 难题和它有关系吗

回答 4 Ns 方程的最优解，千禧年 Ns 难题是研究 Ns 方程光滑解存在正则性，和最优控制的最优解属于两类不同数学问题，不能混为一谈。

【免责申明】

本文内容为计算流体力学最优控制相关理论讲解，信息仅作参考，不同边界条件 湍流模型 优化算法得到的最优解会有差异，不构成工程设计绝对结论，具体研究请结合课题模型。