CN104699993B - 辐射能传播蒙特卡洛算法中的局部网格及光子加密方法 - Google Patents
辐射能传播蒙特卡洛算法中的局部网格及光子加密方法 Download PDFInfo
- Publication number
- CN104699993B CN104699993B CN201510149785.3A CN201510149785A CN104699993B CN 104699993 B CN104699993 B CN 104699993B CN 201510149785 A CN201510149785 A CN 201510149785A CN 104699993 B CN104699993 B CN 104699993B
- Authority
- CN
- China
- Prior art keywords
- grid
- photon
- size
- mesh opening
- medium
- 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.)
- Active
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 11
- 230000003287 optical effect Effects 0.000 claims abstract description 12
- 238000009826 distribution Methods 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 7
- 238000000342 Monte Carlo simulation Methods 0.000 abstract description 18
- 238000004364 calculation method Methods 0.000 description 16
- 238000005516 engineering process Methods 0.000 description 8
- 238000010521 absorption reaction Methods 0.000 description 6
- 210000002615 epidermis Anatomy 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000644 propagated effect Effects 0.000 description 5
- 238000002835 absorbance Methods 0.000 description 4
- 210000004204 blood vessel Anatomy 0.000 description 4
- 210000003491 skin Anatomy 0.000 description 4
- 229940088503 dermalzone Drugs 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 239000013028 medium composition Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 2
- 238000004513 sizing Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 206010048768 Dermatosis Diseases 0.000 description 1
- 241000132179 Eurotium medium Species 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 239000003989 dielectric material Substances 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000012770 industrial material Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 208000017520 skin disease Diseases 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Landscapes
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Optical Communication System (AREA)
Abstract
本发明提供一种辐射能传播蒙特卡洛算法中的局部网格及光子加密方法：首先，对传播介质内不同组成部分依据对应组成部分的光学特性及形状结构选用尺寸不同的网格；其次，根据最稀疏网格的网格尺寸选取所述传播介质内光传播计算所需的总光子数，总光子数的大小应满足按所述最稀疏网格的网格尺寸计算时网格内光子能量分布的统计学要求，当光子由网格尺寸较大的稀疏网格进入网格尺寸较小的细密网格时，光子进行分裂，本发明在现有“基于网格的蒙特卡洛方法”的基础上，依据传播介质的光学特性及形状结构，提出局部网格及光子加密方法，极大地提高了蒙特卡洛方法的计算效率。
Description
技术领域
本发明属于热辐射计算领域，主要涉及针对不同传播介质或具有多种组成的传播介质内辐射能传播蒙特卡洛计算方法的改进算法。
背景技术
辐射传播是热量传递的主要三种形势之一，在工业生产及生活领域具有广泛的应用。例如利用复色光辐射能——太阳能进行发电，或利用单色光辐射能——激光进行工业材料处理或临床激光热疗。由于辐射能在介质内部(如太阳能集热器、工业材料或生物组织内)的传播并不可见，因而数值计算成为了一种有效地“可视化”方法，以增强对辐射能与传播介质间相互作用的理解并以此提高其工程实际应用的效果。在众多的热辐射计算方法中蒙特卡洛方法是目前应用最广的辐射能传播及分布求解方法之一。
蒙特卡洛方法是利用统计学模型求解物理过程的数学方法。由于辐射能传播介质的几何形状往往具有高度复杂性，光子在介质内不同组分界面(以下简称界面)上的折射及散射计算是蒙特卡洛方法实施中的关键问题。Wang等的多层蒙特卡洛模拟方法(Multi-layered Monte Carlo method)是目前应用较广的光传播蒙特卡洛方法(Lihong W,Jacques SL,Liqiong Z.MCML-Monte Carlo modeling of light transport in multi-layered tissues.Computer Methods and Programs in Biomedicine,1995,47(2):131-146.)，然而该方法只能处理与坐标平面垂直的界面。Pfefer等人开发的“基于网格的蒙特卡洛方法”(Voxel-based Monte Carlo method)是目前公认的可以较好处理这类问题的蒙特卡洛方法(Pfefer TJ,Barton JK,Chan EK,et al.A three-dimensional modularadaptable grid numerical model for light propagation during laser irradiationof skin tissue.Ieee Journal of Selected Topics in Quantum Electronics,1996,2(4):934-942.)。与传统蒙特卡洛方法中网格单纯的能量统计功能不同，基于网格的蒙特卡洛方法最大的特点就是光子的运行与网格设置紧密相连，连续的复杂曲面被离散的微小的网格界面所替代，利用光子与网格界面的相互作用替代与曲面间的相互作用。通过该方法极大地简化了复杂界面上折射、反射的求解难度。
Binzoni等人研究结果表明，“基于网格的蒙特卡洛计算方法”中，由于网格尺寸选取不当造成的计算误差可达20％以上(Binzoni T,Leung TS,Giust R,et al.Lighttransport in tissue by 3D Monte Carlo:Influence of boundaryvoxelization.Computer Methods and Programs in Biomedicine,2008,89(1):14-23.)。另一方面，蒙特卡洛方法具有统计学特性，其结果依赖于光子数的选择，光子数越多其计算结果越接近真实结果，网格数越多所需要的光子数也越多。因而，蒙特卡洛方法的计算效率与网格及光子数的选取有着密切的关系。但是，截至目前为止，针对蒙特卡洛方法中网格尺寸及光子数选择判据的研究依然未见报道。蒙特卡洛计算中网格尺寸及光子数的选择多根据经验而无科学依据。
发明内容
本发明的目的在于提供一种辐射能传播蒙特卡洛算法中的局部网格及光子加密方法，可以显著提高基于网络的蒙特卡洛算法计算辐射能传播及分布时的效率。
为达到上述目的，本发明采用了以下技术方案：
首先，对辐射能传播介质(以下简称介质)不同组成部分依据该介质对应组成部分的光学特性及形状结构选用尺寸不同的网格；其次，根据最稀疏网格的网格尺寸选取所述介质内辐射能传播蒙特卡洛计算所需的总光子数，总光子数的大小应满足按所述最稀疏网格的网格尺寸计算时网格内光子能量分布的统计学要求，当光子由网格尺寸较大的稀疏网格进入网格尺寸较小的细密网格时，光子进行分裂，分裂后光子的数目应满足所述细密网格内的光子能量分布的统计学要求。
对所述介质中光吸收较弱的组成部分采用网格尺寸较大的稀疏网格，而对所述介质中光吸收较强的组成部分采用网格尺寸较小的细密网格；对所述介质中形状结构较简单的组成部分采用网格尺寸较大的稀疏网格，而对所述介质中形状结构较复杂的组成部分采用网格尺寸较小的细密网格。
所述介质不同组成部分的网格尺寸的最大允许值依据如下关系式进行选取：
其中，Δw*是无量纲网格尺寸，Δw*＝Δwμt，Δw为实际网格尺寸，μt为辐射能在介质中的衰减系数，是介质中不同组成部分边界上无量纲曲率半径最大值，rmax是介质中不同组成部分边界上曲率半径最大值。
为满足光子能量分布的统计学要求，单个网格内平均光子数应大于等于5。
本发明的有益效果体现在：
本发明在现有基于网格的蒙特卡洛方法的基础上，依据介质的光学特性及形状结构，提出局部网格及光子加密方法，极大地提高了蒙特卡洛方法的计算效率。
附图说明
图1是含有单根圆柱辐射能吸收体的传播介质示意图；
图2是网格加密区示意图；
图3是非加密区网格划分示意图；
图4是加密区网格划分示意图；
图5是光子数加密示意图；
图6是网格局部加密对蒙特卡洛计算时间(Computation time)的影响；
图7是光子增值倍数对蒙特卡洛计算时间的影响；
图8是无光子分裂(a)及光子分裂10倍(b)及100倍(c)后，加密区光密度分布图。
具体实施方式
下面结合附图和实施例对本发明做详细描述。
利用蒙特卡洛方法计算介质辐射能传播及吸收时网格及光子加密的基本思路是：依据介质各组成部分(简称成分)的形状结构及光学参数，针对各个组成部分划分不同尺寸的网格。光学吸收性强，且形状复杂的组成部分要选取较小的网格尺寸。随后，依据计算区域内最大网格选取蒙特卡洛计算所需的光子数，当光子运行过程中进入较细密的网格区域时，光子分裂为若干个相同权值的子光子，以满足细密网格内蒙特卡洛方法对于光子数的需求。
本发明首先从吸收能量的准确统计及形状不规则界面的准确刻画这两方面要求出发，推导出网格尺寸选取的判据：网格大小的选取与介质的光学性质(μt)及介质内界面的曲率半径(r)有关，光在介质内衰减的越快或介质内界面曲率越大，则所需网格尺寸越小，反之则越大。无量纲网格尺寸可依据如下关系式进行选取：
其中，Δw*是无量纲网格尺寸，Δw*＝Δwμt，Δw为实际网格尺寸，rmax是介质中不同组成部分边界上曲率半径最大值，即介质不同组成部分边界上曲率最大值，是介质不同组成部分边界上无量纲曲率半径最大值，μt为辐射能在所述介质内的衰减系数。
另外，利用单相高吸收介质内能量密度分布一维分析解与不同光子数蒙特卡洛方法计算结果进行对比，引出光子数的选择判据：计算中光子数的选择与网格尺寸及数量关系紧密，网格越小，网格数量越大，则计算所需光子数越多。单个网格内平均光子数应大于等于5，即光子使用总数应大于等于5倍网格数。
随后，本发明根据所得判据发展出“网格及光子加密技术”，可显著提高辐射能传输蒙特卡洛算法的计算效率。通过以上判据可以看出：介质各成分间网格尺寸选取值并不相同，辐射能吸收较强及形状较复杂的介质成分所需的网格尺寸小于辐射能吸收较弱、结构简单的介质成分；并且光子数的选取与网格大小相关，为了满足统计学特性，网格越小所需的光子数就越多。在复杂介质内辐射能传播及分布蒙特卡洛计算中，对不同介质成分采用均匀网格时，网格尺寸选取必须以强吸收、形状复杂的成分为依据，但这对于弱吸收、形状简单的成分而言显然是一种计算资源的浪费。因而，本发明在介质内建立非均匀网格，依据介质的辐射能吸收特性及形状对网格进行局部加密，可以提高“基于网格的蒙特卡洛方法”的计算效率。另一方面，依据判据使用非均匀网格时，若光子数的选取以小网格为依据，则对于大网格而言也是一种计算资源的浪费。因而在使用非均匀网格模拟时，应对网格局部加密区域辅以光子分裂技术，以进一步提高“基于网格的蒙特卡洛方法”的计算效率。
以下，本发明以典型单色光辐射能应用——血管性皮肤病激光热疗中辐射能传输及分布模拟为例，说明网格及光子局部加密技术实施的基本思路及其对计算效率的影响。针对图1所示的单根血管皮肤模型，选择对治疗激光吸收性较强的血管区域为网格加密区。在网格划分中，介质各成分间网格尺寸的选择应遵循网格选取准则。含有单根血管皮肤组织的非均匀网格划分的具体过程如下：
首先，确定加密区的范围。加密区应该包含整个血管，可以略大于血管面积，如图2所示。
随后，对吸收特性较弱的表皮及真皮区域采用稀疏网格进行网格划分。网格划分结果如图3所示。
最后，对包含血管的区域采用细密网格进行网格局部加密。整体网格划分结果二维图如图4所示。本发明所采用的算例中，血管轴平行于y轴，且与xz平面垂直，因而在本发明的网格划分时对其他xz平面采用相同操作即可。
光子局部加密：计算中光子数的选择与网格尺寸及数量关系紧密，网格越小，所需的光子数越多，单个网格内平均光子数应大于等于5。光子使用总数应大于等于5倍网格数。在局部网格及光子加密技术中，非加密区光子数的选择，以及加密区光子数的加密倍数亦应遵循各选取准则。
首先，根据网格尺寸较大区域(如图4表皮及真皮区域)计算所需的总光子数，其大小为以当地网格尺寸计算所得的整个计算区域网格个数的五倍。
另一方面，为了满足细密网格区域(图4中血管附近区域)每个网格内的统计值具有统计学意义，光子进入该区域时应进行光子分裂，光子的分裂倍数应使分裂后光子数满足细密网格内网格数需求。
光子加密的具体实施方法是：首先，规定一个光子加密区域，该区域可以大于也可以等于网格加密区，如图5所示。假设，模拟中第k个光子碰到该区域边界，则该光子分裂成n个，同时为了保证能量守恒，每个新光子的能量为原来的1/n。新光子一个接一个进入光子加密区进行传播，直到n个光子走完为止再重新生成第k+1个光子。由皮肤表面进入组织。
本发明首先采用无网格加密及光子加密的原始程序对图1所示组织模型的光传播进行计算，网格尺寸为了满足血液内的吸收统计特性取为10μm，光子数为107个。在HPxw6600Workstation上所需的计算时间为17.5个小时。以此结果作为标准结果，来说明网格及光子加密对于计算效率的提高及计算精度的影响。
采用非均匀网格技术，不同真皮网格尺寸(Δwd)下蒙特卡洛计算所需的时间如图6所示。对直径为120μm的血管而言，衰减系数为191cm-1，对于圆柱形血管来说，rmax＝60μm， 根据公式(1)计算可得网格尺寸Δw＝0.2/μt＝10.47μm。同时对于表皮而言，衰减系数为79.38cm-1，其与真皮组织的边界为平面，因而rmax趋近于无穷大，根据公式(1)，Δw＝0.2/μt＝25.19μm。对于真皮而言，与表皮类似，Δw＝0.2/μt＝0.2/13.0＝143.9μm，即真皮网格尺寸小于等于143.9μm即可。因而，计算中选取表皮层内网格尺寸为20μm，加密区域网格尺寸为10μm，计算采用的光子数为107个。从图中可以看出，当采用非均匀网格时计算时间明显短于均匀网格。并且，随着真皮内网格尺寸的增加，程序运行的时间不断缩短。当真皮中网格选取为80μm时，程序运行时间仅需要5.5个小时。说明了网格加密技术对计算效率的提高程度。
将光子加密技术与网格加密技术相结合，光子加密倍数n其对计算时间的影响如图7所示。从图中可以看出，仅仅采用网格加密可将计算时间由17.5小时缩短至5.5个小时。若采用加密网格的同时，在加密区内将光子增值10倍，可将计算时间缩短为3.5个小时。然而将n提高到100时，计算时间为2.3个小时，比原始计算时间17.5小时明显降低。
当计算中选取表皮层内网格尺寸为20μm，真皮区域网格尺寸为80μm，加密区域网格尺寸为10μm时，从图8中可以看出，当光子加密数n分别为10倍(图8b)及100倍(图8c)时，加密区光密度分布与图8a无光子加密时基本相同。因而可以说明光子加密可以保证蒙特卡洛方法的计算精度。
本发明的特点如下：针对形状简单、光吸收特性较弱的介质区域采用较粗的计算网格，而在形状复杂、光吸收特性较强的组成部分实施网格局部加密。蒙特卡洛计算时光子数仅需满足粗网格计算需求即可，仅对局部网格加密区域计算时实施光子分裂，使局部加密后的光子数满足加密网格区域的计算需求。
总之，针对复杂结构介质内辐射能传播蒙特卡洛计算时，依据本发明提出的蒙特卡洛方法网格及光子数选取判据，对介质不同组成部分采取局部网格及光子加密，满足计算精度的同时降低计算量，提高了计算效率。
Claims (4)
1.一种辐射能传播蒙特卡洛算法中的局部网格及光子加密方法，其特征在于：包括以下步骤：
首先，对辐射能传播介质不同组成部分依据该介质对应组成部分的光学特性及形状结构选用尺寸不同的网格；其次，根据最稀疏网格的网格尺寸选取所述介质内辐射能传播蒙特卡洛计算所需的总光子数，总光子数的大小应满足按所述最稀疏网格的网格尺寸计算时网格内光子能量分布的统计学要求，当光子由网格尺寸较大的稀疏网格进入网格尺寸较小的细密网格时，光子进行分裂，分裂后光子的数目应满足所述细密网格内的光子能量分布的统计学要求。
2.根据权利要求1所述一种辐射能传播蒙特卡洛算法中的局部网格及光子加密方法，其特征在于：对所述介质中光吸收较弱的组成部分采用网格尺寸较大的稀疏网格，而对所述介质中光吸收较强的组成部分采用网格尺寸较小的细密网格；对所述介质中形状结构较简单的组成部分采用网格尺寸较大的稀疏网格，而对所述介质中形状结构较复杂的组成部分采用网格尺寸较小的细密网格。
3.根据权利要求1所述一种辐射能传播蒙特卡洛算法中的局部网格及光子加密方法，其特征在于：所述介质不同组成部分的网格尺寸的最大允许值依据如下关系式进行选取：
<mrow>
<mfenced open='{' close=''>
<mtable>
<mtr>
<mtd>
<msup>
<mi>&Delta;w</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mn>0.2</mn>
<msubsup>
<mi>r</mi>
<mi>max</mi>
<mo>*</mo>
</msubsup>
</mtd>
<mtd>
<msubsup>
<mi>r</mi>
<mi>max</mi>
<mo>*</mo>
</msubsup>
<mo>&le;</mo>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>&Delta;w</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mn>0.2</mn>
</mtd>
<mtd>
<msubsup>
<mi>r</mi>
<mi>max</mi>
<mo>*</mo>
</msubsup>
<mo>&GreaterEqual;</mo>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，Δw*是无量纲网格尺寸，Δw*＝Δwμt，Δw为实际网格尺寸，μt为辐射能在介质中的衰减系数，是介质中不同组成部分边界上无量纲曲率半径最大值，rmax是介质中不同组成部分边界上曲率半径最大值。
4.根据权利要求1所述一种辐射能传播蒙特卡洛算法中的局部网格及光子加密方法，其特征在于：为满足光子能量分布的统计学要求，单个网格内平均光子数应大于等于5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510149785.3A CN104699993B (zh) | 2015-03-31 | 2015-03-31 | 辐射能传播蒙特卡洛算法中的局部网格及光子加密方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510149785.3A CN104699993B (zh) | 2015-03-31 | 2015-03-31 | 辐射能传播蒙特卡洛算法中的局部网格及光子加密方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104699993A CN104699993A (zh) | 2015-06-10 |
| CN104699993B true CN104699993B (zh) | 2017-10-20 |
Family
ID=53347105
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510149785.3A Active CN104699993B (zh) | 2015-03-31 | 2015-03-31 | 辐射能传播蒙特卡洛算法中的局部网格及光子加密方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104699993B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108478192B (zh) * | 2018-04-11 | 2021-01-15 | 西安交通大学 | 一种用于估测皮肤组织微血管深度的测量系统 |
| CN113361032B (zh) * | 2021-06-03 | 2023-11-10 | 西安交通大学 | 一种考虑流场时空演化的网格无关性分析方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6148272A (en) * | 1998-11-12 | 2000-11-14 | The Regents Of The University Of California | System and method for radiation dose calculation within sub-volumes of a monte carlo based particle transport grid |
| CN101739503A (zh) * | 2008-11-12 | 2010-06-16 | 中国科学院自动化研究所 | 一种预测光在生物组织中传播模型的实现方法 |
| CN103295188A (zh) * | 2012-02-28 | 2013-09-11 | 上海联影医疗科技有限公司 | X射线蒙特卡洛模拟的路径积分方法 |
| CN103356170A (zh) * | 2013-05-24 | 2013-10-23 | 天津大学 | 用于带异质体组织光学参数重建的快速蒙特卡罗成像方法 |
-
2015
- 2015-03-31 CN CN201510149785.3A patent/CN104699993B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6148272A (en) * | 1998-11-12 | 2000-11-14 | The Regents Of The University Of California | System and method for radiation dose calculation within sub-volumes of a monte carlo based particle transport grid |
| CN101739503A (zh) * | 2008-11-12 | 2010-06-16 | 中国科学院自动化研究所 | 一种预测光在生物组织中传播模型的实现方法 |
| CN103295188A (zh) * | 2012-02-28 | 2013-09-11 | 上海联影医疗科技有限公司 | X射线蒙特卡洛模拟的路径积分方法 |
| CN103356170A (zh) * | 2013-05-24 | 2013-10-23 | 天津大学 | 用于带异质体组织光学参数重建的快速蒙特卡罗成像方法 |
Non-Patent Citations (3)
| Title |
|---|
| A tetrahedron-based inhomogeneous Monte Carlo optical simulator;H shen et al.;《Phys Med Biol》;20100422;第947-962页 * |
| Influence of voxel size on specific absorbed fractions and S-values in a mouse vexel phantom;A.Mohammadi et al.;《Radiation Protection Dosimetry》;20101230;第258-263页 * |
| 一种模拟生物组织内光传播的三维几何蒙特卡洛方法;张永 等;《中国激光》;20150130;第42卷(第1期);第0104003-1-0104003-8页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104699993A (zh) | 2015-06-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Vassiliev et al. | Validation of a new grid-based Boltzmann equation solver for dose calculation in radiotherapy with photon beams | |
| CA2557888C (en) | Computation of radiating particle and wave distributions using a generalized discrete field constructed from representative ray sets | |
| Jia et al. | GPU-based fast Monte Carlo simulation for radiotherapy dose calculation | |
| Popple et al. | Comprehensive evaluation of a commercial macro Monte Carlo electron dose calculation implementation using a standard verification data set | |
| Wilson et al. | State-of-the-art 3-D radiation transport methods for fusion energy systems | |
| Martínez‐Rovira et al. | Monte Carlo‐based treatment planning system calculation engine for microbeam radiation therapy | |
| CN105740573B (zh) | 一种用于放射线剂量计算的双步蒙特卡洛模拟方法 | |
| Jiang et al. | Implicit, pressure-based incompressible Navier-Stokes equations solver for unstructured meshes | |
| Hauck et al. | A collision-based hybrid method for time-dependent, linear, kinetic transport equations | |
| CN104699993B (zh) | 辐射能传播蒙特卡洛算法中的局部网格及光子加密方法 | |
| Chen et al. | Preliminary study on CAD-based method of characteristics for neutron transport calculation | |
| TW201825050A (zh) | 基於醫學影像的組織元素質量比例解構方法及幾何模型建立方法 | |
| Fonseca et al. | A medical image-based graphical platform—Features, applications and relevance for brachytherapy | |
| Zhang et al. | Comprehensive evaluation and clinical implementation of commercially available Monte Carlo dose calculation algorithm | |
| Gifford et al. | Optimization of deterministic transport parameters for the calculation of the dose distribution around a high dose‐rate brachytherapy source | |
| CN106202613B (zh) | 一种适用于耦合输运计算的大规模并行方法 | |
| Jabbari et al. | A fast Monte Carlo code for proton transport in radiation therapy based on MCNPX | |
| Luu et al. | FPGA-based Monte Carlo computation of light absorption for photodynamic cancer therapy | |
| Yang et al. | A gamma-ray dose rate assessment method for arbitrary shape geometries based on voxelization algorithm | |
| Takada et al. | Development of Monte Carlo based real-time treatment planning system with fast calculation algorithm for boron neutron capture therapy | |
| Berthon et al. | Numerical methods for balance laws with space dependent flux: application to radiotherapy dose calculation | |
| Herty et al. | On optimal treatment planning in radiotherapy governed by transport equations | |
| Bohm et al. | The Impact of Simplifications on 3-D Neutronics Analysis of Blanket Modules in ITER | |
| Yan et al. | A dual-mesh Monte Carlo algorithm using a coarse tetrahedral mesh and voxel output | |
| Fan et al. | An Improved Convolution Algorithm of Cylinder String Based on Unified Coordinate System |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210519 Address after: 200333 room 5091, 5th floor, No.5 Lane 600, Yunling West Road, Putuo District, Shanghai Patentee after: Shanghai Jiding Information Technology Co.,Ltd. Address before: 710049 No. 28, Xianning Road, Xi'an, Shaanxi Patentee before: XI'AN JIAOTONG University |