CN111581876A - 一种用于圆柱坐标系粒子模拟的粒子源求解方法 - Google Patents

一种用于圆柱坐标系粒子模拟的粒子源求解方法 Download PDF

Info

Publication number
CN111581876A
CN111581876A CN202010400284.9A CN202010400284A CN111581876A CN 111581876 A CN111581876 A CN 111581876A CN 202010400284 A CN202010400284 A CN 202010400284A CN 111581876 A CN111581876 A CN 111581876A
Authority
CN
China
Prior art keywords
particle
grid
source
particles
delta
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010400284.9A
Other languages
English (en)
Other versions
CN111581876B (zh
Inventor
金晓林
周明鉴
郭胜龙
黄桃
杨中海
李斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202010400284.9A priority Critical patent/CN111581876B/zh
Publication of CN111581876A publication Critical patent/CN111581876A/zh
Application granted granted Critical
Publication of CN111581876B publication Critical patent/CN111581876B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/25—Design optimisation, verification or simulation using particle-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及圆柱坐标系下粒子模拟粒子源求解领域，具体为一种用于圆柱坐标系粒子模拟的粒子源求解方法。本发明采用了一种新的体积计算公式，通过根据当前求解粒子的位置构造一个以该粒子为中心的网格，然后使用构造所得网格求解体积，进而再计算电荷源；并进一步的提供了当rn<Δr/2时，步骤3使用体积公式(2)进行粒子源求解时，通过添加修正系数χ使得计算结果更加准确。减小了传统粒子源求解中由于体积计算带来的物理噪声，使得粒子源求解结果的精度进一步增加。

Description

一种用于圆柱坐标系粒子模拟的粒子源求解方法
技术领域
本发明涉及圆柱坐标系下粒子模拟粒子源求解领域，具体为一种用于圆柱坐标系粒子模拟的粒子源求解方法。
背景技术
在物理研究领域中，许多物理问题为非线性强耦合问题，仅仅依靠传统的解析方法与分析很难获得准确的结果，而试验方法也常常受限于设备和环境。随着计算机技术的发展，数值计算方法成为等离子体物理研究领域一种重要的研究方法。粒子模拟方法是研究带电粒子运动与电磁场相互作用的数值计算方法之一。
常规粒子模拟方法计算流程的核心步骤为以下三步：电磁场求解、粒子运动求解和粒子源求解。采用时域有限差分的方法求解离散麦克斯韦方程组，得到网格点上的电磁场。粒子运动求解在解得的电磁场分布基础上，根据粒子的位置，计算得到粒子所属网格编号，再按照一定的方式，将该网格点上的电磁场插值求解得到该粒子所在位置处的电磁场，再对带电粒子洛伦兹运动方程进行求解，更新带电粒子的位置和速度等参量。带电粒子更新状态后，根据粒子的新位置计算得到其所属网格编号，然后将该粒子权重分配，得到网格点的电荷源或是电流源。
由于粒子源求解是粒子模拟算法中的核心步骤之一，它的求解精度对整体粒子模拟计算精度有重要的影响。为了求解网格点上的电磁场分布，需要求得网格点上的电荷源或电流源，因此粒子源求解需要按照一定的分配方法(例如线性分配、面积权重分配等)将粒子电荷源或者是电流源分配到相关网格点上。在不同坐标系下，粒子源求解有不同的计算公式。对于某些物理问题，使用圆柱坐标系进行建模表征能带来诸多便利。假设有m个带电粒子，当前求解的粒子编号为n，所属网格编号为(i,j,k)，每个粒子的带电量大小为qn(n＝1,2,......,m)，速度为vn，网格的体积大小为Vi,j,k，径向步长为Δr，轴向步长为Δz，角向步长为
Figure BDA0002489123750000011
分配方式以线性分配为例，假设归一化后粒子位置与格点在i,j,k三个方向上的距离分别为ωn1,ωn2,ωn3，则传统的圆柱坐标系下粒子模拟的粒子源求解方法可用公式(1)表示：
Figure BDA0002489123750000021
由公式(1)可知，电荷源ρi,j,k求解对每个带电粒子逐个进行计算，并且按照一定的分配方式将电荷源分配到网格的对应格点上，同理，电流源Ji,j,k求解先计算得到每个带电粒子运动产生的电流，再按照一定的分配方式将电流源分配到网格的对应格点上。从公式(1)中可以看出，传统的粒子源计算方法使用的体积为带电粒子所在网格的体积，其体积计算公式中体积与粒子所在网格编号的关系如图1所示。可以看出，假设粒子的径向位置r满足 iΔr<r<(i+1)Δr时，无论r取何值，体积均采用了该粒子所在网格的体积即
Figure BDA0002489123750000022
这在图像上表征为一个阶跃函数，该种体积计算方法带来了较大的物理噪声。
发明内容
针对上述存在的问题，为解决粒子源求解因体积计算方法存在较大噪声的问题，本发明专利提出了一种用于圆柱坐标系粒子模拟的粒子源求解方法。
一种用于圆柱坐标系粒子模拟的粒子源求解方法，具体步骤如下:
步骤1、将圆柱坐标系下建立的模拟区域进行网格划分，获取粒子模拟方法粒子运动求解之后粒子的属性，包括粒子的带电量、粒子的速度与位置、粒子所在网格的编号。
步骤2、构造网格。
设有m个带电粒子，当前求解粒子的编号为n，所属网格编号为(i,j,k)，其电量为qn，速度为vn，粒子的径向位置为rn。根据当前求解粒子的位置构造一个以粒子为中心的网格，角度大小不变，使用构造的网格的体积来计算电荷源。如图2所示，构造的网格以当前求解粒子的径向位置rn为中心，分别沿径向内外扩充Δr/2，轴向前后扩充Δz/2，最终得到的网格径向步长为Δr，轴向步长为Δz。
特别地，当rn<Δr/2时，构造的网格是以轴线为中心，径向步长为Δr，轴向步长为Δz，角向步长为
Figure BDA0002489123750000023
的扇形柱体。
步骤3、计算步骤2构造网格的体积Vn，体积Vn计算公式如(2)所示：
Figure BDA0002489123750000031
步骤4、根据步骤3的体积计算公式(2)，将当前求解粒子电荷源或者是电流源分配到相关网格点上进行求解。
进一步的，所述步骤4分配粒子电荷源或电流源的方法为线性分配，假设归一化后当前求解粒子位置与格点在i,j,k三个方向上的距离分别为ωn1,ωn2,ωn3，采用公式(3)求解网格点上粒子源：
Figure BDA0002489123750000032
步骤5、循环执行步骤2至步骤4，直至遍历完所有m个粒子。
进一步的，所述步骤3中使用体积公式(2)进行粒子源求解，当rn<Δr/2时，添加修正系数χ使得计算结果更加准确，修正系数χ和添加修正系数的体积计算公式如公式(4)所示：
Figure BDA0002489123750000033
本发明采用了一种圆柱坐标系下粒子模拟的粒子源求解方法，通过根据当前求解粒子的位置构造一个以粒子为中心的网格，然后使用构造所得网格求解体积，进而再计算电荷源。减小了传统粒子源求解中由于体积计算带来的物理噪声，使得粒子源求解结果的精度进一步增加。
附图说明
图1为传统体积计算方法体积与径向距离的关系；
图2为本发明以粒子为中心构造网格的平面构造示意图；
图3为本发明以粒子为中心构造网格的轴向构造示意图；
图4为实施例体积计算方法中体积与径向距离的关系。
具体实施方式
下面结合具体实施例和附图对本发明做进一步详细说明。
步骤1、在圆柱坐标系下建立高为100cm，半径为50cm的圆柱作为模拟区域，模拟区域内无静场分布。将该模拟区域进行网格划分，划分的网格的轴向步长为Δz＝1cm，径向步长为Δr＝1cm，角向步长为
Figure BDA0002489123750000043
靠近轴线的网格编号较小，网格的起始编号i＝0,j＝0,k＝0。在圆柱区域均匀放置1×105个电子，采用固定权重，权重设置为5×104，则粒子带电量为 q＝8.45×10-15C，粒子的速度为v＝1.0272×106m/s，方向随机。
获取粒子模拟方法粒子运动求解之后粒子的属性，包括粒子的带电量、粒子的速度与位置、粒子所在网格的编号。
步骤2、按照粒子的编号顺序对每个粒子的粒子源进行计算。假设当前求解的粒子编号为n(n＝1,2,......,1×105)，所属网格编号为(i,j,k),根据该粒子的径向位置rn构造网格。构造的网格以粒子的径向位置rn为中心，分别沿径向内外扩充0.5cm，轴向前后扩充0.5cm，角向左右扩充π/12以保证角度大小不变。当rn≤0.5cm时，构造的网格以轴线为中心，半径为1cm，高度为1cm，角度为π/6的扇形柱体。
步骤3、使用体积公式(2)进行粒子源求解时，当所求解的粒子的径向位置rn≤0.5cm，根据电流连续定理，为了保证电流连续性定理，添加修正系数χ使得计算结果更为准确，修正系数χ由公式(5)确定，计算得到χ＝0.8，修正后的体积计算公式曲线图如图3所示。
Figure BDA0002489123750000041
步骤4、使用构造的网格计算体积Vn，公式如式(6)所示：
Figure BDA0002489123750000042
步骤5、根据新的体积计算公式以及归一化后粒子在原网格上距离格点的距离ωn1,ωn2，并以线性分配为例，采用公式(7)求解网格点上粒子源：
Figure BDA0002489123750000051
由图1和图3的对比可知，原来的体积计算方法有较大的物理噪声，而新的体积计算方法减小了这种物理噪声，使得粒子源求解结果更加准确。

Claims (3)

1.一种用于圆柱坐标系粒子模拟的粒子源求解方法，包括以下步骤：
步骤1、将圆柱坐标系下建立的模拟区域进行网格划分，获取粒子模拟方法粒子运动求解之后粒子的属性，包括粒子的带电量、粒子的速度与位置、粒子所在网格的编号；
步骤2、构造网格；
设有m个带电粒子，当前求解粒子的编号为n，所属网格编号为(i,j,k)，其电量为qn，速度为vn，粒子的径向位置为rn；根据当前求解粒子的位置构造一个以粒子为中心的网格，角度大小不变，使用构造的网格的体积来计算电荷源；构造的网格以当前求解粒子的径向位置rn为中心，分别沿径向内外扩充Δr/2，轴向前后扩充Δz/2，最终得到的网格径向步长为Δr，轴向步长为Δz；
当rn<Δr/2时，构造的网格是以轴线为中心，径向步长为Δr，轴向步长为Δz，角向步长为
Figure FDA0002489123740000013
的扇形柱体。
步骤3、计算步骤2构造网格的体积Vn，体积Vn计算公式如(2)所示：
Figure FDA0002489123740000011
步骤4、根据步骤3的体积计算公式(2)，将当前求解粒子电荷源或者是电流源分配到相关网格点上进行求解。
步骤5、循环执行步骤2至步骤4，直至遍历完所有m个粒子。
2.如权利要求1所述用于圆柱坐标系粒子模拟粒子源求解的方法，其特征在于:
所述步骤4分配粒子电荷源或电流源的方法为线性分配，设归一化后当前求解粒子位置与格点在i,j,k三个方向上的距离分别为ωn1,ωn2,ωn3，采用公式(3)求解网格点上粒子源：
Figure FDA0002489123740000012
3.如权利要求1所述用于圆柱坐标系粒子模拟粒子源求解的方法，其特征在于：
所述步骤3中使用体积公式(2)进行粒子源求解，当rn<Δr/2时，添加修正系数χ，修正系数χ和添加修正系数的体积计算公式如公式(4)所示：
Figure FDA0002489123740000021
CN202010400284.9A 2020-05-13 2020-05-13 一种用于圆柱坐标系粒子模拟的粒子源求解方法 Active CN111581876B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010400284.9A CN111581876B (zh) 2020-05-13 2020-05-13 一种用于圆柱坐标系粒子模拟的粒子源求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010400284.9A CN111581876B (zh) 2020-05-13 2020-05-13 一种用于圆柱坐标系粒子模拟的粒子源求解方法

Publications (2)

Publication Number Publication Date
CN111581876A true CN111581876A (zh) 2020-08-25
CN111581876B CN111581876B (zh) 2023-07-07

Family

ID=72123113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010400284.9A Active CN111581876B (zh) 2020-05-13 2020-05-13 一种用于圆柱坐标系粒子模拟的粒子源求解方法

Country Status (1)

Country Link
CN (1) CN111581876B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705969A (zh) * 2021-07-26 2021-11-26 西安交通大学 一种圆柱坐标系下电流电荷分配方法及系统
CN114707379A (zh) * 2022-03-29 2022-07-05 西北工业大学 一种基于圆柱坐标系约束下的冷喷涂颗粒碰撞模拟方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102449728A (zh) * 2009-05-29 2012-05-09 塞莫费雪科学（不来梅）有限公司 带电粒子分析仪以及带电粒子分离方法
US20140257772A1 (en) * 2013-03-06 2014-09-11 Exa Corporation Flow-induced noise source identification
CN108460188A (zh) * 2018-02-05 2018-08-28 电子科技大学 一种应用于pic静电模型的电荷分配有限元fem求解算法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102449728A (zh) * 2009-05-29 2012-05-09 塞莫费雪科学（不来梅）有限公司 带电粒子分析仪以及带电粒子分离方法
US20120138785A1 (en) * 2009-05-29 2012-06-07 Makarov Alexander A Charged Particle Analysers And Methods Of Separating Charged Particles
US20140257772A1 (en) * 2013-03-06 2014-09-11 Exa Corporation Flow-induced noise source identification
CN108460188A (zh) * 2018-02-05 2018-08-28 电子科技大学 一种应用于pic静电模型的电荷分配有限元fem求解算法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113705969A (zh) * 2021-07-26 2021-11-26 西安交通大学 一种圆柱坐标系下电流电荷分配方法及系统
CN113705969B (zh) * 2021-07-26 2022-12-09 西安交通大学 一种圆柱坐标系下电流电荷分配方法及系统
CN114707379A (zh) * 2022-03-29 2022-07-05 西北工业大学 一种基于圆柱坐标系约束下的冷喷涂颗粒碰撞模拟方法

Also Published As

Publication number Publication date
CN111581876B (zh) 2023-07-07

Similar Documents

Publication Publication Date Title
Bagheri et al. Comparison of six simulation codes for positive streamers in air
Wang et al. Implicit and electrostatic particle-in-cell/Monte Carlo model in two-dimensional and axisymmetric geometry: I. Analysis of numerical techniques
An et al. An improved iteration loop for the three dimensional quasi-static particle-in-cell algorithm: QuickPIC
Miyake et al. New electromagnetic particle simulation code for the analysis of spacecraft-plasma interactions
CN111581876A (zh) 一种用于圆柱坐标系粒子模拟的粒子源求解方法
Candlish et al. RAyMOND: an N-body and hydrodynamics code for MOND
CN108416107B (zh) 一种应用于pic的推动粒子运动有限元算法
Cao et al. An iterative immersed finite element method for an electric potential interface problem based on given surface electric quantity
Ong et al. A fast algorithm for three-dimensional potential fields calculation: fast Fourier transform on multipoles
JP6157941B2 (ja) プラズマシミュレーション方法及びプラズマシミュレーションプログラム
Luzum et al. Methods for systematic study of nuclear structure in high-energy collisions
Pfeiffer et al. A Particle-in-Cell solver based on a high-order hybridizable discontinuous Galerkin spectral element method on unstructured curved meshes
Frenod et al. An exponential integrator for a highly oscillatory Vlasov equation
Lund et al. Weak scaling of the parallel immersed-finite-element particle-in-cell (PIFE-PIC) framework with lunar plasma charging simulations
Jo et al. Multi-dimensional electrostatic plasma simulations using the particle-in-cell method for the low-temperature plasmas for materials processing
Garza et al. High-order Chebyshev-based Nyström methods for electromagnetics
Cornet et al. A new algorithm for charge deposition for multiple-grid method for PIC simulations in r–z cylindrical coordinates
CN108460188B (zh) 一种应用于pic静电模型的电荷分配有限元fem求解算法
Roucka et al. Extending PIC models to higher pressures—enhanced model of collisions
Rivara et al. Multithread parallelization of lepp-bisection algorithms
Senichev et al. Storage ring EDM simulation: methods and results
Voropaeva et al. New motion algorithm in the particle-in-cell method
CN113821915A (zh) 一种轴对称电子光学系统的快速计算方法
Rafieian Najaf Abadi et al. Dependence of emittance of the He+ beam extracted from a Penning ion source on the design of the source, a Particle-In-Cell simulation
Fox Parallelization of particle-in-cell simulation modeling Hall-effect thrusters

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant