电离层加热中朗缪尔崩塌的Vlasov方程数值模拟刘默然1周晨1.doc

电离层加热中朗缪尔崩塌的Vlasov方程数值模拟 刘默然1,周晨1* ,赵正予1,张援农1 1 武汉大学电子信息学院,空间物理系,430072 通常的电离层等离子体数值模拟技术可以按照数学模型分为两类，一类是基于等离子体的流体模型，通过等离子体的密度、流速、温度等局部平均量的演化方程进行的数值模拟研究；另一类是基于动理学模型，即在向空间中研究等离子体分布函数的演化。对于等离子体来说，流体力学的连续性方程、动量方程和能量方程实际上是Vlasov方程的一阶、二阶和三阶矩方程，是Vlasov方程对于速度空间做平均的结果，因此必须适当地假设速度空间分布函数以及局域输运系数才能建立。 流体模拟中有关电离层加热中朗缪尔湍动现象主要是基于Zakhrov模型。本文从动理学模拟出发，对Vlasov方程进行数值模拟，探究了电离层加热中朗缪尔崩塌的过程。 Vlasov方程如下所示： (1) 在此，考虑各向同性介质中的等离子体，即忽略地磁场的影响。对于等离子体所受的力，只考虑外加电场和和极化电场对等离子体产生的力。 对于三维的Vlasov方程，事实上是3(空间维度)+3(相空间纬度)维的模拟，这对于现有的计算能力是一个很大的挑战。基于此，我们假设等离子体是水平分层的，即只在一个方向上(x方向)存在梯度，则空间维度简化为一维。对于相空间的简化，我们假设外加电场是沿x方向的,则整个等离子体的受力