Cooperation among enterprises can bring overall and individual performance improvement,and a smooth coordination method is indispensable.However,due to the lack of customized coordination methods,cooperation in the do...Cooperation among enterprises can bring overall and individual performance improvement,and a smooth coordination method is indispensable.However,due to the lack of customized coordination methods,cooperation in the downstream oil supply chain cannot be carried out smoothly.This paper intends to propose a multi-party coordination method to promote cooperation between oil shippers and pipeline operator by optimizing oil transportation,oil substitution and pipeline pricing schemes.An integrated game-theoretic modeling and analysis approach is developed to characterize the operation behaviors of all stakeholders in the downstream oil supply chain.The proposed mixed integer nonlinear programming model constrains supply and demand capacity,transportation routes,oil substitution rules and pipeline freight levels.Logarithm transformation and price discretization are introduced for model linear approximation.Simulation experiments are carried out in the oil distribution system in South China.The results show that compared to the business-as-usual scheme,the new scheme saves transportation cost by 3.48%,increases pipeline turnover by 5.7%,and reduces energy consumption and emissions by 7.66%and 6.77%.It is proved that the proposed method improves the revenue of the whole system,achieves fair revenue distribution,and also improves the energy and environmental benefits of the oil supply chain.展开更多
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.展开更多
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.展开更多
Monte Carlo simulations are frequently utilized in radiation dose assessments.However,many researchers find the prevailing computing platforms to be intricate.This highlights a pressing need for a specialized framewor...Monte Carlo simulations are frequently utilized in radiation dose assessments.However,many researchers find the prevailing computing platforms to be intricate.This highlights a pressing need for a specialized framework for phantom dose evalua-tion.To address this gap,we developed a user-friendly radiation dose assessment platform using the Monte Carlo toolkit,Geant4.The Tsinghua University Phantom Dose(THUDosePD)augments the flexibility of Monte Carlo simulations in dosi-metric research.Originating from THUDose,a code with generic,functional,and application layers,THUDosePD focuses predominantly on anatomical phantom dose assessment.Additionally,it enables medical exposure simulation,intricate geometry creation,and supports both three-dimensional radiation dose analysis and phantom format transformations.The system operates on a multi-threaded parallel CPU architecture,with some modules enhanced for GPU parallel computing.Benchmark tests on the ICRP reference male illustrated the capabilities of THUDosePD in phantom dose assessment,covering the effective dose,three-dimensional dose distribution,and three-dimensional organ dose.We also conducted a voxelization conversion on the polygon mesh phantom,demonstrating the method’s efficiency and consistency.Extended applications based on THUDosePD further underline its broad adaptability.This intuitive,three-dimensional platform stands out as a valuable tool for phantom radiation dosimetry research.展开更多
Prediction models were proposed to estimate the reduced Townsend ionization coefficient and ionization cross-section.A shape function of the reduced Townsend ionization coefficient curves was derived from the ionizati...Prediction models were proposed to estimate the reduced Townsend ionization coefficient and ionization cross-section.A shape function of the reduced Townsend ionization coefficient curves was derived from the ionization collision probability model.The function had three parameters:the first ionization potential energy,A_(α),and B_(α).A_(α)and B_(α)were related to the molecule symmetry and size.The polarization of molecules could characterize the molecule symmetry.The multi-layer molecular cross-section(MMCS)was proposed to describe the contributions of electrons and molecule radius on different molecule surfaces to collisions.A prediction model of the ionization cross-section was also proposed based on Aα.The molecule parameters were calculated by the Becke3–Lee–Yang–Parr(B3LYP)method and the 6–311G**basis set.We used available data of 30 and 23 gases,respectively,to build the prediction models of reduced Townsend ionization coefficients and ionization cross-sections.The relationships between the molecular parameters Aαand Bαand the ionization cross-section were built up via nonlinear fittings.The determination coefficients R^(2)of Aα,Bα,and the ionization cross-section were 0.877,0.887,and 0.838,respectively.The results showed that the accuracy of models was positively correlated with the molecule symmetry and reduced electric field.This was mainly related to the accuracy of the MMCS model in predicting Aα.The MMCS model needed to be improved to describe the collision direction selectivity caused by the molecule asymmetry.Under a high reduced electric field,that error of Aαhad less influence on the prediction results.However,the prediction results for single atoms with high symmetry were poor.This may be due to the absolute error of the model close to single atoms’reduced Townsend ionization coefficients.The models could provide the basis for gas insulation prediction and discharge calculations,especially for symmetric molecules under a high electric field.展开更多
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.展开更多
This paper intends to complete the primary logistics planning of oil products under the imbalance of supply and demand. An integrated mathematical programming model is developed to simultaneously find the balance betw...This paper intends to complete the primary logistics planning of oil products under the imbalance of supply and demand. An integrated mathematical programming model is developed to simultaneously find the balance between supply and demand, and optimize the logistics scheme. The model takes minimum logistics cost and resource adjustment cost as the objective function, and takes supply and demand capacity, transportation capacity, mass balance, and resource adjustment rules as constraints.Three adjustment rules are considered in the model, including resource adjustment within oil suppliers,within oil consumers, and between oil consumers. The model is tested on a large-scale primary logistics of a state-owned petroleum enterprise, involving 37 affiliated refineries, 31 procurement departments,286 market depots and dedicated consumers. After the unified optimization, the supply and demand imbalance is eased by 97% and the total cost is saved by 7%, which proves the effectiveness and applicability of the proposed model.展开更多
The broad-energy germanium(BEGe)detector,with the ability of background discrimination using pulse shape discrimination,is a competitive candidate for neutrinoless double beta decay(ovββ)experiments.In this paper,we...The broad-energy germanium(BEGe)detector,with the ability of background discrimination using pulse shape discrimination,is a competitive candidate for neutrinoless double beta decay(ovββ)experiments.In this paper,we report our measurements of key parameters for detector modeling in a commercial p-type BEGe detector.Point-like sources are used to investigate energy resolution and linearity of the detector.A cylindrical volume source is used for efficiency calibration.With an assembled device for source positioning and a collimated ^(133)Ba source,the detector is scanned to check its active volume.Using an^(241)Am point-like source,the dead layer thicknesses is measured at about 0.17 mm on the front and 1.18 mm on the side.The detector characterization is of importance for BEGe detectors to be used in the ovββ experiments at China JinPing underground Laboratory(CJPL).展开更多
Hydrodesulfurization(HDS)of sour crude oil is an effective way to address the corrosion problems in refineries,and is an economic way to process sour crude oil in an existing refinery built for sweet oil.In the curren...Hydrodesulfurization(HDS)of sour crude oil is an effective way to address the corrosion problems in refineries,and is an economic way to process sour crude oil in an existing refinery built for sweet oil.In the current study,the HDS of Siberian crude oil was carried out in a slurry reactor.The Co-Mo,Ni-Mo,and Ni-W catalysts supported onγ-Al2O3 were compared at the temperature of 340℃ and the pressure of 4.5 MPa.The HDS activity follows the order of Co-Mo N Ni-Mo N Ni-Wat a high concentration of H2S,and the difference between Co-Mo and Ni-Mo becomes insignificant at a low concentration of H2S.The influence of reaction temperature 320-360℃ and reaction pressure 3-5.5 MPa was investigated,and both play a positive role in the HDS reaction.A kinetic model over Ni-Mo/Al2O3 in the slurry reactor was established.The activation energy is estimated as 60.34 kJ·mol−1;the orders of sulfur components and hydrogen partial pressure are 1.43 and 1.30,respectively.The kinetic parameters are compared with those in a trickle-bed reactor,implying that the mass transfer is greatly enhanced in the slurry reactor.The back mixing effect is present in the slurry reactor and can be reduced by a multi-stage design,which would lead to higher reactor efficiency in industrial application.展开更多
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.展开更多
The Monte Carlo(MC)simulation is regarded as the gold standard for dose calculation in brachytherapy,but it consumes a large amount of computing resources.The development of heterogeneous computing makes it possible t...The Monte Carlo(MC)simulation is regarded as the gold standard for dose calculation in brachytherapy,but it consumes a large amount of computing resources.The development of heterogeneous computing makes it possible to substantially accelerate calculations with hardware accelerators.Accordingly,this study develops a fast MC tool,called THUBrachy,which can be accelerated by several types of hardware accelerators.THUBrachy can simulate photons with energy less than 3 MeV and considers all photon interactions in the energy range.It was benchmarked against the American Association of Physicists in Medicine Task Group No.43 Report using a water phantom and validated with Geant4 using a clinical case.A performance test was conducted using the clinical case,showing that a multicore central processing unit,Intel Xeon Phi,and graphics processing unit(GPU)can efficiently accelerate the simulation.GPU-accelerated THUBrachy is the fastest version,which is 200 times faster than the serial version and approximately 500 times faster than Geant4.The proposed tool shows great potential for fast and accurate dose calculations in clinical applications.展开更多
基金partially supported by the Science Foundation of China University of Petroleum,Beijing(2462023XKBH013)the National Natural Science Foundation of China(52202405)。
文摘Cooperation among enterprises can bring overall and individual performance improvement,and a smooth coordination method is indispensable.However,due to the lack of customized coordination methods,cooperation in the downstream oil supply chain cannot be carried out smoothly.This paper intends to propose a multi-party coordination method to promote cooperation between oil shippers and pipeline operator by optimizing oil transportation,oil substitution and pipeline pricing schemes.An integrated game-theoretic modeling and analysis approach is developed to characterize the operation behaviors of all stakeholders in the downstream oil supply chain.The proposed mixed integer nonlinear programming model constrains supply and demand capacity,transportation routes,oil substitution rules and pipeline freight levels.Logarithm transformation and price discretization are introduced for model linear approximation.Simulation experiments are carried out in the oil distribution system in South China.The results show that compared to the business-as-usual scheme,the new scheme saves transportation cost by 3.48%,increases pipeline turnover by 5.7%,and reduces energy consumption and emissions by 7.66%and 6.77%.It is proved that the proposed method improves the revenue of the whole system,achieves fair revenue distribution,and also improves the energy and environmental benefits of the oil supply chain.
基金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.
基金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.
基金This work was supported by the National Natural Science Foundation of China(General Program)(Nos.12175114,U2167209)the Foundation of Key Laboratory of Metrology and Calibration Technology(No.JLKG2022001C001)+2 种基金the Platform Development foundation of China Institute for Radiation Protection(No.YP21030101)the National Key R&D Program of China(No.2021YFF0603600)the Tsinghua University Initiative Scientific Research Program(No.20211080081).
文摘Monte Carlo simulations are frequently utilized in radiation dose assessments.However,many researchers find the prevailing computing platforms to be intricate.This highlights a pressing need for a specialized framework for phantom dose evalua-tion.To address this gap,we developed a user-friendly radiation dose assessment platform using the Monte Carlo toolkit,Geant4.The Tsinghua University Phantom Dose(THUDosePD)augments the flexibility of Monte Carlo simulations in dosi-metric research.Originating from THUDose,a code with generic,functional,and application layers,THUDosePD focuses predominantly on anatomical phantom dose assessment.Additionally,it enables medical exposure simulation,intricate geometry creation,and supports both three-dimensional radiation dose analysis and phantom format transformations.The system operates on a multi-threaded parallel CPU architecture,with some modules enhanced for GPU parallel computing.Benchmark tests on the ICRP reference male illustrated the capabilities of THUDosePD in phantom dose assessment,covering the effective dose,three-dimensional dose distribution,and three-dimensional organ dose.We also conducted a voxelization conversion on the polygon mesh phantom,demonstrating the method’s efficiency and consistency.Extended applications based on THUDosePD further underline its broad adaptability.This intuitive,three-dimensional platform stands out as a valuable tool for phantom radiation dosimetry research.
基金supported by National Natural Science Foundation of China(No.U1966211)National Key R&D Program of China(No.2021YFB2401400)。
文摘Prediction models were proposed to estimate the reduced Townsend ionization coefficient and ionization cross-section.A shape function of the reduced Townsend ionization coefficient curves was derived from the ionization collision probability model.The function had three parameters:the first ionization potential energy,A_(α),and B_(α).A_(α)and B_(α)were related to the molecule symmetry and size.The polarization of molecules could characterize the molecule symmetry.The multi-layer molecular cross-section(MMCS)was proposed to describe the contributions of electrons and molecule radius on different molecule surfaces to collisions.A prediction model of the ionization cross-section was also proposed based on Aα.The molecule parameters were calculated by the Becke3–Lee–Yang–Parr(B3LYP)method and the 6–311G**basis set.We used available data of 30 and 23 gases,respectively,to build the prediction models of reduced Townsend ionization coefficients and ionization cross-sections.The relationships between the molecular parameters Aαand Bαand the ionization cross-section were built up via nonlinear fittings.The determination coefficients R^(2)of Aα,Bα,and the ionization cross-section were 0.877,0.887,and 0.838,respectively.The results showed that the accuracy of models was positively correlated with the molecule symmetry and reduced electric field.This was mainly related to the accuracy of the MMCS model in predicting Aα.The MMCS model needed to be improved to describe the collision direction selectivity caused by the molecule asymmetry.Under a high reduced electric field,that error of Aαhad less influence on the prediction results.However,the prediction results for single atoms with high symmetry were poor.This may be due to the absolute error of the model close to single atoms’reduced Townsend ionization coefficients.The models could provide the basis for gas insulation prediction and discharge calculations,especially for symmetric molecules under a high electric field.
基金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.
基金partially supported by the National Natural Science Foundation of China (51874325)the Science Foundation of China University of PetroleumBeijing (2462021BJRC009)。
文摘This paper intends to complete the primary logistics planning of oil products under the imbalance of supply and demand. An integrated mathematical programming model is developed to simultaneously find the balance between supply and demand, and optimize the logistics scheme. The model takes minimum logistics cost and resource adjustment cost as the objective function, and takes supply and demand capacity, transportation capacity, mass balance, and resource adjustment rules as constraints.Three adjustment rules are considered in the model, including resource adjustment within oil suppliers,within oil consumers, and between oil consumers. The model is tested on a large-scale primary logistics of a state-owned petroleum enterprise, involving 37 affiliated refineries, 31 procurement departments,286 market depots and dedicated consumers. After the unified optimization, the supply and demand imbalance is eased by 97% and the total cost is saved by 7%, which proves the effectiveness and applicability of the proposed model.
基金supported by National Natural Science Foundation of China(Nos.11175099&11355001)Tsinghua University Initiative Scientific Research Program(Nos.20151080354&2014Z21016)
文摘The broad-energy germanium(BEGe)detector,with the ability of background discrimination using pulse shape discrimination,is a competitive candidate for neutrinoless double beta decay(ovββ)experiments.In this paper,we report our measurements of key parameters for detector modeling in a commercial p-type BEGe detector.Point-like sources are used to investigate energy resolution and linearity of the detector.A cylindrical volume source is used for efficiency calibration.With an assembled device for source positioning and a collimated ^(133)Ba source,the detector is scanned to check its active volume.Using an^(241)Am point-like source,the dead layer thicknesses is measured at about 0.17 mm on the front and 1.18 mm on the side.The detector characterization is of importance for BEGe detectors to be used in the ovββ experiments at China JinPing underground Laboratory(CJPL).
文摘Hydrodesulfurization(HDS)of sour crude oil is an effective way to address the corrosion problems in refineries,and is an economic way to process sour crude oil in an existing refinery built for sweet oil.In the current study,the HDS of Siberian crude oil was carried out in a slurry reactor.The Co-Mo,Ni-Mo,and Ni-W catalysts supported onγ-Al2O3 were compared at the temperature of 340℃ and the pressure of 4.5 MPa.The HDS activity follows the order of Co-Mo N Ni-Mo N Ni-Wat a high concentration of H2S,and the difference between Co-Mo and Ni-Mo becomes insignificant at a low concentration of H2S.The influence of reaction temperature 320-360℃ and reaction pressure 3-5.5 MPa was investigated,and both play a positive role in the HDS reaction.A kinetic model over Ni-Mo/Al2O3 in the slurry reactor was established.The activation energy is estimated as 60.34 kJ·mol−1;the orders of sulfur components and hydrogen partial pressure are 1.43 and 1.30,respectively.The kinetic parameters are compared with those in a trickle-bed reactor,implying that the mass transfer is greatly enhanced in the slurry reactor.The back mixing effect is present in the slurry reactor and can be reduced by a multi-stage design,which would lead to higher reactor efficiency in industrial application.
基金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(No.11875036)。
文摘The Monte Carlo(MC)simulation is regarded as the gold standard for dose calculation in brachytherapy,but it consumes a large amount of computing resources.The development of heterogeneous computing makes it possible to substantially accelerate calculations with hardware accelerators.Accordingly,this study develops a fast MC tool,called THUBrachy,which can be accelerated by several types of hardware accelerators.THUBrachy can simulate photons with energy less than 3 MeV and considers all photon interactions in the energy range.It was benchmarked against the American Association of Physicists in Medicine Task Group No.43 Report using a water phantom and validated with Geant4 using a clinical case.A performance test was conducted using the clinical case,showing that a multicore central processing unit,Intel Xeon Phi,and graphics processing unit(GPU)can efficiently accelerate the simulation.GPU-accelerated THUBrachy is the fastest version,which is 200 times faster than the serial version and approximately 500 times faster than Geant4.The proposed tool shows great potential for fast and accurate dose calculations in clinical applications.