The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbit...The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbits. Demonstrating significantly higher efficiency, this advanced method is capable of accomplishing the simulation in less than one-third of the time of directly computing the guiding center motion. In contrast to the CPP-based Boris algorithm, this approach inherits the advantages of the AVF method, yielding stable trajectories even achieved with a tenfold time step and reducing the energy error by two orders of magnitude. By comparing these two CPP algorithms with the traditional RK4 method, the numerical results indicate a remarkable performance in terms of both the computational efficiency and error elimination. Moreover, we verify the properties of slow manifold integrators and successfully observe the bounce on both sides of the limiting slow manifold with deliberately chosen perturbed initial conditions. To evaluate the practical value of the methods, we conduct simulations in non-axisymmetric perturbation magnetic fields as part of the experiments,demonstrating that our CPP-based AVF method can handle simulations under complex magnetic field configurations with high accuracy, which the CPP-based Boris algorithm lacks. Through numerical experiments, we demonstrate that the CPP can replace guiding center dynamics in using energy-preserving algorithms for computations, providing a new, efficient, as well as stable approach for applying structure-preserving algorithms in plasma simulations.展开更多
基金supported by National Natural Science Foundation of China (Nos. 11975068 and 11925501)the National Key R&D Program of China (No. 2022YFE03090000)the Fundamental Research Funds for the Central Universities (No. DUT22ZD215)。
文摘The classical Pauli particle(CPP) serves as a slow manifold, substituting the conventional guiding center dynamics. Based on the CPP, we utilize the averaged vector field(AVF) method in the computations of drift orbits. Demonstrating significantly higher efficiency, this advanced method is capable of accomplishing the simulation in less than one-third of the time of directly computing the guiding center motion. In contrast to the CPP-based Boris algorithm, this approach inherits the advantages of the AVF method, yielding stable trajectories even achieved with a tenfold time step and reducing the energy error by two orders of magnitude. By comparing these two CPP algorithms with the traditional RK4 method, the numerical results indicate a remarkable performance in terms of both the computational efficiency and error elimination. Moreover, we verify the properties of slow manifold integrators and successfully observe the bounce on both sides of the limiting slow manifold with deliberately chosen perturbed initial conditions. To evaluate the practical value of the methods, we conduct simulations in non-axisymmetric perturbation magnetic fields as part of the experiments,demonstrating that our CPP-based AVF method can handle simulations under complex magnetic field configurations with high accuracy, which the CPP-based Boris algorithm lacks. Through numerical experiments, we demonstrate that the CPP can replace guiding center dynamics in using energy-preserving algorithms for computations, providing a new, efficient, as well as stable approach for applying structure-preserving algorithms in plasma simulations.
文摘背景:腰椎小关节炎是引起下腰痛的一个主要原因,目前主要依靠MRI进行初步定性诊断,但仍有一定漏诊、误诊的概率发生,因此MR T2^(*)mapping成像技术有望成为定量检查腰椎小关节炎软骨损伤的重要检测手段。目的:探讨MR T2^(*)mapping成像技术在定量分析腰椎小关节炎软骨损伤退变中的应用价值。方法:收集南京医科大学第四附属医院2020年4月至2022年3月门诊或住院合并下腰痛共110例患者,设为病例组;同时招募无症状志愿者80例,设为对照组。对所有纳入对象L1-S1的小关节行3.0 T MR扫描,获取T2^(*)mapping横断位图像和T2WI图像,分别对所有小关节软骨进行Weishaupt分级及T2^(*)值测量,收集数据并行统计学分析。不同小关节Weishaupt分级之间小关节软骨T2^(*)值比较采用单因素方差分析。结果与结论:①经统计分析发现,病例组腰椎小关节软骨T2^(*)值(17.6±1.5)ms明显较对照组(21.4±1.3)ms降低,差异有显著性意义(P<0.05);②在病例组中,随着腰椎小关节Weishaupt分级增加,小关节软骨T2^(*)值也呈逐渐下降趋势,且这种差异有显著性意义(P<0.05);③提示T2^(*)mapping能够较好地显示腰椎小关节软骨损伤的早期病理变化,腰椎小关节软骨的T2^(*)值能够定量评估腰椎小关节的软骨损伤程度;T2^(*)mapping成像技术能为影像学诊断腰椎小关节炎软骨早期损伤提供很好的理论依据,具有重要的临床应用价值。
