计算天文数值计算lec15高级算法之CIP法lec15.pptVIP

19. 高级算法之CIP法 Yabe (矢部) et al. 1991, Comput. Phys. Commun., 66, 219 2001, JCP, 169, 556 发展了CIP (Cubic Interpolated Profile) 高精度差分法，举例如下： 若 c为常数，则其解为 ，其差分形式为 则新时刻的分布化为老时刻某空间点的分布 其求解之最简方法为一次迎风格式（以c 0为例） 但其缺点是强数值耗散。 CIP方法在两格点间 采用3次多项式拟合 拟合系数满足此函数及其导数值在 i-1, i 格点处连续的条件: 若u0 则 i-1? i+1， ?x ? - ?x 因此，CIP法要求同时求解 为此，对原方程求导 同样， 也满足 综上所述， 初条由 和 同时给出 。 计算实例：单波问题 非线性方程 1. 非对流项 即 上两式求解分两步实现 2. 对流项 1. 非对流项 采用时间前差，空间中心差分 常可用交错网格 2. 对流项 采用CIP方法 若u0 则 i-1? i+1， ?x ? - ?x 1D理想流体Euler方程 为使激波解稳定，可在压力中加人为粘性项 计算实例：激波管问题 20. 高级算法之PIC法 F. H. Harlow. 1964, Methods in Comput. Phys., 3, 319 发展了PIC (Particle-in-cell) 方法。它综合了Euler格式和Lagrange格式的优点，在Euler网格上定义场变量，用质点表征流体。主要思想是在Euler网格中配置离散的Lagrange流体质点。 Euler格式：适用于求解含大扭曲、大滑移问题 Lagrange格式：适用于求解撞击界面的多维流体问题。 要点：各物理量（除质量外）均定义在网格中心，混合网格 的质量和能量分别按物质给出。 初始时质点分布： 以2维均匀网格为例，每个网格内质点数目反映局地密度，且至少要有?100个质点；应使计算末每个网格至少要有3－4个质点（也许应更多）。经验表明，一个网格内质点的分布不宜十分规则，采用菱形分布或不规则分布利于减少质点流过网格边界引起的非物理计算波动。 以二维为例，则网格(i, j)的性质为其中的各质点的总和： 考虑人为 粘性项 q 令总共有K个质点，每个质点具有其 质量 ，坐标 ，速度 其中 计算步骤 1. Euler步 以轴对称为例 动量 忽略输运项， 上述方程组变为 对动量方程采用FTCS (forward time-centered space) 格式求中间速度 能量 网格界面 p 按算术平均值计算，若旁边网格是空网格，置界面 p=0；若界面是固壁，则取界面p等于网格p 。 取人为粘性 对能量方程采用交错（拉链）差分 z向同理 若旁边网格是空网格，则U及V均取本网格中心值，若界面是固壁，则法向速度取为0。 PIC方法的特点就是利用质点的运动来计算通过网格边界的输运量 。 2. 质点步 （1）确定质点速度 在Euler步中得出的速度 是网格中心的速度，而质点一般不在中心，每个质点的速度可取为与其相邻的四个网格速度的加权平均值（第k个质点）: 设四个网格的公共交点是 约定：把覆盖(i, j)网格的部分记为4，对角部分记为1， 沿r方向的部分记为2，沿z方向的部分记为3。 （2）计算质点的新坐标 若覆盖的网格为空网格，空网格的速度取为质点所在网格的速度，或取此网格的面积为0。若覆盖的网格在固壁外，其速度取为相邻内网格的反射速度。 （3）重新分配质点 质点运动后可能留在原网格，也可能移动到新网格，从而引起网格质量、动量和能量的变化 缺点 对下列情况效果不是非常好： 非常不均匀的粒子分布 correlated systems complex geometries 21. 高级算法之SPH法 PIC方法同时使用Eulerian和Lagrangian elements，其欧拉网格的平均密度等于网格中的质点数目。 能否不需网格而得到局地密度？ Lucy, 1977, AJ, 82, 1013 Gingold Monaghan, 1977, MNRAS, 181, 375 发展了SPH (Smoothed Particle Hydrodynamics光滑粒子流体动力学