二维三温流体力学计算中时间步长的自动控制.pdfVIP

。。。。。 计FI 算畀 物例 理理 第 卷第6期 v01．22．N。．6 OF PHYSICS 2005年11月 CHINESEJOURNALCOMPUTATIONAL No，．j2005 【文章编号】1001．246x12005)06-0515—05 二维三温流体力学计算中时间步长的自动控制 符尚武， 戴自换， 邬吉明 (北京应用物理与计算数学研究所计算物理实验室，北京 100088) [摘要】 研究了二维三温流体力学计算中时间步长的控制问题，提出了控制时间步长的多种约束条件：既考 还考虑了Lagrange网格密度(或体积)的相对变化以及三温能量方程的迭代次数等．在计算过程中这些约束条件 可以随时自动地改变时间步长，以最经济和合理的时间步长完成计算．最后给出了数值实验结果． [关键词】 二维三温；流体力学计算；时间步长控制 [中图分类号】 0241 [文献标识码】 A O 引言 在现代科学技术诸多领域的研究中都要涉及二维流体力学，特别是要涉及带有热传导的多介质可压缩 二维流体力学．带有热传导的多介质可压缩二维流体的力学运动是很复杂的运动，它的数值模拟具有相当大 的难度．目前，求解多介质二维流体力学方程组多采用Lagrange方法，它的坐标建立在流体质团上，跟随流体 运动，描述流体运动的计算网格的大小形状随时随地都在改变，而计算过程所使用的时间步长不可能一成不 变，否则既不经济也不可能得到满意的计算结果．时间步长的选取关系到数值计算的稳定性和精确度，如何 控制计算过程中的时间步长，以最经济和合理的时间步长完成计算是二维流体力学计算中值得探讨的一个 重要问题．目前国内外公开发表的关于流体力学，特别是关于辐射流体力学计算中的时间步长控制的专题论 AlamosNational 文很少，1999年美国Los 的论文n]，该文仅涉及到单一的辐射扩散模型问题．本文研究了多介质可压缩二维三温流体力学方程组计算 中时间步长的控制问题，给出了控制时间步长的多种约束条件，并给出了数值实验结果． 1微分方程和定解条件 在平面直角或柱对称Lagrange坐标系中可压缩二维三温流体力学方程组为 訾+lDV·￡，=o， (1) 訾+VP_0’ (2) Fi：叫d(re1)， (3) c¨孥+(Pi+”3a6川i／dV。+m ch孥+(髓等)訾+y V幔一。一耽 初始条件： a=e，i，r， u(菇，Y，to)=0，p(x，Y，to)=po(菇，y)，瓦(z，Y，to)=贮(菇，Y)， 边值条件： [收稿日期】2004—07一Ol；[修回日期]2004—10—25 [基金项目】国家自然科学基金、CAEP科学技术基金和国防科技重点实验室基金资助项目 [作者简介]符尚武(1941一)，男，陕西兴平，研究员，从事计算流体力学方面的研究． 万方数据万方数据 516 计 算 物 理 第22卷 固壁和轴(戈轴)对称条件：U。=0，(F。)。=0； 自由边界条件：(F。)。=(F；)。=0，Tr=T，(省，Y，t)，只=0； 联接条件：在不同介质交界面处，u。，F。，P。连续． e。，e；，e。分别了电子、离子、光