球面上偏微分方程的数值求解的研究.pdf

球面上偏微分方程的数值求解研究 摘 要 本文将立方球体技术引入球面偏微分方程的数值求解之中，该方法既避免了球面 坐标上数值求解偏微分方程的难点，又在数值计算的精度和效能方面有很大的潜力，尤 其是在并行数值计算方面有谱方法不可比拟的优势。 首先本文给出了浅水波方程的六种表达形式，及其一套测试数值求解算例，然后给 出了cube．sphere坐标的构造方法，与球坐标、直角坐标系的转换规则，和该坐标系内 相应微分算子的表达式。在此基础上在立方球体坐标系下推导出浅水波方程的形式，使 达式。为了解决由于边界插值所引起的数值不稳定性，引入了四阶的人工粘性项(一v4l。 通过对测试方案的求解和分析，可以看出该方法有如下的优点：(1)完全显式格式(2) 时间上的四阶精度(3)仅仅通过改变一个参数就可以获得空间上的高阶精度(4)易于 在并行的计算机上有效的执行，而且该数值方案有较强的健壮性。 本文是国内首次应用cube—sphere技术数值求解球面偏微分方程。本文的工作为 数值天气预报的基础工作向前推进了一步。 关键词：浅水波方程、立方球体技术、龙格库塔方法、人工粘性、MDQ方法、数值天气 预报 PartialDifferentialof in Sphere Numerically Equation Study XuYabo Jianzhong Supervisor：Prof．Chang ABSTRACT Inthis methodis tosolve differentialof thesis，cube·sphereadopted partial equation avoidsthe differentialin equation sphere．This difficultyofnumericallysolvingpartial sphere coordinationandhas latent innumericalcalculation’sand great capacity accuracypotency． ithas over methodin numericalcalculation． great spectral parallel Especially advantage Firstthesixformsof ofshallowwater andasuiteofnumerical formulation equations testmethodsare conformationmethod coordinationandtransform given．Then ofcube-sphere between and and rules coordination coordination sphere fight-angle corresponding are inthe t