间断陡削化方程Ⅰ人为刚性源项方法刚性常微分方程的计算方法….PDF

间断陡削化方程Ⅰ:人为刚性源项方法:刚性常微分方程的计算方法 … 第六图书馆 本文针对界面运动的计算提出了适当附加人为刚性源项计算接触间断的方法,数值实验表明这一方法是有效的。本文针对界面运动的 计算提出了适当附加人为刚性源项计算接触间断的方法,数值实验表明这一方法是有效的。人工刚性源项 间断陡削化 刚性 常 微分方程空气动力学学报邬华谟 孙耿中国科学计算数学与科学工程计算所1998第六图书馆 Vd．16．No．1 Mar．，1998 f、、， 厂 卷年 间断陡削化方法 I：人为刚性源项方法 一 刚性常微分方程的计算方法的应用★ 盟 中国科学院计算数学与科学工程计算所 ·北京 100080 2 孙 耿 王厘尔 中国孚iii再所，北京100080 02 · 摘要 本文针对界面运动的计算提出了适当附加人为刚性擐项计算接触间断的方法。数值 实验表明了这一方法是有效的。 田龃姓 关键词 人工刚性振项} 中图分类号 O241．3，O351．2 0 引 言 雉 揪 间断在 自然界是普遍存在的现象，计算流体力学中常常要赴理分隔不同流动状态的澈渡 问断，切向间断或分隔不同介质的界面接触阃断。这类间断的精确计算是非常困难的，同时对 于正确模拟力学过程又是十分重要的。自从 50年代NeuIlⅧ 和Richt~yer在文献[1]中引人 人为粘性项计算激波 以来，各类阃断的计算方法研究一直相当活跃，成功地发展了众多的计算 同断解的方法与理论，促进了计算流体力学这门新学科的诞生与发展。虽然计算同断特别是激 弓 渡强同断的计算方法 已日趋成熟，但人们对改进接触同断的计算方法的研究仍很活跃。计算实 践与理论分析表明线性方程的强间断(接触间断)的计算 比拟线性方程强间断 (澈波)解的计算 复杂，因此提高计算线性方程接触间断的分辨率是很有意义的。 计算强间断问题的有效算法很多。但总的来说可以分为阃断跟踪算法及捕捉算法两大类。 问断跟踪法有明显跟踪与隐含跟踪两种。属于明显跟踪方法的有特征线方法，奇性分离方法嘲 是一种隐式特征线方法 在胯含跟踪方法里阃断不明显出现，但有关阅断性质的量，如阋断位 置和传播速度等都记录在相应阿格区同里 文献[3]的方法就属于这一类。界于明显与隐含跟 踪之间的有广义 Lasrange方法H]。在这里同断传播线已作为坐标线 ，从而可以很精确地处理接 触问断。在间断捕捉算法中，文献[1]的人为粘性方法是这一方面的先驱．此后这一思想得到了 很大的发展。在方程中显式附加物理牯性璜的思想发展为在算法中附加格式牯性。使得捕捉 ★国隶 自然科学蓦盒和国家摹登资助项 目 奉丈于1997年3月l0日收割，9月27日收到修改稿。 第 1期 邬华谟、孙 耿等： 间断陡削化方法 1：人为刚性源项方法