基于伴随网格自适应研究.pdfVIP

第六届两北地区计算物理学术会议论文集 2008 基于伴随的网格自适应研究 杨振虎 (中国航空工业第六三一研究所气动数值模拟航空科技重点实验室) 摘要网格自适应是一种有效消除计算结果对网格依赖性，提高计算精度的有效办法，本文给出了一种基于伴 随方法的后验误差估计网格自适应方法，其传感器设计为控制方程的加权残差，权值取为与大多数工程关心的整 体量如升阻力系数等对应的伴随变量。数值实践表明，该基于该方法的网格自适应结果能显著改善整体量的精准 度，提高计算的可信度。 关键词网格自适应：伴随方法：误差分析：传感器 中图分类号0241 文献标识码A 1．引言 根据误差准则，网格自适应在局部进行细化或粗化网格。人们针对各种网格类型开发了大量 一阶1--阶导数，这些简单的传感器有各种问题。第一，选择的传感器对流场或计算量可能不合 适，例如，速度差分传感器可能在边界层梯度值较大而值较大，即使分辨率已经足够大，但是在 更感兴趣的涡分离结构处却不能很好的指示出来。第二，这些传感器没有权值，或者只是限于粗 燥的依赖于几何，导致可能在不必要的地方细化网格，浪费网格点。例如，从翼型上脱落的剪切 层在尾流下游会发生对流现象，但是在1．2个弦长的尾流下游，即使很差的剪切层分辨率也不会 对对翼型周围的流动产生很大的影响而不必要具有很高的分辨率。第三，对不连续问题，这些传 感器可能不收敛，例如激波随着细化不断的变得陡峭，但是这些地方不一定是最需要细化的。这 将导致不正确的网格收敛解。 近年来，基于伴随方法的后验误差估计研究取得了很大进展。伴随方法基于这样的事实一 大量计算工程最关心的是整体量如阻力／升力的误差。每个网格点上的伴随解可以看作是单位源 项对关心的整体量的线性影响。在误差分析中，源项可以在一定程度上看作是数值解与原始偏微 分方程的逼近程度的度量即残差，因此，伴随解给出了计算关心量残差的权影响。在伴随变量值 很小的地方，即使残差很大，其对关心的整体量影响也是很小的，因此在这里是不需要进行网格 加密的。另一方面，伴随变量值很大的地方，即使残差很小，其对关心的整体量影响也很显著， 因此在这里需要进行局部网格加密。 对给定函数的误差，自适应的策略就是细化那些对误差贡献最显著的网格单元。文中首先给 出了该方法在3D非结构网格Euler计算中的基础理论；然后给出了三角翼模型的计算情况，数 据实验表明，该方法是一种非常有效提高计算精准的方法。 2．误差分析与自适应原理 设U为非线性微分方程的解 Ⅳ(“)=0 “∈Q 边界条件为 D(“)=0 材∈aQ 线性微分算子L。，B。定义为 啦=蛳坐掣 -62- 第六届两北地区计算物理学术会议论文集 2008 D(u+藏)一D(掰) B。fi-=lim “ 占—0 非线性函数，(甜)定义为， J(u+蔚)一，(“) lim =(g(甜)，历)o+(办，C。万)∞ 口—’0 对应的伴随问题定义为： C’，=g(“) ’，∈Q 废v=h v∈孢 伴随等价定义为： (tv，万)+(B：v，C。万)m=(v，L。历)+(c|：v，B。万)鼬 上式对所有的1，，历成立。 设甜矗，V||1分别为离散主控偏微分方程及其边界条件与其离散伴随偏微分方程及其边界条件 的渐进数值解，假设足够可微，定义g^，hh为： 己17h=gA，瓦vh=魂 一B似 ) lI B ‰ lu+O(u^叫) lrjO C(U,lih)= ，k． I扯+口(Ⅳ．一H) O l g(蛐^)2k(“+口(甜。一材))dp， O 毗脚