一种考虑随机和认知不确定性的CFD模型修正方法

一种考虑随机和认知不确定性的cfd模型修正方法
技术领域
1.本发明涉及cfd模型修正技术领域，具体涉及一种考虑随机和认知不确定性的cfd模型修正方法。


背景技术：

2.近年来随着高性能计算机技术的迅速发展，计算流体力学(computational fluid dynamics,cfd)由于具有成本低、周期短、可模拟各种工况等众多优势，在航空航天领域产品设计和性能评估中已经成为不可或缺的重要手段。在实际的cfd数值模拟中，普遍存在着诸如来流条件、几何偏差等客观存在的随机不确定性(aleatory uncertainty)以及湍流模型及其系数等认知不确定性(epistemic uncertainty)，使得cfd数值模拟结果与试验测量结果或多或少地存在偏差。
3.认知不确定性本质上由于知识缺乏或数据不足而导致，可加以控制减小其对系统输出的影响，在系统受随机和认知混合不确定性影响下，确定性情况下确认好的cfd仿真模型极有可能对不确定性非常敏感，对于新的流动情况极有可能带来很大误差，而将认知不确定性当作随机不确定性必然存在分布及分布参数的假设偏差，可见传统基于经验和试凑的模型修正方法存在不足，因此需要一种混合不确定性影响下的cfd数值模拟模型修正方法，提高cfd的可信度和预测能力。


技术实现要素：

