一种基于动态因子驱动递进式多场耦合算法.pdfVIP

一种基于动态因子驱动的递进式多场耦合算法∗ 何再兴，张树有，伊国栋 浙江大学机械与能源工程学院，杭州(310027) E-mail ：2251hezaix@163.com 摘 要：以热——结构耦合为模型，提出了一种基于动态因子驱动的递进式多场耦合算法。 该算法首先根据实际物理量变化的大小将模型初步分裂成显式积分区域和隐式积分区域，对 不同区域分别进行物理场分割计算，然后检验显式积分区域边界场量的变化大小是否超过临 界值，对超过临界值的显式积分区域，根据其移动增量控制动态因子对该区域其进行二次分 裂，并对划分后的隐式积分子区域进行进一步的积分运算，从而得到符合精度要求的运算结 果。本文的算法在低压断路器过载保护元件——双金属片在通电流情况下的分析中得到实际 应用验证。 关键词：多场耦合，分裂，分割，动态因子，二次划分 中图分类号：TH702 1. 引言 由于工程实际中存在大量的多物理场的耦合作用，多场耦合(multiphysics problem) 已经 成为CAE 的研究热点之一[1,2] 。耦合场仿真使用的数值方法大多数以牛顿法和松弛法为基础， 多场耦合算法在构造上主要有两类方法：区域分裂法与物理场分割算法。 区域分裂法主要是根据物理模型不同区域的场变化大小，分别采用隐式和显式积分，以 达到在保证计算精度的前提下减小计算量的目的。直接积分方法是多场耦合常用的数值方 法，该方法主要有显式积分法和隐式积分法两种。隐式方法的优点在于数值稳定性和能量耗 散性能较好，精度高；显式积分法最大的优点是计算效率高，但是它是一种条件稳定的算法， 精度低。为了把隐式算法的数值稳定性和高精度特性与显式算法的计算效率结合起来，许多 学者在这方面做了大量的研究。Hughes[3,4]基于单元分离法提出了一种显隐式混合算法， Chiang[5]通过设计方程组并行求解器来实现显隐式混合并行计算的功能。 物理场分割算法( partitioned solution)主要是将不同物理场分割开来，在一个计算步中只 含有一个物理场量。用预测值替代某一场变量的真实值,时间积分在仅考虑一个场变量的子 系统中进行,而相互作用项作为外力来考虑,从而实现不同物理场的解耦。分割算法在不同物 理场的求解分开进行,降低了计算复杂度，提高了计算效率。主要的算法有Huang[6],[7]提出的 无条件稳定的交叉迭代法，以及多次迭代方法[8]等分割算法。 分裂算法和分割算法分别通过不同的方式来减小计算量，分裂算法通过显隐式积分的有 效混合来提高运算效率，但是求解耦合场方程组的运算量仍然很大；分割算法若是全部采用 隐式示积分，则运算精度虽然可以保证，但是计算效率会比较低。随着工程实际中越来越复 杂的耦合场问题的出现，迫切需要一种高效与高精度的多场耦合算法。鉴于此，本文提出了 一种综合区域分裂与物理场分割算法优点的多场耦合算法，并通过实例进行了验证。 2. 热—结构耦合模型 在一般三维问题中，瞬态温度的场变量φ(x , y , z ,t)在直角坐标中应满足的微分方程是 （边界条件应该满足：Γ+Γ +Γ =Γ，Γ是Ω域的全部边界）： 1 2 3 ∗ 本课题由国家863 项目（2006AA04Z114 ）和高等学校博士学科点专项科研基金（20040335060 ）资助。 - 1 - ∂φ ∂ ∂φ ∂⎛ ⎞ ⎛ ∂φ ∂ ∂φ⎞ ⎛ ⎞ ρc − k − k − k −ρQ 0 (1) ⎜ x ⎟ ⎜ x ⎟ ⎜ x ⎟ t t x y y z z ∂ ∂ ⎝ ∂ ⎠ ∂ ⎝ ∂ ⎠ ∂ ⎝ ∂