简介:研制了二维多介质流体程序,主要包括单介质内高精度流体力学计算,多介质混合网格内各种介质输运过程和压力驰豫平衡过程计算、实际状态方程的黎曼解计算。流体计算分别采用高分辨两步PPM(ParabolicPiecewiseMethod)算法、TVD(TotalVariationDiminishing)算法和FCT(FluxCorrected—Transport)算法,流体界面追踪采用VOF(Volume-of-Fluid)。数值求解可压缩多流体方程组和可压缩VOF方程。二维界面追踪分别采用一阶精度Youngs方法和二阶精度Elivira方法,三维界面追踪采用一阶精度Youngs方法,
简介:重点研究了不同预热对烧蚀RT弱非线性模耦合的影响。首先不同预热情况烧蚀RT不稳定性的线性增长率曲线明显不同。不同预热情况的模拟结果都表明:烧蚀RT的模耦合系数是k/kc的函数,二次和三次谐波的产生系数C(k)和D(k)遵从定标关系C(k):c1(1-c2k/kc),D(k)=d1-d2k/kc+d3(k/kc)^2。弱预热和中等预热情况,二次谐波的产生系数约只有经典RT的50%,三次谐波的产生系数约只有经典RT的20%。两体模耦合起主要作用,n次高次谐波可近似看作是由(”一1)次两体模耦合产生的。由于烧蚀RT不稳定性两体模耦合系数小于1,所以高次谐波的产生系数很快衰减。由于非线性作用变弱,单模扰动的非线性饱和阈值明显增大,意味着基模有较长时间的线性增长,另外,向长波长方向的模耦合系数明显大于经典尺丁不稳定性数值,因此,烧蚀RT非线性模耦合容易形成长的尖顶,从而对点火构成严重威胁。
简介:采用PPM方法数值求解Euler方程;采用Shyue提出的考虑压力平衡的混合网格状态方程的处理方法,完成R-M不稳定性问题后期混合的数值模拟。界面不稳定性后期混合具有明显的三维特征,二维计算不能分辨后期混合流体团之间三维扭曲拉伸现象,因此要求三维数值模拟。另一方面,界面不稳定性后期,通过非线性作用,小尺度运动被充分激发,必须模拟从大尺度到小尺度的级串现象,因此数值模拟要求很高的空间分辨率,要求大规模数值计算。由此我们采用MPI、应用区域分解方法完成程序并行化,并行程序具有较好的可扩展性。
简介:强激光在冕区等离子体中传播到临界面附近生成相对论电子和相对论电子束流在随后较长一段稠密等离子体区的能量传输是快点火中的关键问题。对快点火条件下的激光等离子体参数,临界面附近产生的前向快电子电流往往超过阿尔芬极限电流,必须在稠密等离子体中产生中和回流,快电子流才能在稠密等离子体中向前输运。横向电磁不稳定性(类Weibel不稳定性,WI)和纵向静电双流不稳定性(TSI)很容易在这种电子双流体系中激发,前向电子束会被调制或成丝状结构,同时激发电磁场,粒子部分动能会转化为电磁场能量。不稳定性在非线性饱和后,发生电流丝的合并、磁场重联等过程,部分电磁场能量会再转化为粒子能量,表现为对离子体的横向加热。Weibel不稳定性的作用可能形成围绕传播电子束的磁通道,对快电子的定向和准直传播是重要的。TSI激发的纵向静电场对磁场通道会有明显的调制甚至破坏作用,直接影响高能电子流从激光吸收区到燃料压缩区的准直传播。