4.本发明提供了一种考虑随机和认知混合不确定性的cfd模型修正方法，能够有效量化混合不确定性下cfd输出的不确定度，并修正具有认知不确定性的模型参数，使得修正后的模型参数更好地反映真实流动状况和流场流动特性，构建高保真的cfd仿真模型，提高cfd数值预测的可信度。
5.本发明的一种考虑随机和认知混合不确定性的cfd模型修正方法，包括如下步骤：将混合不确定性下的不确定度评估、模型确认、灵敏度分析以及参数修正的闭环流程反复迭代，直至模型满足精度要求且迭代修正满足收敛条件；
6.其中，所述不确定度评估过程为：对cfd数值模拟中的认知不确定性和随机不确定性进行梳理，基于证据理论统一表征；在证据理论框架下完成混合不确定性的综合传播和系统响应的不确定度量化；所述模型确认过程为：建立基于距离准则的模型确认度量指标，判断cfd的预测结果和试验结果的吻合度，若cfd预测精度不满足要求，则开展混合不确定性下的全局灵敏度分析；通过所述灵敏度分析发掘出随机不确定性存在时对cfd响应影响较大的认知不确定性参数，所述灵敏度分析为基于概率包络面积变化率的全局灵敏度分析；所述参数修正具体方式为：基于似然样本策略对发掘出的参数开展模型参数修正。
7.其中，所述不确定度评估过程中，构建cfd数值模拟的代理模型，预测输出响应。
8.其中，所述不确定性量化具体过程为：
9.步骤1-1：收集认知不确定性变量信息，构建识别框架根据工程经验或权
威专家预测确定基本可信度分配其中表示第i个认知不确定性变量的第k个焦元；
10.步骤1-2：收集随机不确定性变量信息，并通过概率分布转化为证据理论表示：对于有界分布，将变量全部取值[a,b]等分为n个子区间[ak,bk](1≤k≤n)，由概率密度函数计算各个子区间的bpa；对于无界分布，将变量取值截断后按照有界分布作相同处理；
[0011]
步骤1-3：假设各维变量均相互独立，通过单个证据变量的焦元和bpa确定联合焦元a
l
和联合bpa为：联合焦元的总数n
total
为：
[0012]
式中，1≤l≤n
total
，其中和分别表示认知和随机不确定性变量的焦元个数；分别是m+n维证据变量第l组焦元的下界和上界向量；
[0013]
步骤1-4：在每个联合焦元a
l
上进行极值分析，获得在联合焦元a
l
上输出响应的最大值和最小值，即响应焦元y
l
，进一步确定在所有随机和认知不确定性变量变化范围内的y的全局最大值和最小值；其中，构建cfd数值模拟模型的代理模型，并结合寻优算法获取极大和极小值；
[0014]
步骤1-5：取阈值为v，且v满足g
min
≤v≤g
max
；令阈值区间gv＝[g
min
,v]，判断响应焦元y
l
和阈值区间gv两集合间的包含关系，计算出信任函数bel和似然函数pl，分别对应了系统响应不确定性度量的下界和上界，具体地：
[0015][0016][0017]
步骤1-6：设阈值v递增的步长为h，令v＝v+h，重复执行步骤1-5，在v从g
min
递增至g
max
的过程中，得到一系列分析结果[v,bel(gv)]和[v,pl(gv)]，由此获得累积信任函数cbf和累积似然函数cpf。
[0018]
其中，基于概率包络面积变化率的全局灵敏度分析具体过程为：
[0019]
步骤2-1：进行包含所有不确定性变量的混合不确定性量化，得到系统输出响应的总概率包络，并根据pl和bel的离散差值采用梯形数值积分的方法计算总的概率包络面积s
t
；
[0020]
步骤2-2:选择待评估的不确定性变量进行剔除，再次进行混合不确定性量化，并计算概率包络面积；
[0021]
步骤2-3：计算剔除不确定性变量前后的概率包络面积变化率，按照概率包络面积变化率的灵敏度分析指标得出变量的灵敏度系数；
[0022]
步骤2-4：重复执行步骤2-2和步骤2-3，将计算出来的所有随机和认知不确定性变量的灵敏度系数分别排序。
[0023]
其中，所述步骤2-2中，对于类型不同的随机和认知不确定性变量采用以下方法进
行剔除：
[0024]
1)随机不确定性变量：用均值代替；
[0025]
2)认知不确定性变量：用概率加权中值代替。
[0026]
其中，基于似然样本策略的模型参数修正的具体过程为：
[0027]
步骤3-1：根据灵敏度分析结果，选择灵敏度系数最大的认知不确定性参数作为待修正的模型参数，对排序靠后的认知不确定性参数做忽略处理；
[0028]
步骤3-2：对随机不确定性参数按其概率分布进行随机抽样，对认知不确定性参数按其基本可信度分配进行重要性抽样；
[0029]
步骤3-3：将步骤3-2产生的输入样本一一代入cfd数值模拟，计算其响应预测数据；
[0030]
步骤3-4：建立cfd数值模拟结果和试验结果间相对误差的距离指标，来量化二者之间的差异；
[0031]
步骤3-5：对所有cfd预测数据点的相对误差从小到大进行排列，若相对误差满足rank(rε)≤λns，说明该预测数据与试验数据吻合较好，则根据cfd预测值反推其对应的输入样本作为似然样本，并予以保留；式中，rank表示排序；λ为截断比例；ns是所有样本数量；
[0032]
步骤3-6：根据上一步的似然样本，确定待修正参数的上、下边界，以此作为其新的证据区间，相应的bpa与初始各焦元概率质量保持一致；对于无需修正的认知不确定性参数以及随机不确定参数，其不确定性表征模型均保持不变，进行下一轮的模型修正迭代；
[0033]
步骤3-7：当最大的相对误差不超过设定阈值，认为模型精度满足要求；当待修正认知不确定性参数的证据区间上下界不再发生明显变化，认为修正结果已经收敛。
[0034]
有益效果：
[0035]
本发明提出的混合不确定性影响下的cfd模型修正方法，可快速有效评估随机和认知混合不确定性对cfd数值模拟系统响应的影响，科学地开展模型修正，所提出的涵盖不确定性量化、灵敏度分析、参数修正的模型修正方案和流程，打破了传统基于经验和试凑的模型修正方法的不足，为cfd数值模拟的模型修正提供了系统有效的解决思路。本发明具体涉及一种融合混合不确定性量化、全局灵敏度分析、参数修正策略的模型修正流程。采用基于证据理论的随机和认知混合不确定性量化方法、基于概率包络面积变化率的全局灵敏度分析方法以及基于似然样本策略的模型参数修正方法，利用证据理论对随机和认知混合不确定性进行综合量化，在证据理论框架下建立概率包络面积变化率的混合不确定性灵敏度分析指标，采用似然样本策略对模型参数进行修正。通过混合不确定度量化得到cfd输出的不确定性区间表示，基于此开展全局灵敏度分析发掘对输出影响较大的模型参数，并根据似然样本策略对其加以修正，使得cfd预测结果与试验结果尽可能吻合。
[0036]
本发明构建cfd数值模拟模型的代理模型，在极值分析过程中，不需要不断调用cfd数值模拟计算，并结合寻优算法获取极大和极小值，计算成本很低。
附图说明
[0037]
图1为本发明考虑混合不确定性的cfd模型修正方法技术路线示意图。
[0038]
图2为本发明基于证据理论的随机和认知混合不确定性量化的流程图。
[0039]
图3为本发明基于概率包络面积变化率的灵敏度分析流程图。
[0040]
图4为本发明基于似然样本策略的模型参数修正流程图。
[0041]
图5为本发明证据理论框架下系统响应的概率分布示意图。
[0042]
图6为本发明概率包络面积变化率指标示意图。
[0043]
图7为本发明参数修正的似然样本策略原理图。
具体实施方式
[0044]
本发明考虑混合不确定性的cfd模型修正方法其技术路线如图1。首先对cfd数值模拟中的认知不确定性(如：湍流模型封闭系数)和随机不确定性(如：来流马赫数和攻角、几何偏差)进行梳理，基于证据理论统一表征；其次，在证据理论框架下完成混合不确定性的综合传播和系统响应的不确定度量化，同时为有效降低计算量，考虑构建cfd数值模拟的代理模型，快速预测输出响应；然后，建立基于距离准则的模型确认度量指标，判断cfd的预测结果和试验结果的吻合度，以决定是否需要对模型参数开展修正；若cfd预测精度不满足要求，则开展混合不确定性下的全局灵敏度分析，发掘出随机不确定性存在时对cfd响应影响较大的认知不确定性参数，基于似然样本策略对该参数开展修正。上述混合不确定性下的不确定度评估、模型确认、灵敏度分析、参数修正的闭环流程将反复迭代，直至模型满足精度要求且迭代修正满足收敛条件，最后将修正好的cfd数值模拟模型用于工程产品的性能分析和优化设计。
[0045]
设响应函数y＝g(xe,xa)，其中表示m维认知不确定性向量，表示n维随机不确定性向量。本发明提供的一种混合不确定性影响下的cfd模型修正方法包括三个核心部分——基于证据理论的随机和认知混合不确定性量化、基于概率包络面积变化率的全局灵敏度分析以及基于似然样本策略的模型参数修正。
[0046]
其中，图2为本发明基于证据理论的随机和认知混合不确定性量化的流程图。基于证据理论的随机和认知混合不确定性量化具体过程为：
[0047]
步骤1-1：收集认知不确定性变量信息，构建识别框架，根据工程经验或权威专家预测确定基本可信度分配表示第i个认知不确定性变量的第k个焦元)。
[0048]
步骤1-2：收集随机不确定性变量信息，并通过概率分布转化为证据理论表示。对于有界分布，将变量全部取值[a,b]等分为n个子区间[ak,bk](1≤k≤n)，由概率密度函数计算各个子区间的bpa，如式(1)和式(2)；对于无界分布，将变量取值截断后按照有界分布作相同处理。
[0049][0050][0051]
式中，表示第j个随机不确定性变量的第k个焦元，pdf(xa)表示随机变量xa的概率密度函数。
[0052]
步骤1-3：假设各维变量均相互独立，通过单个证据变量(由证据理论统一表示的随机或认知不确定性变量)的焦元和bpa确定联合焦元a
l
和联合bpa如式(3)和式(4)，其中联合焦元的总数n
total
如式(5)。
[0053][0054][0055][0056]
式中，1≤l≤n
total
，其中和分别表示认知和随机不确定性变量的焦元个数；分别是m+n维证据变量第l组焦元的下界和上界向量。
[0057]
步骤1-4：在每个联合焦元a
l
上进行极值分析，获得在联合焦元a
l
上输出响应的最大值和最小值，即响应焦元y
l
，如式(6)。进一步确定在所有随机和认知不确定性变量变化范围内的y的全局最大值和最小值[g
min
,g
max
]。
[0058][0059]
由上式可知，极值分析需要不断调用cfd数值模拟计算，计算成本难以接受，因此构建cfd数值模拟模型的代理模型，并结合寻优算法获取极大和极小值。
[0060]
步骤1-5：取阈值为v，且v满足g
min
≤v≤g
max
。令阈值区间gv＝[g
min
,v]，判断响应焦元y
l
和阈值区间gv两集合间的包含关系，根据式(7)和式(8)计算出信任函数bel和似然函数pl，分别对应了系统响应不确定性度量的下界和上界。
[0061][0062][0063]
步骤1-6：设阈值v递增的步长为h，令v＝v+h，重复执行步骤1-5，在v从g
min
递增至g
max
的过程中，得到一系列分析结果[v,bel(gv)]和[v,pl(gv)]，由此获得累积信任函数cbf和累积似然函数cpf，其示意如图5。
[0064]
图3为本发明基于概率包络面积变化率的灵敏度分析流程图。基于概率包络面积变化率的全局灵敏度分析具体过程为：
[0065]
步骤2-1：按照基于证据理论的随机和认知混合不确定性量化方法，进行包含所有不确定性变量的混合不确定性量化，得到系统输出响应的总概率包络，并根据pl和bel的离散差值采用梯形数值积分的方法计算总的概率包络面积s
t
。
[0066]
步骤2-2:选择待评估的不确定性变量xk(1≤k≤m+n)进行剔除(或裁剪)，再次进行混合不确定性量化，并计算概率包络面积对于类型不同的随机和认知不确定性变量采用以下方法进行剔除(或裁剪)：
[0067]
1)随机不确定性变量：用均值代替。
[0068]
2)认知不确定性变量：用概率加权中值代替。
[0069]
概率加权中值的计算方式如下：
[0070][0071]
式中，n(xk)表示认知不确定性变量的焦元个数；m(a
ik
)、分别是该变量第i个焦元的概率质量、下界和上界。
[0072]
步骤2-3：计算剔除(或裁剪)不确定性变量xk前后的概率包络面积变化率，按照下式所示的概率包络面积变化率的灵敏度分析指标得出变量xk的灵敏度系数sk，直观示意如图6。
[0073][0074]
式中，sk表示第k(1≤k≤m+n)个不确定性变量的灵敏度系数；s
t
表示考虑所有不确定性变量的情况下，对系统输出响应进行不确定性量化后所得的总概率包络面积；表示剔除(或裁剪)变量xk的不确定性后，再次对系统输出响应进行混合不确定性量化所得的新概率包络面积。
[0075]
步骤2-4：重复执行步骤2-2和步骤2-3，将计算出来的所有随机和认知不确定性变量的灵敏度系数分别排序。
[0076]
图4为本发明基于似然样本策略的模型参数修正流程图。基于似然样本策略的模型参数修正的具体过程为：
[0077]
步骤3-1：根据灵敏度分析结果，选择灵敏度系数最大的认知不确定性参数作为待修正的模型参数，对排序靠后的认知不确定性参数做忽略处理以实现降维。
[0078]
步骤3-2：对随机不确定性参数按其概率分布进行随机抽样，对认知不确定性参数按其基本可信度分配进行重要性抽样，基本可信度分配高的子区间会生成更多的样本。
[0079]
步骤3-3：将步骤3-2产生的ns个输入样本一一代入cfd数值模拟，计算其响应预测数据y
p
。
[0080]
步骤3-4：考虑到实际工程中由于成本和时间的限制，往往无法进行完整的确认试验，试验数据较为稀缺，因此无法对cfd数值模拟开展不确定性下的精细化模型确认。如式(11)，建立cfd数值模拟结果和试验结果间相对误差(relativeerror,re)的距离指标，来量化二者之间的差异，从而对cfd数值模拟进行简单确认。
[0081][0082]
式中，ye是试验数据，n
′
是试验数据个数，ns是所有样本数量。
[0083]
步骤3-5：对所有cfd预测数据点的相对误差从小到大进行排列，若相对误差满足式(12)，说明该预测数据与试验数据吻合较好，则根据cfd预测值反推其对应的输入样本作为似然样本，并予以保留，其原理如图7。
[0084]
rank(rε)≤λns(12)
[0085]
式中，rank表示排序；λ为截断比例，决定了似然样本的总数λns，根据初始预测数
据与试验数据的差异程度的不同，λ取值在[1％,10％]范围内为宜。
[0086]
步骤3-6：根据上一步的似然样本，确定待修正参数的上、下边界，以此作为其新的证据区间，相应的bpa与初始各焦元概率质量保持一致。对于无需修正的认知不确定性参数以及随机不确定参数，其不确定性表征模型均保持不变，进行下一轮的模型修正迭代。
[0087]
步骤3-7：当最大的相对误差不超过ε，认为模型精度满足要求；同时考虑到迭代修正的收敛性，当待修正认知不确定性参数的证据区间上下界不再发生明显变化，认为修正结果已经收敛。
[0088][0089]
式中，ε和δ为设定的收敛阈值，j为迭代数；a
low
、a
up
分别是待修正参数xk证据区间的下界和上界。
[0090]
为了验证本发明所提出方法的有效性，特将本发明所提出的模型修正方法应用于跨音速三维机翼onera m6的cfd气动特性预测问题，考虑sa湍流模型中的6个封闭系数c
b1
、c
b2
、c
v1
、c
w2
、c
w3
、σ具有认知不确定性，来流条件中马赫数ma和攻角α具有随机不确定性，综合分析混合不确定性对机翼升力系数c
l
的影响。m6机翼的几何建模参考nasa官网公布的数据，cfd计算由fluent软件实现，将升力系数cfd仿真结果与文献参考试验结果对比。
[0091]
通过对spalart allmaras(sa)模型湍流系数进行灵敏度分析，发现对升力系数影响最大的是c
b1
，因此对其进行修正。表1展示了模型修正前后c
b1
的迭代结果，通过对比参考试验数据与进行模型修正前与模型修正后的升力系数的不确定性区间，可以发现修正后升力系数的不确定度明显减小，sa模型系数c
b1
经过修正后取值位于0.09附近，优于fluent软件的推荐值0.1355。取式(13)中ε＝1％、δ＝0.001，两次迭代即满足收敛条件，cfd的预测结果十分接近参考试验数据，提出的模型修正方法具有很快的收敛速度。依据该实施例，可以证明本专利公开的一种混合不确定性影响下的cfd模型修正方法是有效的。
[0092]
表1模型修正的迭代结果
[0093][0094]
以上为本发明的较佳实施实施例，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。