Paul K. Romano , Samuel Pasmann , Patrick C. Shriwise , Charles P.S. Swanson
{"title":"应用于蒙特卡罗粒子输运模拟的叠加网格计算材料体积分数","authors":"Paul K. Romano , Samuel Pasmann , Patrick C. Shriwise , Charles P.S. Swanson","doi":"10.1016/j.fusengdes.2025.115364","DOIUrl":null,"url":null,"abstract":"<div><div>We present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Storing the accumulated track-length data using an open-addressed hash table minimizes the potential memory footprint and, when paired with lock-free atomic compare-and-swap probing, removes thread contention to enable efficient parallel execution. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/<span><math><mi>R</mi></math></span>, where <span><math><mi>R</mi></math></span> is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.</div></div>","PeriodicalId":55133,"journal":{"name":"Fusion Engineering and Design","volume":"220 ","pages":"Article 115364"},"PeriodicalIF":2.0000,"publicationDate":"2025-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations\",\"authors\":\"Paul K. Romano , Samuel Pasmann , Patrick C. Shriwise , Charles P.S. Swanson\",\"doi\":\"10.1016/j.fusengdes.2025.115364\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><div>We present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Storing the accumulated track-length data using an open-addressed hash table minimizes the potential memory footprint and, when paired with lock-free atomic compare-and-swap probing, removes thread contention to enable efficient parallel execution. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/<span><math><mi>R</mi></math></span>, where <span><math><mi>R</mi></math></span> is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.</div></div>\",\"PeriodicalId\":55133,\"journal\":{\"name\":\"Fusion Engineering and Design\",\"volume\":\"220 \",\"pages\":\"Article 115364\"},\"PeriodicalIF\":2.0000,\"publicationDate\":\"2025-08-01\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Fusion Engineering and Design\",\"FirstCategoryId\":\"5\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0920379625005605\",\"RegionNum\":3,\"RegionCategory\":\"工程技术\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"NUCLEAR SCIENCE & TECHNOLOGY\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Fusion Engineering and Design","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0920379625005605","RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"NUCLEAR SCIENCE & TECHNOLOGY","Score":null,"Total":0}
Computing material volume fractions on a superimposed mesh as applied to Monte Carlo particle transport simulations
We present a newly implemented ray tracing algorithm in OpenMC for efficiently computing material volume fractions on superimposed meshes in complex geometries. By firing rays along each coordinate direction through the geometry, the approach accumulates track-length data in each mesh element, thereby determining the fractional composition of each material. Storing the accumulated track-length data using an open-addressed hash table minimizes the potential memory footprint and, when paired with lock-free atomic compare-and-swap probing, removes thread contention to enable efficient parallel execution. Scaling studies on three different models—a random tetrahedra configuration, the Frascati Neutron Generator ITER dose rate benchmark, and a stellarator design—show excellent parallel performance, with nearly linear speedup on modern multi-threaded and distributed-memory systems. An analysis of the residual error relative to high-resolution reference solutions demonstrated that under optimal conditions it decreases as 1/, where is the number of rays fired, making it straightforward to achieve user-prescribed accuracy. This new functionality enables practical, mesh-based approaches for detailed nuclear analyses in production Monte Carlo workflows without resorting to expensive, fully conformal or unstructured meshing.
期刊介绍:
The journal accepts papers about experiments (both plasma and technology), theory, models, methods, and designs in areas relating to technology, engineering, and applied science aspects of magnetic and inertial fusion energy. Specific areas of interest include: MFE and IFE design studies for experiments and reactors; fusion nuclear technologies and materials, including blankets and shields; analysis of reactor plasmas; plasma heating, fuelling, and vacuum systems; drivers, targets, and special technologies for IFE, controls and diagnostics; fuel cycle analysis and tritium reprocessing and handling; operations and remote maintenance of reactors; safety, decommissioning, and waste management; economic and environmental analysis of components and systems.