有限时间奇性问题在JASMIN上开发研究.pdfVIP

郭红等：有限时间奇竹n日题在JAsMT～卜的开发研究 I 有限时间奇性问题在JASMN上的开发研究 郭 红’，李 若2 100094；2．北京大学数学科学学院，北京100871) (1-北京应用物理与计算数学研究所，北京 【摘要】 三维不可压N—s方程是否会产生出有限时问奇性是流体力学中最具挑战性的问题之一．在数值计算 中，随着网格不断加密，会因对计算资源需求过大导致无法深入地开展研究工作。我们将建立并行自适应算法，解 决计算需求过大和资源不足的矛盾。并行自适应算法的研究要求同时具备多方面的知识与技能，开发相应的程序具 有一定难度。我们将在JASMIN框架的支持下，研制可有效使用现代高性能计算机的并行自适应计算应用程序，来帮 助深入研究有限时间奇性问题。 【关键词】 有限时间奇性：并行计算；自适应网格加密；JASMIN 0 引 言 开问题之一，也是七个克莱千禧问题中唯一一个和应用数学有关的问题。关于这个问题的理论分析 长期以来处于困顿状态，在数值计算方面，有两个算例尤为受到关注，包括反对称涡管结构和 Kida-Pelz高对称初值。在2005年，侯和李的计算衷明，反对称涡管结构不会发展出有限时间奇性 …，这个结果受到广泛关注。 的涡量增长，可能产生出有限时间奇性；1997年，Pelz基于其1994年的计算结果，提出可能发展 出奇性的解的涡线模型，并使用数值计算表明，这个涡线模型确实具有稳定的奇性怕1；2006年，侯 和李对1994年的计算结果进行重复，发现计算过程中的舍入误差是不会破坏解的结构的，从而为大 规模并行计算对此问题的模拟提供了可能：进一步的，侯和李构造了具有周期结构的涡线模型，发 展了相应的数值算法，对周期结构的涡线模型进行模拟，完全重复了周期涡线模型的稳定的奇性现 象。 上述工作均是采用拟谱方法实现的。拟谱方法具有精度高等优点，因此一直以来成为计算有限 时间奇性问题的酋选方法。然而拟谱方法无法通过局部加密网格来捕捉物理敏感区域，随着网格的 不断加密，对计算资源要求过大导致无法深入地开展研究工作。如何解决这一问题，能对Kida-Pelz 高对称初值进行比较深入的研究，成为有限时间奇性研究领域的迫切任务。 相对于串行计算，并行计算可以充分利用当前的高性能计算环境，扩大问题规模，提高求解精 度。相对于单层网格计算，多层自适应网格计算Ⅲ采用细化局部网格区域的方法，在相同数值精度 的情况下，可以减少网格单元，节省计算资源。并行自适应计算融合了这两种方法，已经成为高性 能、大规模科学与工程数值模拟的最要手段。但是并行自适应算法的研究要求同时具备物理、计算 数学、并行计算和自适应方法等多方面的知识。同时，并行自适应算法的实现涉及复杂的数据结构 和并行程序设计及优化等。这些对于从事物理研究、计算方法研究或者岛性能计算研究的人都是一 个巨大的挑战。 由北京应用物理与计算数学研究所研制的JASMIN框架哺1，针对科学计算中的结构网格应用，通 过封装高性能的数据结构、集成成熟的数值算法、屏蔽大规模并行和网格自适应 1郭红(1978～)、女、汉族、河北深泽、助理研究员、理学博士、主要研究方向：计算流体力学、并行计算，通讯地址： 北京8009信箱17分箱，100088 270 第!。口q隔中囝窄气动力学物理气休动力学学术交漓会论文集 III 的计算技术，支撑物理建模、数值方法、离性能算法的创新研究，加速研制可有效使用现代高性能 计算机的并行自适应计算应用程序。在JASMIN框架的支持下，用户只需提供求解偏微分方程组的离 散格式、初边值条件和误差评估方法，就可以获得并行自适应计算应用程序。 本工作采用有限差分离散方法结合h一自适应结构网格加密技术，应用JASMIN框架对不可压三 维N—S方程进行并行求解，希望能对Kida—Pelz高对称初值进行比较深入的研究。 1 控制方程 考虑涡量一流函数型的三维不可压N-S方程 一△痧=历 (1) 厅=V×痧