计算天文数值计算lec2解域、网格的划分与差分的构成lec2分析.pptVIP

2. 解域的离散 3. 网格的分类 不均匀的原则 自适应网格 4. 规则的欧拉网格（一维） 4. 规则的欧拉网格（二维） 5. 网格的生成 6. 差分处理概貌 7. 空间导数的差分常用形式 8. 空间导数的差分一般方法 * * /~chenpf/c/courses/compastro/lec1.ppt 1可逐渐用2－15代替 http://www.cfca.nao.ac.jp/hpc/misc/EANAM2004 超新星中子星黑洞/ 恒星行星形成/ 恒星大气太阳系/ 吸积盘喷流/ 星际介质物理/ 宇宙学星系/ 解域、网格的划分与差分的构成 1.解域 有限的空间区域 （突出问题的重点） 选取原则 物理边界：根据观测或 物理原理、定律 几何边界：周期性或对称性 推荐使用 计算边界：无穷域（如耀斑 模拟的上、左右边界） 数值求解的空间区域 边界 将解域分为许多小格子，即网格 grid或grid mesh 网格的交点称为格点或节点 grid point 格子可为三角形、四边形甚至六边形 计算机只能处理有限的信息，故须对连续区域进行有限分割 欧拉网格 （空间上固定） 拉格朗日网格（随流体一起运动） 混合型网格（任意运动） 规则网格：网格直接由坐标线组成 不规则网格： 可不均匀 原则1：保证足够多的格点分辨物理结构（如激波、陡变区、电流片等） 如：一般方法中至少5个格点分辨激波，10个左右格点分辨电流片 原则2：尽量缓慢变化 如：采用等比递增，等比系数不大 原理： 开始一个粗略的网格 估计局地截断误差 确定需细化的区域 叠加上更精细网格 继续确定下一级细化网格直至满足给定的截断误差 对任一格点，时间步长和空间步长同比缩小 ***: 其中仅初始粗网格的间距是固定的 /users/dagh/Tutorial/jcb_mitra/tut_jcb_mitra.html 针对局部陡变区域 AMR: Adaptive Mesh Refinement 星风 耀斑环内动力学 自然坐标系 太阳风的纬度分布 耀斑中磁重联过程 直角 或柱坐标中r-z 极坐标 喷流 需要专门的网格生成技术和工具 近年发展了无网格mesh free模拟 以一维对流扩散方程为例 初条： 边条： 令时间步长为 则 格点分布 x tn tn+1 以差分代替微分 代入方程 f x i-1 i i+1 x tn tn+1 微分方程的差分形式 微分方程的差分形式 当步长 ? 和间距 ? ?0，原则上可获得精确解 时空步长须满足一定的约束条件才使得差分方程与微分方程相容，从而保证计算稳定； 小心源于中心差分的棋盘误差 chess error ； 将来可知上述差分处理不理想；设计理想的（如可靠性、精确度）差分格式是有限差分的核心。 注意 以二维均匀网格为例 一阶偏微分 前向差分（forward difference） 后向差分（backward difference） 中心差分（central difference） 三点加权平均差分 三点中心差分 ………… 由此可见，形式不唯一 理想（磁）流体方程组为一阶，考虑耗散等项后为二阶 二阶偏微分 混合偏微分 以二维非均匀网格为例 其中 1. 待定系数法 试构造 的差分表达式，要求截断误差为 且只应用三个网格点 （均匀间距） 令 对 在 附近按Taylor级数展开 令 中 及 的系数为0， 而 的系数为1 代入 即 作业 试构造 的差分表达式，要求截断误差为 且只应用三个网格点 其中 与 不相等， 同样利用待定系数法求得 其中