11. 初始条件和边界条件.pptVIP

11. 初始条件和边界条件 举例 初始条件 （一）定态解析解 等温静止平衡大气 磁静平衡大气 （二）半解析方法 （三）观测值作初条 （四）数值方法 11.2 边界条件 一、 物理边界条件 二、 几何边界条件 三、 计算边界条件 四、 奇异边界的处理 * * 它们也是数值方法成败的关键 纷繁复杂的天体现象满足几乎同一方程， 这很大程度上归因于定解条件：初边值 Chen Shibata (2000) 解释日冕物质抛射 Chen et al. (2002) 解释EIT波 非定常问题： 定态解作初始态，再引入物理扰动 如电阻、足点运动、压力脉冲等等 定 常 问 题： 迭代 迭代初值 势场 直角坐标系下的拱形场 柱坐标系下线偶极子场 球坐标系子午面偶极子场 势场 球坐标系赤道面内 开场 闭场 螺旋场 无力场 直角坐标系 柱坐标系下一维无力场 (1) 均匀引力场 (2) 万有引力场 (1) 一维电流片 (2) 位于 无限长 半圆柱面电流片 注意：解析解作初值只是解析上的平衡，而非数值上的平衡，有时会导致扰动的传播，有时则很快演化到数值平衡态。 方法一：在方程（如运动方程）中减去初值 方法二：修正部分系数 通常是部分观测量， 如太阳大气温度的分布。 选定某物理量分布后由差分方法求出其它量的分布，如 静力平衡 磁静平衡 对定常和非定常问题都涉及边界条件，它比初条更为复杂，为数模之弱点和难点。它不同于内点，通常不能用差分形式。 物理边界条件：根据观测或 物理原理、定律 几何边界条件：周期性或对称性 计算边界条件：无穷域（如耀斑 模拟的上、左右边界） 给定某量或其导数的值 Neumann第一、二、三边值问题 周期性、对称性(推荐使用) 1. 无穷域 有限边界：人为边界 2. 有时出现在物理边界上，因变量不能全部给定，只能部分给定，其余由相容关系给出 例：新浮磁流 例：对称边界 开拓网格 正常网格 有时会导致非物理结果 推荐使用 参见P. F. Chen, C. Fang, Y. Q. Hu, 2000, Chinese Science Bulletin, 45, 798 以密度为例 1. 过给定 0 1 2 3 稳定，但有强反射 2. 等值外推（0次外推） 精度及稳定性均好 3. 线性外推（1次外推） 可视为无反射条件 （某种数值外推） 4. 增量等值外推 单边差分 它等效于n时间层的线性外推 v0时宜采用，它具有迎风特征，v0不宜 经验表明(1)一般情况下外推公式精度越高，计算越不稳定；(2)单边差分和等值外推效果较好；(3)有重力时，可采用压力梯度力与重力平衡；(4)有些量不能等值外推，如MHD方程组中的磁通量函数（见Chen et al. 1999, ApJ, 513, 516）；(5)磁场的边条总是一个麻烦，可在边界采用势场近似。 1. 轴对称情况 作傅立叶展开 源于坐标系的非单值性。既不属于计算边界，又不能标定物理边界，如柱坐标中的r=0或球坐标的θ=0 ，但?可取0－2? 采用开拓网格 2. 非轴对称情况（采用非开拓网格） 即极轴，与?无关，故 （1）. 对p, ?,Br,vr （2）. 对B?、 B?、v?、v? 但较复杂，提倡极轴上的点的值由周围点的值取平均，只是注意对v?、v?须将之化为直角坐标系下分量后在平均 代入原方程组，得出Fm的方程组 以球坐标为例 (1) 等截面管道中一元不定常等熵流动 五、 特征线方法 连续性方程 （1） 运动方程 （2） 能量方程 （3） 沿着C+: 沿着C－: 沿着C0: 黎曼不变量 (5.2) 一维不定常线性奇次一阶偏微分方程组的特征形式 其特征方程为 特征线 若矩阵A的m个本征值全为实数，则原方程为双曲型。设A的m个本征值为， 相应的本征矢量为 以 分别左乘原方程 原方程的特征形式 对一维不定常线性奇次方程 为不变量，对非线性情况，变化缓慢 (5.3) 近特征线和近特征形式 多维MHD问题特征几何及方程的特征形式十分复杂，故往往计算量太大以致于无法进行，为此人们引入近特征线或投影特征线的概念，它是特征线在某一时空坐标平面内的投影。一般来说，只要所选时空坐标平面基本与流动方向和磁场方向一致，近特征线与特征线相差不远。 其中已选适当单位消去?0 其中v为二重根 阿尔芬波速 快波速 慢波速 近特征线 为