关于期末考试新版7.pptVIP

《计算天文》期末汇报暨 天文系第四届数值模拟研讨会 时间：2009年1月9日下午2：00－6：00 地点：天文系202会议室 报告形式：口头（均利用ppt文件），每人（组）10－15分钟 联系人：P. F. Chen 联系电话电子邮件：chenpf@nju.edu.cn * jgyj 《计算天文》期末汇报暨 天文系第五届数值模拟研讨会 时间：2011年7月10日下午2：00－6：00 地点：天文系202会议室 报告形式：口头（均利用ppt文件），每人（组）10－15分钟 联系人：P. F. Chen 联系电话电子邮件：chenpf@nju.edu.cn * jgyj 《计算天文》期末汇报暨 天文系第六届数值模拟研讨会 时间：2013年9月1日下午2：00－6：00 地点：天文系202会议室 报告形式：口头（均利用ppt文件），每人（组）10－15分钟 联系人：P. F. Chen 联系电话电子邮件：chenpf@nju.edu.cn * jgyj 注意事项： 请大家好好准备； 可以单独做个题目，也可2人为一组合作做题目； 题目可以是我提供的，也可结合你们各自的研究自行选择； 除作报告外，每人（组）请将模拟工作整理出1－3页纸（形式如下），以便于我编辑成文集后发送给与会人士： 太阳耀斑环动力学数值模拟 张三 （南京大学天文系，南京 210093） 1 引言 他人工作等. 2 数值方法 ….. 3 结果 ……含图 4 讨论/结论 参考文献 * jgyj 注意事项： ppt报告中请多用些图甚至动画来显示你的结果； 期间有任何问题请及时与我(苏富特楼203室)联系，讨论； 请记住几个时间表： 6月20日前：请发信件将题目和成员告诉我：我的电子邮箱是 chenpf@nju.edu.cn; BBS帐号是MHD；两者任一。 8月31日前：所有同学将论文(pdf)及报告(ppt文件)发给我 9月1日下午：天文系202室研讨会 云台及个别同学在7月5日上午10点@天文系202室研讨会 * jgyj 19. 高级算法之CIP法 Yabe (矢部) et al. 1991, Comput. Phys. Commun., 66, 219 2001, JCP, 169, 556 发展了CIP (Cubic Interpolated Profile) 高精度差分法，举例如下： 若 c为常数，则其解为 ，其差分形式为 则新时刻的分布化为老时刻某空间点的分布 * jgyj 其求解之最简方法为一次迎风格式（以c 0为例） 但其缺点是强数值耗散。 CIP方法在两格点间 采用3次多项式拟合 拟合系数满足此函数及其导数值在 i-1, i 格点处连续的条件: 若u0 则 i-1? i+1， ?x ? - ?x * jgyj 因此，CIP法要求同时求解 为此，对原方程求导 同样， 也满足 综上所述， 初条由 和 同时给出 。 * jgyj 计算实例：单波问题 * jgyj 非线性方程 1. 非对流项 即 上两式求解分两步实现 2. 对流项 * jgyj 1. 非对流项 采用时间前差，空间中心差分 常可用交错网格 * jgyj 2. 对流项 采用CIP方法 若u0 则 i-1? i+1， ?x ? - ?x * jgyj 1D理想流体Euler方程 为使激波解稳定，可在压力中加人为粘性项 计算实例：激波管问题 * jgyj 20. 高级算法之PIC法 F. H. Harlow. 1964, Methods in Comput. Phys., 3, 319 发展了PIC (Particle-in-cell) 方法。它综合了Euler格式和Lagrange格式的优点，在Euler网格上定义场变量，用质点表征流体。主要思想是在Euler网格中配置离散的Lagrange流体质点。 Euler格式：适用于求解含大扭曲、大滑移问题 Lagrange格式：适用于求解撞击界面的多维流体问题。 * jgyj 要点：各物理量（除质量外）均定义在网格中心，混合网格 的质量和能量分别按物质给出。 初