MHD的一些困难问题的数值求解.pdfVIP

MHD的一些困难问题的数值求解 牟宗泽 龙永兴彭点云黄林 王中天 (核工业西南物理研究院 成都 61004约 用精确局部基构造出的新差分格式求解一些困难的等离子体物理问题，避 免了传统方法引起解的非物理振荡和解的畸变。理论和数值实验表明新方 法能解决一些用现有方法无法满意计算的数值问题 . 引 言 等离子体丰富的物理现象提出了一些至今求解仍然很困难的数值问题.如奇扰动(刚性)问题、 高振荡间匙、反向点及多尺度问题等。求解这些问题时，如用色散型离散格式，其数值解出现非物 理振荡:如数值格式呈耗散型、则数值耗散淹没物理耗散，其解失去真实性。这些问的题研究一直 为许多学者关汪乞含刊.相继发展了一些迎风格式、高阶格式、Ak调格式(FCT.TVD),它们总体上更 服了一些非物理振荡 但在解的高梯度区仍然出现非物理提荡。高阶格式需要更多的计值点，边界 的低阶处理也会影响内部精度。另外函数限制器会影响到问题收敛到稳态解川。对高振荡问题，至 今尚无有效方法求解，大多数方法因是求其某种平均量及其光滑部分，而放弃难于求解的点的函数 值 传统方法质疑 为行么通常传统方法不能满意地求解上述儿类有主富物理背景的典型问题?传统方法的基础是 Favlor展开式前面几项，红 u._1一u_, h一生二h0 (1) 2h up6s,120 U,为 亡扮 式可近似表式成 u = ，L二u._,十O(h2) (2) 2h 鱿，二u，一2u+ ‘，__-u-=th厂=_-u=seh* (3) h} 2 360 (3)式近似表为 u一2u+u, 十0(h-) (4) 力2 9它各种迎风格式、高阶格式也都是建立在这个荃础上的.这就要求解EE数是连续、光滑外 或者解函数是缓慢变化的 或者离散网格是极其精细的。总之要保证解函数相对于网格而言是缓变 的。这就涉及到通常所说的解的稳定性和精确性问题。色散型格式引起解的非物理振荡:耗散型格 式引起解数值的失夏。对于对流扩散微分方程表现为要求局部 Reynolds数小于临界值。当物理解 有极A的边界层效应或高振荡现象时 ’这个条件是相当苛刻的。如当解函数局部性质具有实的或虚 的指数幂次时。它的iavlor展开式为 万乃2 刃h eu二I+肋 + + — 十 (5) 2 6 fih, 才h* 10儿6 Cos(肠)=1- 岛 十 — 一 — + (6) ‘盆 24 720 — 298 一 当司;hiRi^h.寸物理上li边界层效应或有高孤荡角和勺清形)日手豁然 ·5少·之,)成仅保留已 的前儿项不是 一个好的近似式、忽略掉的黄的fftEl能比保留的量的值还要大 这5.不能认为忽略掉 h的n次is以后的贫，这个公式的近似是O(h)的 苦则解的特性或数值v失a 正确的计算应保 证解的行为特性的正确和数值的精确