The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female populat...The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female population to replace the traditional standard that utilizes simplified breast models,for the accurate estimation of the mean glandular dose of a patient undergoing digital breast tomosynthesis(DBT).The first set of detailed Chinese female breast models and representative breast parameters was constructed.Considering backscatter radiation and computational efficiency,we improved the combination of these models and the Chinese reference adult female whole-body voxel phantom.Image acquisition for four commercial DBT systems that are widely employed in China was simulated using the Monte Carlo method to obtain the normalized glandular dose coefficients of DBT(D_(gN)^(DBT))and the glandular depth dose(D_(g)^(dep)(z))for different breast characteristics and X-ray spectra.We calculated a series of D_(gN)^(DBT) values for breasts with different percentage mass glandularities(5%,25%,50%,75%,and 100%)and compressed breast thicknesses(2,3,4,5,6,and 7 cm)at various tube potentials(25,28,30,32,35,and 49 kV)and target/filter combinations(W/Rh,W/Al,Mo/Mo,Rh/Rh,and Rh/Ag).The parameter dependence of the breast characteristics and beam conditions on D_(gN)^(DBT) in detailed breast models was investigated.The D_(gN)^(DBT) results were 14.6-51.0%lower than those of the traditional dosimetry standard in China.The difference in D_(gN)^(DBT) was mainly due to a decrease in the depth of the main energy deposition area caused by the glandular distribution along the depth direction.The results obtained in this study may be used to improve breast dosimetry in China and provide more detailed information on risk assessment during DBT.展开更多
High-purity germanium(HPGe)detectors,which are used for direct dark matter detection,have the advantages of a low threshold and excellent energy resolution.The surface passivation of HPGe has become crucial for achiev...High-purity germanium(HPGe)detectors,which are used for direct dark matter detection,have the advantages of a low threshold and excellent energy resolution.The surface passivation of HPGe has become crucial for achieving an extremely low energy threshold.In this study,first-principles simulations,passivation film preparation,and metal oxide semiconductor(MOS)capacitor characterization were combined to study surface passivation.Theoretical calculations of the energy band structure of the -H,-OH,and -NH_(2) passivation groups on the surface of Ge were performed,and the interface state density and potential with five different passivation groups with N/O atomic ratios were accurately analyzed to obtain a stable surface state.Based on the theoretical calculation results,the surface passivation layers of the Ge_(2)ON_(2) film were prepared via magnetron sputtering in accordance with the optimum atomic ratio structure.The microstructure,C-V,and I-V electrical properties of the layers,and the passivation effect of the Al/Ge_(2)ON_(2)/Ge MOS were characterized to test the interface state density.The mean interface state density obtained by the Terman method was 8.4×10^(11) cm^(-2) eV^(-1).The processing of germanium oxynitrogen passivation films is expected to be used in direct dark matter detection of the HPGe detector surface passivation technology to reduce the detector leakage currents.展开更多
GPU-based Monte Carlo(MC)simulations are highly valued for their potential to improve both the computational efficiency and accuracy of radiotherapy.However,in proton therapy,these methods often simplify human tissues...GPU-based Monte Carlo(MC)simulations are highly valued for their potential to improve both the computational efficiency and accuracy of radiotherapy.However,in proton therapy,these methods often simplify human tissues as water for nuclear reactions,disregarding their true elemental composition and thereby potentially compromising calculation accuracy.Consequently,this study developed the program g MCAP(GPU-based proton MC Algorithm for Proton therapy),incorporating precise discrete interactions,and established a refined nuclear reaction model(REFINED)that considers the actual materials of the human body.Compared to the approximate water model(APPROX),the REFINED model demonstrated an improvement in calculation accuracy of 3%.In particular,in high-density tissue regions,the maximum dose deviation between the REFINED and APPROX models was up to 15%.In summary,the g MCAP program can efficiently simulate 1 million protons within 1 s while significantly enhancing dose calculation accuracy in high-density tissues,thus providing a more precise and efficient engine for proton radiotherapy dose calculations in clinical practice.展开更多
A Compton camera prototype has been developed using a pixelated CZT detector with 4-by-4 pixels.Signals of the detector are read out by a VASTAT ASIC that is controlled by a self-developed DAQ board. The DAQ software ...A Compton camera prototype has been developed using a pixelated CZT detector with 4-by-4 pixels.Signals of the detector are read out by a VASTAT ASIC that is controlled by a self-developed DAQ board. The DAQ software is developed using LabVIEW, and the offline Compton imaging codes are written in C++. The prototype has been successfully calibrated, and its capabilities for source detection, spectroscopy, and Compton imaging have been demonstrated using a Cs-137 source.The angular resolution of the 662 keV line is 36° FWHM for the simple back-projection method and 9.6° FWHM for the MLEM reconstruction method. The system is ready to be extended to 11-by-11 pixels in the future, and a better imaging quality can be expected due to the better relative position resolution.展开更多
Automatic conversion from a computer-aided design(CAD) model to Monte Carlo geometry is one of the most effective methods for large-scale and detailed Monte Carlo modeling. The CAD to Monte Carlo geometry converter(CM...Automatic conversion from a computer-aided design(CAD) model to Monte Carlo geometry is one of the most effective methods for large-scale and detailed Monte Carlo modeling. The CAD to Monte Carlo geometry converter(CMGC) is a newly developed conversion code based on the boundary representation to constructive solid geometry(BRep→CSG) conversion method. The goal of the conversion process in the CMGC is to generate an appropriate CSG representation to achieve highly efficient Monte Carlo simulations. We designed a complete solid decomposition scheme to split a complex solid into as few nonoverlapping simple sub-solids as possible. In the complete solid decomposition scheme, the complex solid is successively split by so-called direct, indirect, and auxiliary splitting surfaces. We defined the splitting edge and designed a method for determining the direct splitting surface based on the splitting edge, then provided a method for determining indirect and auxiliary splitting surfaces based on solid vertices. Only the sub-solids that contain concave boundary faces need to be supplemented with auxiliary surfaces because the solid is completely decomposed, which will reduce the redundancy in the CSG expression. After decomposition, these sub-solids are located on only one side of their natural and auxiliary surfaces;thus, each sub-solid can be described by the intersections of a series of half-spaces or geometrical primitives. The CMGC has a friendly graphical user interface and can convert a CAD model into geometry input files for several Monte Carlo codes. The reliability of the CMGC was evaluated by converting several complex models and calculating the relative volume errors. Moreover, JMCT was used to test the efficiency of the Monte Carlo simulation. The results showed that the converted models performed well in particle transport calculations.展开更多
Polytetrafluoroethylene(PTFE)is a low-background polymer that is applied to several applications in rare-event detection and underground low-background experiments.PTFE-based electronic substrates are important for re...Polytetrafluoroethylene(PTFE)is a low-background polymer that is applied to several applications in rare-event detection and underground low-background experiments.PTFE-based electronic substrates are important for reducing the detection limit of high-purity germanium detectors and scintillator calorimeters,which are widely applied in dark matter and 0υββdetection experiments.The traditional adhesive bonding method between PTFE and copper is not conducive to working in liquid nitrogen and extremely low-temperature environments.To avoid adhesive bonding,PTFE must be processed for surface metallization owing to the mismatch between the PTFE and copper conductive layer.Low-background PTFE matrix composites(m-PTFE)were selected to improve the electrical and mechanical properties of PTFE by introducing SiO_(2)/TiO_(2) particles.The microstructures,surface elements,and electrical properties of PTFE and m-PTFE were characterized and analyzed following ion implantation.PTFE and m-PTFE surfaces were found to be broken,degraded,and cross-linked by ion implantation,resulting in C=C conjugated double bonds,increased surface energy,and increased surface roughness.Comparably,the surface roughness,bond strength,and conjugated double bonds of m-PTFE were significantly more intense than those of PTFE.Moreover,the interface bonding theory between PTFE and the metal copper foil was analyzed using the direct metallization principle.Therefore,the peel strength of the optimized electronic substrates was higher than that of the industrial standard at extremely low temperatures,while maintaining excellent electrical properties.展开更多
We present the design, fabrication, and characterization of two new types of terahertz photoconductive emitters. One has an asymmetric four-contact electrode structure and the other has an arc-shaped electrode structu...We present the design, fabrication, and characterization of two new types of terahertz photoconductive emitters. One has an asymmetric four-contact electrode structure and the other has an arc-shaped electrode structure, which are all modified from a traditional strip line antenna. Numerical simulations and real experiments confirm the good performance of the proposed antennas. An amplitude increase of about 40% is experimentally observed for the terahertz signals generated from the new structures. The special electrode structure and its induced local bias field enhancement are responsible for this radiation efficiency improvement. Our work demonstrates the feasibility of developing highly efficient terahertz photoconductive emitters by optimizing the electrode structure.展开更多
A four-vane radio-frequency quadrupole(RFQ)accelerator is under construction for the Compact Pulsed Hadron Source(CPHS)project at Tsinghua University.The 3 m-long RFQ will accelerate a 50 keV proton beam from the ECR ...A four-vane radio-frequency quadrupole(RFQ)accelerator is under construction for the Compact Pulsed Hadron Source(CPHS)project at Tsinghua University.The 3 m-long RFQ will accelerate a 50 keV proton beam from the ECR source to 3 MeV,and deliver it to the downstream drift tube linac(DTL)with a peak current of 50 mA,pulse length of 0.5 ms and beam duty factor of 2.5%.The inter-vane voltage is designed to increase with the longitudinal position to produce a short RFQ.Coupling plates are therefore not necessary.The cavity cross section and vane-tip geometry are tailored as a function of the longitudinal position,while limiting the peak surface electric field to 1.8 Kilpatrick.The RFQ is designed,manufactured,and installed at Tsinghua University.We also present the tuning and cold test results of the RFQ accelerator.After final tuning,the relative error of the quadrupole field is within 2%,and the admixture of the two dipole modes are less than 2%of the quadrupole mode.展开更多
This study proposes a source distribution inversion convolutional neural network (SDICNN), which is deep neural network model for the inversion of complex source distributions, to solve inversion problems involving fi...This study proposes a source distribution inversion convolutional neural network (SDICNN), which is deep neural network model for the inversion of complex source distributions, to solve inversion problems involving fixed-source distributions. A function is developed to obtain the distribution information of complex source terms from radiation parameters at individual sampling points in space. The SDICNN comprises two components:a fully connected network and a convolutional neural network. The fully connected network mainly extracts the parameter measurement information from the sampling points,whereas the convolutional neural network mainly completes the fine inversion of the source-term distribution. Finally, the SDICNN obtains a high-resolution source-term distribution image. In this study, the proposed source-term inversion method is evaluated based on typical geometric scenarios. The results show that, unlike the conventional fully connected neural network, the SDICNN model can extract the two-dimensional distribution features of the source terms, and its inversion results are better. In addition, the effects of the shielding mechanism and number of sampling points on the inversion process are examined. In summary, the result of this study can facilitate the accurate assessment of dose distributions in nuclear facilities.展开更多
In most coherent imaging modality,speckle noise is a major cause that blurs the boundary of tissues and degrades the image contrast.By utilizing the unique properties of supercontinuum(SC)generated by noise-like pulse...In most coherent imaging modality,speckle noise is a major cause that blurs the boundary of tissues and degrades the image contrast.By utilizing the unique properties of supercontinuum(SC)generated by noise-like pulses(NLPs)and a simple multi-frame averaging technique,we achieved significant speckle reduction in spectral domain optical coherence tomography(SD-OCT).We quantitatively compared the speckle of our proposed method with those of conventional swept source OCT(SS-OCT)and SD-OCT based on commercial light sources.The experimental results show that SC pumped by NLPs combined with noncoherent averaging method achieves better denoising performance in terms of contrast to noise ratio(CNR).展开更多
Global variance reduction is a bottleneck in Monte Carlo shielding calculations.The global variance reduction problem requires that the statistical error of the entire space is uniform.This study proposed a grid-AIS m...Global variance reduction is a bottleneck in Monte Carlo shielding calculations.The global variance reduction problem requires that the statistical error of the entire space is uniform.This study proposed a grid-AIS method for the global variance reduction problem based on the AIS method,which was implemented in the Monte Carlo program MCShield.The proposed method was validated using the VENUS-Ⅲ international benchmark problem and a self-shielding calculation example.The results from the VENUS-Ⅲ benchmark problem showed that the grid-AIS method achieved a significant reduction in the variance of the statistical errors of the MESH grids,decreasing from 1.08×10^(-2) to 3.84×10^(-3),representing a 64.00% reduction.This demonstrates that the grid-AIS method is effective in addressing global issues.The results of the selfshielding calculation demonstrate that the grid-AIS method produced accurate computational results.Moreover,the grid-AIS method exhibited a computational efficiency approximately one order of magnitude higher than that of the AIS method and approximately two orders of magnitude higher than that of the conventional Monte Carlo method.展开更多
基金supported by the National Natural Science Foundation of China(Nos.U2167209 and 12175114)the National Key R&D Program of China(No.2021YFF0603600).
文摘The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female population to replace the traditional standard that utilizes simplified breast models,for the accurate estimation of the mean glandular dose of a patient undergoing digital breast tomosynthesis(DBT).The first set of detailed Chinese female breast models and representative breast parameters was constructed.Considering backscatter radiation and computational efficiency,we improved the combination of these models and the Chinese reference adult female whole-body voxel phantom.Image acquisition for four commercial DBT systems that are widely employed in China was simulated using the Monte Carlo method to obtain the normalized glandular dose coefficients of DBT(D_(gN)^(DBT))and the glandular depth dose(D_(g)^(dep)(z))for different breast characteristics and X-ray spectra.We calculated a series of D_(gN)^(DBT) values for breasts with different percentage mass glandularities(5%,25%,50%,75%,and 100%)and compressed breast thicknesses(2,3,4,5,6,and 7 cm)at various tube potentials(25,28,30,32,35,and 49 kV)and target/filter combinations(W/Rh,W/Al,Mo/Mo,Rh/Rh,and Rh/Ag).The parameter dependence of the breast characteristics and beam conditions on D_(gN)^(DBT) in detailed breast models was investigated.The D_(gN)^(DBT) results were 14.6-51.0%lower than those of the traditional dosimetry standard in China.The difference in D_(gN)^(DBT) was mainly due to a decrease in the depth of the main energy deposition area caused by the glandular distribution along the depth direction.The results obtained in this study may be used to improve breast dosimetry in China and provide more detailed information on risk assessment during DBT.
基金supported by the National Natural Science Foundation of China(No.12005017).
文摘High-purity germanium(HPGe)detectors,which are used for direct dark matter detection,have the advantages of a low threshold and excellent energy resolution.The surface passivation of HPGe has become crucial for achieving an extremely low energy threshold.In this study,first-principles simulations,passivation film preparation,and metal oxide semiconductor(MOS)capacitor characterization were combined to study surface passivation.Theoretical calculations of the energy band structure of the -H,-OH,and -NH_(2) passivation groups on the surface of Ge were performed,and the interface state density and potential with five different passivation groups with N/O atomic ratios were accurately analyzed to obtain a stable surface state.Based on the theoretical calculation results,the surface passivation layers of the Ge_(2)ON_(2) film were prepared via magnetron sputtering in accordance with the optimum atomic ratio structure.The microstructure,C-V,and I-V electrical properties of the layers,and the passivation effect of the Al/Ge_(2)ON_(2)/Ge MOS were characterized to test the interface state density.The mean interface state density obtained by the Terman method was 8.4×10^(11) cm^(-2) eV^(-1).The processing of germanium oxynitrogen passivation films is expected to be used in direct dark matter detection of the HPGe detector surface passivation technology to reduce the detector leakage currents.
文摘GPU-based Monte Carlo(MC)simulations are highly valued for their potential to improve both the computational efficiency and accuracy of radiotherapy.However,in proton therapy,these methods often simplify human tissues as water for nuclear reactions,disregarding their true elemental composition and thereby potentially compromising calculation accuracy.Consequently,this study developed the program g MCAP(GPU-based proton MC Algorithm for Proton therapy),incorporating precise discrete interactions,and established a refined nuclear reaction model(REFINED)that considers the actual materials of the human body.Compared to the approximate water model(APPROX),the REFINED model demonstrated an improvement in calculation accuracy of 3%.In particular,in high-density tissue regions,the maximum dose deviation between the REFINED and APPROX models was up to 15%.In summary,the g MCAP program can efficiently simulate 1 million protons within 1 s while significantly enhancing dose calculation accuracy in high-density tissues,thus providing a more precise and efficient engine for proton radiotherapy dose calculations in clinical practice.
文摘A Compton camera prototype has been developed using a pixelated CZT detector with 4-by-4 pixels.Signals of the detector are read out by a VASTAT ASIC that is controlled by a self-developed DAQ board. The DAQ software is developed using LabVIEW, and the offline Compton imaging codes are written in C++. The prototype has been successfully calibrated, and its capabilities for source detection, spectroscopy, and Compton imaging have been demonstrated using a Cs-137 source.The angular resolution of the 662 keV line is 36° FWHM for the simple back-projection method and 9.6° FWHM for the MLEM reconstruction method. The system is ready to be extended to 11-by-11 pixels in the future, and a better imaging quality can be expected due to the better relative position resolution.
基金Supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China （No. 2013ZX06002001- 007）, the National Key Scientific Instrument and Equipment Development Projects, China （No. 2012YQ180118） and the National Natural Science Foundation of China （Nos. 11275110, 11075091 and 11105081）.
基金the National Natural Science Foundation of China(No.11805017)。
文摘Automatic conversion from a computer-aided design(CAD) model to Monte Carlo geometry is one of the most effective methods for large-scale and detailed Monte Carlo modeling. The CAD to Monte Carlo geometry converter(CMGC) is a newly developed conversion code based on the boundary representation to constructive solid geometry(BRep→CSG) conversion method. The goal of the conversion process in the CMGC is to generate an appropriate CSG representation to achieve highly efficient Monte Carlo simulations. We designed a complete solid decomposition scheme to split a complex solid into as few nonoverlapping simple sub-solids as possible. In the complete solid decomposition scheme, the complex solid is successively split by so-called direct, indirect, and auxiliary splitting surfaces. We defined the splitting edge and designed a method for determining the direct splitting surface based on the splitting edge, then provided a method for determining indirect and auxiliary splitting surfaces based on solid vertices. Only the sub-solids that contain concave boundary faces need to be supplemented with auxiliary surfaces because the solid is completely decomposed, which will reduce the redundancy in the CSG expression. After decomposition, these sub-solids are located on only one side of their natural and auxiliary surfaces;thus, each sub-solid can be described by the intersections of a series of half-spaces or geometrical primitives. The CMGC has a friendly graphical user interface and can convert a CAD model into geometry input files for several Monte Carlo codes. The reliability of the CMGC was evaluated by converting several complex models and calculating the relative volume errors. Moreover, JMCT was used to test the efficiency of the Monte Carlo simulation. The results showed that the converted models performed well in particle transport calculations.
基金supported by the National Natural Science Foundation of China(Nos.12141502 and 12005017).
文摘Polytetrafluoroethylene(PTFE)is a low-background polymer that is applied to several applications in rare-event detection and underground low-background experiments.PTFE-based electronic substrates are important for reducing the detection limit of high-purity germanium detectors and scintillator calorimeters,which are widely applied in dark matter and 0υββdetection experiments.The traditional adhesive bonding method between PTFE and copper is not conducive to working in liquid nitrogen and extremely low-temperature environments.To avoid adhesive bonding,PTFE must be processed for surface metallization owing to the mismatch between the PTFE and copper conductive layer.Low-background PTFE matrix composites(m-PTFE)were selected to improve the electrical and mechanical properties of PTFE by introducing SiO_(2)/TiO_(2) particles.The microstructures,surface elements,and electrical properties of PTFE and m-PTFE were characterized and analyzed following ion implantation.PTFE and m-PTFE surfaces were found to be broken,degraded,and cross-linked by ion implantation,resulting in C=C conjugated double bonds,increased surface energy,and increased surface roughness.Comparably,the surface roughness,bond strength,and conjugated double bonds of m-PTFE were significantly more intense than those of PTFE.Moreover,the interface bonding theory between PTFE and the metal copper foil was analyzed using the direct metallization principle.Therefore,the peel strength of the optimized electronic substrates was higher than that of the industrial standard at extremely low temperatures,while maintaining excellent electrical properties.
基金supported by the National Science and Technology Support Program of China under Grant No.2013BAK14B03
文摘We present the design, fabrication, and characterization of two new types of terahertz photoconductive emitters. One has an asymmetric four-contact electrode structure and the other has an arc-shaped electrode structure, which are all modified from a traditional strip line antenna. Numerical simulations and real experiments confirm the good performance of the proposed antennas. An amplitude increase of about 40% is experimentally observed for the terahertz signals generated from the new structures. The special electrode structure and its induced local bias field enhancement are responsible for this radiation efficiency improvement. Our work demonstrates the feasibility of developing highly efficient terahertz photoconductive emitters by optimizing the electrode structure.
基金the National Natural Science Foundation of China under Grant Nos 91126003 and 11175096.
文摘A four-vane radio-frequency quadrupole(RFQ)accelerator is under construction for the Compact Pulsed Hadron Source(CPHS)project at Tsinghua University.The 3 m-long RFQ will accelerate a 50 keV proton beam from the ECR source to 3 MeV,and deliver it to the downstream drift tube linac(DTL)with a peak current of 50 mA,pulse length of 0.5 ms and beam duty factor of 2.5%.The inter-vane voltage is designed to increase with the longitudinal position to produce a short RFQ.Coupling plates are therefore not necessary.The cavity cross section and vane-tip geometry are tailored as a function of the longitudinal position,while limiting the peak surface electric field to 1.8 Kilpatrick.The RFQ is designed,manufactured,and installed at Tsinghua University.We also present the tuning and cold test results of the RFQ accelerator.After final tuning,the relative error of the quadrupole field is within 2%,and the admixture of the two dipole modes are less than 2%of the quadrupole mode.
基金supported by the Platform Development Foundation of the China Institute for Radiation Protection (No. YP21030101)the National Natural Science Foundation of China (General Program)(Nos. 12175114, U2167209)+1 种基金the National Key R&D Program of China (No. 2021YFF0603600)the Tsinghua University Initiative Scientific Research Program (No. 20211080081)。
文摘This study proposes a source distribution inversion convolutional neural network (SDICNN), which is deep neural network model for the inversion of complex source distributions, to solve inversion problems involving fixed-source distributions. A function is developed to obtain the distribution information of complex source terms from radiation parameters at individual sampling points in space. The SDICNN comprises two components:a fully connected network and a convolutional neural network. The fully connected network mainly extracts the parameter measurement information from the sampling points,whereas the convolutional neural network mainly completes the fine inversion of the source-term distribution. Finally, the SDICNN obtains a high-resolution source-term distribution image. In this study, the proposed source-term inversion method is evaluated based on typical geometric scenarios. The results show that, unlike the conventional fully connected neural network, the SDICNN model can extract the two-dimensional distribution features of the source terms, and its inversion results are better. In addition, the effects of the shielding mechanism and number of sampling points on the inversion process are examined. In summary, the result of this study can facilitate the accurate assessment of dose distributions in nuclear facilities.
基金This research was supported by the National Natural Science Foundation of China(61227807,61575108 and 61505034)Tsinghua Initiative Scientific Research Program(2013THZ02-3)C.L.Pan and Y.J.You were supported in part by the Ministry of Science and Technology,Taiwan,under Grant 103-2622-E-007-006-CC2 and by the National Tsing Hua University Research Program Grant 104N2711E1.C.L.Pan is also supported by the Air Force Office of Scientific Research FA2386-13-1-4086.The first three authors contribute equally to this paper。
文摘In most coherent imaging modality,speckle noise is a major cause that blurs the boundary of tissues and degrades the image contrast.By utilizing the unique properties of supercontinuum(SC)generated by noise-like pulses(NLPs)and a simple multi-frame averaging technique,we achieved significant speckle reduction in spectral domain optical coherence tomography(SD-OCT).We quantitatively compared the speckle of our proposed method with those of conventional swept source OCT(SS-OCT)and SD-OCT based on commercial light sources.The experimental results show that SC pumped by NLPs combined with noncoherent averaging method achieves better denoising performance in terms of contrast to noise ratio(CNR).
基金supported by the Platform Development Foundation of the China Institute for Radiation Protection(No.YP21030101)the National Natural Science Foundation of China(General Program)(Nos.12175114,U2167209)+1 种基金the National Key R&D Program of China(No.2021YFF0603600)the Tsinghua University Initiative Scientific Research Program(No.20211080081).
文摘Global variance reduction is a bottleneck in Monte Carlo shielding calculations.The global variance reduction problem requires that the statistical error of the entire space is uniform.This study proposed a grid-AIS method for the global variance reduction problem based on the AIS method,which was implemented in the Monte Carlo program MCShield.The proposed method was validated using the VENUS-Ⅲ international benchmark problem and a self-shielding calculation example.The results from the VENUS-Ⅲ benchmark problem showed that the grid-AIS method achieved a significant reduction in the variance of the statistical errors of the MESH grids,decreasing from 1.08×10^(-2) to 3.84×10^(-3),representing a 64.00% reduction.This demonstrates that the grid-AIS method is effective in addressing global issues.The results of the selfshielding calculation demonstrate that the grid-AIS method produced accurate computational results.Moreover,the grid-AIS method exhibited a computational efficiency approximately one order of magnitude higher than that of the AIS method and approximately two orders of magnitude higher than that of the conventional Monte Carlo method.
