磁螺度的三维数值模拟_20110607.ppt.pptVIP

* * 磁螺度的三维数值模拟 杨尚斌 合作者：J?rg Büchner Jean Santos Carol * 磁螺度的定义 Magnetic helicity 规范不变性 * 相对磁螺度 (Ref. Berger 1984) E=mgh h=0 h=z E=0 E=mgz Gravitational potential energy * 重力势能 相对磁螺度 The relative magnetic helicity is gauge invariant * 相对磁螺度的传输 Magnetic helicity change rate in a volume : The first term represents helicity dissipation inside the volume The second term represents the helicity flux across the boundary * Ideal conducting fluid 1 2 1. The first term is due to footpoint motion of the flux. 2. The second term is due to emergence of new magnetic flux. * Deducing V Local Correlation Tracking can be used to calculate the first term The two terms in the formula can combine together under the ideal magneto hydrodynamics if we use LCT method to deduce V Demoulin et al. (2003) Chae et al. (2001) * Magnetic helicity at a certain time t We use the two formulas to calculate the magnetic helicity in the observation. V S * Counterclockwise Vortex H0 V S V S Clockwise Vortex H0 * 磁螺度与数值模拟 磁螺度的传输从观测上我们可以直接得到，并取得了一系列 的观测结果；但是得到的磁场螺度是否真正包含了所有传输的 螺度? 这需要用数值模拟的方法来检验 （Demoulin and Pariat 2009)； (Devore 2000) * Fan and Gibson 2004 * (Amari 2003) Cheung et al. 2005 * Ravindra et al. 2008 * x y z * 六个边界上: 以底部边界为例 边界条件 * 求解泊松方程 (a) (b) We use a fast Poisson solver based on the HODIE finite-difference scheme (Boisvert, R. 1984) to resolve of the two Poisson equation . We also clean the divergence of and * Magnetic helicity information from the data Magnetic helicity in the simulation box Magnetic helicity dissipation in the simulation box Magnetic helicity flux across the boundary For one data with grid points 131×131×131. It takes about 5 minutes to get one magnetic helicity result. * 数据驱动（Data-Driven）的MHD模型 Buchner et al. 2004ab; Buchner 2006; Santos and Buchner 2007 Santos et al. 2008 欧姆定律： 安培定律： 状态方程： 质量守恒方程 动量方程 磁感应方程 能量方程 * * 理想磁场流体中 Helicity in the simulation box Integral of