计算天文讲义 定常问题的跨音速流.pptVIP

14. 定常问题的跨音速流 15. 非定常跨音速流（拟一维） 边界条件 初始条件 结果 16. 可视化及几个技巧 * * 尽管Biemann (1951)指出太阳持续发射corpuscles (微粒)从而导致彗尾总背着太阳，Chapman (1957)仍然考虑含热传导的静止太阳大气 Parker的太阳风解 假设日冕等温 Parker (1958, ApJ, 128, 664)指出日冕不能处于静力平衡，而会膨胀，从而形成太阳风。当时极具争议。1959年Lunik III及Venus I卫星证实。 不同的C导致不同的解，其中跨声速解V(对应C=-3)符合太阳风解，VI符合吸积盘解。I和II为非物理解。III为超声速解，IV为亚声速解。 在v=vc处存在临界点，对应 模型之一： 考虑沿磁流管的流动的一维多方过程 原则上给定S(r)分布及r0处的边界值即可求出各物理量的分布。 对球对称流，引入无量纲参数： 用此式调节U0寻找跨声速解十分不方便 对其求导，还原为微分方程 其中 由临界点条件：P=Q=0 为判断该临界点的性质，在其附近作代换 一级近似 下页 其通解为 定义 结点 鞍点 中心点 焦点 对本问题， 无太阳风解 由罗必塔法则求得 存在正解的条件 将U*和Z*代入原无量纲方程 此为边界参数U0和H须满足的关系，即在r0只能给定其二者之一。 另，可由临界点处Q=0证明该临界点为声速点 (1) 选定 (2) 选取边条T0和U0二者之一（可取在太阳表面或地球处） (3) 计算临界点坐标 和 ，以及由罗必塔法则求临界点的 (4) 从临界点出发，取合适步长，采用龙格库塔法向两边对 微分方程进行积分，其起始斜率为 (5) 给定参考密度 ，由连续性方程求 ，由多方关系给出 (6) 调整 使得地球处或太阳表面处的参数与观测值接近。 综上所述，采用多方模型求解太阳风的步骤如下： 若选定T0 ，则由 定H，然后由 定U0 若选定U0 ，则由上式求H 其它模型： 考虑热传导甚至磁场等 在多方模型中，我们须给定两个自由参数T0和?0（ 除外）。而热传导模型中只能给定一个自由参数。就和观测值比较而言，热传导模型比多方模型更为理想。 考虑磁场后将会出现3个临界点（慢磁声速、阿尔芬速、快磁声速）。可采用相同方法求临界点，然后由临界点出发对微分方程进行积分。 此外，有时可采用非定常方程组，取某初始条件，计算在一定边界条件下系统演化的收敛解。 采用MacCormack格式 predictor corrector 无量纲 化方程 左：2条右行特征线，1条左行线，故需指定2个变量，允许1个变量浮动 右：3条右行特征线，故允许3个变量浮动 对u采用线性外推，其它量固定，即u1=2u2-u3; ρ1=1; T1=1 对所有量采用线性外推 取截面为 理论上知道随坐标增加 密度温度降低，而速度增加 取线性分布 原则上任意，但不接近真实的话收敛慢，甚至发散 三维数据可视化效果较好 IRAF：一般天文资料处理 Pgplot：利用Fortran程序调用 Origin： Matlab： AVS： IDL： Vis5D：免费软件/~billh/vis5d.html 生成VRML文件 MAXIMA 可以在 /projects/maxima 免费下载。 ？如何制作动画 方法1：直接生成MPEG或先生成一系列GIF图像， 然后用Windows软件或Linux的命令whirlgif 即可将此系列图合成为一个动画GIF图。 方法2：利用Java script语言在IE或Netscape中控制一 系列图像的显示（见附件）。 ？如何将两幅灰度图叠加 方法1：直接将两个数组相加 方法2：新数组的奇数点取第一个数组的值，偶数 点取第二个数组的值 ？如何将GIF等图像中的曲线变成数据 方法：IDL aa(where(aa eq 123))=0 ？如何表现3维线条（如磁力线）