[发明专利]一种基于伴随同化算法的流场压力反演方法

说明书

本发明涉及流场压力计算技术领域，涉及一种基于伴随同化算法的流场压力反演方法，包括以下步骤：S1、建立工况流动问题的有限体积模型；S2、开展稳态流场计算；S3、在特征表面布置一定数目的压力测量节点；S4、获得与流场数值计算边界条件相一致的实验壁面压力测量结果；S5、建立基于壁面压力驱动的连续伴随优化方法。本发明利用数据同化试图结合实验和数值模拟的优点，通过调整数值模拟中的不确定参数来匹配实验不确定范围内的测量数据，产生既符合流动控制方程又与实验数据非常相似的解；同时，通过流场的合理外推，可以在更大尺度上测量到研究对象的流动形态。

技术领域

本发明涉及流场压力计算技术领域，尤其涉及一种基于伴随同化算法的流场压力反演方法。

背景技术

流动数值模拟在工程界得到了广泛的应用。然而，由于求解过程中存在的诸多问题：如不准确的湍流模型、错误的数值设置(包括几何和边界条件)、不准确的来流湍流强度以及数值耗散等等，导致数值模拟方法的准确应用大打折扣。为弥补这一缺陷，通常采用实验方法验证数值模拟结果的准确性。然而，由于测量误差、物理模型或风洞参数上的缺陷(变形、振动、上游条件不均匀等)，实验结果也容易因此而产生不确定性。同时，受测量视野限制，实验值在时间或空间上往往是欠解析的，或者只包含关于整个流场的部分信息。

总的来说，已有专利主要集中在通过速度测量重建流场，只有有限的研究提到了使用壁压作为测量值的可能性。由于壁压在实际应用中比速度更容易获得，它是提高流场重建性能的一个很有前景的信息源。然而，当前采用壁压测量来进行数据同化时，算法即使准确地恢复了壁面压力系数，也很难足以重建分离区域的流动。这一缺陷主要由两个方面引起，1、壁压对远离壁面的流动影响有限，这可能是由于局部流动和壁面测量之间的相关性随着距离的增加而衰减，2、由于不合适的雷诺应力分解形式，平均压力场不能在基于稳态涡流粘度的动量方程中准确恢复；这些问题将导致对平均流场压力的预测不准确，从而导致速度的预测不准确。

发明内容

本发明所要解决的技术问题是：利用数据同化试图结合实验和数值模拟的优点，通过调整数值模拟中的不确定参数来匹配实验不确定范围内的测量数据，产生既符合流动控制方程又与实验数据非常相似的解；重建后的流场能够准确恢复实验中无法或难以测量的量，如局部压力脉动、涡量和雷诺应力；同时，通过流场的合理外推，可以在更大尺度上测量到研究对象的流动形态。

本发明所采用的技术方案：一种基于伴随同化算法的流场压力反演方法，包括以下步骤：

S1、建立工况流动问题的有限体积模型，采用网格划分软件对工况计算域划分网格；

S2、在给定的流场边界条件下，结合涡粘模型，基于稳态RANS方程开展稳态流场计算；

具体的，稳态流场计算的结果包括速度，压力以及涡粘性系数的全场分布；

S3、开展工况流动问题的实验测量，在工况特征表面布置一定数目的压力测量节点，将通过压力传感器测量得到的压力值记为p_exp；

进一步的，压力测量节点数范围为10-20个，

S4、获得与流场数值计算边界条件相一致的实验壁面压力测量结果；其中，必须确保实验测量和流场计算的边界条件相一致；

S5、建立基于壁面压力驱动的连续伴随优化方法，以S2计算得到的稳态流场作为输入，开展迭代求解，反演得到最终的流场分布，具体包括如下步骤：

S51、流场计算的原始方程具体如下：

1)、首先给出基于涡粘模型的数值计算的原始控制方程：

其中，U_i是平均速度矢量，υ是运动粘度，ρ是流体密度，p是平均压力，f_i是雷诺应力矢量；