求解热传导正问题及反问题的数值方法研究

来源 :太原理工大学 | 被引量 : 4次 | 上传用户:song52535
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
传热学是一门与人类生活息息相关的学科,需要有效的方法来探索和解决实际工程及现实生活中的热传导问题。由于实际问题的复杂性,利用各种数值方法来求解实际的热传导问题成为传热学研究的重点之一。本文主要基于网格方法(有限元)、无网格方法(RBF)以及有背景网格的“类”无网格法(光滑有限元)来构建高效、稳定的数值模型,解决以下复杂而重要的热传导问题:(1)问题域边界随时间变化的热传导问题;(2)逆向边值定解问题;(3)物体边界重构的热传导问题;(4)Cauchy不适定热传导反问题。首先,作者对大型设备内部边界随时间变化的热传导问题进行研究,如流动的钢水对炼钢炉炉壁侵蚀,使得炉壁内部边界随时间变化。在标准FEM解决非稳态热传导问题时,时间项采用常规的有限差分法,所以在时间推进的过程中需要上一时刻的节点温度。但是因为边界的变化导致不同时刻网格的变化,从而使得部分节点温度缺失。因此,本文提出了两种不同的FEM模型来计算有移动边界的非稳态热传导问题。第一种模型是通过扩展的插值模型来计算上一时刻缺失的节点温度,然后再利用标准FEM算法计算当前时刻的温度分布。该模型可以获得较精确的解,但是因为需要在单元中进行额外的插值,运算时间较长。在第二种模型中,考虑到边界随时间变化,使得节点温度不仅与时间相关,而且与边界变化率相关,因此基于全微分公式,通过数学推导获得一个修正项公式。修正项消除了因为直接采用标准FEM模型求解移动边界非稳态热传导问题带来的误差,从而作者建立了一个统一的处理边界固定和边界移动问题的FEM算法。并通过大量数值实例,说明第二种模型更精确,而且运算速度也比第一种模型提高了10%左右。作者研究的第二个比较复杂的热传导问题是具有不适定性的Cauchy反问题。首先从标准FEM系统方程中分块抽取出与Cauchy边界热流相关的局部节点温度方程组,然后建立Cauchy边界上已知热流与局部节点温度方程组的关系,最终建立未知热流方程组,将该问题转换成一个正问题求解。为了消除反问题的不适定性,获得稳定、可靠的数值结果,作者采用SVD方法来消除“不稳定模态”。即依靠截断一部分小特征值来防止输入误差的扩大。为了能够自动选择截断特征值的个数,本文提出了一个自动选择算法:通过选择近似热流值和给定热流差值的最小值来确定最合适的特征值个数。接下来,作者探索利用径向基函数建立一个独特的时空统一模型来求解混合介质墙的非稳态Cauchy反问题。并深入研究混合介质墙内侧的边界重构问题。实际工业生产中,通过传感器获取炼钢炉外表面的温度和热流,本文通过将时间和空间统一处理,建立时空统一的RBF反问题模型来求解混合介质内侧的温度。在时空统一模型中,设置“影响因子”来调整时间和空间对于“距离”的影响。此外,该模型没有采用常用的逐层求解的办法,而是形成多层统一的整体方程组,避免了逐层求解中带来的误差积累。为了获得更加精确的解,本文提出了一种自动算法来选择每层的形状参数c:最合适的c可以使得近似温度值和给定温度值的误差最小。为消除反问题的不适定性,本章利用Tikhonov正则化方法来保证RBF统一时空模型能获得稳定、精确的解。对于边界重构问题,需要在上述时空统一的RBF反问题模型基础上设置一个虚拟边界,通过查找满足内部边界温度要求的节点位置来模拟动态变化的内部边界。在本论文研究的第四个复杂问题中,作者建立一个通用、高效的S-FEM求解器来解决2D和3D的稳态及非稳态的热传导问题。S-FEM是近年来G.R.Liu提出的一种新的数值方法,具有许多FEM不具备的独特性质。本文提出的高效、精确的算法能实现不同的S-FEM模型,与当前研究人员普遍采用的面积/体积平均的方法不同,本文实现的光滑域线段/光滑域子面积分的方法更通用、效率更高。因为该方法不仅可以应用到线性单元中,而且可以直接运用到高阶元,因此本文提出的求解器更通用、性能更高。这也是第一篇严格按照G.R.Liu弱弱理论实现S-FEM方法的文章。此外,作者还提出了建立不同2D和3D光滑域的高性能算法,并在建立光滑域的过程中获得所有关系矩阵以及计算光滑温度梯度矩阵所需要的数据。在这个过程中,一个非常简洁的算法被提出来计算高斯点处单位法向量的方向。为了提高求解器的运算速度,所有数据都保存到数据库中,减少查询和读取的时间。在后处理部分,利用光滑域光滑温度梯度重构了背景网格节点以及单元的温度梯度。在通用S-FEM稳态热问题模型的基础上,采用不同的差分格式建立S-FEM非稳态热问题求解算法,并通过大量数值实验,讨论不同S-FEM模型在不同差分格式中表现出的数值特性。在作者研究的第五个热传导问题中,S-FEM方法被用来求解不适定的Cauchy热传导反问题。利用矩阵分块技术,从S-FEM系统方程组中抽取和未知热流相关的局部节点温度方程组,结合Cauchy边界热流和局部节点温度方程组关系,建立基于光滑域的热流方程组。虽然和第二个问题的模型相似,但是S-FEM反问题模型基于光滑域,可以获得更精确的解。为了比较不同正则化方法,作者分别采用SVD和Tikhonov方法用来消除热流方程组系数矩阵的奇异性。同时,本文还通过数值实例结果,讨论了ES-FEM和NS-FEM在求解Cauchy热传导反问题时表现出的不同数值特性。
其他文献
采用计算机三维造型得到汽车仪表板的实体模型,利用MoldFlow软件对模型进行了有限元分析,完成了浇注系统、冷却系统及成形工艺的设定。在对熔体充模、保压、冷却进行模拟的过
随着全球经济的快速发展,电网大规模互联成为电力系统发展的必然趋势,电网的复杂程度也随之增加。在提高系统运行经济性的同时,对电网的供电可靠性提出了更为苛刻的要求。局
《拍案惊奇》和《二刻拍案惊奇》是由明代作家凌濛初创作而成。由于小说的受众大部分都是文化程度较低的下层人士,其创作时大致采用了通俗易懂的口语,因此“二拍”基本上反映了
我们知道,一个国家经济的增长,总是来自于劳动分工与专业化引起的劳动生产率的提高以及市场的深化和扩展。换句话说人类社会的经济增长说到底和市场的深化和扩展是一回事。从
背景和目的:钙化性主动脉瓣疾病是全世界范围内发病率最高的心脏瓣膜疾病之一。而随着老龄化社会的来临,这一数字仍将不断增加,造成极大的社会医疗负担。目前对于该病,除了手
为分析碳纤维木质复合材料(CFRW)的电学性能,采用木质纤维作为基体材料,加入碳纤维得到具有导电性的新型复合材料。通过对CFRW表面电阻率测试实验数据分析可知,碳纤维赋予CFR
以人为本是当代社会的本质和核心,也是中国成人教育发展的动力和方向。文章在分析了人本化成人教育内涵和我国成人高等教育价值取向偏差之后.对新时期成人教育如何体现人本价值
背景:在经皮椎体成形治疗骨质疏松椎体压缩性骨折时,较高的骨水泥渗漏率限制了其应用,如何避免或减少骨水泥渗漏成为目前本领域研究热点。目的:比较应用高粘度及低粘度骨水泥
目的通过生物力学测试比较上下平行2枚和倒三角形3枚A0空心加压螺钉固定Pauwels Ⅱ型股骨颈骨折的强度,为股骨颈骨折选择更合适的临床治疗方法提供参考。方法 将12根入选的股
目的评估微信公众平台在北京市区级结核病防治机构健康教育工作中的应用现状。方法对北京市区级结核病防治机构进行微信公众平台应用基本情况问卷调查,同时检索认证公众号,对