非匹配网格间数据传递

分类:耦合分析·多物理场 | 2026-10-01 改写
CAE visualization for mapping dissimilar meshes theory - technical simulation diagram
非整合网格间数据转移

理论:传递数值还是传递总量

概述

🙋

把流体分析的热流传给结构分析后,温度对不上了。是网格不同的原因吗?

🎓

可能性很大。流体侧为分辨边界层通常壁面附近网格很细,结构侧较粗,所以边界面上两套网格不对齐。热流怎样传递会改变结果。在接收方点的位置读取值的方法(一致型)很简单,但热流分布中有尖峰时会漏掉峰值,传递的总热量就变了。热是应该守恒的量,总量一变,结构整体温度就偏了。按重叠面积平均的守恒型能保持总量,但粗网格会把峰抹平。用哪种要按传递量的性质决定。

一致型与守恒型映射

$$ \text{consistent: } q_j^{T} = \sum_i N_i(x_j)\,q_i^{S},\qquad \text{conservative: } q_j^{T} = \frac{1}{|\Gamma_j|}\sum_i |\Gamma_j \cap \Gamma_i|\,q_i^{S} $$

S为发送方(流体),T为接收方(结构),N为发送方形函数(插值),|Γⱼ ∩ Γᵢ|为两单元重叠面积(一维为长度)。

总量守恒

$$ P = \int_\Gamma q\,d\Gamma,\qquad \sum_j q_j^T |\Gamma_j| = \sum_i q_i^S |\Gamma_i|\ \ (\text{conservative}) $$

算例设定:长1 m的边界。发送方97单元,热流q(x) = 1000·(1 + 0.8 sin 2πx + 3 exp(−((x − 0.3)/0.02)²)) W/m²(x = 0.3处有尖峰)。接收方为10单元和40单元。

Coffee Break 随笔闲谈

连接不对齐的网格

1990年代耦合分析开始普及时,让边界网格对齐很费工夫,很多时候几乎不可能。于是研究了在不对齐网格间传值的映射方法,气动弹性领域提出了让力与位移映射组合以保持虚功的思路。现在的耦合库(preCICE等)可以按量选择一致型或守恒型。

算例

例1:粗接收方(10单元)

方法传递总热量总量误差最大热流
最近点值1,003.5 W/m−9.30%1,813 W/m²
线性插值(一致型)1,002.3 W/m−9.41%1,810 W/m²
按重叠平均(守恒型)1,106.4 W/m0.00%2,317 W/m²

例2:较细接收方(40单元)

方法总量误差最大热流
最近点值−1.52%3,518 W/m²
线性插值(一致型)−0.16%3,753 W/m²
按重叠平均(守恒型)0.00%3,699 W/m²

发送方总量为1,106.35 W/m,最大单元值为4,571 W/m²。

🙋

10单元时,一致型会丢掉近一成的热量。

🎓

如果接收方单元中心没有落在峰上,按点读取的方法就会漏掉峰值。峰宽(约0.04 m)比接收方单元长度(0.1 m)窄,峰的热量几乎传不过去。热流、力这类合计有意义的量原则上用守恒型传递,温度、位移这类场值用一致型传递。流固耦合中,从结构到流体用一致型传位移、从流体到结构用守恒型传力,边界上的能量交换也是正确的。但守恒型下峰高也会在接收方网格上被抹平,要正确求局部最高温度,接收方网格也要细到能分辨峰宽。加密到40单元后,各方法的峰高都恢复了很多。

映射的选择与确认

  1. 区分传递的量是场值(温度、位移)还是合计有意义的量(热流、力)。
  2. 场值用一致型,热流和力用守恒型传递。
  3. 比较传递前后的总热量、合力和力矩。
  4. 确认两侧边界形状一致(无间隙或嵌入)。
  5. 局部值重要时,按分布的细致程度加密接收方网格。
Coffee Break 随笔闲谈

“被过度冷却的部件”

某部件热分析中,把流体分析求得的传热结果传到结构粗网格后,部件温度比实测低10℃以上。调查发现传过去的热流合计比流体侧少约8%,高温气体冲击的窄小范围的热量丢失了。改用守恒型映射并加密该范围的结构网格后,热量差消失,温度与实测相差在几℃以内。此后把比较映射前后热量合计作为标准确认步骤。

常见错误

错误与对策

错误后果对策
用一致型传热流热量增减改用守恒型
不确认总量察觉不到误差比较前后合计
接收方过粗峰被抹平局部加密
边界形状偏差数值被外推使形状一致
力与位移分别选择能量不守恒保持虚功的组合
🙋

我想了解相关内容。

关联模拟器

在本领域的交互式模拟器中体验理论

模拟器列表

相关分析领域

结构分析电磁分析热分析
本文评价
感谢您的反馈!
有帮助
更详细
报错
有帮助
0
更详细
0
报错
0